Investigation of the Tank Model for Urban Storm Water Management

Abstract
The effect of impervious surfaces in urbanizing watersheds on stormwater runoff volumes and peak flood flows is investigated using a lumped hydrologic model, known as the Tank model, for the urbanized Mimico Creek watershed in Ontario. This model includes a separate algorithm for the simulation of stormwater runoff from impervious surfaces and an algorithm for the simulation of runoff and baseflow generated from pervious areas of the watershed. The simulated daily flow records of the calibrated model matched well with observed historic flow records. Since Mimico Creek is currently at 50% imperviousness; the model to simulate the sensitivity to change in streamflow for both 40% and 60% imperviousness scenarios is demonstrated. This study presents a simple novel management tool for assessing the effect of change in imperviousness of an urban watershed on streamflow volumes and peak flood flows.
1 Introduction
Urbanization is an ongoing trend in virtually all nations. High-density urbanization results in significant increase in impermeable surface areas, with resultant flooding problems. However, new urban stormwater management systems have been introduced over the last few decades, known as low impact development (LID) systems. The primary focus of LID systems is to promote infiltration and to decrease the volume of surface runoff and peak flood flows by changing the effective imperviousness.
Various hydrologic models, including the BC Water Balance Model; PC-SWMM, EPA National Stormwater Calculator, WinSLAMM, MIDS GUI Credit Calculator, and GIS-Calculator have been used to investigate the application of LID systems (Yates 1996; Elliott and Trowsdale 2007; Dietz 2007; Dietz and Clausen 2008; Finney et al. 2010; Basri 2013; Saini & Singh 2016). This study investigates the application of a Tank model to assess the effects of imperviousness on urban streamflows. Although the Tank model has been widely used in Asian countries, no evidence of its application in Canada was found.
2 Study Area
2.1 General Information
Mimico Creek begins in Brampton, then flows through part of Mississauga and past Toronto Pearson International Airport, and finally ends in Lake Ontario. The length of Mimico Creek is 33 km, and the area of the basin is 77 km2 (Nirupama et al. 2014). The shaded area in Figure 1 shows the area of the Mimico Creek watershed. Mimico Creek watershed is one of the most highly developed watersheds in the Greater Toronto Area (TRCA 2012).
Figure 1 Mimico Creek watershed area (Source: TRCA 2012).
The hydrometric station of Mimico Creek is located in Islington, 5 km away from its mouth, namely 02HC033.
2.2 Imperviousness Calculation
UDFCD (2001) recommends percentage imperviousness values for different types of land use; the selected values are shown in Table 1.
Table 1 Urban Drainage and Flood Control District recommended percentage imperviousness values (UDFCD, 2001).
Land use type | % Watershed area | % Impervious cover for this land use |
Commercial | 5.4 | 85% |
Industrial | 20.4 | 85% |
Institutional | 5.7 | 50% |
Open | 18.3 | 0% |
Roads | 11.9 | 100% |
Residential | 38.3 | 36% |
However, based on the breakdown of the land uses and the percentage imperviousness of each land use category, the overall imperviousness of the Mimico Creek watershed is about 50%.
2.3 Observed Streamflow Data
The daily flow data in Mimico Creek hydrometric station at Islington (02HC033) are available from 1965 to 2016, and water level records from 2002 to 2016. In this paper, three years of continuous data were used to calibrate the Tank model, for the period 1990–1992. The observed streamflow data of 2001 were used for Tank model validation.
2.4 Precipitation
Precipitation is the only input data for the Tank model. For this research, the daily precipitation data were obtained from Environment Canada. The selected meteorological station is Toronto Lester B. Pearson International Airport, which is near the centre of the study area and 9.18 km away from the hydrometric station, see Figure 2.
Figure 2 Drainage area for the Mimico Creek watershed.
The input to the Tank model is effective precipitation, not actual precipitation; however, the lack of evaporation data means that the effective precipitation cannot be calculated directly from the actual precipitation; hence the effective precipitation for the permeable Tank model were obtained using Equation 1, which is the actual precipitation multiplied by the weight of total observed discharge and total precipitation.
![]() |
(1) |
where:
EPi | = | effective precipitation (mm) on the ith day, |
Pi | = | total precipitation (mm) on the ith day, |
Obi | = | the ith day observed discharge (mm), and |
n | = | total number of days. |
Following the Urban Storm Drainage Criteria Manual (UDFCD 2001) recommendation, 2.54 mm is assumed as the depression loss for impervious land cover.
The average annual precipitation for the Mimico Creek watershed is 808 mm. The daily precipitation data for this meteorological station are available from 1953 to 2013, and for this research, the period 1990–1992 was selected to calculate the Tank model input for calibration and sensitivity analysis. The observed precipitation of 2001 was used for the model validation.
3 Tank Model
3.1 Structure
The Tank model structure depends on the type of land use. A watershed does not consist of a single type of land use, but consists typically of two or more types. Basri (2013) summarizes how the Tank model structure varies with the land uses, such as urbanized area with only one tank, while garden areas and forest area are with three tanks and four tanks in vertical series, respectively.
To analyse how the imperviousness affects a city, a combination tank model was used for the Mimico Creek watershed. The study watershed was divided into permeable and impermeable areas. The impermeable tank model, which is to simulate the runoff from the 50% impermeable urban land use, is built like a settlement tank with only one tank having one side outlet. For the permeable tank model, the two layer tank model (Figure 3) was considered as an acceptable option to simulate the 50% permeable urban land use. For the impermeable tank model, effective precipitation as the single input of this type of model, was calculated by subtracting 0.1 in. (2.54 mm) from the daily precipitation. The output from the side outlet is surface runoff. For the permeable tank model, precipitation is the single input into the top tank, followed by subtraction of the evaporation from the same tank. Both of these two processes can be calculated using Equation 1, known as the effective precipitation.
There are two side outlets in the top tank, and one in the bottom tank for the permeable model. The outputs from the side outlets are termed runoffs. The output from the bottom outlet of the top tank is the infiltration, which would occur with the onset of precipitation. The principle of the two types of Tank model is to simulate the runoff base flow on the water balance. Chow et al. (1988) summarized the water balance equation for a watershed as shown in Equation 2. The water balance equations and its components of the Tank model for investigation are shown in Equations 2 through 10. Figure 3 presents the schematic diagram of the structure of the Tank model. These two Tank models were built independently, but their simulation results were assigned different weights according to the imperviousness areas of the case study watershed.
Figure 3 Schematic of the Tank model.
The interpretation of the parameters of the combination Tank model is as shown in Table 2. A spreadsheet was applied to conduct these equations to simulate the runoff.
Table 2 Tank model parameters used for calibration.
Parameter | Unit | Description |
a1 | 1/d | Surface runoff coefficient of permeable Tank model |
a2 | 1/d | Intermediate runoff coefficient of permeable Tank model |
a3 | 1/d | Baseflow coefficient of permeable Tank model |
a4 | 1/d | Surface runoff coefficient of impermeable Tank model |
b1 | 1/d | Infiltration coefficient of permeable Tank model |
h1 | mm | Height of surface outlet of permeable Tank model |
h2 | mm | Height of intermediate outlet of permeable Tank model |
H1 | mm | Initial storage of top tank of permeable Tank model |
H2 | mm | Initial storage of bottom tank of permeable Tank model |
H3 | mm | Initial storage of impermeable Tank model |
Q1 | mm/d | Surface runoff of permeable Tank model |
Q2 | mm/d | Intermediate runoff of permeable Tank model |
Q3 | mm/d | Baseflow of permeable Tank model |
Q4 | mm/d | Surface runoff of impermeable Tank model |
Q | mm/d | Total simulated discharge of the modified Tank model |
Z1 | mm/d | Infiltration of permeable Tank model |
X1(t) | mm | Storage of top tank of permeable Tank model at t day |
X2(t) | mm | Storage of bottom tank of permeable Tank model at t day |
X3(t) | mm | Storage of impermeable Tank model at t day |
t | d | Time lag |
EP | mm/d | Effective precipitation |
The total simulation discharge of this type of Tank model consists of the sum of 50% of permeable tank model simulation discharge and 50% of permeable tank model simulation discharge.
Q1(t) = a1 × [X1(t) – h1] | (2) |
Q2(t) = a2 × [X1(t) – h2] | (3) |
Q3(t) = a3 × X2(t) | (4) |
Q4(t) = a4 × X3(t) | (5) |
Z(t) = b1 × X1(t) | (6) |
X1(t) = X1(t – 1) + EP − Q1(t) – Q2(t) − Z(t) | (7) |
X2(t) = X2(t – 1) + Z(t) – Q3(t) | (8) |
X3(t) = X3(t – 1) + EP – Q4(t) | (9) |
Q(t) = 50% × [Q1(t) + Q2(t) + Q3(t)] + 50% × Q4(t) | (10) |
At the initial time (t = 0), the initial conditions of water level storage in both the permeable tank model and the impermeable tank model are given as H1, H2 and H3.
3.2 Calibration
The parameters of the Tank model were calibrated using the Solver tool of Microsoft Excel by minimizing the root mean square error (RMSE) between observed and simulated daily flow records. Table 3 summarizes the results of the calibrated parameters.
Table 3 Calibrated parameters of the Tank model.
a1 | a2 | a3 | a4 | b1 | h1 | h2 | H1 | H2 | H3 |
0.07 | 0.01 | 0.001 | 1.62 | 0.90 | 0.04 | 0.02 | 0 | 256 | 1 |
Figure 4 shows the observed and simulated daily streamflows using the calibrated Tank model that shows a good match for a range of large and small storm events that occurred in the year 1991.
Figure 4 Comparison of simulated and observed hydrograph of the Tank model in 1991.
However, comparison of the simulation results for 1992 (Figure 5) indicates that unmatched peak events occur in the winter and spring periods. This means the Tank model may not perform well in the winter period in Canada. To improve the performance of the Tank model in the future, the snowmelt algorithm of the Tank model mentioned in Sugawara (1985) should be taken into consideration.
Figure 5 Comparison of simulated and observed hydrograph of the Tank model in 1992.
3.3 Tank Model Application
Volume Analysis
After the Tank model was calibrated, it was applied to analyse the relationship between the runoff and imperviousness of the urban area. Total runoff volume is one of the important hydrologic indicators, especially for design and retrofit of an urban drainage system with low impact development (LID) technology to increase permeable surface coverage. Table 4 indicates that total runoff volume will decrease from 347 mm to 303 mm if the percentage impervious cover is reduced from 50% to 40%.
Table 4 Runoff volume for various imperviousness percentages.
Runoff Volume (mm) | % ImperviousCover |
362 | 60% |
347 | 50% |
303 | 40% |
Analysis of Peak Discharge
Another important hydrologic indicator is peak discharge. The accurate prediction of peak discharge can effectively prevent property damage and protect the lives threatened. The values in Table 5 are the simulation results of the mean daily flood peaks varying with imperviousness during the period 1990–1992. The mean daily flood peaks will decrease significantly if thepercentage impervious cover is reduced from 50% today down to 40% through the implementation of permeable pavers and other LID systems to increase permeable surface coverage.
Table 5 Mean daily flood peaks variances with imperviousness from 1990 to 1992.
Water Year | Imperviousness | ||
60% | 50% | 40% | |
(m3/s) | (m3/s) | (m3/s) | |
1990 | 9.50 | 8.74 | 6.59 |
1991 | 12.31 | 11.38 | 8.58 |
1992 | 16.86 | 15.59 | 11.75 |
Average | 12.89 | 11.90 | 8.97 |
% change | +8% | Base Case | -25% |
4 Conclusion
A hydrologic model, known as the Tank model was used to investigate the effect of imperviousness on stormwater runoff volumes and peak flood flows. The Tank model includes algorithms for the simulation of stormwater runoff from impervious surfaces and both runoff and baseflow generated from pervious areas of the watershed. The model was successfully calibrated and tested, using historic records, for the case study watershed of Mimico Creek, Ontario for simulating daily flows.
Mimico Creek watershed is currently at 50% imperviousness; however, the calibrated Tank model was used to investigate two very different future imperviousness scenarios. First, the implementation of extensive LID systems (e.g. porous pavements, green roofs and infiltration galleries) to reduce imperviousness down to 40%, which would reduce mean annual streamflow volume by 12% and drop the peak flood flows by 25%. Second, continue urbanization under a business-as-usual scenario to 60% imperviousness, which will increase runoff volumes by 4% and peak flood flows by 8%.
In conclusion, the Tank model is a simple management tool for preliminary investigation of the potential effect of an increase or a decrease in the imperviousness of urbanizing watersheds on both runoff volumes and peak flood flows, in support of the implementation of low impact development technologies to protect the health of urban streams and their aquatic life.
Acknowledgment
The authors would like to acknowledge the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) for their support of this research.
References
- Basri, H. 2013. “Development of Rainfall–Runoff Model Using Tank Model: Problems and Challenges in Province of Aceh, Indonesia.” Aceh International Journal of Science and Technology 2 (1).
- Chow, V. T., D. R. Maidment and L. W. Mays 1988. Applied Hydrology. New York: McGraw-Hill.
- Dietz, M. E. 2007. “Low Impact Development Practices: A Review of Current Research and Recommendations for Future Directions.” Water, Air, & Soil Pollution 186 (1-4): 351–63.
- Dietz, M. E. and J. C. Clausen. 2008. “Stormwater Runoff and Export Changes with Development in a Traditional and Low Impact Subdivision.” Journal of Environmental Management 87 (4): 560–6.
- Elliott, A. H. and S. A. Trowsdale. 2007. “A Review of Models for Low Impact Urban Stormwater Drainage.” Environmental Modelling & Software 22 (3): 394–405.
- Finney, K., B. Gharabaghi, E. A. McBean, R. P. Rudra and G. MacMillan. 2010. “Compost Biofilters for Highway Stormwater Runoff Treatment.” Water Quality Research Journal of Canada 45 (4): 391–402.
- Nirupama, N., C. Armenakis and M. Montpetit. 2014. “Is Flooding in Toronto a Concern?” Natural Hazards 72 (2): 1259–64.
- Saini, S. and A. Singh. 2016. “Modeling LID using SWMM5 and MIDS Credit Calculator: Credit Valley Conservation’s Elm Drive Case Study.” Journal of Water Management Modeling 24: C403. doi: 10.14796/JWMM.C403.
- Sugawara, M. 1985. “Tank model.” Journal of Geography 94 (4): 209–21.
- TRCA (Toronto and Region Conservation Authority). 2012. “Watershed Features: Etobicoke Creek and Mimico Creek.” Toronto: TRCA. https://trca.ca/conservation/watershed-management/etobicoke-mimico-creek/watershed-features/.
- UDFCD (Urban Drainage and Flood Control District). 2001. In Urban Storm Drainage Criteria Manual, vol. 3, Stormwater Quality, 143-56. Denver, CO: Urban Drainage and Flood Control District. http://udfcd.org/criteria-manual.
- Yates, D. N. 1996. “WatBal: An Integrated Water Balance Model for Climate Impact Assessment of River Basin Runoff.” International Journal of Water Resources Development 12 (2): 121–40.