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.     

FLOOD MODELLING WITH GIS-DERIVED DISTRIBUTED UNIT HYDROGRAPHS

The concept of a spatially distributed unit hydrograph is based on the fact that the unit hydrograph can be derived from the time–area curve of a watershed by the S-curve method. The time–area diagram is a graph of cumulative drainage area contributing to discharge at the watershed outlet within a specified time of travel. Accurate determination of the time–area diagram is made possible by using a GIS. The GIS is used to describe the connectivity of the links in the watershed flow network and to calculate distances and travel times to the watershed outlet for various points within the watershed. Overland flow travel times are calculated by the kinematic wave equation for time to equilibrium; channel flow times are based on the Manning and continuity equations. To account for channel storage, travel times for channel reaches are increased by a percentage depending on the channel reach length and geometry. With GIS capability for rainfall mapping, the assumption of a uniform spatial rainfall distribution is no longer necessary; hence the term, spatially distributed unit hydrograph. An example of the application for the Waiparous Creek in the Alberta Foothills is given. IDRISI is used to develop a simple digital elevation model of the 229 km² watershed, using 1 km × 1 km grid cells. A grid of flow directions is developed and used to create an equivalent channel network. Excess rainfall for each 1 km × 1 km cell is individually computed by the Soil Conservation Service (SCS) runoff curve method and routed through the equivalent channel network to obtain the time–area curve. The derived unit hydrograph gave excellent results in simulating an observed flood hydrograph. The distributed unit hydrograph is no longer a lumped model, since it accounts for internal distribution of rainfall and runoff. It is derived for a watershed without the need for observed rainfall and discharge data, because it is essentially a geomorphoclimatic approach. As such, it allows the derivation of watershed responses (hydrographs) to inputs of various magnitudes, thus eliminating the assumption of proportionality of input and output if needed. The superposition of outputs is retained in simulating flood hydrographs by convolution, since it has been shown that some non-linear systems satisfy the principle of superposition. The distributed unit hydrograph appears to be a very promising rainfall runoff model based on GIS technology.     

Review of field observations of runoff generation in the tropics

Over the last thirty years, French hydrologists have observed flow characteristics at the outlets of small tropical watersheds, characteristics which vary between arid zones, humid zones and equatorial rainforest zones. The main physical phenomena related to the various types of runoff generation have been observed and analysed. The generation of runoff and subsurface runoff following tropical rainstorms depends on soil surface formations, texture and structure as well as vegetation cover, geological substrata, weathered material and climate.

Analysis is also made of experimental rainfall-runoff relationships, the shape of flood hydrographs and the conditions necessary for the occurrence of steady base flow. Scale effects between experimental plot and watershed have been studied to identify any effect this may have upon runoff generation and flow. Tables detail watershed characteristics and relationships between these and rainfall and runoff.     

Geospatially based distributed rainfall‐runoff modelling for simulation of internal and outlet responses in a semi‐forested lower Himalayan watershed

This study presents a Geographic Information System (GIS)‐based distributed rainfall‐runoff model for simulating surface flows in small to large watersheds during isolated storm events. The model takes into account the amount of interception storage to be filled using a modified Merriam ( 1960 ) approach before estimating infiltration by the Smith and Parlange ( 1978 ) method. The mechanics of overland and channel flow are modelled by the kinematic wave approximation of the Saint Venant equations which are then numerically solved by the weighted four‐point implicit finite difference method. In this modelling the watershed was discretized into overland planes and channels using the algorithms proposed by Garbrecht and Martz ( 1999 ). The model code was first validated by comparing the model output with an analytical solution for a hypothetical plane. Then the model was tested in a medium‐sized semi‐forested watershed of Pathri Rao located in the Shivalik ranges of the Garhwal Himalayas, India. Initially, a local sensitivity analysis was performed to identify the parameters to which the model outputs like runoff volume, peak flow and time to peak flow are sensitive. Before going for model validation, calibration was performed using the Ordered‐Physics‐based Parameter Adjustment (OPPA) method. The proposed Physically Based Distributed (PBD) model was then evaluated both at the watershed outlet as well as at the internal gauging station, making this study a first of its kind in Indian watersheds. The results of performance evaluation indicate that the model has simulated the runoff hydrographs reasonably well within the watershed as well as at the watershed outlet with the same set of calibrated parameters. The model also simulates, realistically, the temporal variation of the spatial distribution of runoff over the watershed and the same has been illustrated graphically. Copyright © 2011 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.     

