Development of design flood hydrographs using probability density functions

Probability density functions (PDFs) are used to fit the shape of hydrographs and have been popularly used for the development of synthetic unit hydrographs by many hydrologists. Nevertheless, modelling the shapes of continuous stream flow hydrographs, which are probabilistic in nature, is rare. In the present study, a novel approach was followed to model the shape of stream flow hydrographs using PDF and subsequently to develop design flood hydrographs for various return periods. Four continuous PDFs, namely, two parameter Beta, Weibull, Gamma and Lognormal, were employed to fit the shape of the hydrographs of 22 years at a site of Brahmani River in eastern India. The shapes of the observed and PDF fitted hydrographs were compared and root mean square errors, error of peak discharge (EQ_P) and error of time to peak (ET_P) were computed. The best‐fitted shape and scale parameters of all PDFs were subjected to frequency analysis and the quartiles corresponding to 20‐, 50‐, 100‐ and 200‐year were estimated. The estimated parameters of each return period were used to develop the flood hydrographs for 20‐, 50‐, 100‐ and 200‐year return periods. The peak discharges of the developed design flood hydrographs were compared with the design discharges estimated from the frequency analysis of 22 years of annual peak discharges at that site. Lognormal‐produced peak discharge was very close to the estimated design discharge in case of 20‐year flood hydrograph. On the other hand, peak discharge obtained using the Weibull PDF had close agreement with the estimated design discharge obtained from frequency analysis in case of 50‐, 100‐ and 200‐year return periods. The ranking of the PDFs based on estimation of peak of design flood hydrograph for 50‐, 100‐ and 200‐year return periods was found to have the following order: Weibull > Beta > Lognormal > Gamma. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Effects of hillslope dynamics and network geometry on the scaling properties of the hydrologic response

This paper investigates the specific contributions of river network geomorphology, hillslope flow dynamics and channel routing to the scaling behavior of the hydrologic response as function of drainage area. Scaling relationships emerged from the observations of geomorphological and hydrological data and were reproduced in previous works through mathematical models, for both idealized self-similar networks and natural basins. Recent literature highlighted that scale invariance of hydrological quantities depends not only on the metrics of the drainage catchment but also on effective flow routing. In this study we employ a geomorphological width function scheme to test the simple scaling hypothesis adopting more realistic dynamic conditions than in previous approaches, specifically taking into account the role of hillslopes. The analysis is based on the derivation of the characteristic distributions of path lengths and travel times, inferred from DEM processing and measurements of rainfall and runoff data. The study area is located in the Tiber River region (central Italy).Results indicate that, while scaling properties clearly emerge when the hydrologic response is defined on the basis of the sole geomorphology, scale invariance is broken when less idealized flow dynamics are taken into account. Lack of scaling appears in particular as a consequence of the catchment to catchment variability of hillslope velocities.     

Effect of growing watershed imperviousness on hydrograph parameters and peak discharge

An increasing impervious area is quickly extending over the Wu‐Tu watershed due to the endless demands of the people. Generally, impervious paving is a major result of urbanization and more recently has had the potential to produce more enormous flood disasters than those of the past. In this study, 40 available rainfall–runoff events were chosen to calibrate the applicable parameters of the models and to determine the relationships between the impervious surfaces and the calibrated parameters. Model inputs came from the outcomes of the block kriging method and the non‐linear programming method. In the optimal process, the shuffled complex evolution method and three criteria were applied to compare the observed and simulated hydrographs. The tendencies of the variations of the parameters with their corresponding imperviousness were established through regression analysis. Ten cases were used to examine the established equations of the parameters and impervious covers. Finally, the design flood routines of various return periods were furnished through use of approaches containing a design storm, block kriging, the SCS model, and a rainfall‐runoff model with established functional relationships. These simulated flood hydrographs were used to compare and understand the past, present, and future hydrological conditions of the watershed studied. In the research results, the time to peak of flood hydrographs for various storms was diminished approximately from 11 h to 6 h in different decrements, whereas peak flow increased respectively from 127 m³ s^?1 to 629 m³ s^?1 for different storm intensities. In addition, this study provides a design diagram for the peak flow ratio to help engineers and designers to construct hydraulic structures efficiently and prevent possible damage to human life and property. Copyright © 2007 John Wiley & Sons, Ltd.     

