One‐ and two‐dimensional modelling of overland flow in semiarid shrubland,Jornada basin,New Mexico

Two distributed parameter models, a one‐dimensional (1D) model and a two‐dimensional (2D) model, are developed to simulate overland flow in two small semiarid shrubland watersheds in the Jornada basin, southern New Mexico. The models are event‐based and represent each watershed by an array of 1‐m² cells, in which the cell size is approximately equal to the average area of the shrubs. Each model uses only six parameters, for which values are obtained from field surveys and rainfall simulation experiments. In the 1D model, flow volumes through a fixed network are computed by a simple finite‐difference solution to the 1D kinematic wave equation. In the 2D model, flow directions and volumes are computed by a second‐order predictor–corrector finite‐difference solution to the 2D kinematic wave equation, in which flow routing is implicit and may vary in response to flow conditions. The models are compared in terms of the runoff hydrograph and the spatial distribution of runoff. The simulation results suggest that both the 1D and the 2D models have much to offer as tools for the large‐scale study of overland flow. Because it is based on a fixed flow network, the 1D model is better suited to the study of runoff due to individual rainfall events, whereas the 2D model may, with further development, be used to study both runoff and erosion during multiple rainfall events in which the dynamic nature of the terrain becomes an important consideration. Copyright © 2006 John Wiley & Sons, Ltd.     

Comparison of catchment and network delineation approaches in complex suburban environments: application to the Chaudanne catchment,France

Suburban areas are subject to strong anthropogenic modifications, which can influence hydrological processes. Sewer systems, ditches, sewer overflow devices and retention basins are introduced and large surface areas are sealed off. The knowledge of accurate flow paths and watershed boundaries in these suburban areas is important for storm water management, hydrological modelling and hydrological data analysis. This study proposes a new method for the determination of the drainage network based on time efficient field investigations and integration of sewer system maps into the drainage network for small catchments of up to 10 km². A new method is also proposed for the delineation of subcatchments and thus the catchment area. The subcatchments are delineated using a combination of an object‐oriented approach in the urban zone and geographical information system–based terrain analysis with flow direction forcing in the rural zone. The method is applied to the Chaudanne catchment, which belongs to the Yzeron river network and is located in the suburban area of Lyon, France. The resulting subcatchment map gives information about subcatchment response and contribution. The method is compared with six other automatic catchment delineation methods based on stream burning, flow direction forcing and calculation of subcatchments for inlet points. None of the automatic methods could correctly represent the catchment area and flow paths observed in the field. The watershed area calculated with these methods differs by as much as 25% from the area computed with the new method. Copyright © 2012 John Wiley & Sons, Ltd.     

Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Two‐dimensional continuous simulation of spatiotemporally varied hydrological processes using the time‐area method

Distributed, continuous hydrologic models promote better understanding of hydrology and enable integrated hydrologic analyses by providing a more detailed picture of water transport processes across the varying landscape. However, such models are not widely used in routine modelling practices, due in part to the extensive data input requirements, computational demands, and complexity of routing algorithms. We developed a two‐dimensional continuous hydrologic model, HYSTAR, using a time‐area method within a grid‐based spatial data model with the goal of providing an alternative way to simulate spatiotemporally varied watershed‐scale hydrologic processes. The model calculates the direct runoff hydrograph by coupling a time‐area routing scheme with a dynamic rainfall excess sub‐model implemented here using a modified curve number method with an hourly time step, explicitly considering downstream ‘reinfiltration’ of routed surface runoff. Soil moisture content is determined at each time interval based on a water balance equation, and overland and channel runoff is routed on time‐area maps, representing spatial variation in hydraulic characteristics for each time interval in a storm event. Simulating runoff hydrographs does not depend on unit hydrograph theory or on solution of the Saint Venant equation, yet retains the simplicity of a unit hydrograph approach and the capability of explicitly simulating two‐dimensional flow routing. The model provided acceptable performance in predicting daily and monthly runoff for a 6‐year period for a watershed in Virginia (USA) using readily available geographic information about the watershed landscape. Spatial and temporal variability in simulated effective runoff depth and time area maps dynamically show the areas of the watershed contributing to the direct runoff hydrograph at the outlet over time, consistent with the variable source area overland flow generation mechanism. The model offers a way to simulate watershed processes and runoff hydrographs using the time‐area method, providing a simple, efficient, and sound framework that explicitly represents mechanisms of spatially and temporally varied hydrologic processes. Copyright © 2015 John Wiley & Sons, Ltd.     

Integration of transfer function model and back propagation neural network for forecasting storm sewer flow in Taipei metropolis