UPLAND EROSION MODELING WITH CASC2D-SED

1 INTRODUCTION Erosion and sedimentation embody the processes of detachment, transport, and deposition of soil particles. Erosion and subsequent deposition can cause major problems. Erosion reduces productivity of cropland, sediment degrades water quality and may carry soil adsorbed polluting chemicals. Deposition in irrigation canals, stream channels and reservoirs reduces structural capacity and requires costly removal. Ideally, an erosion model should describe the physical processes…     

Reliability analysis of hydraulic structures considering unit hydrograph uncertainty

Unit hydrographs (UHs), along with design rainfalls, are frequently used to determine the discharge hydrograph for design and evaluation of hydraulic structures. Due to the presence of various uncertainties in its derivation, the resulting UH is inevitably subject to uncertainty. Consequently, the performance of hydraulic structures under the design storm condition is uncertain. This paper integrates the linearly constrained Monte-Carlo simulation with the UH theory and routing techniques to evaluate the reliability of hydraulic structures. The linear constraint is considered because the water volume of each generated design direct runoff hydrograph should be equal to that of the design effective rainfall hyetograph or the water volume of each generated UH must be equal to one inch (or cm) over the watershed. For illustration, the proposed methodology is applied to evaluate the overtopping risk of a hypothetical flood detention reservoir downstream of Tong-Tou watershed in Taiwan.     

Impact assessment of climatic and land-use changes on flood runoff in southeast Queensland

Abstract

Over the past century, land-use has changed in southeast Queensland, and when coupled with climatic change, the risk of flooding has increased. This research aims to examine impacts of climate and land-use changes on flood runoff in southeast Queensland, Australia. A rainfall–runoff model, RORB, was calibrated and validated using observed flood hydrographs for one rural and one urbanized catchment, for 1961–1990. The validated model was then used to generate flood hydrographs using projected rainfall based on two climate models: the Geophysical Fluid Dynamics Laboratory Climate Model 2.1 (GFDL CM2.1) and the Conformal-Cubic Atmospheric Model (CCAM), for 2016–2045. Projected daily rainfall for the two contrasting periods was used to derive adjustment factors for a given frequency of occurrence. Two land-use change scenarios were used to evaluate likely impacts. Based on the projected rainfall, the results showed that, in both catchments, future flood magnitudes are unlikely to increase for large flood events. Extreme land-use change would significantly impact flooding in the rural catchment, but not the urbanized catchment.

Editor Z.W. Kundzewicz; Associate editor Y. Gyasi-Agyei     

Reliability analysis of hydraulic structures considering unit hydrograph uncertainty

Unit hydrographs (UHs), along with design rainfalls, are frequently used to determine the discharge hydrograph for design and evaluation of hydraulic structures. Due to the presence of various uncertainties in its derivation, the resulting UH is inevitably subject to uncertainty. Consequently, the performance of hydraulic structures under the design storm condition is uncertain. This paper integrates the linearly constrained Monte-Carlo simulation with the UH theory and routing techniques to evaluate the reliability of hydraulic structures. The linear constraint is considered because the water volume of each generated design direct runoff hydrograph should be equal to that of the design effective rainfall hyetograph or the water volume of each generated UH must be equal to one inch (or cm) over the watershed. For illustration, the proposed methodology is applied to evaluate the overtopping risk of a hypothetical flood detention reservoir downstream of Tong-Tou watershed in Taiwan.     

CHARACTERISTICS OF SEDIMENT TRANSPORT ALONG A RIVER REACH WITH A RESERVOIR

