Effects of the catchment runoff coefficient on the performance of TOPMODEL in rainfall–runoff modelling

Effects of the catchment runoff coefficient on the performance of TOPMODEL in simulating catchment rainfall–runoff relationships are investigated in this paper, with an aim to improve TOPMODEL's simulation efficiency in catchments with a low runoff coefficient. Application of TOPMODEL in the semi‐arid Yihe catchment, with an area of 2623 km² in the Yellow River basin of China, produced a Nash–Sutcliffe model efficiency of about 80%. To investigate how the catchment runoff coefficient affects the performance of TOPMODEL, the whole observed discharge series of the Yihe catchment is multiplied with a larger‐than‐unity scale factor to obtain an amplified discharge series. Then TOPMODEL is used to simulate the amplified discharge series given the original rainfall and evaporation data. For a set of different scale factors, TOPMODEL efficiency is plotted against the corresponding catchment runoff coefficient and it is found that the efficiency of TOPMODEL increases with the increasing catchment runoff coefficient before reaching a peak (e.g. about 90%); after the peak, however, the efficiency of TOPMODEL decreases with the increasing catchment runoff coefficient. Based on this finding, an approach called the discharge amplification method is proposed to enhance the simulation efficiency of TOPMODEL in rainfall–runoff modelling in catchments with a low runoff coefficient. Copyright © 2004 John Wiley & Sons, Ltd.     

Development of a continuous soil moisture accounting procedure for curve number methodology and its behaviour with different evapotranspiration methods

The curve number (CN) method is widely used for rainfall–runoff modelling in continuous hydrologic simulation models. A sound continuous soil moisture accounting procedure is necessary for models using the CN method. For shallow soils and soils with low storage, the existing methods have limitations in their ability to reproduce the observed runoff. Therefore, a simple one‐parameter model based on the Soil Conservation Society CN procedure is developed for use in continuous hydrologic simulation. The sensitivity of the model parameter to runoff predictions was also analysed. In addition, the behaviour of the procedure developed and the existing continuous soil moisture accounting procedure used in hydrologic models, in combination with Penman–Monteith and Hargreaves evapotranspiration (ET) methods was also analysed. The new CN methodology, its behaviour and the sensitivity of the depletion coefficient (model parameter) were tested in four United States Geological Survey defined eight‐digit watersheds in different water resources regions of the USA using the SWAT model. In addition to easy parameterization for calibration, the one‐parameter model developed performed adequately in predicting runoff. When tested for shallow soils, the parameter is found to be very sensitive to surface runoff and subsurface flow and less sensitive to ET. Copyright © 2007 John Wiley & Sons, Ltd.     

A distributed monthly hydrological model for integrating spatial variations of basin topography and rainfall

Hydrological models at a monthly time‐scale are important tools for hydrological analysis, such as in impact assessment of climate change and regional water resources planning. Traditionally, monthly models adopt a conceptual, lumped‐parameter approach and cannot account for spatial variations of basin characteristics and climatic inputs. A large requirement for data often severely limits the utility of physically based, distributed‐parameter models. Based on the variable‐source‐area concept, we considered basin topography and rainfall to be two major factors whose spatial variations play a dominant role in runoff generation and developed a monthly model that is able to account for their influences in the spatial and temporal dynamics of water balance. As a hybrid of the Xinanjiang model and TOPMODEL, the new model is constructed by innovatively making use of the highly acclaimed simulation techniques in the two existing models. A major contribution of this model development study is to adopt the technique of implicit representation of soil moisture characteristics in the Xinanjiang model and use the TOPMODEL concept to integrate terrain variations into runoff simulation. Specifically, the TOPMODEL topographic index ln(a/tanβ) is converted into an index of relative difficulty in runoff generation (IRDG) and then the cumulative frequency distribution of IRDG is used to substitute the parabolic curve, which represents the spatial variation of soil storage capacity in the Xinanjiang model. Digital elevation model data play a key role in the modelling procedures on a geographical information system platform, including basin segmentation, estimation of rainfall for each sub‐basin and computation of terrain characteristics. Other monthly data for model calibration and validation are rainfall, pan evaporation and runoff. The new model has only three parameters to be estimated, i.e. watershed‐average field capacity WM, pan coefficient η and runoff generation coefficient α. Sensitivity analysis demonstrates that runoff is least sensitive to WM and, therefore, it can be determined by a prior estimation based on the climate and soil properties of the study basin. The other two parameters can be determined using optimization methods. Model testing was carried out in a number of nested sub‐basins of two watersheds (Yuanjiang River and Dongjiang River) in the humid region in central and southern China. Simulation results show that the model is capable of describing spatial and temporal variations of water balance components, including soil moisture content, evapotranspiration and runoff, over the watershed. With a minimal requirement for input data and parameterization, this terrain‐based distributed model is a valuable contribution to the ever‐advancing technology of hydrological modelling. Copyright © 2006 John Wiley & Sons, Ltd.     