Nitrate removal processes in the riverbed during a single-peak flood event

The effect of the interplay between unsteady flow and bedform in a flood event on nitrogen cycling in the hyporheic zone (HZ) remains poorly understood. In this study, a reactive transport groundwater model with different flood hydrographs was proposed to investigate the effect of modified hyporheic flow on nitrate dynamics in the HZ, including nitrate source-sink function, response to the single-peak flood event and removal efficiency. The results demonstrate that there exists an optimal range of river channel gradients that could enhance the biogeochemical reactions (respiration, nitrification and denitrification) in a flood event. The HZ acts as a nitrate sink especially after the flood event, and its source-sink function is independent of the unsteady discharge/stage conditions. The nitrate in the HZ has a hysteretic response to peak stage/discharge, and its removal efficiency is decreased by up to 70% compared to steady flow conditions. These findings not only provide a better understanding of nitrogen dynamics under the effect of unsteady channel flow, but also can be applied for river restoration to efficiently remove nitrate in the HZ by modifying river channel gradients.     

Improvement of surface run‐off in the hydrological model ParFlow by a scale‐consistent river parameterization

We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels.     

A simulation experiment for optimal design hyetograph selection

The aim of this work is to assess the accuracy of literature design hyetographs for the evaluation of peak discharges during flood events. Five design hyetographs are examined in a set of simulations, based upon the following steps: (i) an ideal river basin is defined, characterized by a Beta distribution shaped unit hydrograph (UH); (ii) 1000 years of synthetic rainfall are artificially generated; (iii) a discharge time‐series is obtained from the convolution of the rainfall time‐series and the UH, and the reference T‐years flood is computed from this series; (iv) for the same return period T, the parameters of the intensity–duration–frequency (IDF) curve are estimated from the 1000 years of synthetic rainfall; (v) five design hyetographs are determined from the IDF curves and are convolved with the discrete UH to find the corresponding design hydrographs; (vi) the hydrograph peaks are compared with the reference T‐years flood and the advantages and drawbacks of each of the five approaches are evaluated. The rainfall and UH parameters are varied, and the whole procedure is repeated to assess the sensitivity of results to the system configuration. We found that all design hyetographs produce flood peak estimates that are consistently biased in most of the climatic and hydrologic conditions considered. In particular, significant underestimation of the design flood results from the adoption of any rectangular hyetograph used in the context of the rational formula. In contrast, the Chicago hyetograph tends to overestimate peak flows. In two cases it is sufficient to multiply the result by a constant scaling factor to obtain robust and nearly unbiased estimates of the design floods. Copyright © 2007 John Wiley & Sons, Ltd.     

A New Method for Estimating Recharge to Unconfined Aquifers Using Differential River Gauging

In semiarid and arid environments, leakage from rivers is a major source of recharge to underlying unconfined aquifers. Differential river gauging is widely used to estimate the recharge. However, the methods commonly applied are limited in that the temporal resolution is event‐scale or longer. In this paper, a novel method is presented for quantifying both the total recharge volume for an event, and variation in recharge rate during an event from hydrographs recorded at the upstream and downstream ends of a river reach. The proposed method is applied to river hydrographs to illustrate the method steps and investigate recharge processes occurring in a sub‐catchment of the Murray Darling Basin (Australia). Interestingly, although it is the large flood events which are commonly assumed to be the main source of recharge to an aquifer, our analysis revealed that the smaller flow events were more important in providing recharge.     

Bayesian estimation of inflow hydrographs in ungauged sites of multiple reach systems

A Bayesian Geostatistical Approach to evaluate unknown upstream flow hydrographs in multiple reach systems is implemented. The methodology was, firstly, tested through three synthetic examples of river confluences, that differ in the available data, boundary conditions and number of the estimated inflow time series. Input discharge hydrographs were routed downstream by means of the widely known HEC-RAS river analysis system to obtain the downstream stage hydrographs used as known observations for the reverse procedure. In almost all cases, the observed water levels were corrupted with random errors to highlight the reliability of the methodology in preventing instabilities and overfitting. Then the procedure was applied to the real case study of the Parma–Baganza river confluence located at the city of Parma (Italy) to assess the tributary Baganza River inflow hydrograph (supposed completely ungauged) using water level data collected downstream on the main reach. The results show that the methodology properly reproduces the unknown inflows even in presence of errors affecting the downstream water levels. The practical applicability of the proposed approach is also demonstrated in complex river systems.     