Based on long-term measurements at three gauging stations, Toudaoguai, Fugu and Hequ, and one meteorological station, this article discusses the features of discharge (Q) and sediment concentration (Cs) of a river reach of the Yellow River with a reservoir located in the Loess Plateau. The impacts of the local sub-watershed between Toudaoguai and Fugu gauging stations on sediment budget to the Yellow River have been analyzed. In addition, the deposition processes in the Tianqiao Reservoir have been investigated. Results show over 80% of the precipitation that falls in the local subwatershed is unable to contribute to the Yellow River runoff process. It is found that the annualmaximum sediment concentration is usually less than 30 kg/m^3 during flood seasons at Toudaoguai Gauging Station, but the sediment concentration varies dramatically at Fugu Gauging Station. About 35% of the sediment eroded in the sub-watersheds between Toudaoguai and Fugu gauging stationswas produced from the Huangfuchuan sub-watershed which has a drainage area accounting only for 10% of the drainage area between Toudaoguai and Fugu gauging stations. The Tianqiao Reservoir generally has deposition during the summer flood season, and scouring during the non-flood season.On average, over 85% of deposited sediment in the reservoir occurs in the 12 km long lower reservoir reach. The volume of annual deposition in the reservoir mainly depends on the volume of water from the local region between Hequ and Fugu gauging stations.     

Assessment on the function of reservoirs for flood control during typhoon seasons based on a distributed hydrological model

A grid‐based distributed hydrological model, PDTank model, is used to simulate hydrological processes in the upper Tone River catchment. The Tone River catchment often suffers from heavy rainfall events during the typhoon seasons. The reservoirs located in the catchment play an important role in flood regulation. Through the coupling of the PDTank model and a reservoir module that combines the storage function and operation function, the PDTank model is used for flood forecasting in this study. By comparing the hydrographs simulated using gauging and radar rainfall data, it is found that the spatial variability of rainfall is an important factor for flood simulation and the accuracy of the hydrographs simulated using radar rainfall data is slightly improved. The simulation of the typhoon flood event numbered No. 9 shows that the reservoirs in the catchment attenuate the peak flood discharge by 423·3 m³/s and validates the potential applicability of the distributed hydrological model on the assessment of function of reservoirs for flood control during typhoon seasons. Copyright © 2011 John Wiley & Sons, Ltd.     

Development of a thresholds approach for real‐time flash flood prediction in complex geomorphological river basins

Flash floods represent one of the deadliest and costliest natural disasters worldwide. The hydrological analysis of a flash flood event contributes in the understanding of the runoff creation process. This study presents the analysis of some flash flood events that took place in a complex geomorphological Mediterranean River basin. The objective of the present work is to develop the thresholds for a real‐time flash flood forecasting model in a complex geomorphological watershed, based on high‐frequency data from strategically located hydrological and meteorological telemetric stations. These stations provide hourly real‐time data which were used to determine hydrological and meteorological parameters. The main characteristics of various hydrographs specified in this study were the runoff coefficients, lag time, time to peak, and the maximum potential retention. The estimation of these hydrometeorological parameters provides the necessary information in order to successfully manage flash floods events. Especially, the time to peak is the most significant hydrological parameter that affects the response time of an oncoming flash flood event. A study of the above parameters is essential for the specification of thresholds which are related to the geomorphological characteristics of the river basin, the rainfall accumulation of an event, the rainfall intensity, the threshold runoff through karstic area, the season during which the rainfall takes place and the time intervals between the rainstorms that affect the soil moisture conditions. All these factors are combined into a real‐time‐threshold flash flood prediction model. Historical flash flood events at the downstream are also used for the validation of the model. An application of the proposed model is presented for the Koiliaris River basin in Crete, Greece. Copyright © 2011 John Wiley & Sons, Ltd.     

On the relationship between hydrographs and chemographs

The spatial representativeness of gauging stations was investigated in two low‐mountainous river basins near the city of Trier, southwest Germany. Longitudinal profiles during low and high flow conditions were sampled in order to identify sources of solutes and to characterize the alteration of flood wave properties during its travel downstream. Numerous hydrographs and chemographs of natural flood events were analysed in detail. Additionally, artificial flood events were investigated to study in‐channel transport processes. During dry weather conditions the gauging station was only representative for a short river segment upstream, owing to discharge and solute concentrations of sources contiguous to the measurement site. During artificial flood events the kinematic wave velocity was considerably faster than the movement of water body and solutes, refuting the idea of a simple mixing process of individual runoff components. Depending on hydrological boundary conditions, the wave at a specific gauge could be entirely composed of old in‐channel water, which notably reduces the spatial representativeness of a sampling site. Natural flood events were characterized by a superimposition of local overland flow, riparian water and the kinematic wave process comprising the downstream conveyance of solutes. Summer floods in particular were marked by a chronological occurrence of distinct individual runoff components originating only from a few contributing areas adjacent to the stream and gauge. Thus, the representativeness of a gauge for processes in the whole basin depends on the distance of the nearest significant source to the station. The consequence of our study is that the assumptions of mixing models are not satisfied in river basins larger than 3 km². Copyright © 2006 John Wiley & Sons, Ltd.     

