SIMHYD with routing is a conceptual rainfall-runoff model that estimates daily stream flow from daily rainfall and areal potential evapotranspiration data. SIMHYD with routing is similar to the standard SIMHYD model but has an additional store to simulate routing of surface flows.

The model contains stores for interception loss, soil moisture, groundwater and routing. The model has nine parameters.

Scale

SIMHYD with routing operates at a functional unit scale and daily time-step.

Principal developer

Cooperative Research Centre for Catchment Hydrology. SIMHYD is a simplified version of the daily conceptual rainfall-runoff model, HYDROLOG, that was developed in 1972 (see Porter 1972; and Porter & McMahon 1975) and the more recent MODHYDROLOG (Chiew & McMahon 1991).

The SIMHYD with routing model has nine parameters as compared to the 17 parameters required for HYDROLOG and the 19 for MODHYDROLOG.

The SIMHYD model was developed by Francis Chiew in Fortran and converted to C# in TIME by Jean-Michel Perraud. Routing was added by Jai Vaze and Phillip Jordan.

Scientific provenance

SIMHYD has been widely applied to a large number of Australian catchments by several hydrologists (Peel et al., 2000; Chiew et al., 2008). The extent of its use outside of Australia is unknown but the conceptual structure is not particularly limited to Australian catchments and with appropriate calibration and testing it is likely that it could be successfully applied in other countries.

Version

Source v4.3

Dependencies

None

Availability/conditions

SIMHYD is automatically installed with Source. SIMHYD is also available through the Rainfall Runoff Library on eWater Toolkit http://www.toolkit.net.au/Tools/RRL.

Flow phase

The structure of the simple lumped conceptual daily rainfall-runoff model, SIMHYD, is shown in Figure 1.

Figure 1. Structure of the SIMHYD with routing rainfall-runoff model

In SIMHYD, daily rainfall first fills the interception store, which is emptied each day by evaporation. The excess rainfall is then subjected to an infiltration function that determines the infiltration capacity. The excess rainfall that exceeds the infiltration capacity becomes infiltration excess runoff.

Moisture that infiltrates is subjected to a soil moisture function that diverts the water to the stream (interflow), groundwater store (recharge) and soil moisture store. Interflow is first estimated as a linear function of the soil wetness (soil moisture level divided by soil moisture capacity). The equation used to simulate interflow therefore attempts to mimic both the interflow and saturation excess runoff processes (with the soil wetness used to reflect parts of the catchment that are saturated from which saturation excess runoff can occur). Groundwater recharge is then estimated, also as a linear function of the soil wetness. The remaining moisture flows into the soil moisture store.

Evapotranspiration from the soil moisture store is estimated as a linear function of the soil wetness, but cannot exceed the atmospherically-controlled rate of areal potential evapotranspiration. The soil moisture store has a finite capacity and overflows into the groundwater store. Base flow from the groundwater store is simulated as a linear recession from the store.

The model therefore estimates runoff generation from three sources - infiltration excess runoff, interflow (and saturation excess runoff) and base flow. Surface flows are routed through a storage.

The fundamental equations of the model are shown in Equation 1 through Equation 10.

Equation 1

Equation 2

Equation 3

Equation 4

Equation 5

Equation 6

Equation 7

Equation 8

Equation 9

Equation 10

Input data

The model requires daily rainfall and potential evapotranspiration data. The rainfall and evaporation data sets need to be continuous (no gaps) and overlapping. Catchment area in km² is required to provide flow output volumes.

Daily rainfall data may be obtained from rain gauges or rainfall surfaces but will need to be converted to a time series record that is spatially representative of the whole catchment. Note that the time that rainfall data is collected may be important. Very often rainfall data is collected in the morning, the usual time is 9 am, and may be more representative of the previous day’s rainfall.

Daily evaporation is an estimate of the spatially averaged evaporation rate of the catchment being modelled. This estimate is subject to the types of land uses that are in the catchment. This may be estimated by applying a crop/land use factor to daily pan or potential evapotranspiration surface data.

