Unit hydrograph estimation with multiple events and prior information: I. Theory and a computer program

Abstract

The method of regularization for estimating unit hydrographs is expanded to allow the inclusion of prior information about the unit hydrograph shape. This may give smooth estimates without any loss in volume. The method is illustrated with prior information from a regression on catchment characteristics and with catchment lag determined from the data. A computer program to implement the method is given together with a sample calculation.     

Rainfall—runoff response,event-based runoff coefficients and hydrograph separation

Abstract

Event-based runoff coefficients can provide information on watershed response. They are useful for catchment comparison to understand how different landscapes “filter” rainfall into event-based runoff and to explain the observed differences with catchment characteristics and related runoff mechanisms. However, the big drawback of this important parameter is the lack of a standard hydrograph separation method preceding its calculation. Event-based runoff coefficients determined with four well-established separation methods, as well as a newly developed separation method, are compared and are shown to differ considerably. This signifies that runoff coefficients reported in the literature often convey less information than required to allow for catchment classification. The new separation technique (constant-k method) is based on the theory of linear storage. Its advantages are that it is theoretically based in determining the end point of an event and that it can also be applied to events with multiple peaks. Furthermore, it is shown that event-based runoff coefficients in combination with simple statistical models improve our understanding of rainfall—runoff response of catchments with sparse data.     

Estimation of suspended sediment loads in a seasonal stream in the wet‐dry tropics,Northern Territory,Australia

One year of instantaneous suspended sediment concentration, C, and instantaneous discharge, Q, data collected at Ngarradj downstream of the Jabiluka mine site indicate that the use of a simple C–Q rating curve is not a reliable method for estimating suspended sediment loads from the Ngarradj catchment. The C–Q data are not only complicated by hysteresis effects within the rising and falling stages of individual events, but also by variable depletion of available suspended sediment through multipeaked runoff events. Parameter values were fitted to an event‐based suspended sediment load–Q relationship as an alternative to the C–Q relationship. Total suspended sediment load and Q data for 10 observed events in the Ngarradj stream catchment were used to fit parameter values to a suspended sediment load–Q relationship, using (a) log–log regression and (b) iterative parameter fitting techniques. A more reliable and statistically significant prediction of suspended sediment load from the Ngarradj catchment is obtained using an event‐based suspended sediment load–Q relationship. Fitting parameters to the event‐based suspended sediment load–Q relationship using iterative techniques better predicts long‐term suspended sediment loads compared with log–log regression techniques. Copyright © 2004 John Wiley & Sons, Ltd.     

Estimation of temporal variation of sediment yield using GIS / Estimation de la variation temporelle de l'exportation sédimentaire grâce à un SIG

Abstract

A GIS-based method is proposed for computation of temporal variation of sediment yield during isolated storm events. Data from three Indian catchments, namely Karso and Nagwa in Jharkhand and Kharkari in Rajasthan, have been used. The Integrated Land and Water Information System (ILWIS) GIS package was used for (a) catchment discretization into cell areas using grid networks, (b) evaluation of the spatial variation in catchment topographical characteristics and land use, and (c) presentation of the results obtained. The process of sediment delivery from grid cells to the catchment outlet is represented by the topographical characteristics of the cells. Unit sediment graphs for the catchments are derived by translation of the sediment yield from the grid cells and routing through a linear storage reservoir. The proposed method is found to provide satisfactory estimates of the temporal variation of sediment yield during isolated storm events. The total sediment yield of a storm event may also be computed using the proposed method.     

Evaluation of satellite-based rainfall products for hydrological modelling in Morocco

ABSTRACT

Several satellite-based precipitation estimates are becoming available at a global scale, providing new possibilities for water resources modelling, particularly in data-sparse regions and developing countries. This work provides a first validation of five different satellite-based precipitation products (TRMM-3B42 v6 and v7, RFE 2.0, PERSIANN-CDR, CMORPH1.0 version 0.x) in the 1785 km² Makhazine catchment (Morocco). Precipitation products are first compared against ground observations. Ten raingauges and four different interpolation methods (inverse distance, nearest neighbour, ordinary kriging and residual kriging with altitude) were used to compute a set of interpolated precipitation reference fields. Second, a parsimonious conceptual hydrological model is considered, with a simulation approach based on the random generation of model parameters drawn from existing parameter set libraries, to compare the different precipitation inputs. The results indicate that (1) all four interpolation methods, except the nearest neighbour approach, give similar and valid precipitation estimates at the catchment scale; (2) among the different satellite-based precipitation estimates verified, the TRMM-3B42 v7 product is the closest to observed precipitation, and (3) despite poor performance at the daily time step when used in the hydrological model, TRMM-3B42 v7 estimates are found adequate to reproduce monthly dynamics of discharge in the catchment. The results provide valuable perspectives for water resources modelling of data-scarce catchments with satellite-based rainfall data in this region.