Impact of imperfect rainfall knowledge on the efficiency and the parameters of watershed models

It is crucial to analyze the sensitivity of watershed (rainfall-runoff) models to imperfect knowledge of rainfall input, in order to judge whether or not they are reliable and robust, especially if they are to be used for operational purposes. In this paper, a new approach to sensitivity analysis is proposed, based on a comparison between the efficiency ratings and parameter values of the models and the quality of rainfall input estimate (GORE and BALANCE indexes, assessing the quality of rainfall time distribution and the total depth respectively). Data from three watersheds of increasing size (71, 1120, and 10700 km²), are used to test three watershed models of varying complexity (three-parameter GR3J model and six-parameter modified versions of TOPMODEL and IHACRES).

These models are able to cope with imperfect rainfall input estimates, and react to improvements in rainfall input accuracy by better performance and reduced variability of efficiency. Two different types of model behavior were identified: the models either benefit from improved rainfall data by producing more consistent parameter values, or they are unable to take advantage of the improvements. Although the watershed size seems to be immaterial, the smaller watersheds appear to need more precise areal rainfall estimates (a higher concentration of raingages) to ensure good modeling results.     

A stochastic framework to assess the performance of flood warning systems based on rainfall‐runoff modeling

Nowadays, Flood Forecasting and Warning Systems (FFWSs) are known as the most inexpensive and efficient non‐structural measures for flood damage mitigation in the world. Benefit to cost of the FFWSs has been reported to be several times of other flood mitigation measures. Beside these advantages, uncertainty in flood predictions is a subject that may affect FFWS's reliability and the benefits of these systems. Determining the reliability of advanced flood warning systems based on the rainfall–runoff models is a challenge in assessment of the FFWS performance which is the subject of this study. In this paper, a stochastic methodology is proposed to provide the uncertainty band of the rainfall–runoff model and to calculate the probability of acceptable forecasts. The proposed method is based on Monte Carlo simulation and multivariate analysis of the predicted time and discharge error data sets. For this purpose, after the calibration of the rainfall–runoff model, the probability distributions of input calibration parameters and uncertainty band of the model are estimated through the Bayesian inference. Then, data sets of the time and discharge errors are calculated using the Monte Carlo simulation, and the probability of acceptable model forecasts is calculated by multivariate analysis of data using copula functions. The proposed approach was applied for a small watershed in Iran as a case study. The results showed using rainfall–runoff modeling based on real‐time precipitation is not enough to attain high performance for FFWSs in small watersheds, and it seems using weather forecasts as the inputs of rainfall–runoff models is essential to increase lead times and the reliability of FFWSs in small watersheds. Copyright © 2013 John Wiley & Sons, Ltd.     

Using a topographic index to distribute variable source area runoff predicted with the SCS curve‐number equation

Because the traditional Soil Conservation Service curve‐number (SCS‐CN) approach continues to be used ubiquitously in water quality models, new application methods are needed that are consistent with variable source area (VSA) hydrological processes in the landscape. We developed and tested a distributed approach for applying the traditional SCS‐CN equation to watersheds where VSA hydrology is a dominant process. Predicting the location of source areas is important for watershed planning because restricting potentially polluting activities from runoff source areas is fundamental to controlling non‐point‐source pollution. The method presented here used the traditional SCS‐CN approach to predict runoff volume and spatial extent of saturated areas and a topographic index, like that used in TOPMODEL, to distribute runoff source areas through watersheds. The resulting distributed CN–VSA method was applied to two subwatersheds of the Delaware basin in the Catskill Mountains region of New York State and one watershed in south‐eastern Australia to produce runoff‐probability maps. Observed saturated area locations in the watersheds agreed with the distributed CN–VSA method. Results showed good agreement with those obtained from the previously validated soil moisture routing (SMR) model. When compared with the traditional SCS‐CN method, the distributed CN–VSA method predicted a similar total volume of runoff, but vastly different locations of runoff generation. Thus, the distributed CN–VSA approach provides a physically based method that is simple enough to be incorporated into water quality models, and other tools that currently use the traditional SCS–CN method, while still adhering to the principles of VSA hydrology. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Application of a data‐based mechanistic modelling (DBM) approach for predicting runoff generation in semi‐arid regions