Reservoir attenuation of floods from ungauged basins

Abstract

In a typical reservoir routing problem, the givens are the inflow hydrograph and reservoir characteristic functions. Flood attenuation investigations can be easily accomplished using a hydrological or hydraulic routing of the inflow hydrograph to obtain the reservoir outflow hydrograph, unless the inflow hydrograph is unavailable. Although attempts for runoff simulation have been made in ungauged basins, there is only a limited degree of success in special cases. Those approaches are, in general, not suitable for basins with a reservoir. The objective of this study is to propose a procedure for flood attenuation estimation in ungauged reservoir basins. In this study, a kinematic-wave based geomorphic IUH model was adopted. The reservoir inflow hydrograph was generated through convolution integration using the rainfall excess and basin geomorphic information. Consequently, a fourth-order Runge-Kutta method was used to route the inflow hydrograph to obtain the reservoir outflow hydrograph without the aid of recorded flow data. Flood attenuation was estimated through the analysis of the inflow and outflow hydrographs of the reservoir. An ungauged reservoir basin in southern Taiwan is presented as an example to show the applicability of the proposed analytical procedure. The analytical results provide valuable information for downstream flood control work for different return periods.     

Synthetic hydrograph generation by hydrological donors

A method to build synthetic hydrographs is introduced, based on 1300 gauging stations in France and Switzerland, covering a wide range of size and climatology. For each station, an average of two floods per year are selected by a peak-over-threshold method, providing about 69?000 hydrographs. For a given catchment, some “donor stations” are selected with criteria of proximity in space, size and runoff production. These donors provide hundreds of hydrographs which can complement the ones recorded locally, or replace them if no hydrograph is available. With this subset of hydrographs, one can estimate the catchment’s average peak-to-volume ratio of floods, and build the corresponding median hydrograph. Another application is, for a given daily discharge sequence (being observed or simulated), to generate a relevant synthetic hydrograph by combining appropriate hydrographs of the subset. These methods are assessed by performing a jack-knife validation on a wide dataset of stations.     

An approach for evaluating the hydrological effects of urbanization and its application

This study focuses mainly on observing urban development in Taiwan's Wu‐Tu watershed from the perspective of urban hydrological theory. An approach is proposed for developing a method for incorporating available meteorological data to define the degree of change in a runoff hydrograph for urbanizing basins. The mean rainfall was estimated using the Kriging method. For calibration, two methods of calculating the effective rainfall (the Φ‐index method and the non‐linear‐programming (NLP) method) were used as model inputs, and the optimal global parameters of the linear reservoir model were then obtained from the shuffled complex evolution (SCE) algorithm. Twenty‐six (1966–1991) and eight (1994–1997) rainfall–runoff events were used for calibration and verification, respectively. The NLP method yielded better results than the Φ‐index method, especially for multipeak rainfall–runoff events. The regression equation determined the relationship between the parameters of the model and impervious areas. A comparison based on the results of the instantaneous unit hydrograph of the study area revealed that three decades of urbanization had increased the peak flow by 27%, and the time to peak was decreased by 4 h. The study simply describes the results of the impact of imperviousness on hydrological modelling. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of moving storms on attainment of equilibrium discharge