Editor M.C. Acreman; Associate editor N. Verhoest     

CURRENT RESEARCH IN BRITISH SNOWMELT RIVER FLOODING

ABSTRACT

An adaptation of the degree-day method has been used to analyse a number of snowmelt events on two catchments as a first step in a programme of research on snowmelt river flooding in Britain. The analysis indicates that the degree-day factor varies during events and between events on the same catchment. A snowmelt event is seen as consisting of three phases, an initial lag phase, a phase of nearly constant degree-day factor and a recession phase. The degree-day factor in the constant part of each event has significant correlation with the total flood volume on both catchments and with liquid precipitation during the snowmelt on one catchment only. Separate procedures are considered necessary for forecasting the initial lag phase and runoff during the recession.     

Improving the calibration strategy of the physically-based model WaSiM-ETH using critical events

Abstract

The use of a physically-based hydrological model for streamflow forecasting is limited by the complexity in the model structure and the data requirements for model calibration. The calibration of such models is a difficult task, and running a complex model for a single simulation can take up to several days, depending on the simulation period and model complexity. The information contained in a time series is not uniformly distributed. Therefore, if we can find the critical events that are important for identification of model parameters, we can facilitate the calibration process. The aim of this study is to test the applicability of the Identification of Critical Events (ICE) algorithm for physically-based models and to test whether ICE algorithm-based calibration depends on any optimization algorithm. The ICE algorithm, which uses the data depth function, was used herein to identify the critical events from a time series. Low depth in multivariate data is an unusual combination and this concept was used to identify the critical events on which the model was then calibrated. The concept is demonstrated by applying the physically-based hydrological model WaSiM-ETH on the Rems catchment, Germany. The model was calibrated on the whole available data, and on critical events selected by the ICE algorithm. In both calibration cases, three different optimization algorithms, shuffled complex evolution (SCE-UA), parameter estimation (PEST) and robust parameter estimation (ROPE), were used. It was found that, for all the optimization algorithms, calibration using only critical events gave very similar performance to that using the whole time series. Hence, the ICE algorithm-based calibration is suitable for physically-based models; it does not depend much on the kind of optimization algorithm. These findings may be useful for calibrating physically-based models on much fewer data.

Editor D. Koutsoyiannis; Associate editor A. Montanari

Citation Singh, S.K., Liang, J.Y., and Bárdossy, A., 2012. Improving calibration strategy of physically-based model WaSiM-ETH using critical events. Hydrological Sciences Journal, 57 (8), 1487–1505.     

Book reviews

Abstract

Statistical and deterministic modelling estimates of flood magnitudes and frequencies that can affect flood-plain ecology in the upper Ahuriri River catchment, a mountainous high country catchment in the New Zealand Southern Alps, were evaluated. Statistical analysis of 46 years of historical data showed that floods are best modelled by the generalized extreme value and lognormal distributions. We evaluated application of the HEC-HMS model to this environment by modelling flood events of various frequencies. Model results were validated and compared with the statistical estimates. The SCS curve number method was used for losses and runoff generation, and the model was very sensitive to curve number. The HEC-HMS flood estimates matched the statistical estimates reasonably well, and, over all return periods, were on average approximately 1% greater. However, the model generally underestimated flood peaks up to the 25-year event and overestimated magnitudes above this. The results compared well with other regional estimates, including studies based on L-moments, and showed that this catchment has smaller floods than other similarly-sized catchments in the Southern Alps.

Editor D. Koutsoyiannis; Associate editor H. Aksoy

Citation Caruso, B.S., Rademaker, M., Balme, A., and Cochrane, T.A., 2013. Flood modelling in a high country mountain catchment, New Zealand: comparing statistical and deterministic model estimates for ecological flows. Hydrological Sciences Journal, 58 (2), 328–341.     

Rainfall–runoff response and event-based runoff coefficients in a humid area (northwest Spain)

Abstract