This paper addresses the application of a data‐based mechanistic (DBM) modelling approach using transfer function models (TFMs) with non‐linear rainfall filtering to predict runoff generation from a semi‐arid catchment (795 km²) in Tanzania. With DBM modelling, time series of rainfall and streamflow were allowed to suggest an appropriate model structure compatible with the data available. The model structures were evaluated by looking at how well the model fitted the data, and how well the parameters of the model were estimated. The results indicated that a parallel model structure is appropriate with a proportion of the runoff being routed through a fast flow pathway and the remainder through a slow flow pathway. Finally, the study employed a Generalized Likelihood Uncertainty Estimation (GLUE) methodology to evaluate the parameter sensitivity and predictive uncertainty based on the feasible parameter ranges chosen from the initial analysis of recession curves and calibration of the TFM. Results showed that parameters that control the slow flow pathway are relatively more sensitive than those that control the fast flow pathway of the hydrograph. Within the GLUE framework, it was found that multiple acceptable parameter sets give a range of predictions. This was found to be an advantage, since it allows the possibility of assessing the uncertainty in predictions as conditioned on the calibration data and then using that uncertainty as part of the decision‐making process arising from any rainfall‐runoff modelling project. Copyright © 2001 John Wiley & Sons, Ltd.     

Refined conceptualization of TOPMODEL for shallow subsurface flows

The TOPMODEL framework was used to derive expressions that account for saturated and unsaturated flow through shallow soil on a hillslope. The resulting equations were the basis for a shallow‐soil TOPMODEL (STOPMODEL). The common TOPMODEL theory implicitly assumes a water table below the entire watershed and this does not conceptually apply to systems hydrologically controlled by shallow interflow of perched groundwater. STOPMODEL provides an approach for extending TOPMODEL's conceptualization to apply to shallow, interflow‐driven watersheds by using soil moisture deficit instead of water table depth as the state variable. Deriving STOPMODEL by using a hydraulic conductivity function that changes exponentially with soil moisture content results in equations that look very similar to those commonly associated with TOPMODEL. This alternative way of conceptualizing TOPMODEL makes the modelling approach available to researchers, planners, and engineers who work in areas where TOPMODEL was previously believed to be unsuited, such as the New York City Watershed in the Catskills region of New York State. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

A spatially distributed Clark’s unit hydrograph based hybrid hydrologic model (Distributed-Clark)

ABSTRACT

A hybrid hydrologic model (Distributed-Clark), which is a lumped conceptual and distributed feature model, was developed based on the combined concept of Clark’s unit hydrograph and its spatial decomposition methods, incorporating refined spatially variable flow dynamics to implement hydrological simulation for spatially distributed rainfall–runoff flow. In Distributed-Clark, the Soil Conservation Service (SCS) curve number method is utilized to estimate spatially distributed runoff depth and a set of separated unit hydrographs is used for runoff routing to obtain a direct runoff flow hydrograph. Case studies (four watersheds in the central part of the USA) using spatially distributed (Thiessen polygon-based) rainfall data of storm events were used to evaluate the model performance. Results demonstrate relatively good fit to observed streamflow, with a Nash-Sutcliffe efficiency (E_NS) of 0.84 and coefficient of determination (R²) of 0.86, as well as a better fit in comparison with outputs of spatially averaged rainfall data simulations for two models including HEC-HMS.     

Power function decay of hydraulic conductivity for a TOPMODEL‐based infiltration routine