Although rainfall is assumed spatially uniform in conventional hydrological modelling for rainfall–runoff simulations, moving storms have been shown to have substantial influence on flow hydrographs. In this study, criteria for attainment of the equilibrium discharge from watersheds subjected to moving storms were examined. Non-linear numerical kinematic-wave models were developed to simulate runoff from an overland plane and from a V-shaped catchment. Dimensional analysis was applied to obtain the independent variables to be used as control factors in performing a series of numerical tests. The results indicate that, for storms moving downstream, runoff can attain equilibrium discharge even though the storm length is shorter than the watershed length and the rainfall duration is less than the time to equilibrium of the watershed for stationary uniform storms. The phenomenon of attainment of equilibrium discharge from watersheds subjected to moving storms is contradictory to conventional hydrologic design, which assumes the storm duration must equal the time to equilibrium to attain the maximum discharge. Copyright © 2007 John Wiley & Sons, Ltd.     

SWAT application to estimate design runoff curve number for South Korean conditions

Runoff estimations based on the standard USDA–NRCS curve number (CN) table without calibration have a tendency to give inaccurate results when the CN values are applied in South Korea which has many high slope watersheds and that has a continental monsoon climate. Particularly for the design flood estimation, accurately calibrated CN values are required because the estimated peak flow is very sensitive to the selection of CN. However, the lack of flood data makes it difficult to calibrate and assign runoff CNs to Korean watersheds. Even if sufficient data are available to estimate CN values, it is also difficult to obtain the direct flows by separating base flows from total runoff hydrographs due to the temporal irregularity of rainfall events and the resulting complex pattern of runoff. Therefore, an alternative method for estimating CNs needs to be developed to overcome these issues. The purpose of this study is to present a method for estimating runoff CNs using the soil and water assessment tool (SWAT) model which can take into account watershed heterogeneities such as climate conditions, land use and soil types. The proposed CN estimation method uses the simulated flow data by SWAT instead of using measured flow data. This method has advantages in estimating CN values spatially for each subbasin division considering watershed characteristics. The use of daily data can reduce the sensitivity to the abnormality that is commonly involved in flow data with a small time scale. The SWAT‐based CN estimation method, combined with the asymptotic CN method, was applied to the Chungju dam watershed in South Korea. A regression equation was then developed from this approach, which was used to estimate CN values that decrease exponentially as rainfall amounts increase and that converge to 60·6 and 79·4 without and with considering subsurface lateral flow, respectively. Furthermore, the CN values for the antecedent moisture conditions were determined using the probabilistic approach. The CN associated with the 50% probability for the Chungju dam watershed is 87·8 which can be taken to be representative of antecedent moisture condition (AMC) II. The CN_I and CN_III associated with 90% and 10% probabilities are 78·9 and 94·1, respectively. The estimated CN_II = 87·8 differs markedly from the geographic information system (GIS)‐based CN 65·0, which implies that the standard USDA–NRCS CN method should be calibrated to the studied area of interest. Copyright © 2010 John Wiley & Sons, Ltd.     

Rainfall/runoff processes in a small peri‐urban catchment in the Andes mountains. The Rumihurcu Quebrada,Quito (Ecuador)

Situated at the foot of the Pichincha volcano, the city of Quito is frequently subjected to hydroclimatic hazards. In 1995 an 11·2 km² watershed, located in the vicinity of the city, was equipped with eight rain gauges and two flow gauges to better understand the local rainfall/runoff transformation processes. Rainfall simulation experiments were carried out on more than 40 one‐square‐metre plots to measure infiltration point‐processes. The high density of measurement devices allowed us to identify the origin and nature of the various contributions to runoff for the different physiographic units of the watershed: urban area from an altitude of 2800 to 3200 m; farmland, pasture and forested land, and finally páramo above 3900 m. Runoff occurs mainly in the lower part of the basin and is caused by urbanization; however, the natural soils of this area can also produce Hortonian runoff, which is predominant in a few events. This contribution can be studied through rainfall simulation experiments. In the upper natural zone, the younger and more permeable soils generate less runoff on the slopes. However, almost permanently saturated contributing areas, which are located in the bottom of the quebradas, may generate flood events, the size of which depends on the extent of the area concerned. Variations in the runoff coefficients are related first to the baseflow and second to the amount of rainfall in the previous 24 h. This analysis, which underlines the complexity of a small, peri‐urban, volcanic catchment, is a necessary preliminary to runoff modelling in an area where very few experiments have been carried out on small catchments. Copyright © 2001 John Wiley & Sons, Ltd.