Daily flow data in ML/day, m³/s or mm/day may be required to calibrate the model.

Selecting stream flow data to use in a river-basin-scale simulation study needs information about the reliability of the data. It is best to use data which are most representative of the stream flow from the catchment. Observed data would normally be selected, except where the data are of poor quality or of unknown reliability.

Parameters or settings

Model parameters are summarised in Table 1.

Table 1. Model Parameters

Parameter	Description	Unitys	Default	Min	Max
Baseflow coeff.	Base flow Coefficient		n.a.	0.0	1.0
Impervious Threshold	Impervious Threshold	mm	n.a.	0.0	5.0
Infiltration Coeff.	Infiltration Coefficient		n.a.	0.0	400
Infiltration shape	Infiltration Shape		n.a.	0.0	10.0
Interflow Coeff.	Interflow Coefficient		n.a.	0.0	1.0
Perv. Fraction	Pervious Fraction		n.a.	0.0	1.0
RISC	Rainfall Interception Store Capacity	mm	n.a.	0.0	5.0
Recharge coefficient	Recharge Coefficient		n.a.	0.0	1.0
SMSC	Soil Moisture Store Capacity	mm	n.a.	1.0	500
K	Surface delay parameter	hours	n.a.	0.5	200
x	inflow to outflow bias parameter		n.a.	0.0	0.5

The relative sensitivity of parameters will vary between catchments but generally the model is most sensitive to the soil moisture store capacity, pervious fraction and base flow index

Output data

The model outputs daily surface and base flow. This may be saved in ML/day, m³/s or mm/day.

Configuration

This model requires calibration and validation.

Reference list

Porter, JW 1972, The synthesis of continuous streamflow Department of Civil Engineering, Monash University, Melbourne, p. 222.

Porter, JW & McMahon, TA 1975, ‘Application of a catchment model in southeastern Australia’, Journal of Hydrology, vol. 24, pp. 121-134.

Chiew, FHS & McMahon, TA 1991, ‘Improved modelling of the groundwater processes in HYDROLOG’, Proceedings of the 20th Hydrology and Water Resources Symposium, Perth, October 1991, Institute of Engineers Australia, pp. 492-497.

Bibliography

Rainfall Runoff Library v1.0.5, June 25, 2004 (http://www.toolkit.net.au/Tools/RRL).

Chiew, F. H. S., M. J. Stewardson and T. A. Mc Mahon (1993). "Comparison of six rainfall-runoff modelling approaches." Journal of Hydrology 147: 1-36.

Chiew, F. H. S. and T. A. McMahon (1994). "Application of the daily-runoff model MODHYDROLOG to 28 Australian catchments." Journal of Hydrology 153: 383-416.

Chiew, FHS, Peel, MC & Western, AW 2002, ‘Application and testing of the simple rainfall-runoff model SIMHYD’, in VP Singh & DK Frevert (eds), Mathematical models of small watershed hydrology and Applications, Water Resources Publication, Colorado, pp.335-367.

Chiew FHS, Vaze J, Viney NR, Jordan PW, Perraud J-M, Zhang L, Teng J, Young WJ, Penaarancibia J, Morden RA, Freebairn A, Austin J, Hill PI, Wiesenfeld CR and Murphy R (2008) Rainfall-runoff modelling across the Murray-Darling Basin. A report to the Australian Government from the CSIRO Murray-Darling Basin Sustainable Yields Project. CSIRO, Australia. 62pp

Peel, M. C., F. H. S. Chiew, A. W. Western and T. A. McMahon (2000). Extension of Unimpaired Monthly Streamflow Data and Regionalisation of Parameter Values to Estimate Streamflow in Ungauged Catchments. Canberra, National Land and Water Resources Audit.

Vaze, J., J.-M. Perraud, N. R. Viney, D. A. Post, J. Teng, B. Wang, J. Lerat and M. Goswami (2010). "Rainfall-runoff modelling across southeast Australia: Datasets, models and results." Australian Journal of Water Resources 14(2): 101-116.

Browser not supported