TOPMODEL rainfall‐runoff hydrologic concepts are based on soil saturation processes, where soil controls on hydrograph recession have been represented by linear, exponential, and power function decay with soil depth. Although these decay formulations have been incorporated into baseflow decay and topographic index computations, only the linear and exponential forms have been incorporated into infiltration subroutines. This study develops a power function formulation of the Green and Ampt infiltration equation for the case where the power n = 1 and 2. This new function was created to represent field measurements in the New York City, USA, Ward Pound Ridge drinking water supply area, and provide support for similar sites reported by other researchers. Derivation of the power‐function‐based Green and Ampt model begins with the Green and Ampt formulation used by Beven in deriving an exponential decay model. Differences between the linear, exponential, and power function infiltration scenarios are sensitive to the relative difference between rainfall rates and hydraulic conductivity. Using a low‐frequency 30 min design storm with 4·8 cm h^?1 rain, the n = 2 power function formulation allows for a faster decay of infiltration and more rapid generation of runoff. Infiltration excess runoff is rare in most forested watersheds, and advantages of the power function infiltration routine may primarily include replication of field‐observed processes in urbanized areas and numerical consistency with power function decay of baseflow and topographic index distributions. Equation development is presented within a TOPMODEL‐based Ward Pound Ridge rainfall‐runoff simulation. Copyright © 2006 John Wiley & Sons, Ltd.     

Multi‐criterial validation of TOPMODEL in a mountainous catchment

The need for powerful validation methods for hydrological models including the evaluation of internal stages and spatially distributed simulations has often been emphasized. In this study a multi‐criterial validation scheme was used for validation of TOPMODEL, a conceptual semi‐distributed rainfall–runoff model. The objective was to test TOPMODEL's capability of adequately representing dominant hydrological processes by simple conceptual approaches. Validation methods differed in the type of data used, in their target and in mode. The model was applied in the humid and mountainous Brugga catchment (40 km²) in south‐west Germany. It was calibrated by a Monte Carlo method based on hourly runoff data. Additional information for validation was derived from a recession analysis, hydrograph separation with environmental tracers and from field surveys, including the mapping of saturated areas. Although runoff simulations were satisfying, inadequacies of the model structure compared with the real situation with regard to hydrological processes in the study area were found. These belong mainly to the concept of variable contributing areas for saturation excess overland flow and their dynamics, which were overestimated by the model. The simple TOPMODEL approach of two flow components was found to be insufficient. The multi‐criterial validation scheme enables not only to demonstrate limitations with regard to process representation, but also to specify where and why these limitations occur. It may serve as a valuable tool for the development of physically sound model modifications. Copyright © 1999 John Wiley & Sons, Ltd.     

Suitability of a lumped rainfall–runoff model for flashy tropical watersheds in New Caledonia

ABSTRACT

The GR4H lumped hourly rainfall–runoff model was assessed for its integration in a ridge-to-reef modelling framework. Particular attention was paid to rainfall representation, robustness of parameter estimates and ability to reproduce the main runoff features. The study was conducted in four tropical mountainous watersheds in New Caledonia, which are exposed to intense rainfall events, large annual climatic variations triggered by El Niño oscillation, and wildfires. The inverse distance and elevation weighting algorithm outperformed other classical rainfall interpolation methods under data-limited conditions. The time span of data needed for robust calibration was site specific and varied from 6–7 years to 10 years, which may be linked to El Niño events and to wildfires. With sufficient data, simulation quality was equivalent during the calibration and validation periods. The GR4H model was generally able to simulate both flash floods and large annual variations. The model was more reliable when simulating wet years and watersheds not subject to land-cover changes.     

HydRun: A MATLAB toolbox for rainfall–runoff analysis

Understanding the nature of streamflow response to precipitation inputs is at the core of hydrological applications and water resource management. Indices such as the base flow index, recession constant, and response lag of a watershed retain an important place in hydrology as metrics to compare watersheds and understand the impact of human activity, geology, geomorphology, soils, and climate on precipitation–runoff relations. Extracting characteristics of the hyetograph–hydrograph relationship is often done manually, which is time consuming and may result in subjective and potentially inconsistent outcomes. Here, we present a MATLAB‐based toolbox, called HydRun, for rapid and flexible rainfall–runoff analysis. HydRun uses a series of flexible routines to extract base flow from the hydrograph and then computes commonly used time instants of the rainfall–runoff relationship. HydRun provides users the flexibility to decide thresholds and limits of analysis, but objectively computes hydrometric indices. The toolkit includes a graphical user interface and example files. In this paper, we apply HydRun to 4 watersheds, 3 in Scotland and 1 in Canada, to demonstrate the software functions and highlight important decisions the user must make in its application.     

Multiple objective evaluation of a simple phosphorus transfer model