The hydrological response of a small agroforestry catchment in northwest Spain (Corbeira catchment, 16 km²) is analysed, with particular focus on rainfall events. Fifty-four rainfall–runoff events, from December 2004 to September 2007, were used to analyse the principal hydrological patterns and show which factors best explain the hydrological response. The nonlinearity between rainfall and runoff showed that the variability in the hydrological response of the catchment was linked to the seasonal dynamics of the rainfall and, to a lesser extent, to evapotranspiration. The runoff coefficient, estimated as the ratio between direct runoff and rainfall volume, on an event basis, was analysed as a function of rainfall characteristics (amount and intensity) and the initial catchment state conditions prior to an event, such as pre-event baseflow and antecedent rainfall index. The results revealed that the hydrological response depends both on the soil humidity conditions at the start of the event and on rainfall amount, whereas rainfall intensity presented only a significant correlation with discharge increment. The antecedent conditions seem to be a key point in runoff production, and they explain much of the response. The hydrographs are characterized by a steep rising limb, a relatively narrow peak discharge and slow recession limb. These data and the observations suggest that the subsurface flow is the dominant runoff process.

Editor Z.W. Kundzewicz; Associate editor T. Wagener

Citation Rodríguez-Blanco, M.L., Taboada-Castro, M.M. and Taboada-Castro, M.T., 2012. Rainfall–runoff response and event-based runoff coefficients in a humid area (northwest Spain). Hydrological Sciences Journal, 57 (3), 445–459.     

Selecting best mapping strategies for storm runoff modeling in a mountainous semi‐arid area

Accurate runoff and soil erosion modeling is constrained by data availability, particularly for physically based models such as OpenLISEM that are data demanding, as the processes are calculated on a cell‐by‐cell basis. The first decision when using such models is to select mapping units that best reflect the spatial variability of the soil and hydraulic properties in the catchment. In environments with limited data, available maps are usually generic, with large units that may lump together the values of the soil properties, affecting the spatial patterns of the predictions and output values in the outlet. Conversely, the output results may be equally acceptable, following the principle of equifinality. To studyhow the mapping method selected affects the model outputs, four types of input maps with different degrees of complexity were created: average values allocated to general soil map units (ASG1), average values allocated to detailed map units (ASG2), values interpolated by ordinary kriging (OK) and interpolated by kriging with external drift (KED). The study area was Ribeira Seca, a 90 km² catchment located in Santiago Island, Cape Verde (West Africa), a semi‐arid country subject to scarce but extreme rainfall during the short tropical summer monsoon. To evaluate the influence of rainfall on runoff and erosion, two storm events with different intensity and duration were considered. OK and KED inputs produced similar results, with the latter being closer to the observed hydrographs. The highest soil losses were obtained with KED (43 ton ha^{? 1} for the strongest event). To improve the results of soil loss predictions, higher accurate spatial information on the processes is needed; however, spatial information of input soil properties alone is not enough in complex landscapes. The results demonstrate the importance of selecting the appropriate mapping strategy to obtain reliable runoff and erosion estimates. Copyright © 2013 John Wiley & Sons, Ltd.     

Sensitivity of hydrological models to uncertainty in rainfall input

Abstract

Different approaches used in hydrological modelling are compared in terms of the way each one takes the rainfall data into account. We examine the errors associated with accounting for rainfall variability, whether in hydrological modelling (distributed vs lumped models) or in computing catchment rainfall, as well as the impact of each approach on the representativeness of the parameters it uses. The database consists of 1859 rainfall events, distributed on 500 basins, located in the southeast of France with areas ranging from 6.2 to 2851 km². The study uses as reference the hydrographs computed by a distributed hydrological model from radar rainfall. This allows us to compare and to test the effects of various simplifications to the process when taking rainfall information (complete rain field vs sampled rainfall) and rainfall–runoff modelling (lumped vs distributed) into account. The results appear to show that, in general, the sampling effect can lead to errors in discharge at the outlet that are as great as, or even greater than, those one would get with a fully lumped approach. We found that small catchments are more sensitive to the uncertainties in catchment rainfall input generated by sampling rainfall data as seen through a raingauge network. Conversely, the larger catchments are more sensitive to uncertainties generated when the spatial variability of rainfall events is not taken into account. These uncertainties can be compensated for relatively easily by recalibrating the parameters of the hydrological model, although such recalibrations cause the parameter in question to completely lose physical meaning.