Typhoons and storms have often brought heavy rainfalls and induced floods that have frequently caused severe damage and loss of life in Taiwan. Our ability to predict sewer discharge and forecast floods in advance during storm seasons plays an important role in flood warning and flood hazard mitigation. In this paper, we develop an integrated model (TFMBPN) for forecasting sewer discharge that combines two traditional models: a transfer function model and a back propagation neural network. We evaluated the integrated model and the two traditional models by applying them to a sewer system of Taipei metropolis during three past typhoon events (NARI, SINLAKU, and NAKR). The performances of the models were evaluated by using predictions of a total of 6 h of sewer flow stages, and six different evaluation indices of the predictions. Finally, an overall performance index was determined to assess the overall performance of each model. Based on these evaluation indices, our analysis shows that TFMBNP yields accurate results that surpass the two traditional models. Thus, TFMBNP appears to be a promising tool for flood forecasting for the Taipei metropolis sewer system. For publication in Stochastic Environmental Research and Risk Analysis.     

Analysis and integrated modelling of groundwater infiltration to sewer networks

Infiltration of groundwater to sewer systems is a problem for the capacity of the system as well as for treatment processes at waste water treatment plants. This paper quantifies the infiltration of groundwater to a sewer system in Frederikshavn Municipality, Denmark, by measurements of sewer flow and novel model set‐up, which simulates the interaction between groundwater and sewer flow. The study area has a separate waste water sewer system, but the discharged volumes from the system are approximately twice the volumes from a tight system without infiltration. The model set‐up makes use of two commercial models: mike she for simulation of groundwater transport and mike urban (mouse ) [DHI, Hørsholm, Denmark] for simulation of sewer flow. By simulating the groundwater level and calibrating infiltration coefficients against sewer flow measurements, it has been possible to estimate the average infiltration to the sewer system with satisfying results. The infiltration processes are indeed complicated and to a large degree heterogeneous throughout the sewer system. The paper shows contribution from both saturated and unsaturated groundwater zones, which makes the modelling process complex. Copyright © 2016 John Wiley & Sons, Ltd.     

Cross-Media Models of the Chesapeake Bay Watershed and Airshed

A continuous, deterministic watershed model of the Chesapeake Bay watershed, linked to an atmospheric deposition model is used to examine nutrient loads to the Chesapeake Bay under different management scenarios. The Hydrologic Simulation Program - Fortran, Version 11 simulation code is used at an hourly time-step for ten years of simulation in the watershed. The Regional Acid Deposition Model simulates management options in reducing atmospheric deposition of nitrogen. Nutrient loads are summed over daily periods and used for loading a simulation of the Chesapeake estuary employing the Chesapeake Bay Estuary Model Package. Averaged over the ten-year simulation, loads are compared for scenarios under 1985 conditions, forecasted conditions in the year 2000, and estimated conditions under a limit of technology scenario. Limit of technology loads are a 50%, 64%, and 42% reduction from the 1985 loads in total nitrogen, total phosphorus, and total suspended solids, respectively. Urban loads, which include point source, on-site wastewater disposal systems, combined sewer overflows, and nonpoint source loads have the highest flux of nutrient loads to the Chesapeake, followed by crop land uses.on assignment from NOAA Air Resources Laboratory     

UNIT HYDROGRAPH DERIVED FROM A SPATIALLY DISTRIBUTED VELOCITY FIELD

A unit hydrograph model is proposed in which the watershed is decomposed into subareas which are individual cells or zones of neighbouring cells. The unit hydrograph is found for each subarea and the response at the outlet to excess rainfall on each subarea is summed to produce the watershed runoff hydrograph. The cell to cell flow path to the watershed outlet is determined from a digital elevation model. A constant flow velocity is assigned to each cell and the time lag between subarea input and response at the watershed outlet is found by integrating the flow time along the path from the subarea to the outlet. The response function for a subarea is modelled as a lagged linear reservoir in which the flow time is equal to the sum of a time of translation and an average residence time in the reservoir. It is shown that the assumption of a spatially varying, but time-invariant, velocity field underlying this model produces a linear system model for all subareas whose outputs can be summed in the manner indicated. An example application is presented for the 8.70 km² Severn watershed at Plynlimon in Wales using a 50 m digital elevation model in which the cell velocity is calculated by modifying an average watershed velocity according to the terrain slope and the drainage area of each cell. The resulting model reasonably reproduces the observed unit hydrograph.     

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.     

Comparing quantitative precipitation forecast methods for prediction of sewer flows in a small urban area

Abstract

Due to the relatively small spatial scale, as well as rapid response, of urban drainage systems, the use of quantitative rainfall forecasts for providing quantitative flow and depth predictions is a challenging task. Such predictions are important when consideration is given to urban pluvial flooding and receiving water quality, and it is worthwhile to investigate the potential for improved forecasting. In this study, three quantitative precipitation forecast methods of increasing complexity were compared and used to create quantitative forecasts of sewer flows 0–3 h ahead in the centre of a small town in the north of England. The HyRaTrac radar nowcast model was employed, as well as two different versions of the more complex STEPS model. The STEPS model was used as a deterministic nowcasting system, and was also blended with the Numerical Weather Prediction (NWP) model MM5 to investigate the potential of increasing forecast lead-times (LTs) using high-resolution NWP. Predictive LTs between 15 and 90 min gave acceptable results, but were a function of the event type. It was concluded that higher resolution rainfall estimation as well as nowcasts are needed for prediction of both local pluvial flooding and combined sewer overflow spill events.