Role of hypsometry and planform in basin hydrologic response

The effects of variations of drainage basin morphometry and relief characteristics on flood peak magnitude and time-to-peak are investigated using simulated stream networks. The networks are produced by three models: headward growth, systematic capture, and minimum power relaxation. Translational and kinematic wave flood routing were used to generate synthetic hydrographs. Peak discharge and time-to-peak are predictable to a high degree by five different sets of morphometric-relief parameters. In order of decreasing order of importance in predictive ability the parameters characterize basin size, relative relief, basin concavity, and basin shape. Both simulated and natural stream networks exhibit strong dependence of planimetric morphometry upon basin concavity. The effect of this dependency is to increase the effect of basin concavity upon flood hydrographs.     

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.     

Simulation of flood hydrographs based on lumped and semi-distributed models for two tropical catchments in Kenya

Abstract

The design and operation of flood management systems require computation of flood hydrographs for both design floods and flood forecasting purposes, since observed data are usually inadequate for these tasks. This is particularly relevant for most developing countries, i.e. mainly for tropical catchments. One possible way of obtaining information about flood hydrographs is through the use of rainfall—runoff models. Two such models, namely the Bochum model and the Nash Cascade—Diskin Infiltration model, which are semi-distributed and lumped models, respectively, were used in the present study. These models were applied to two catchments in Kenya with drainage areas of 6.71 km² and 26.03 km². A set of 13 selected rainfall—runoff events was used to calibrate and validate the models. The physical parameters required by the models were derived from catchment characteristics using GIS and remote sensing data while the conceptual parameters were obtained by optimization. The flood hydrographs simulated using the parameters so derived indicated that it is possible to use the two models in this tropical environment.     

A real-time stage Muskingum forecasting model for a site without rating curve

Abstract

An adaptive model for on-line stage forecasting is proposed for river reaches where significant lateral inflow contributions occur. The model is based on the Muskingum method and requires the estimation of four parameters if the downstream rating curve is unknown; otherwise only two parameters have to be determined. As the choice of the forecast lead time is linked to wave travel time along the reach, to increase the lead time, a schematization of two connected river reaches is also investigated. The variability of lateral inflow is accounted for through an on-line adaptive procedure. Calibration and validation of the model were carried out by applying it to different flood events observed in two equipped river reaches of the upper-middle Tiber basin in central Italy, characterized by a significant contributing drainage area. Even if the rating curve is unknown at the downstream section, the forecast stage hydrographs were found in good agreement with those observed. Errors in peak stage and time to peak along with the persistence coefficient values show that the model has potential as a practical tool for on-line flood risk management.     

Dissecting the effect of rainfall variability on the statistical structure of peak flows

This study examines the role of rainfall variability on the spatial scaling structure of peak flows using the Whitewater River basin in Kansas as an illustration. Specifically, we investigate the effect of rainfall on the scatter, the scale break and the power law (peak flows vs. upstream areas) regression exponent. We illustrate why considering individual hydrographs at the outlet of a basin can lead to misleading interpretations of the effects of rainfall variability. We begin with the simple scenario of a basin receiving spatially uniform rainfall of varying intensities and durations and subsequently investigate the role of storm advection velocity, storm variability characterized by variance, spatial correlation and intermittency. Finally, we use a realistic space–time rainfall field obtained from a popular rainfall model that combines the aforementioned features. For each of these scenarios, we employ a recent formulation of flow velocity for a network of channels, assume idealized conditions of runoff generation and flow dynamics and calculate peak flow scaling exponents, which are then compared to the scaling exponent of the width function maxima. Our results show that the peak flow scaling exponent is always larger than the width function scaling exponent. The simulation scenarios are used to identify the smaller scale basins, whose response is dominated by the rainfall variability and the larger scale basins, which are driven by rainfall volume, river network aggregation and flow dynamics. The rainfall variability has a greater impact on peak flows at smaller scales. The effect of rainfall variability is reduced for larger scale basins as the river network aggregates and smoothes out the storm variability. The results obtained from simple scenarios are used to make rigorous interpretations of the peak flow scaling structure that is obtained from rainfall generated with the space–time rainfall model and realistic rainfall fields derived from NEXRAD radar data.     