Citation Arnaud, P., Lavabre, J., Fouchier, C., Diss, S. & Javelle, P. (2011) Sensitivity of hydrological models to uncertainty of rainfall input. Hydrol. Sci. J. 56(3), 397–410.     

Hydrological processes in macrocatchment water harvesting in the arid region of Tunisia: the traditional system of tabias/Processus hydrologiques au sein d’un aménagement de collecte des eaux dans la région aride tunisienne: le système traditionnel des tabias

Abstract

Abstract In arid Tunisia, a tabia system is a traditional macrocatchment water harvesting system. It consists of a runoff area, which occupies two thirds of the slope and is traditionally used for grazing; and one to five cropped plots within U-shaped soil banks arranged in a cascade in the third downstream area. These ?run-on? areas accumulate and store the occasional runoff. Each soil bank is constructed with a discharge weir that allows modification of the flooded area and discharge of excess water towards downstream plots. Such a harvesting system, located in an area with 140 mm annual rainfall, was instrumented during four hydrological years (1995–1999) and 45 rainfall events were recorded. Eleven of these events gave a measurable inflow to at least one of the four plots. The observations showed that the traditional tabia system reduced total surface runoff from the catchment to essentially zero. The harvesting system significantly reduced peaks of surface runoff within the catchment, which also reduced erosion hazards. The cultivated area of about 5% of the total catchment could be supplied by a harvested water amount corresponding to about seven times the amount of each rainfall event larger than 20 mm.     

The use of maximum entropy to increase the informational content of hydrological networks by additional gauges

ABSTRACT

The optimization and extension of existing gauging networks are a challenging task, which can be done under consideration of many different aspects. One possibility is to maximize the obtained information on regional hydrological characteristics by new gauges compared to existing ones. For this, information theory approaches are most suitable. Here, the principle of maximum entropy is applied to calculate the probability of non-similarity of catchments to determine locations of new gauges according to the catchment characteristics that are most relevant for the hydrological conditions. The realization in an interactive application, provided online, makes use easy for practitioners and scientists. Goodness-of-fit measures are applied to investigate the explanatory power of the model and the contribution of each characteristic to the model, which gives information on the most influential properties of the catchment. The relevance of the proposed approach is proven by comparing hydrological signatures between similar and non-similar catchment.     

Integration of remote sensing data into hydrological models for reservoir management

Abstract

The purpose of this paper is to present the methodology set up to derive catchment soil moisture from Earth Observation (EO) data using microwave spaceborne Synthetic Aperture Radar (SAR) images from ERS satellites and to study the improvements brought about by an assimilation of this information into hydrological models. The methodology used to derive EO data is based on the appropriate selection of land cover types for which the radar signal is mainly sensitive to soil moisture variations. Then a hydrological model is chosen, which can take advantage of the new information brought by remote sensing. The assimilation of soil moisture deduced from EO data into hydrological models is based principally on model parameter updating. The main assumption of this method is that the better the model simulates the current hydrological system, the better the following forecast will be. Another methodology used is a sequential one based on Kalman filtering. These methods have been put forward for use in the European AIMWATER project on the Seine catchment upstream of Paris (France) where dams are operated to alleviate floods in the Paris area.     

Estimation of monthly runoff from small catchments in India

Abstract

Data from 31 non-snowfed catchments in India having catchment areas less than 1515 km² have been analysed to develop a simple method for the estimation of monthly runoff for the monsoon months of June to October. One of the parameters of this method was found to vary with the catchment area, the percentage of forest cover in the catchment and the monthly average temperature. The value of another parameter of the proposed method was found to be constant during any one month in a hydrologically homogeneous region. The method proposed herein is useful for estimating the monthly runoff during the monsoon period from catchments having scarce data.     

Changes in flood regime by use of the modified curve number method

Abstract

This paper describes the use of a simple two stage rainfall-runoff model in which a curve number (CN) principle is used to calculate the soil water content and, subsequently, the rainfall contribution to direct runoff and groundwater flow. The maximum soil water retention, S, is used to express various characteristics of a catchment (infiltration rate, soil cover and land use, as in the CN method) relevant to flood formation. Using historical flood events, the model is calibrated, and the statistical distribution parameters of peak flows determined. With the same historical input data scenarios (rainfall), sets of flood hydrographs are simulated for various values of the parameter S, and corresponding distribution parameters of peak flows are determined. This procedure is used to demonstrate possible changes in flood regime to be expected due to changes of the catchment soil properties and its vegetation cover. A case study is presented for the River Hron catchment, area 582 km², in the mountainous region of central Slovakia.     