Editor D. Koutsoyiannis; Guest editor R.J. Moore     

Application of distributed hydrological models for predictions in ungauged basins: a method to quantify predictive uncertainty

Stream flow predictions in ungauged basins are one of the most challenging tasks in surface water hydrology because of nonavailability of data and system heterogeneity. This study proposes a method to quantify stream flow predictive uncertainty of distributed hydrologic models for ungauged basins. The method is based on the concepts of deriving probability distribution of model's sensitive parameters by using measured data from a gauged basin and transferring the distribution to hydrologically similar ungauged basins for stream flow predictions. A Monte Carlo simulation of the hydrologic model using sampled parameter sets with assumed probability distribution is conducted. The posterior probability distributions of the sensitive parameters are then computed using a Bayesian approach. In addition, preselected threshold values of likelihood measure of simulations are employed for sizing the parameter range, which helps reduce the predictive uncertainty. The proposed method is illustrated through two case studies using two hydrologically independent sub‐basins in the Cedar Creek watershed located in Texas, USA, using the Soil and Water Assessment Tool (SWAT) model. The probability distribution of the SWAT parameters is derived from the data from one of the sub‐basins and is applied for simulation in the other sub‐basin considered as pseudo‐ungauged. In order to assess the robustness of the method, the numerical exercise is repeated by reversing the gauged and pseudo‐ungauged basins. The results are subsequently compared with the measured stream flow from the sub‐basins. It is observed that the measured stream flow in the pseudo‐ungauged basin lies well within the estimated confidence band of predicted stream flow. Copyright © 2013 John Wiley & Sons, Ltd.     

Application of a watershed runoff model to North-east Pond River,Newfoundland: to study water balance and hydrological characteristics owing to atmospheric change

The hydrological sensitivities to long-term climate change of a watershed in Eastern Canada were analysed using a deterministic watershed runoff model developed to simulate watershed acidification. This model was modified to study atmospheric change effects in the watershed. Water balance modelling techniques, modified for assessing climate effects, were developed and tested for a watershed using atmospheric change scenarios from both state of the art general circulation models and a series of hypothetical scenarios. The model computed daily surface, inter- and groundwater flows from the watershed. The moisture, infiltration and recharge rate are also computed in the soil reservoirs. The thirty years of simulated data can be used to evaluate the effects of climatic change on soil moisture, recharge rate and surface and subsurface flow systems. The interaction between surface and subsurface water is discussed in relation to climate change. These hydrological results raise the possibility of major environmental and socioeconomic difficulties and have significant implications for future water resource planning and management. © 1997 John Wiley & Sons, Ltd.     

Development of a polder module in the SWAT model: SWATpld for simulating polder areas in south‐eastern China

Polders are one of the most common artificial hydrological entities in the plain river network regions of China. Due to enclosed dikes, manual drainage, and irrigation intake operations, polders have had a significant impact on the hydrological processes of these areas. Distributed hydrological models are effective tools to understand and reproduce the hydrological processes of a watershed. To date, however, few models are able to simulate the drainage and irrigation intake interactions of polders at a watershed scale. This study develops a modified version of the Soil and Water Assessment Tool (SWAT) model, which is designed to better represent polders (SWATpld). The SWATpld model simulates drainage and irrigation intake processes by calculating the excess‐water storage in the inner rivers and irrigation schedule for paddy rice in the polder. Both SWAT and SWATpld models were tested for the Liyang watershed. SWATpld outperformed SWAT in simulating the daily discharge and intake of the experimental polder and predicting the monthly peak flow at the outlet of the Liyang watershed, which suggests that the modified model simulates the hydrological responses of the study watershed with polder operations more realistically than the original SWAT model does. Further evaluation at various locations and in various climate conditions would increase the confidence of this model.     

Inference model derivation with a pattern analysis for predicting the risk of microbial pollution in a sewer system

