The Laurenson Runoff Routing procedure can be applied to rural and urban catchments in both design storm and continuous simulation situations. It allows for the non-linear response from catchments over a large range of event magnitudes.
The Laurenson procedure is also known as the RAFTS model, and is originally from a paper Aitken, A.P. (1975) Hydrologic investigation and design of urban stormwater drainage systems. Aust. Water Resources Council Tech. Paper No. 10, Dept of the Environment and Conservation, A.G.P.S., Canberra.
However the equation is repeated in: ARR Volume 1 on page 187.
See also the reference Sobinoff, P. Pola, J.P. and O’Loughlin, G.C. (1983) Runoff routing parameters for the Newcastle-Sydney-Wollongong region. Hydrol. And Water Resources Symposium 1983, inst. Engrs Aust., Natl Conf. publ. No. 83/13, pp. 28-32.
The equation is only suitable for an n value of 0.025.
Further modifications were made by Willing and Partners following calibration runs on project catchments between 1976 and 1985.
The hydrological data requirements are: catchment area, slope, degree of urbanisation, loss rates, observed or design rainfall.
These data are used to compute the storage delay coefficient for each of the sub-catchments and hence to develop the non-linear runoff hydrograph. A default exponent is adopted, although the user may override this value with either a different non-linear exponent or a rating table of flow vs an exponent to define different degrees of catchment non-linear response.
Each sub-catchment is divided into 10 sub-areas. Each of the sub-areas is treated as a cascading non-linear storage obeying the relationship:
S = volume of storage (hrs x m3/s)
B = storage delay time coefficient
n = storage non-linearity exponent by default is set to -0.285 (AR&R, 1997)
The rainfall is applied to each sub-area, an excess computed and the excess converted into an instantaneous inflow. This instantaneous flow is then routed through the sub-area storages to develop an individual sub-catchment outlet hydrograph.
Initial Rainfall Loss
Specifies the depth of rainfall that infiltrates before any runoff occurs.
Continuous Rainfall Loss
Specifies the loss that occurs after the initial loss has been satisfied.
B can be entered or calculated from:
B = 0.285 A^0.52 (1 + U)^-1.97Sc^-0.50
A = catchment area (km2)
U = fraction of catchment urbanised
Sc = modified equal area slope (%)
Storage non-linearity exponent by default is set to -0.285
Can be used as a multiplication factor used during calibration to modify the calculated B. It will default to value of 1.
PERN is an adaption factor used to introduce the effect of pervious catchment roughness as follows:
Fraction of catchment urbanised.
Modified equal area slope (%).
Laurenson E.M. (1964) A Catchment Storage Model for Runoff Routing, Journal of Hydrology, Volume 2, pp 141-163
Institution of Engineers, Australia. (1987) Australian Rainfall and Runoff - A Guide to Flood Estimation (edited by D.H. Pilgrim and R.P. Canterford), Volumes 1 and 2, Revised Edition, November, 1987, Canberra