Basin-scale analysis of rainfall and runoff in Peru (1969–2004): Pacific,Titicaca and Amazonas drainages

Abstract

We present a procedure for estimating Q95 low flows in both gauged and ungauged catchments where Q95 is the flow that is exceeded 95% of the time. For each step of the estimation procedure, a number of alternative methods was tested on the Austrian data set by leave-one-out cross-validation, and the method that performed best was used in the final procedure. To maximise the accuracy of the estimates, we combined relevant sources of information including long streamflow records, short streamflow records, and catchment characteristics, according to data availability. Rather than deriving a single low flow estimate for each catchment, we estimated lower and upper confidence limits to allow local information to be incorporated in a practical application of the procedure. The components of the procedure consist of temporal (climate) adjustments for short records; grouping catchments into eight seasonality-based regions; regional regressions of low flows with catchment characteristics; spatial adjustments for exploiting local streamflow data; and uncertainty assessment. The results are maps of lower and upper confidence limits of low flow discharges for 21 000 sub-catchments in Austria.     

Scale impacts on spatial variability in reference evapotranspiration

ABSTRACT

Evapotranspiration (ET) is one of the most important components in the hydrological cycle, and a key variable in hydrological modelling and water resources management. However, understanding the impacts of spatial variability in ET and the appropriate scale at which ET data should be incorporated into hydrological models, particularly at the regional scale, is often overlooked. This is in contrast to dealing with the spatial variability in rainfall data where existing guidance is widely available. This paper assesses the impacts of scale on the estimation of reference ET (ETo) by comparing data from individual weather stations against values derived from three national datasets, at varying resolutions. These include the UK Climate Impacts Programme 50 km climatology (UKCP₅₀), the UK Met Office 5 km climatology (UKMO₅) and the regional values published in the Agricultural Climate of England and Wales (ACEW). The national datasets were compared against the individual weather station data and the UKMO₅ was shown to provide the best estimate of ETo at a given site. The potential impacts on catchment modelling were then considered by mapping variance in ETo to show how geographical location and catchment size can have a major impact, with small lowland catchments having much higher variance than those with much larger areas or in the uplands. Some important implications for catchment hydrological modelling are highlighted.

Editor D. Koutsoyiannis; Associate editor L. Ruiz     

COMMON UNITGRAPHS FOR SETS OF RUNOFF EVENTS. PART 1: UNITGRAPH IDENTIFICATION FROM STREAMFLOW DATA

It is demonstrated that a unitgraph can be obtained without using rainfall data, provided data from at least two runoff events are available. A numerical method has been devised which calculates this common unitgraph for a set of surface runoff events and at the same time determines the input (rainfall excess) for each event. The method has been successfully tested on catchments ranging in size from 0.4 to 600 km²; it requires streamflow observations taken at intervals which retain all significant frequencies in the runoff hydrographs. The method also requires application of a baseflow separation procedure which is consistent for all events. The new approach has the potential to lead to more objective studies of the effects of catchment changes on the unitgraph and provides scope for comparisons of the common unitgraphs with geomorphological instantaneous unitgraphs.     

Infiltration in the Corgo River basin (northern Portugal): coupling water balances with rainfall—runoff regressions on a monthly basis

Abstract

A combination of water balances and rainfall—runoff regressions is used to calculate infiltration, overland flow, baseflow and change to the surface water reservoir, on a monthly basis; evapotranspiration from the underground reservoir, on an annual basis; and a lag phase of maximum infiltration and maximum baseflow within a hydrological year. The water balance equations are written for catchment areas formed on crystalline rocks and located in temperate climates. The regression lines are fitted to precipitations and river flows. In a first run, the model is tested with the Corgo River hydrographic basin, a small watershed in the Trás-os-Montes and Alto Douro province, northern Portugal. The results compare favourably with results of other groups, working under similar environmental conditions. The sensitivity of the model to changes in the basin characteristics and climate is tested by a second run using data from the Terva River basin, a nearby catchment that is much smaller than the Corgo basin and has a much lower effective precipitation, defined here as a difference between precipitation and potential evapotranspiration. As a consequence of having a lower effective precipitation, the river dry-out starts earlier in the Terva (May) than in the Corgo (June).