Developing a mathematical model for predicting fecal coliform bacteria concentration is very important because it can provide a basis for water quality management decisions that can minimize microbial pollution risk to the public. This paper introduces a hybrid modeling methodology which is a combined use of a neural network-based pattern analysis and an evolutionary process model induction system. The neural network-based pattern analysis technique is applied to extract knowledge on inter-relationships between fecal coliform concentrations and other measurable variables in a sewer system. Based on the result of neural network-based pattern analysis, an evolutionary process model induction system is used to derive mathematical inference models that can predict fecal coliform bacteria concentration from easily measurable variables instead of directly measuring fecal coliform bacteria concentration in a sewer system. The neural network-based pattern analysis extracts that temperature and ammonia concentration are the most important driving forces leading to an increase in fecal coliform bacteria concentration in the sewer system at Paraparaumu City, New Zealand. Fecal coliform bacteria concentration is also positively correlated with dissolved phosphorus and inversely with flow rate. The multivariate inference models that are able to predict fecal coliform bacteria concentration are successfully derived as functions of flow rate, temperature, ammonia, and dissolved phosphorus in the form of understandable mathematical formulae using the evolutionary process model induction system, even if a priori mathematical knowledge of the dynamic nature of fecal coliform bacteria is poor. The multivariate inference models evolved by the evolutionary process model induction system produce a slightly better performance than the multi-layer perceptron neural network model.     

A kinematic‐wave‐based distributed watershed model using FEM,GIS and remotely sensed data

For the appropriate management of water resources in a watershed, it is essential to calculate the time distribution of runoff for the given rainfall event. In this paper, a kinematic‐wave‐based distributed watershed model using finite element method (FEM), geographical information systems (GIS) and remote‐sensing‐based approach is presented for the runoff simulation of small watersheds. The kinematic wave equations are solved using FEM for overland and channel flow to generate runoff at the outlet of the watershed concerned. The interception loss is calculated by an empirical model based on leaf area index (LAI). The Green‐Ampt Mein Larson (GAML) model is used for the estimation of infiltration. Remotely sensed data has been used to extract land use (LU)/land cover (LC). GIS have been used to prepare finite element grid and input files such as Manning's roughness and slope. The developed overland flow model has been checked with an analytical solution for a hypothetical watershed. The model has been applied to a gauged watershed and an ungauged watershed. From the results, it is seen that the model is able to simulate the hydrographs reasonably well. A sensitivity analysis of the model is carried out with the calibrated infiltration parameters, overland flow Manning's roughness, channel flow Manning's roughness, time step and grid size. The present model is useful in predicting the hydrograph in small, ungauged watersheds. Copyright © 2007 John Wiley & Sons, Ltd.     

Effectiveness of BIOECODS for peak flow attenuation: an appraisal using InfoWorks SD

The Bio-Ecological Drainage System (BIOECODS) is a sustainable drainage system, which adopts a “control at source” approach for urban storm water management in Malaysia. This study attempts to model a small-scale BIOECODS using InfoWorks SD. New modelling techniques are used to fully integrate the surface and on-line subsurface conveyance system, in which overland flow routing is described by a storm water management model that uses a nonlinear reservoir method and the kinematic wave approximation of the St Venant equation, and subsurface flow is described by the Horton method in conjunction with the Soil Conservation Service (SCS) curve number (CN) method. The observed water levels at primary outlets are compared with those obtained from model simulation. The modelling approach has been proven successful as the hydrographs (predicted and observed) match each other closely, with a mean error in the range of 4.58–7.32%. Results from the model showed that the BIOECODS is able to attenuate peak flow by 60–75%, and increase the lag time by 20 min within an area of <28?300 m² when compared with a traditional drainage system.     

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.     

Watershed modeling of rainfall excess transformation into runoff

In this paper an attempt is made to present a distributed physiographic conceptual model that uses the principles of flow continuity and momentum. For this purpose, the watershed under study is divided into subwatersheds keeping in view the drainage patterns and characteristics. Then the main tributaries are identified and their drainage areas are delineated to form tributary subwatersheds. The main channel subwatersheds have taken care of the remaining area in the vicinity of the main channel. The kinematic wave theory is applied for the overland runoff computations from these subwatersheds. Further, the overland flows are superimposed onto the main channel. The dynamic wave theory is used to route the flows through the main channel to compute the watershed responses at the outlet. The proposed model is tested onto a natural watershed. The computations were performed for few storm events. Comparison of the significant parameters of the computed and the observed hydrographs shows that the maximum relative error in prediction is 5.8%. Thus, the results are satisfactory. Better results can be obtained when measured rainfall-excess data are available or a more realistic loss index is adopted for rainfall-excess separation.     

湖北牛山湖小型鱼类的群落结构及多样性

流域数据模型是流域特征的语义、行为和规则的表达,是进行集成流域模拟和管理的空间数据组织的重要内容．本文以长江三角洲太湖流域上游的西苕溪流域为例,将汇流单元分为自然流域、“大包围”、圩区三种形式,分别建立了各自内部的河湖网络关系,即山区由自然流域单元,内包含水库(或湖泊)、河流、水工点的树状河湖网络关系;平原区由人工汇流单元,内包含大包围、圩区、湖泊、湿地、河道、水工点的网状河湖网络关系,为建立适合我国的流域数据模型做了有益的探索．