A simple phosphorus (P) transfer model of the Welland catchment, UK, is evaluated against multiple objective functions using a Monte Carlo approach that combines calibration, identifiability, sensitivity and uncertainty analysis. The model is based on simple conceptual rainfall‐runoff and river routing components, combined with estimates of the daily non‐point source load derived from annual landuse‐based export coefficients, disaggregated as a function of the runoff. The model has limited data requirements, consistent with data availability, and is parsimoneous with respect to the number of parameters identified through inverse modelling. The best performing parameter sets capture the main aspects of the observed flow and total P (TP) concentrations and provide a suitable basis for a decision‐support tool. However, a trade‐off is evident between matching the observed flow peaks, flow recessions and TP concentrations simultaneously, highlighting some limitations of the model structure and/or calibration data. Model analysis indicates that daily non‐point source load cannot be described as a function of near‐surface runoff and land use alone, but that other influences, including seasonality, are important. However, further model development to improve performance is likely to introduce additional complexity (in terms of parameter numbers), and hence additional problems of parameter identifiability and output uncertainty, which in turn raises issues of the information content of the available data. Copyright © 2004 John Wiley & Sons, Ltd.     

Hydrological model sensitivity to parameter and radar rainfall estimation uncertainty

Radar estimates of rainfall are being increasingly applied to flood forecasting applications. Errors are inherent both in the process of estimating rainfall from radar and in the modelling of the rainfall–runoff transformation. The study aims at building a framework for the assessment of uncertainty that is consistent with the limitations of the model and data available and that allows a direct quantitative comparison between model predictions obtained by using radar and raingauge rainfall inputs. The study uses radar data from a mountainous region in northern Italy where complex topography amplifies radar errors due to radar beam occlusion and variability of precipitation with height. These errors, together with other error sources, are adjusted by applying a radar rainfall estimation algorithm. Radar rainfall estimates, adjusted and not, are used as an input to TOPMODEL for flood simulation over the Posina catchment (116 km²). Hydrological model parameter uncertainty is explicitly accounted for by use of the GLUE (Generalized Likelihood Uncertainty Estimation). Statistics are proposed to evaluate both the wideness of the uncertainty limits and the percentage of observations which fall within the uncertainty bounds. Results show the critical importance of proper adjustment of radar estimates and the use of radar estimates as close to ground as possible. Uncertainties affecting runoff predictions from adjusted radar data are close to those obtained by using a dense raingauge network, at least for the lowest radar observations available. Copyright © 2004 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.     

Incorporating rainfall intensity into daily rainfall records for simulating runoff and infiltration into soil profiles

Hydrological modelling often implies the use of rainfall data. Its quality and resolution directly affect the accuracy of the simulation results. This study illustrates that a simple approach of incorporating rainfall intensity information in daily rainfall records significantly improves the simulation of surface runoff and rainfall infiltration into soil profiles. The procedure is developed using a frequency analysis on rainfall data of the Royal Meteorological Institute of Belgium, collected with a resolution of 10 min and for a consecutive period of 61 years. The frequency analysis of the data allowed the incorporation of rainfall intensity information into daily rainfall records. To test the effect of this approach the surface runoff and water flow into three different soil types was simulated using the HYDRUS‐1D model for a typical dry, normal and wet year. The simulation results whereby the observed 10‐min rainfall data was used as input were considered as the reference. Comparative analysis revealed that the simulations using the 10 min rainfall data deducted from the incorporation of rainfall intensity into daily rainfall records, deviate a maximum 1·2% from the reference and produce much better results than the Soil Conservation Service (SCS) runoff curve‐number method because rainfall intensity is considered in the procedure presented. The SCS curve‐number method typical overestimates surface runoff during periods of low rainfall intensity (winter) and underestimate runoff during periods of high rainfall intensities (summer). Copyright © 2002 John Wiley & Sons, Ltd.     

APPROXIMATING RAINFALL-RUNOFF MODELLING RESPONSE USING A STOCHASTIC INTEGRAL EQUATION

Rainfall-runoff modelling uncertainty can be analysed by the use of a stochastic integral formulation. The stochastic integral equation can be based on the rainfall–runoff model input of model rainfall or model rainfall excess. Similarly, the stochastic integral equation can be based on the rainfall–runoff model output of the modelled runoff hydrograph. The residual between actual measured runoff data and modelled runoff (from the rainfall–runoff model) is analysed here by the use of a stochastic integral equation. This approach is used to develop a set of convolution integral transfer function realizations that represent the chosen rainfall–runoff modelling error. The resulting stochastic integral component is a distribution of possible residual outcomes that may be directly added to the rainfall–runoff model's deterministic outcome, to develop a distribution of probable runoff hydrograph realizations from the chosen rainfall–runoff model.