Evidence of Self‐Organized Criticality in riverbank mass failures: a matter of perspective?

A growing body of field, theoretical and numerical modelling studies suggests that predicting river response to even major changes in input variables is difficult. Rivers are seen to adjust rapidly and variably through time and space as well as changing independently of major driving variables. Concepts such as Self‐Organized Criticality (SOC) are considered to better reflect the complex interactions and adjustments occurring in systems than traditional approaches of cause and effect. This study tests the hypothesis that riverbank mass failures which occurred both prior to, and during, an extreme flood event in southeast Queensland (SEQ) in 2011 are a manifestation of SOC. Each wet‐flow failure is somewhat analogous to the ‘avalanche’ described in the initial sand‐pile experiments of Bak et al. (Physical Review Letters, 1987, 59(4), 381–384) and, due to the use of multitemporal LiDAR, the time period of instability can be effectively constrained to that surrounding the flood event. The data is examined with respect to the key factors thought to be significant in evaluating the existence of SOC including; non‐linear temporal dynamics in the occurrence of disturbance events within the system; an inverse power‐law relation between the magnitude and frequency of the events; the existence of a critical state to which the system readjusts after a disturbance; the existence of a cascading processes mechanism by which the same process can initiate both low‐magnitude and high‐magnitude events. While there was a significant change in the frequency of mass failures pre‐ and post‐flood, suggesting non‐linear temporal dynamics in the occurrence of disturbance events, the data did not fit an inverse power‐law within acceptable probability and other models were found to fit the data better. Likewise, determining a single ‘critical’ state is problematic when a variety of feedbacks and multiple modes of adjustment are likely to have operated throughout this high magnitude event. Overall, the extent to which the data supports a self‐organized critical state is variable and highly dependent upon inferential arguments. Investigating the existence of SOC, however, provided results and insights that are useful to the management and future prediction of these features. Copyright © 2014 John Wiley & Sons, Ltd.     

Nonlinear transformation of unit hydrograph

Unit hydrograph (UH) and its numerous derivatives have been popular for estimation of flood hydrographs. Two major assumptions still overshadow UH applications. One is the linearity and the other is time invariance. In theory, only peak discharge of an equilibrium hydrograph follows linear proportionality to excess rainfall intensity. In trying to relax the linearity constraint, this paper aims to propose a nonlinear way of transforming a given UH to other general hydrographs. The transformation or mapping technique relies on a simple rainfall ratio raised to a power less than unity. The case of nonlinear transformation is illustrated for a number of watershed geometries with either known kinematic wave analytic solutions or observed data. The nonlinear UH approach also relaxes the assumption of constant time base of the UH. The proposed nonlinear UH transformation may thus be viewed as a major step in closing the gap between physically based and traditional UH-based surface runoff simulation approaches.     

A hydro‐economic modelling framework for flood damage estimation and the role of riparian vegetation

A modelling framework for the quick estimate of flood inundation and the resultant damages is developed in this paper. The model, called the flood economic impact analysis system (FEIAS), can be applied to a river reach of any hydrogeological river basin. For the development of the integrated modelling framework, three models were employed: (1) a modelling scheme based on the Hydrological Simulation Program FORTRAN model that was developed for any geomorphological river basin, (2) a river flow/floodplain model, and (3) a flood loss estimation model. The first sub‐model of the flood economic impact analysis system simulates the hydrological processes for extended periods of time, and its output is used as input to a second component, the river/floodplain model. The hydraulic model MIKE 11 (quasi‐2D) is the river/floodplain model employed in this study. The simulated flood parameters from the hydraulic model MIKE 11 (quasi‐2D) are passed, at the end of each time step, to a third component, the flood loss model for the estimation of flood damage. In the present work, emphasis was given to the seasonal variation of Manning's coefficient (n), which is an important parameter for the determination of the flood inundation in hydraulic modelling. High values of Manning's coefficient for a channel indicate high flow resistance. The riparian vegetation can have a large impact on channel resistance. The modelling framework developed in this paper was used to investigate the role of riparian vegetation in reducing flood damage. Moreover, it was used to investigate the influence of cutting riparian vegetation scenarios on the flow characteristics. The proposed framework was applied to the downstream part of the Koiliaris River basin in Crete, Greece, and was tested and validated with historical data. Copyright © 2012 John Wiley & Sons, Ltd.     

Extension of observed flood series by combining a distributed hydro-meteorological model and a copula-based model

Long flood series are required to accurately estimate flood quantiles associated with high return periods, in order to design and assess the risk in hydraulic structures such as dams. However, observed flood series are commonly short. Flood series can be extended through hydro-meteorological modelling, yet the computational effort can be very demanding in case of a distributed model with a short time step is considered to obtain an accurate flood hydrograph characterisation. Statistical models can also be used, where the copula approach is spreading for performing multivariate flood frequency analyses. Nevertheless, the selection of the copula to characterise the dependence structure of short data series involves a large uncertainty. In the present study, a methodology to extend flood series by combining both approaches is introduced. First, the minimum number of flood hydrographs required to be simulated by a spatially distributed hydro-meteorological model is identified in terms of the uncertainty of quantile estimates obtained by both copula and marginal distributions. Second, a large synthetic sample is generated by a bivariate copula-based model, reducing the computation time required by the hydro-meteorological model. The hydro-meteorological modelling chain consists of the RainSim stochastic rainfall generator and the Real-time Interactive Basin Simulator (RIBS) rainfall-runoff model. The proposed procedure is applied to a case study in Spain. As a result, a large synthetic sample of peak-volume pairs is stochastically generated, keeping the statistical properties of the simulated series generated by the hydro-meteorological model. This method reduces the computation time consumed. The extended sample, consisting of the joint simulated and synthetic sample, can be used for improving flood risk assessment studies.     

Flood hydrograph reconstruction from the peak flow value in ephemeral streams using a simplified robust single‐parameter model

Flood hydrographs from ephemeral streams in arid areas provide valuable information for assessing run‐off and groundwater recharge. However, such data are often scarce or incomplete, especially in hyper‐arid regions. The hypothesis of this study was that it is possible to reconstruct a hydrograph of a specific point along an ephemeral stream with the knowledge of only the peak flow rate of a flood event at that point and that this can be done at almost every point along the stream. The feasibility of this approach lies in the shape of the recession stage of the flood hydrograph, which is known to be a repeating phenomenon. The recession stage comes immediately after the peak flow rate, when it begins its decline, and lasts until the flood is extinguished. A general shape of the flood recession stage can be provided. Because the recession stage represents ~80% of the duration of a flood event, it can provide a general idea of the flood hydrograph's shape. A simple model based on geometric progression is suggested to describe the repeating recession stage of a flood. The advantage of the proposed model is that it requires only one parameter: the recession characteristic at a fixed point along the ephemeral stream, termed recession coefficient q. By knowing the recession coefficient of a fixed point and the peak flow rate of a flood event at that point, one can plot the flood hydrograph. A good agreement is shown between the observed and computed values of the recession stage. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling time varying scouring at bed sills

In this paper a modelling approach is presented to predict local scour under time varying flow conditions. The approach is validated using experimental data of unsteady scour at bed sills. The model is based on a number of hypotheses concerning the characteristics of the flow hydrograph, the temporal evolution of the scour and the geometry of the scour hole. A key assumption is that, at any time, the scour depth evolves at the same rate as in an equivalent steady flow. The assumption is supported by existing evidence of geometrical affinity and similarity of scour holes formed under different steady hydraulic conditions. Experimental data are presented that show the scour hole development downstream of bed sills due to flood hydrographs follow a predictable pattern. Numerical simulations are performed with the same input parameters used in the experimental tests but with no post‐simulation calibration. Comparison between the experimental and model results indicates good correspondence, especially in the rising limb of the flow hydrograph. This suggests that the underlying assumptions used in the modelling approach are appropriate. In principle, the approach is general and can be applied to a wide range of environments (e.g. bed sills, step‐pool systems) in which scouring at rapid bed elevation changes caused by time varying flows occurs, provided appropriate scaling information is available, and the scour response to steady flow conditions can be estimated. Copyright © 2011 John Wiley & Sons, Ltd.