Printer Friendly

Modelling stormwater runoff, quality, and pollutant loads in a large urban catchment.

Abstract. Identification of stormwater runoff, its pollution load, and their implications for land use is essential in implementing stormwater management strategies. Hydrologic modelling provides an opportunity to assess them at limited data resources. In this study, the stormwater management model SWMM5 is applied for model development for a large basin in Tallinn. A geographic information system tool is used for subcatchment delineation, identification of directly connected impervious areas (DCIAs), and preparation of catchment input parameters. The model is calibrated and verified using sampled storm events to estimate event mean concentrations and annual loads. The predictive capability of the model for quantity is good and for quality moderate. The findings from the model show the percentage of the impervious area in the large catchment to be low at 19.7%. Although DCIAs, in particular roads and roofs, have relatively smaller areas they significantly impact runoff production (up to 75%) and loads (up to 66% total phosphorus and 71% total suspended solids). The first flush at the beginning of runoff is less important in case of a low intensity of rainfall, but heavy rain and snowmelt generate substantial runoff and pollution loads. When grab sampling is applied, it should focus on the medium and large events within 6 hours of storm commencement in order to achieve better mass estimations.

Key words: hydrologic modelling, impervious area, event mean concentrations, mass loads, first flush.

Abbreviations and symbols
ADD--antecedent dry days            SA1--time weighted sampling
oxygen demand                       SA2--random grab sampling
CC--correlation                     SA3--grab sampling within 6 h
coefficient                         irrespective of storm size
Cr--commercial roofs                SA4--grab sampling within
                                    6 h of medium and large storms
connected impervious area           SWMM--stormwater management model
DEE--Department of
Environmental Engineering
of Tallinn                          TIA--total impervious area
University of Technology            TN--total nitrogen
DNCIA--directly not
connected impervious area           TP--total phosphorus
EIA--effective impervious area      TSS--total suspended solids
ELLE--Estonian, Latvian, and
Lithuanian Environment group        A--catchment area
EMC--event mean concentration       depression storage
GIS--geographic information system  depression storage
NOF--normalized objective function  surface roughness
NSC--Nash-Sutcliffe coefficient     surface roughness
                                    % imp--percentage
R--residential area                 of impervious area
Rd--roads                           [S.sub.o]--average slope
RMSD--root mean square deviation    Sc--sensitivity coefficient
RE--relative error                  W--catchment width
Rr--residential roofs


Stormwater runoff has induced flooding and degraded the receiving water bodies in many urbanized areas. Runoff reduction and the control of pollution loads at coastal areas has been a goal for the Baltic Sea states for a number of past decades (HELCOM, 2002, 2007). This is reflected in national strategies, acts, and regulations such as, in Estonia, the Estonian Water Act (Riigikogu, 2016), Tallinn stormwater strategy until 2030 (Tallinna Linnavolikogu, 2012), and Tallinn City development plan (Tallinna Linnavolikogu, 2013), among others. The estimation of actual stormwater quality and quantity has been a requirement in evaluating the compliance of stormwater management regulations and in implementing effective control measures, and it is an inherently difficult task when there is a lack of sufficient information (Chiew and McMahon, 1999; Park and Stenstrom, 2008). Although effective, extensive monitoring campaigns are not always feasible due to resource availability and associated uncertainties. In such a situation, stormwater modelling is a helpful tool, which uses limited data resources and can simulate time intervals beyond the monitoring period (Vezzaro and Mikkelsen, 2012).

The constituents of stormwater runoff mostly depend on the rainfall and catchment characteristics (Nazahiyah et al., 2007; Liu et al., 2012, 2013) as well as on the upstream human activities (Brezonik and Stadelmann, 2002; Park et al., 2009). Models are primarily underpinned by rainfall (rainfall intensity, rainfall duration, antecedent dry condition) and catchment characteristics such as land use type and surface imperviousness (Park et al., 2009; Liu et al., 2013). However, it is also necessary to consider urban forms such as the road layout and spatial distribution of urban areas, which have significant influences on pollutant generation (Liu et al., 2012, 2013). Moreover, mixed land use, typical in large catchment basins, produces a higher variability in stormwater quality, as there is a high interspersion of various land use types (Lee et al., 2009). Therefore, a large catchment behaves in a different way than a small urban catchment in runoff and pollutant loading.

One of the important catchment characteristics is surface imperviousness. A large proportion of impervious surface accelerates the runoff rate and produces shorter times of concentration or lag times as well as an increased peak flow and runoff volume; however, the proximity of impervious surfaces to a drainage system plays a significant role in these interactions (Shuster et al., 2005). Estimation based on total impervious area overestimates the quantity and quality of stormwater (Shuster et al., 2005; Ebrahimian et al., 2016) and a larger area of the catchment is more crucial to add up the amount. Effective impervious area (EIA) is used by many researches instead of total impervious area to effect closer results. For EIA directly connected impervious areas (DCIAs) need to be accounted because the part that is not directly connected to storm drainage has less impact on the output (Lee and Heaney, 2003; Ebrahimian et al., 2016). Identification of the potential DCIAs is essential in implementing the possible solution for reducing runoff and pollution load.

The reviews of the stormwater models by Elliott and Trowsdale (2007) and Jayasooriya and Ng (2014) found that the US Environmental Protection Agency's Storm Water Management Model (EPA SWMM) has an effective performance in simulating stormwater quality and quantity. It is a widely used water quality model that has the capacity for both single-event and continuous simulation in the prediction of flows and pollutant concentrations. Many studies have indicated that SWMM has reasonable accuracy when the model outcomes are calibrated and validated (Jayasooriya and Ng, 2014). There is broad applicability of SWMM for simulating runoff and quality dynamics: hydrology assessment for pre- and post-development condition (Jang et al., 2007), pollutant wash-off water quality analyses (Tsihrintzis and Hamid, 1998; Temprano et al., 2006; Lee et al., 2010), combined sewer overflow modelling and assessment (Zhang and Li, 2015), flood forecasting (Han et al., 2014), and stormwater treatment facilities modelling and assessment (Burszta-Adamiak and Mrowiec, 2013; Zhang et al., 2013).

Many researchers applied the model for assessing runoff and loads with calibration and verification (Tsihrintzis and Hamid, 1998; Temprano et al., 2006; Nazahiyah et al., 2007; Tan et al., 2008; Lee et al., 2010; Chow et al., 2012; Mancipe-Munoz et al., 2014; Rosa et al., 2015). These studies are mostly for small catchment basins. Tan et al. (2008) and Mancipe-Munoz et al. (2014) worked on comparatively large urban catchments but mainly focused on runoff calibration, although they worked also for continuous events. There are few studies on quality dynamics for large catchment basins. In Estonia, stormwater modelling can rarely be found. Hood et al. (2007) used SWMM to estimate the flow and pollution load of a moderate size subcatchment in Tallinn, but they did not evaluate its predictive capabilities adequately. Six small subcatchments with distinctive land use (transportation, residential, and commercial) in Tallinn were analysed for runoff and suspended solids by Koppel et al. (2014). The runoff dynamics of mixed land use can be different from single land use and can provide a resultant runoff coefficient that is different from the coefficient of single land use (Lee et al., 2009). The studies by Hood et al. (2007) and Koppel et al. (2014) used rainfall data from Harku station, which is almost 5-7 km away from the studied sites, but no stations nearer the sites were available at that time. However, an analysis using distant stations can lead to unconvincing results. There is also room to look at the impact of directly connected impervious areas (DCIAs). Lee and Heaney (2003) modelled the hydrologic performance of DCIAs and reported that DCIAs are the main areas contributing to runoff and have a most pronounced effect on urban hydrology. The DCIA or the connectivity to the urban area at the catchment scale influences the hydrologic response (Yang et al., 2011; Burns et al., 2015).

This study is mainly focused on model development for the Mustoja basin in Tallinn. It aims to identify the sensitive input parameters and potential DCIAs that influence the runoff and loads in a large catchment basin. It will estimate event mean concentrations (EMCs) and mass loads and finally use these results to evaluate the practiced sampling campaigns.


2.1. Description of the studied catchment area

The study was conducted in the Mustoja basin, a large catchment basin in Tallinn City (Fig. 1). Tallinn is the capital of Estonia where approximately 32% of the population is centred. It is in the north-eastern part of Europe on the shores of the Gulf of Finland. It has a humid continental climate with moderately warm summers and cold snowy winters. Normally, average monthly temperatures range from -5.7 [degrees]C in February to 16.3 [degrees]C in July (EWS, 2015). Total annual rainfall in Tallinn is 704 mm and average monthly rainfall varies from 32 mm in April to 86 mm in August. Snow cover usually lasts from mid-December to late March.

The study catchment comprises approximately 10.24 k[m.sup.2] (30% of the Tallinn area), which for the most part spreads over the Kristiine and Mustamae districts in the upstream side. Runoff mainly flows through underground pipe networks and ditches to the downstream natural channel where three pipes under Marja, Haabersti, and Mustjoe streets intersect. Finally this water is discharged to Kopli Bay and the Baltic Sea. The land use within the catchment is mainly residential with detached houses and apartment buildings in the upstream side. Industrial and commercial areas are dominant features in the downstream side within the catchment.


In this study, EPA SWMM version 5.1 (SWMM5) was used to simulate stormwater quantity and quality. It is a comprehensive hydrological and water quality model used for single event or long-term events in urban areas (Rossman, 2010). The model is common throughout the world for planning, analysis, and design related to stormwater drainage systems in urban areas. The US EPA provides this model and related graphical user interface free of charge. The SWMM conceptual model comprises four major environmental compartments: (i) the atmosphere compartment, e.g. rain gauge objects, accounts for precipitation and pollutants from the air; (ii) the land surface compartment, which is represented through subcatchment objects, models areas that receive precipitation and generate runoff; (iii) the groundwater compartment, e.g. aquifer objects, receives infiltration from the land surface and provides input to the transport compartment; and (iv) the transport compartment, whose routes flow from runoff source areas through a network of pipes, channels, etc. (Rossman, 2010). SWMM offers a selection of three different built-in infiltration models and flow routing methods. The infiltration loss on the pervious area was estimated using Horton's equation (Horton, 1933) and the runoff transport was computed using the dynamic wave routing method under the complete Saint-Venant equations that consider the backwater effects with 30 s time steps.

2.2. Model development

Model development can be divided into two steps. In the first step, runoff modelling is constructed, which requires information on catchment characteristics, conveyance system, rainfall, and infiltration. In the second step, the quality model is developed using the runoff model, which requires build-up and wash-off components.

Catchment characteristics are catchment area (A), catchment width (W), average slope ([S.sub.o]), percentage of the impervious area (% imp), surface depression storage, and surface roughness. Surface depression storage includes impervious ([D.sub.imp]) and pervious depression storage ([D.sub.per]) while surface roughness includes impervious ([N.sub.imp]) and pervious surface roughness ([N.sub.per]).

Subcatchments were delineated using the drainage networks and surface slope. The drainage network was provided by AS Tallinna Vesi in digital format, which includes the location/elevation of stormwater pipes and manholes, pipe diameter, pipe material, and year of installation (Table 1).

The Mustoja catchment basin is mainly served with a separate storm drainage system of approximately 51 km of pipes with the diameter varying from 0.15 to 2 m and 4 km of drain ditches constructed within 79 years. The surface slope was determined using a 1-m resolution digital elevation map (DEM) provided by the Estonian Land Board. A geographic information system (GIS) toolbox was used to prepare the catchment properties. The Mustoja basin was divided into 378 sub-catchments (see in Fig. 1) ranging in size from 0.06 to 23.8 ha with the slope varying in the range 0-11.1%. The width of each catchment was computed applying the longest flow path method within the catchment. The width ranged from 25.8 to 3002.8 m.

For estimating the percentage of the impervious area, the land use map for the year 2014 was obtained from the Estonian Land Board. According to the Estonian National Topographic Database, the Mustoja basin had impervious surfaces of almost 50%, most of which are roads, fields of production sites, and buildings (Table 2). Pervious surfaces were mainly formed by the green areas, which made up a quarter of the catchment area, and by forest and private yards with 1/5 of the catchment area. The total impervious area (TIA) does not contribute to the actual urban runoff. Instead, EIA or the portion of TIA that is hydraulically connected to the storm sewer system is an important parameter in determining this runoff (Lee and Heaney, 2003; Ebrahimian et al., 2016). The DCIAs were identified and separated from the unconnected parts following the procedure explained by Lee and Heaney (2003). Field investigation, areal maps, online maps, and storm pipe details were used, and imperviousness was spatially analysed using GIS. Numerous land uses were grouped together to form a simple classification including residential, commercial, industrial, forest, water bodies, roads, and roofs (Table 3). Roads and roofs are potential urban land uses for runoff and pollutants (Ballo et al., 2009). Approximately 55.7% of the areas were determined as residential (R), 7% as industrial (I), 10.2% as roads (Rd), 23% as residential roofs (Rr), 3.9% as commercial roofs (Cr), and 0.2% as forest and water (Table 3). Of the total catchment area DCIAs were found to form nearly 27% (278 ha) of which TIA made up 79% (217 ha) and EIA 56% (156 ha). We estimated EIA because all the runoff does not enter the inlets despite the fact that DCIA only represents the land use connected to the drainage system. Nevertheless, when the directly not connected impervious area (DNCIA) that is not connected to the drainage system was considered, TIA was 44% (448.9 ha) and EIA was 19% (198.4 ha). The major part of EIA was DCIA roads, followed by DCIA commercial roof, industrial areas, and residential areas.

2.2.1. Rainfall details

Stormwater runoff estimation in the Tallinn region has usually been based on the Tallinn-Harku station, about 6 km from the pilot basin, where the meteorological observations started in 1805. The rainfall varies in time, space, and altitude among stations even within a catchment. Rainfall measuring equipment must be close to or within the catchment in order to take account of the climatic factors, characterization and modelling of the stormwater drainage system in a statistically reliable way (Lowe et al., 2014). Since 2013, AS Tallinna Vesi has installed nine other rainfall measuring stations. Within the Mustoja basin, in the downstream area near Saarma Street, the Department of Environmental Engineering of Tallinn University of Technology installed a tipping bucket rain gauge 'Saarma' in May 2014. The influential stations were identified as Tondi 90 and Saarma (see Fig. 1) applying the Thiessen polygon method (Fetter, 2001) using a GIS toolbox. The Tondi 90 station was within the basin in the upstream area near Tondi Street. One minute rainfall data from these two stations were used for the simulation so that the stations close to the subcatchments could feed them the corresponding data. The average monthly temperatures (in Harku station) and rainfalls recorded in three stations (Tondi 90, Saarma, and Harku) during Jan. 2014 to Feb. 2016 are presented in Fig. 2. August 2014 and July 2015 were the wettest months while February to April were dry months in 2014 and 2015, and October 2015 was the driest month. Compared with the normal rainfall pattern in the Harku station over 1981-2010 (EWS, 2015), these two years were dry years with approximately 34% and 26% less rain in 2014 and 2015, respectively.

The storm events recorded from May 2014 to Dec. 2015 in the Tondi 90 and Saarma stations were similar in rainfall characteristics (Table 4). Over 2015, 144 to 156 events occurred with average rainfall intensity 0.1 to 4.8 mm/h, producing 523 to 550 mm of rainfall. When the intensities were averaged over the years, the intensity amounted to 0.5-0.7 mm/h. Most of the storms (90%) during these two years were below 1.4 mm/h. An extreme event during the record period reached the peak of 19.4 mm/h on 22 Aug. 2014. The average total event rains were small in range, 1.7-2.2 mm for a duration of 6.1-8.1 h and inter-event time of 52.8-59.8 h.


2.2.2. Storm events and sampling details

Three storm events were used for model calibration and four events for validation. The calibrated events were during 9 September 2014 (event 1), 22-23 September 2014 (event 2), and 4-5 December 2015 (event 3). The rainfall characteristics (intensity, duration, antecedent dry days (ADD), and total rain) during these events are presented in Table 5. Event 1 and event 2 were captured after a long dry period of 12.1 and 12.6 days; the former had intense rain with 2.2 mm/h of rain of short duration (2.3 h) and the latter had medium rain with 0.8 mm/h of a long duration (26 h). Event 3 was recorded during wet weather conditions with 1.5 ADD and was below 1.4 mm/h; this falls between the medium and 90 percentile across all recorded years. Event 1 did not have quality data; therefore, it was used only for runoff calibration.

The automatic flow measurement unit installed by AS Tallinna Vesi provided flow rate for each 10 min for event 1 and event 2. At the same time, Tallinn University of Technology installed a water level measurement gauge to form a discharge-rating curve. The time interval sampling approach called sampling approach 1 (SA1) was used to collect samples during the events. Grab samples were also taken two times a week and a total of 104 samples were taken during the period from 6/11/2014 to 16/12/2015, which is called random sampling approach or sampling approach 2 (SA2). On each occasion, the flow was measured along with instantaneous water level using a flow tracker. Samples were analysed by competent water quality laboratories certified by the Estonian Accreditation Centre (EAK), which consistently followed the standard EVS-EN ISO/IEC 17025:2006 (EVS, 2015). Many parameters, e.g. temperature, dissolved oxygen, conductivity, pH, biological oxygen demand (BOD), total suspended solids (TSS), ammonia, total nitrogen (TN), total phosphorus (TP), chloride, oils, microbiological parameters, and heavy metals, were measured. Heavy metal samples were taken every week, and a total of 30 samples were determined for analysing the content of Cd, Cr, Cu, Pb, Ni, Fe, and Zn. In this study, TSS and TP parameters with a strong linear relationship (regression coefficient 0.72-0.86) were simulated and the Water Quality Laboratory at Tallinn University of Technology was involved in analysing these parameters.

Other input parameters for catchment properties were adopted from the range provided in SWMM User's manual (Huber et al., 1988), books (Bedient and Huber, 1988; Wanielista, 1990), and published papers (Temprano et al., 2006; Nazahiyah et al., 2007), e.g. 0.3 mm for [D.sub.imp], 2.5 mm for [D.sub.per], 25% for zero depression storage, 0.01 for surface roughness coefficient for overland flow on impervious portion, 0.3 for surface roughness coefficient for overland flow on pervious portion, etc. (Table 6). The infiltration parameters (maximum infiltration rate, minimum infiltration rate, and decay constant) were often difficult to identify and were adopted from the range provided by Bedient and Huber (1988) and Nazahiyah et al. (2007) based on the soil type used. Pollutant buildup depends on the land use, dirt and dust accumulation, and ADD, whereas pollutant wash-off depends on land use, runoff rate, and removal efficiency. The build-up and wash-off parameters provided by Chow et al. (2012) for different land uses were similar to the values used by Koppel et al. (2014) in the local context, and those values were adopted to the associated land use types.

2.2.3. Sensitivity and model performance test

The robustness of the results from the model was analysed for sensitivity by varying the values of the input parameters. For example, % imp and width were changed within [+ or -]10% and other parameters within the range as presented in Table 6. Sensitivity coefficient (Sc), which measures the effect of change in one factor on another (James and Burges, 1982; Chow et al., 2012), was the indicator for this analysis. The model goodness of fit was tested with four types of indicators as used by Chow et al. (2012) and Koppel et al. (2014): correlation coefficient, relative error (RE), normalized objective function (NOF), and Nash--Sutcliffe coefficient (NSC). Correlation coefficient (CC) is a general statistical measure of the linear relationship between observed and predicted values. Generally, the best calibration requires CC to be as close to 1 as possible. The RE for flow and peak flow measures the average tendency of the simulated data to be larger or smaller than observed (Chow et al., 2012). The optimal value of RE is 0.0, but within [+ or -]10% can be acceptable with low magnitude values indicating accurate model simulation. Let O denote the observed value, S the simulated value, [O.sub.i] the ith observed and [S.sub.i] the ith simulated value, and n the total number of observations. Then RE is calculated as follows:


NOF, which is the root mean square deviation (RMSD) normalized with the mean of the observed values, is expressed as follows:


The optimal value of NOF is 0 but 0 to 1 can be acceptable for calibration (Kornecki et al., 1999). NSC is calculated by Eq. (3) (Nash and Sutcliffe, 1970).


NSC ranges between 0 and 1.0, with NSC = 1 being the optimal value. Values between 0.0 and 1.0 are generally viewed as acceptable levels of performance.


3.1. Results of sensitivity analysis

In the studied basin, the sensitive input parameters were found to be % imp, W, [D.sub.imp], and [N.sub.imp] (Table 7). Among them, % imp has a significant influence on runoff flow rate at Sc = 0.82 and on peak flow at Sc = 0.68. Catchment width has a moderate influence with Sc of 0.44 for flow rate and 0.36 for peak flow. Both [D.sub.imp] and [N.sub.imp] have a negative coefficient, which indicates that the output values will increase with a decrease in these input parameters. In this large basin, they are not highly significant being in the range from -0.019 to -0.029 for [D.sub.imp] and from -0.008 to -0.018 for [N.sub.imp]. However, [D.sub.imp] is comparatively more influential for initial peak flow, having Sc of -0.029, which indicates that the depression storages regulate the first peak in the hydrograph. Tan et al. (2008) and Chow et al. (2012) found a similar influence for % imp. However, the width of their subcatchments was more sensitive to peak flow than runoff depth, which was not the case in this study. Instead, the influence there was nearly equal on peak flow and flow rate. Unlike in their study, in our case the basin was large and the sensitivity to peak flow was not strong due to width and [N.sub.imp].

3.2. Calibration and verification results

The model is evaluated for the predictive capabilities through four indicators using Eqs (1)-(3). The results in Table 8 and Table 9 show that the CC, NOF, and NSC are within the acceptable range for runoff quantity for all events, indicating the modelled and measured runoff to be in a good relationship. The model for runoff simulation is acceptable at three important indicators: CC is 0.87-1 i.e. close to 1; NOF is 0.1-0.3, i.e. between 0 to 1; and NSC is 0.3-1, i.e. close to 1. However, RE is beyond [+ or -]10% at flow error less than 20% and peak error is between -27.4% and +21.2%. Similarly, the indicators such as CC and NOF for quality simulation (Table 9) suggest that the predictive capability of the model is in an acceptable range being 0.4-1 for CC and 0.1-0.5 for NOF. Here NSC and RE are not in a good range; NSC is poor for TSS but it is good for TP at 0.4-0.7, while RE can go beyond the acceptable range for both quality parameters. Overall, the model is acceptable for runoff quantity, but it provides less accurate estimation for quality performance. Calibration through increasing the number of events can reduce the error to some extent.

The calibrated results for runoff quantity are presented in Table 6 and for quality in Table 10. These are the best-fit values for the input parameters after testing the goodness of fit. The output fits at 0.9 factor of % imp (Table 6); as a result the ultimate value is 19.7%, where the runoff coefficients for different land use vary as follows: for R from 0.05 to 0.2, Rr from 0.15 to 0.20, RC from 0.2 to 0.26, I from 0.09 to 0.45, Cr from 0 to 0.56, and Rd from 0.03 to 0.57, depending on the proximity of the connection to the drainage system. Temprano et al. (2006) in residential areas of Spain and Hood et al. (2007) in the highly residential Tallinn area found lower % imp values of 15.9% and 7.2%, respectively. A relatively large pervious area within the sub-catchments and the effect of mixed land use have resulted in lower imperviousness. The depression storage used in this model is the average of dry and wet weather, which is determined as 0.7 mm for the basin. It is higher than the values obtained by Hood et al. (2007) and Koppel et al. (2014), where they proposed 0.15 to 0.29 mm for the Tallinn area. Depression storage for impervious area is sensitive to the initial peak flow, and it varied from 0.3 to 1 in this study depending on the weather conditions. Chow et al. (2012) provided this storage for residential area 0.2 in wet weather and 0.8 in dry weather, while it can increase to 0.75 in wet and 1.05 in dry conditions for commercial area. Temprano et al. (2006) suggested for residential area a higher value of 2.5 mm. Therefore, the value obtained for this basin can be justified. Manning's roughness for impervious area is 0.0135, since the impervious surface characteristics within this area are concrete/asphalt paving and/or a gravelled surface (Huber et al., 1988). A similar value is also used by Koppel et al. (2014) and Chow et al. (2012) in their studies.

Build-up and wash-off parameters determine the amount of pollutants discharged to drains. The more impervious the area, the more wash-off coefficients were found. The build-up is slightly higher in the commercial area than in the residential area and the industrial area has a slightly higher build-up than the commercial area (Table 10). Our analysis is different from the previous studies listed in Table 10 in the aspect of land use separation. In our study, roads and roofs were separated to investigate their effect on build-up and wash-off components. Wash-off coefficients for DCIA roads and roofs are higher than for DNCIA, which indicates that they are crucial for pollutant load. Maximum build-up coefficients are relatively small compared to findings by Temprano et al. (2006) and Hood et al. (2007) and the model fit values are similar to the findings of Chow et al. (2012), suggesting that there are cleaning activities on the basin upstream. Street cleaning is active in the basin as it is one of the action plans of the stormwater strategy in Tallinn (Tallinna Linnavolikogu, 2012). In the model, the cleaning efficiency used is 30% to 50% for TSS and up to 90% for TP in DCIA roads at an interval of 7 to 14 days.

Four storm events--08/06/2014, 06/11/2014, 21/05/2015, and 06/08/2015 - were used for model verification. Evaluation was performed by comparing the observed and simulated runoff volume and peak flow for quantity and observed/simulated event load and peak load for TSS and TP using their RE (Table 11). The storm events for this verification have few observations and are not suitable for use with other indicators. A single event of 06/08/2015 was recorded to verify the quality performance. For volume RE% is [+ or -]10% ranging from -9.6% to 6% and for peak flow RE% is nearly [+ or -]10% in a range from -11.0% to 7.8%. Such RE% ranges indicate that the model can sufficiently predict stormwater runoff. However, the quality prediction is moderate for TSS and weak for TP having RE% beyond [+ or -]10%. Therefore, the verification results also suggest that more events of quality observations are needed to calibrate and verify quality performance. Leecaster et al. (2002) proposed seven storms or ~50% of the storms to get annual concentrations within 10% uncertainty whereas Bertrand-Krajewski (2007) recommended at least ten storms to calibrate regression models with a smaller confidence interval. The higher the number of storms captured and calibrated, the lesser model error can be expected.

3.3. Stormwater pollution load: TSS and TP

The simulated event-based concentrations and loads are presented in Table 12 and in Fig. 3. The EMCs of TSS in event 2 and event 3 exceed the national stormwater limit value of 40 mg/L (Vabariigi valitsus, 2013). The literature review by Gobel et al. (2007) showed that EMCs variation is large and depends on urban forms or land use. In their study the TSS loadings ranged from 0.2 to 937 mg/L with roofs (13-120 mg/L) and high-density traffic areas (99-937 mg/L) having high TSS discharges. Loadings of TP were in the range from 0.01 to 0.5 mg/L; 0.06-0.5 mg/L came from roofs and 0.23-0.34 mg/L from traffic areas. In five small catchments with a total area of 116 hectares in a Polish city the concentrations of TSS varied in the range 1.8-736 mg/L (mean ~31 mg/L) (Baralkiewicz et al., 2014). Compared to Estonian regulations, the concentrations of TP in our study were below the national stormwater limit value of 1 mg/L (Vabariigi valitsus, 2013) but higher than the limit of poor status of river quality levels of class IV (i.e. 0.1 mg/L in annex 4, Regulation No. 59) (Keskkonna-minister, 2010). A study of three stormwater catchments in Paris resulted in the TP content range of 0.3-3.52 mg/L (Zgheib et al., 2012) and in a Polish city the TP content varied from 0.02 to 0.57 mg/L (mean ~0.17 mg/L) (Baralkiewicz et al., 2014). The EMC results of the Mustoja basin fell in the range provided by different studies, although simulated TP loads are higher than the usual range. Illegal discharges into the sewerage system can be suspected in the basin, which probably attributed to the higher TP loads. The peak concentrations were higher than the national limits: ~3 to 7 times as high as the limit for TSS and 1.4 to 2.7 times the national TP limits.

The effect on runoff can be observed for up to 19.7 h for event 1 where a single rainfall occurred. The baseflow in this event has a great influence on increasing the volume of runoff (Fig. 3), accounting for approximately 67% of the total runoff. Therefore, the duration of runoff is crucial when estimating event total volume. The storm-flow volume showed an increasing tendency proportional to total rainfall whereas stormflow mass load is likely to increase in proportion to rainfall intensity. In all three events, stormflow was more polluted than baseflow, making up nearly 90% of the total load for TSS and TP, which do not respond to the volume of runoff whether it is higher or lower than the baseflow.


When simulating continuous rainfall for the years 2014 and 2015, the evaporation loss was also considered according to the daily average temperature obtained from the Harku station. Water from the surface is continuously vaporized during the dry period. Annual evaporation was estimated as 0.6 ML (0.014%) in both years. The annual outfall loadings were found to be 97.8 tonne TSS and 1.5 tonne TP from 4400 ML of runoff in 2014 and 110.7 tonne TSS and 1.7 tonne TP from 4500 ML of runoff in 2015 (Fig. 4). The simulated SWMM results were compared with the results from three monitoring programmes. The first one is the monitoring programme during 2012--2014 conducted by the Department of Environmental Engineering of Tallinn University of Technology (DEE) and AS Tallinna Vesi at the outlet of the Mustoja basin approximately 500 m downstream from the studied outlet. This programme was commissioned by Tallinn Environmental Board and the samples were taken six times a year. The same methodology was used in the second monitoring programme but it was conducted by the Estonian, Latvian, and Lithuanian Environment Group (ELLE) in 2015. The third monitoring programme was SA2 conducted by DEE, which was different in methodology in terms of sampling interval as the samples were taken twice a week. The total annual runoff was calculated based on daily average flow, and the total annual loads were calculated based on mean flow and mean concentrations of pollutants (Fig. 4).

In the simulation, the system rainfall was a combination of rainfalls from the Tondi 90 and Saarma stations while DEE, ELLE, and SA2 used rainfall from the Harku station (Fig. 4). Compared to rainfall from the Harku station, system rainfall is by 25% lower for 2014 and 11% lower for 2015. However, the modelled runoffs differ from the DEE and SA2 runoffs to some extent. The modelled runoff is nearly 10% higher than DEE in 2014 and 8% lower than SA2 in 2015. On the contrary, the ELLE runoff is significantly higher at 54% than the simulated runoff, although the annual rainfall is just 11% higher. One reason could be the method of calculation and another reason could be errors in measured values, which can make a difference in the annual runoff and loads. In the ELLE measurement, higher measured flow rates and concentrations have resulted in higher runoff and loads. ELLE TSS is enormously high being about 400% the simulated runoff. However, TSS loads from DEE and SA2 are 35% and 60% lower than the SWMM results for 2014 and 2015. There are quite significant differences in TP loads: 60% between DEE and SWMM, 59% between SA2 and SWMM, and 43% between ELLE and SWMM. These deviations can be attributed to the load calculation method and further to the mean concentration and mean flow because the error in the means could magnify the estimates. Overall, the modelled runoff and loads are higher than estimations from DEE and SA2, but they are significantly lower than ELLE measurements.


3.4. Main contributing impervious surfaces

In the Mustoja basin, the directly connected impervious area or DCIA was determined to be 26.8% of the total land use area (see Table 3). Such percentage of DCIAs was found to be quite effective for the runoff and quality output at 80.1% for peak flow, 75.1% for total runoff volume, 70.5% for TSS, and 66.1% for TP (Table 13). The road area and commercial roof occupy nearly 77% of the total EIA. Thus the DCIA roads and roofs have a higher contribution to the runoff and pollution load, even though the overall runoff coefficient was found to be less at 0.18. Most of the areas in the basin are either pervious or not connected directly to storm drainage. These areas constitute about 73.2% of the total area.

3.5. First flush effect

It would be interesting to examine the first flush phenomenon in the basin if it exists, in order to control the initial contaminated portion through isolation or diverting the stormwater from the road or roof surfaces to treatment facilities. The presence of the first flush phenomenon was assessed by plotting the cumulative fraction of the pollutant load against the cumulative fraction of the runoff volume (Bertrand-Krajewski et al., 1998) as in Fig. 5. The higher loading during storm runoff suggests the presence of first flush (Lee et al., 2004). To account for the intensity of the first flush, the pollutant load swept along by the first 30% of the volume was measured (Temprano et al., 2006; Nazahiyah et al., 2007). In Fig. 5, four events (event 2 to event 5) were used for TSS and TP loadings. One event was during the snowmelt period (event 4: 26/01/2016 to 28/01/2016) and another immediately after the snowmelt (event 5: 29/01/2016 to 01/02/2016). Event 2 and event 3 had no influence from snowmelt, as the former had antecedent dry condition of 12.6 days and the latter had antecedent wet conditions of 1.5 days. The pollutant loadings were 38%, 28%, 50%, and 39% of TSS and 45%, 36%, 45%, and 33% of TP swept by 30% of the runoff volume in events 2, 3, 4, and 5, respectively. The degree of the first flush for event 2 is not high, being 38% for TSS and 45% for TP, but the deviation from the diagonal line in Fig. 5 is evident after 40% of the runoff volume. It suggests that the flushing of the pollutant load is higher at a later stage within the last 60% of the runoff. The number of ADD was almost 13 for this event but the intensity of the rainfall played a greater role than ADD, as a similar result was obtained for TP in event 2 but the deviation is smaller. The first flushes at the initial portion of runoff for events 2, 3, and 5 are not significant. The maximum intensities of rainfall within the 30% of runoff for these events are identified as small at 0.7 mm/h, 1.7 mm/h, and 0.9 mm/h. Instead, the latter part within the last 60% of the volume is more important for the pollutant load. Similar findings were obtained by McCarthy (2009), who found the first flush at the end of the event. Moreover, in a relatively pervious area, Maestre and Pitt (2005) did not observe any first flush. However, Lee et al. (2004) suggested that a seasonal first flush can occur in most cases. Our study also showed that the snowmelt event had a greater influence on the first flush at 50% for TSS and 45% for TP, indicating the effect of a seasonal first flush. Before event 4, there was an extensive period of temperatures below zero, and pollutants accumulated with the snow packing, which was washed off after the temperature became positive during this event.


3.6. Comparison of loadings from four sampling approaches

Finally, the sampling approaches are evaluated for the mass estimations, which provide information for choosing the appropriate sampling option. Four sampling approaches: time weighted sampling (SA1), random grab sampling (SA2), grab sampling within 6 h irrespective of storm size (SA3), and grab sampling within 6 h of medium and large storms (SA4) were analysed for annual average flow, annual average pollutants concentration, and annual pollutants load. Samples from SA3 and SA4 are formed from the data set of SA2 using the recorded time of sampling and rainfall time. Volume weighted mass load estimation from time weighted sampling is considered as the base load because the estimation has a smaller error (Leecaster et al., 2002) compared to other sampling approaches. SA2 is a random sampling that does not respond to the time of corresponding storms. SA3 does not take account of the influencing storm size but is taken within 6 h of the commencement of rain. SA4 is the sampling approach recommended for grab sampling by Lee et al. (2007) and Khan et al. (2006), in which grab samples were taken within 6 h of medium and large storm events. The volume of total runoff is calculated based on the runoff coefficient, catchment area, and runoff depth. Comparison of these four approaches with simulated output is presented in Fig. 6.

The sampling method SA1 with volume weighted estimation is close to the simulated flow rate and annual loads, which indicates that time weighted sampling with volume weighted calculation as specified by Leecaster et al. (2002) provides a smaller error. It is considered as the actual load in the analysis of sampling approaches. The uncertainty of the model in quality estimation has probably built such an error and difference in TSS and TP from SA1. In the case of flow rates and TSS loads obtained from the sampling approach SA4, they are close to the estimations from SA1 and SWMM simulation at around 350 L/s and 60 t/yr, although TP loads have double the deviation. Random sampling or SA2 has often estimated low values of less than half the flow and 2/3 of loads. This is due to the fact that small rainfall or baseflows were mainly sampled in this approach. SA3, which is taken irrespective of the storm size as usual practice, has all of its mean outputs in the middle range of SA2 and SA4. SA4 with mean estimation, on the other hand, is nearer to the actual flow, TSS, and TP. This approach has limitations in identifying medium and large storms. It is difficult to predict the storm size before the rainfall ends. Alternatively, SA3 samples can be used to overcome the limitation after the storm details are retrieved. Therefore, SA3 is a practical and suitable sampling method, which assists in determining the samples of medium and large storms within 6 h after the start of rainfall, although it requires rainfall information and a greater number of samples than SA4.


The model was developed for a large basin in Tallinn and was calibrated and verified for the observed storm events. Additionally, sensitivity analyses for finding sensitive parameters, imperviousness identification, first flush study to search for a possibility of controlling the initial portion of runoff, and evaluation of sampling approaches were carried out. The following conclusions were drawn.

* Sensitivity analysis showed the model to be sensitive to the percentage of the impervious area for predicting both flow rate and peak flow. Impervious depression storage regulates the initial peak flow. Impervious surface roughness and width of catchment have weak connections to the model predictions.

* For the studied large basin, the percentage of the impervious area was found to be 19.7% where the runoff coefficients for different land use vary from 0.05 to 0.57 depending on the proximity of the connection to the drainage system. The overall imperviousness is relatively low, depicting a large basin with mixed land use, which has a high impact on reducing the runoff coefficient. The average depression storage for dry and wet weather was used and the effective value was found to be 0.7 mm.

* The duration of runoff is crucial when estimating event total volume. Stormflow volume has an increasing tendency with total rainfall, whereas stormflow mass load is related to intensity. In the events, pollution is attributed more to stormflow than baseflow.

* The DCIA is an important factor for impervious estimation. Roads and roofs, which are directly connected to storm drainage, are crucial elements of DCIA and these impervious areas can contribute up to 75% of runoff and 66-71% of load.

* During low-intensity rainfalls, the first flush was found to be less effective. The first flush can occur later during 60% of the runoff volume. Therefore, the implementation of treatment facilities to control initial runoff and pollution load may not be effective and more research on the first flush of large basins is required. It was found that the first flush effect during the snowmelt period existed, and its impact cannot be ignored. It will be interesting to simulate snow-melt once there are sufficient input parameters and observations.

* A reliable sampling approach at limited resources is grab sampling within 6 h of storm events, as it provides results close to the simulated results and those obtained by the time weighted sampling method. This approach should focus on capturing medium and large storm events.


Overall, the model development provided information about catchment's stormwater dynamics even at limited resources. The model developed in SWMM had a good performance for quantity estimation, but the quality results were of moderate accuracy. Calibration and verification using more storm events will increase the accuracy. Its applicability will increase even during wintertime once snow period simulation is included. Nevertheless, this study has found crucial impervious areas, percentage of the impervious area, runoff, EMCs, and pollution loads, and evaluated the sampling approaches used. The information will be helpful in planning pollution reduction strategies, implementing pollution control facilities, and designing monitoring programmes in the field of storm-water management.


The authors wish to acknowledge the support of Tallinn Environment Board, the Estonian Ministry of the Environment, and the Estonian Environmental Investment Centre for providing finances for the research. The publication costs of this article were covered by the Estonian Academy of Sciences.


Ballo, S., Liu, M., Hou, L., and Chang, J. 2009. Pollutants in stormwater runoff in Shanghai (China): implications for management of urban runoff pollution. Prog. Nat. Sci., 19(7), 873-880.

Baralkiewicz, D., Chudzinska, M., Szpakowska, B., Swierk, D., Goldyn, R., and Dondajewska, R. 2014. Storm water contamination and its effect on the quality of urban surface waters. Environ. Monit. Assess., 186(10), 6789-6803.

Bedient, P. B. and Huber, W. C. 1988. Hydrology and Floodplain Analysis. Addison-Wesley, Reading, MA.

Bertrand-Krajewski, J.-L. 2007. Stormwater pollutant loads modelling: epistemological aspects and case studies on the influence of field data sets on calibration and verification. Water Sci. Technol., 55(4), 1-17.

Bertrand-Krajewski, J.-L., Chebbo, G., and Saget, A. 1998. Distribution of pollutant mass vs volume in stormwater discharges and the first flush phenomenon. Water Res., 32(8), 2341-2356.

Brezonik, P. L. and Stadelmann, T. H. 2002. Analysis and predictive models of stormwater runoff volumes, loads, and pollutant concentrations from watersheds in the Twin Cities metropolitan area, Minnesota, USA. Water Res., 36(7), 1743-1757.

Burns, M. J., Walsh, C. J., Fletcher, T. D., Ladson, A. R., and Hatt B. E. 2015. A landscape measure of urban stormwater runoff effects is a better predictor of stream condition than a suite of hydrologic factors. Ecohydrology, 8(1), 160-171.

Burszta-Adamiak, E. and Mrowiec, M. 2013. Modelling of green roofs' hydrologic performance using EPA's SWMM. Water Sci. Technol., 68(1), 36-42.

Chiew, F. H. S. and McMahon, T. A. 1999. Modelling runoff and diffuse pollution loads in urban areas. Water Sci. Technol., 39(12), 241-248.

Chow, M. F., Yusop, Z., and Toriman, M. E. 2012. Modelling runoff quantity and quality in tropical urban catchments using Storm Water Management Model. Int. J. Environ. Sci. Technol., 9(4), 737-748.

Ebrahimian, A., Gulliver, J. S., and Wilson, B. N. 2016. Effective impervious area for runoff in urban watersheds. Hydrol. Process., 30(20), 3717-3729.

Elliott, A. H. and Trowsdale, S. A. 2007. A review of models for low impact urban stormwater drainage. Environ.

Modell. Softw., 22(3), 394-405.

EVS. 2015. EVS-EN ISO/IEC 17025:2006. General requirements for the competence of testing and calibration laboratories. Estonian Centre for Standardisation.

[EWS] Estonian Weather Service. 2015. Climate Normal, 2015. (accessed 2016-06-28).

Fetter, C. W. 2001. Applied Hydrogeology. Prentice Hall, Upper Saddle River, NJ.

Gobel, P., Dierkes, C., and Coldewey, W. G. 2007. Storm water runoff concentration matrix for urban areas. J. Contam. Hydrol., 91(1-2), 26-42.

Han, K., Kim, Y., Kim, B., Famiglietti, J. S., and Sanders, B. F. 2014. Calibration of stormwater management model using flood extent data. Proceedings of the ICE - Water Management, 167(1), 17-29.

HELCOM. 2002. HELCOM Recommendation 23/5, Reduction of Discharges from Urban Areas by the Proper Management of Storm Water Systems. 23rd Meeting, Helsinki, Finland, 5-7 March 2002, Minutes of the Meeting. Annex 7. (accessed 2016-06-15).

HELCOM. 2007. HELCOM Baltic Sea Action Plan. HELCOM Ministerial Meeting. Krakow, Poland, 15 November 2007. Sea action plan/BSAP_Final.pdf (accessed 2016-07-12).

Hood, M., Reihan, A., and Loigu, E. 2007. Modeling urban stormwater runoff pollution in Tallinn, Estonia. UWM, international symposium on new directions in urban water management, September. Paris, France.

Horton, R. E. 1933. The role of infiltration in the hydrologic

cycle. Eos Trans. AGU, 14(1), 446-460.

Huber, W. C., Dickinson, R. E., Barnwell, T. O. Jr., and Branch, A. 1988. Storm Water Management Model, Version 4: User's Manual. EPA/600/3-88-001a. (accessed 2016-06-12).

James, L. and Burges, S. 1982. Selection, calibration, and testing of hydrologic models. In Hydrologic Modeling of Small Watersheds (Haan, C. T., Johnson, H. P., and Brakensiek, D. L., eds), 1st edition, pp. 437-472. American Society of Agricultural Engineers, St. Joseph, Michigan.

Jang, S., Cho, M., Yoon, J., Yoon, Y., Kim, S., Kim, G., et al. 2007. Using SWMM as a tool for hydrologic impact

assessment. Desalination, 212(1), 344-356.

Jayasooriya, V. M. and Ng, A. W. M. 2014. Tools for modeling of stormwater management and economics of green infrastructure practices: a review. Water Air Soil Pollut., 225(8), 1-20.

Keskkonnaminister. 2010. Keskkonnaministri 28. juuli 2009. a maaruse nr 44 "Pinnaveekogumite moodustamise kord ja nende pinnaveekogumite nimestik, mille seisundiklass tuleb maarata, pinnaveekogumite seisundiklassid ja seisundiklassidele vastavad kvaliteedinaitajate vaartused ning seisundiklasside maaramise kord" muutmine [Amendment of the Minister of the Environment Regulation No. 44 of 28 July 2009 "Surface water body formation order and list of surface water bodies whose status class has to be determined, status classes of surface water bodies and quality indicator values corresponding to class statuses and order for determining status classes"]. RT I, 25.11.2010 (in Estonian). (accessed 2016-08-20).

Khan, S., Lau, S.-L., Kayhanian, M., and Stenstrom, M. K. 2006. Oil and grease measurement in highway runoff--sampling time and event mean concentrations. J. Environ.

Eng., 132(3), 415-422.

Koppel, T., Vassiljev, A., Puust, R., and Laanearu, J. 2014. Modelling of stormwater discharge and quality in urban

area. Int. J. Ecol. Sci. Environ. Eng., 1(3), 80-90.

Kornecki, T. S., Sabbagh, G. J., and Storm, D. E. 1999. Evaluation of runoff, erosion and phosphorus modeling system--SIMPLE. JAWRA, 35(4), 807-820.

Lee, J. G. and Heaney, J. P. 2003. Estimation of urban imperviousness and its impacts on storm water systems. J. Water Res. Pl.-ASCE, 129(5), 419-426.

Lee, H., Lau, S.-L., Kayhanian, M., and Stenstrom, M. K. 2004. Seasonal first flush phenomenon of urban stormwater discharges. Water Res., 38(19), 4153-4163.

Lee, H., Swamikannu, X., Radulescu, D., Kim, S.-J., and Stenstrom, M. K. 2007. Design of stormwater monitoring programs. Water Res., 41(18), 4186-4196.

Lee, S.-W., Hwang, S.-J., Lee, S.-B., Hwang, H.-S., and Sung, H.-C. 2009. Landscape ecological approach to the relationships of land use patterns in watersheds to water quality characteristics. Landscape Urban Plan., 92(2), 80-89.

Lee, S. C., Park, I. H., Lee, J. I., Kim, H. M., and Ha, S. R. 2010. Application of SWMM for evaluating NPS reduction performance of BMPs. Desalination and Water Treatment, 19(1-3), 173-183.

Leecaster, M. K., Schiff, K., and Tiefenthaler, L. L. 2002. Assessment of efficient sampling designs for urban stormwater monitoring. Water Res., 36(6), 1556-1564.

Liu, A., Goonetilleke, A., and Egodawatta, P. 2012. Inadequacy of land use and impervious area fraction for determining urban stormwater quality. Water Resour. Manag., 26(8), 2259-2265.

Liu, A., Egodawatta, P., Guan, Y., and Goonetilleke, A. 2013. Influence of rainfall and catchment characteristics on urban stormwater quality. Sci. Total Environ., 444, 255-262.

Lowe, R., Thorndahl, S., Mikkelsen, P. S., Rasmussen, M. R., and Madsen, H. 2014. Probabilistic online runoff forecasting for urban catchments using inputs from rain gauges as well as statically and dynamically adjusted weather radar. J. Hydrol., 512, 397-407.

Maestre, A. and Pitt, R. 2005. The National Stormwater Quality Database, Version 1.1. A compilation and analysis of NPDES stormwater monitoring information, 2005. l (accessed 2016-08-10).

Mancipe-Munoz, N. A., Buchberger, S. G., Suidan, M. T., and Lu, T. 2014. Calibration of rainfall--runoff model in urban watersheds for stormwater management assessment. J. Water Res. Pl.-ASCE, 140(6), 05014001.

McCarthy, D. T. 2009. A traditional first flush assessment of E. coli in urban stormwater runoff. Water Sci. Technol., 60(11), 2749-2757.

Nash, J. E. and Sutcliffe, J. V. 1970. River flow forecasting through conceptual models part I--A discussion of principles. J. Hydrol., 10(3), 282-290.

Nazahiyah, R., Yusop, Z., and Abustan, I. 2007. Stormwater quality and pollution loading from an urban residential catchment in Johor, Malaysia. Water Sci. Technol., 56(7), 1-9.

Park, M.-H. and Stenstrom, M. K. 2008. Identification of roads for urban runoff pollution management. In IEEE International Geoscience and Remote Sensing Symposium, IGARSS 2008. July 7-11, 2008, Boston, Massachusetts, USA. Proceedings.

Park, M.-H., Swamikannu, X., and Stenstrom, M. K. 2009. Accuracy and precision of the volume--concentration method for urban stormwater modeling. Water Res., 43(11), 2773-2786.

Riigikogu. 2016. Veeseadus [Water Act]. RT I, 27.12.2016 (in Estonian). (accessed 2016-08-10).

Rosa, D. J., Clausen, J. C., and Dietz, M. E. 2015. Calibration and verification of SWMM for low impact development. JAWRA, 51(3), 746-757.

Rossman, L. A. 2010. Storm Water Management Model User's Manual, Version 5.0. National Risk Management Research Laboratory, Office of Research and Development, US Environmental Protection Agency, Cincinnati, OH.

Shuster, W. D., Bonta, J., Thurston, H., Warnemuende, E., and Smith, D. R. 2005. Impacts of impervious surface on watershed hydrology: a review. Urban Water Journal, 2(4), 263-275.

Tallinna Linnavolikogu. 2012. Tallinna sademevee strateegia aastani 2030 [Tallinn Stormwater Strategy until 2030]. RT IV, 09.03.2013, 41 (in Estonian), (accessed 2016-03-16).

Tallinna Linnavolikogu. 2013. Tallinna arengukava 2014--2020 [Tallinn Development Plan 2014-2020]. RT IV, 25.06.2013, 41 (in Estonian). (accessed 2016-03-16).

Tan, S. B., Chua, L. H., Shuy, E. B., Lo, E. Y.-M., and Lim, L. W. 2008. Performances of rainfall-runoff models calibrated over single and continuous storm flow events. J. Hydrol. Eng., 13(7), 597-607.

Temprano, J., Arango, O., Cagiao, J., Suarez, J., and Tejero, I. 2006. Stormwater quality calibration by SWMM: a case study in Northern Spain. Water SA, 32(1), 55-63.

Tsihrintzis, V. A. and Hamid, R. 1998. Runoff quality prediction from small urban catchments using SWMM. Hydrol. Process., 12(2), 311-329.

Vabariigi valitsus. 2013. Reovee puhastamise ning heit- ja sademevee suublasse juhtimise kohta esitatavad nouded, heit- ja sademevee reostusnaitajate piirmaarad ning nende nouete taitmise kontrollimise meetmed [Regulation on Wastewater and Stormwater Management Requirements, Pollution Parameters and Compliance Limits with the Control Measures]. RT I, 13.06.2013, 13 (in Estonian). (accessed 2016-08-10).

Vezzaro, L. and Mikkelsen, P. S. 2012. Application of global sensitivity analysis and uncertainty quantification in dynamic modelling of micropollutants in stormwater runoff. Environ. Modell. Softw., 27-28, 40-51.

Wanielista, M. P. 1990. Hydrology and Water Quality Control. John Wiley and Sons, New York.

Wanielista, M. P., Kersten, R. D., and Eaglin, R. 1997. Hydrology: Water Quantity and Quality Control. John Wiley and Sons, New York.

Yang, G., Bowling, L. C., Cherkauer, K. A., and Pijanowski, B. C. 2011. The impact of urban development on hydrologic regime from catchment to basin scales. Landscape Urban Plan., 103(2), 237-247.

Zgheib, S., Moilleron, R., and Chebbo, G. 2012. Priority pollutants in urban stormwater: part 1--case of separate storm sewers. Water Res., 46(20), 6683-6692.

Zhang, W. and Li, T. 2015. The influence of objective function and acceptability threshold on uncertainty assessment of an urban drainage hydraulic model with generalized likelihood uncertainty estimation methodology. Water Resour. Manag., 29(6), 2059-2072.

Zhang, G., Hamlett, J. M., Reed, P., and Tang, Y. 2013. Multi-objective optimization of low impact development designs in an urbanizing watershed. OJOp, 2(4), 95-108.

Suurlinna valgala aravoolu, veekvaliteedi ja koormuse modelleerimine

Bharat Maharjan, Karin Pachel ja Enn Loigu

Sademevee aravoolu, saastekoormuse ja maakasutuse maaratlemine on sademevee kaitlemise strateegia elluviimiseks vaga oluline. Modelleerimine voimaldab hinnata neid naitajaid ka piiratud andmehulga pohjal. Kaesolevas uurimuses on rakendatud sademeveekaitlemise mudelit SWMM5 ja arendatud seda uhe Tallinna valgala tarbeks. Geograafilist infosusteemi GIS kasutati alamvalgala ja aravoolu otseselt mojutavate vettpidavate pinnakatetega alade piiritlemiseks ning valgala sisendi parameetrite ettevalmistamiseks. Mudelit kalibreeriti ja kontrolliti juhuslike valing-vihmasadude pohjal, et hinnata saasteainete keskmisi kontsentratsioone ning aastasi koormusi. Mudel voimaldab paris hasti prognoosida arajuhtimist vajava sademevee hulka, kuid selle koostist ennustab ainult moodukal maaral. Mudeli rakendamisel leiti, et kovakattega alade suure osakaalu tottu uuritaval valgalal jaab veepidavus alla 19,7%. Otseselt aravoolule moju avaldavate vettpidavate pindade, eriti teede ja katuste osakaal on suhteliselt madalam, kuid avaldab olulist moju aravoolu kujunemisele (kuni 75%) ning koormustele (kuni 66% ja 71%). Valingvihmade ja lume sulamise aegne aravool ning koormus moodustavad olulise osa aastasest kogukoormusest. Kui rakendatakse pistelist proovivotumetoodikat, peab see suunatud olema keskmistele ja tugevatele sajuhoogudele ning toimuma kuue tunni jooksul peale saju algust, et saada paremaid mootmistulemusi.

Bharat Maharjan*, Karin Pachel, and Enn Loigu

Department of Civil Engineering and Architecture, School of Engineering, Tallinn University of Technology, Ehitajate tee 5, 19086 Tallinn, Estonia

Received 28 September 2016, accepted 29 October 2016, available online 29 May 2017

(*) Corresponding author,
Table 1. Pipe network details in the Mustoja catchment basin

Diameter, m   Number   Length, m  Installation year
             of pipes

0.1-3.0       124       3 335.2   1966--2015
3.0-5.0       652      18 941.3   1956--2015
5.0-7.0       544      17 319.4   1936--2015
7.0-9.0        50       1 506.9   1956--2005
9.0-1.5       197       7 295.1   1956--2005
1.5-2.0        41       1 852.1   1966--2005
Grand total  1608      50 250.0     79

 Diameter, m Typical installation year  Typical material

0.1-3.0      2008/1968            PP, ABS
3.0-5.0      2008/1978/1972       ABS, PP, PVC, concrete
5.0-7.0      1968/1978            Concrete
7.0-9.0      1968/1965            Concrete
9.0-1.5      2003/1972/2000/1967  Concrete
1.5-2.0      1974/1998/1972       Concrete
Grand total

Table 2. Land use details within the Mustoja catchment

Land use                     Area, ha  % of
                                       land use

PERVIOUS                       512.4    50.0
Pond                             0.2     0.0
Stagnant water body/unknown      1.7     0.2
Forest                          39.0     3.8
Shrubbery                        0.8     0.1
Cellar                           0.0     0.0
Green area                     251.0    24.5
Other open area                 60.3     5.9
Private yard                   159.4    15.6
IMPERVIOUS                     511.5    50.0
Ruins                            0.5     0.1
Horticultural land               1.5     0.2
Field of production            104.2    10.2
Greenhouse                       4.0     0.4
Basement                         0.2     0.0
Building under construction      0.0     0.0
Other facility                   0.5     0.0
Road                           235.3    23.0
Building                        99.4     9.7
Garage                           8.8     0.9
Manufacturing building          54.2     5.3
Shelter                          2.5     0.2
Bridge                           0.3     0.0
Grand total                   1024.0   100.0

Table 3. Simplified land use classified into DCIA and DNCIA of

Simplified         Area          TIA,   EIA,
land use           ha      %     ha     ha

Commercial roof     55.9   5.5   50.3   33.8
Industrial          32.9   3.2   14.0   17.6
Residential roof    38.6   3.8   34.7   10.4

Road               147.5  14.4  118.0   94.2

Total              274.8  26.8  217.0  156.0
Commercial roof     13.6   1.3   12.2    1.4
Industrial          71.3   7.0   26.4    3.6
Mixed residential    4.5   0.4    2.1    1.4
and commercial
Residential        481.8  47.1   75.3   24.7

Residential roof    48.1   4.7   43.3    4.8
Road                88.1   8.6   70.5    6.1

Forest              39.8   3.9    2.0    0.4
Water                1.9   0.2    0.0    0.0
Total              749.2  73.2  231.9   42.4


Simplified         Land use details
land use

Commercial roof    Commercial building, manufacturing building
Industrial         Field of production
Residential roof   Building under construction, commercial building,
                   building, shelter
Road               Bridge, roads (road as drain, DCIA road, feeder
                   road, one
                   side vegetated street)
Commercial roof    DNCIA manufacturing building
Industrial         Field of production
Mixed residential  Residential and commercial building
and commercial
Residential        Basement, cellar, garage, green area,
                   horticulture land, other
                   open area, private yard, Private yard, ruins
Residential roof   DNCIA building, manufacturing building
Road               Side vegetated street, road (road along swale,
                   road along
                   swale channel)
Forest             Forest, shrubbery
Water              Pond, stagnant water body/unknown

Table 4. Rainfall characteristics in 2014 and 2015

Station   Year  No. of  Total  Mean rainfall intensity,  Peak,
(months)        events  rain,  mm/h                      mm/h
                        mm     [Range] mean, median,     [Range]
                               0.9 percentile            mean

Tondi     2014  107     363.4  [0.1-4.65] 0.6, 0.4, 1.3   [0.1-7.1]
(V-XII)                                                        1.3
Saarma    2014  98      423.2    [0.2-18.8] 1, 0.5, 1.4  [0.2-19.4]
(V-XII)                                                        1.9
Tondi     2015  156     522.7   [0.1-5.3] 0.5, 0.3, 1.1  [0.1-10.3]
(I-XII)                                                        1.1
Saarma    2015  144     549.6   [0.2-4.8] 0.7, 0.4, 1.3  [0.2-16.4]
(I-XII)                                                        1.3

Station   Total event   Duration,   Inter-event
(months)  rain, mm      h           time, h
          [Range]       [Range]     [Range]
          mean          mean        mean

Tondi     [0.05-10.25]  [0.5-45.5]    [5.75-328]
(V-XII)            1.7           7          52.8
Saarma        [0.1-14]  [0.5-52.5]  [5.75-313.5]
(V-XII)            2.2         6.5          57.7
Tondi      [0.05-17.1]  [0.5-60.5]   [8.5-558.5]
(I-XII)            1.7         8.1          55.2
Saarma      [0.1-15.7]  [0.5-54.5]    [10-352.5]
(I-XII)            1.9         7.7          59.8

Table 5. Characteristics of sampled storm events

Event              Antecedent  Duration,  Mean rainfall    Peak, mm/h
                   dry days    h          intensity, mm/h

Event 1            12.1         2.3       2.2              8.7
(9 Sep. 2014)
Event 2            12.6        26.0       0.8              3.3
(22-23 Sep. 2014)
Event 3             1.5        13.5       1.2              3.1
(4-5 Dec. 2015)

Event              Total event   Rainfall intensity,
                   rainfall, mm  mm/day

Event 1            5.1           4.3
(9 Sep. 2014)
Event 2            6.2           4.3
(22-23 Sep. 2014)
Event 3            9.7           7.6
(4-5 Dec. 2015)

Table 6. Calibrated values of input parameters for simulating runoff

Parameter                   Range (Reference)

% Imp factor                [+ or -]10%
W factor                    [+ or -]10%
[D.sub.imp]                 0.3 to 2.3 (Huber et al., 1988)
[D.sub.per]                 2.5 to 5.1 (Huber et al., 1988)
[N.sub.imp]                 0.01 to 0.03 (Wanielista et al., 1997)
[N.sub.per]                 0.02 to 0.45 (Huber et al., 1988)
Maximum infiltration, mm/h  50 to 200 (Bedient and Huber, 1988)
Minimum infiltration, mm/h  0.5 to 12 (Nazahiyah et al., 2007)
Decay constant, L/h         0.000389 to 0.0039 L/s i.e. 1.4 to 14 L/h
                            (Nazahiyah et al., 2007)

Parameter                    Calibrated values

% Imp factor                 0.9
W factor                     1
[D.sub.imp]                  0.7
[D.sub.per]                  3
[N.sub.imp]                  0.0135
[N.sub.per]                  0.2
Maximum infiltration, mm/h  50
Minimum infiltration, mm/h   0.5
Decay constant, L/h

Table 7. Sensitivity coefficient for input parameters

Parameter        Sensitivity coefficient
               Runoff  Runoff             Peak flow
               volume  flow rate

% Imp factor   0.610   0.816              0.687
W factor       0.309   0.439              0.357
[D.sub.imp]   -0.029  -0.019             -0.029
[D.sub.per]    0.000   0.000              0.000
[N.sub.imp]   -0.001  -0.018             -0.008
[N.sub.per]    0.000   0.000              0.000

NA--not available.

Table 8. Goodness of fit for runoff quantity simulation

Goodness of     Acceptable         Event 1               Event 2
fit indicators  range        Flow rate  Peak flow  Flow rate

CC              Close to 1   0.98         NA       0.87
RE %            [+ or -]10%   4.1       -27.4      5.6
NOF             0 to 1        0.3         0.3      0.2
NSC             Close to 1    1.0          NA      0.6

Goodness of     Event 2          Event 3
fit indicators  Peak flow  Flow rate  Peak flow

CC               0.99       0.94       1.00
RE %             0.6       19.4       21.2
NOF              0.1        0.3        0.2
NSC              0.8        0.7        0.3

Table 9. Goodness of fit for quality simulation

Goodness of     Acceptable       Event 2            Event 3
fit indicators  range        TSS      Peak TSS  TSS      Peak TSS

CC              Close to 1     0.8       0.8     0.4      -1.0
RE%             [+ or -]10%  -1.12     -14.5     4.3     -13.4
NOF             0 to 1         0.4       0.3     0.3       0.3
NSC             Close to 1     0.7     -48.8    -0.5      -1.0

Goodness of          Event 2           Event 3
fit indicators   TP       Peak TP  TP       Peak TP

CC                0.9      1.0     0.9      0.9
RE%             -33.0    -15.6    11.0     -7.9
NOF               0.5      0.2     0.2      0.1
NSC               0.4      0.4     0.7      0.5

Table 10. Calibrated build-up and wash-off input parameters

Land use                              Total phosphorus
                       Build-up                  Wash-off
                       Coeff.      Exponent      Coeff.      Exponent

This study
R                      0.2         0.02          0.2         0.7
Mixed R and C          0.2         0.03          0.3         1
I                      0.25        0.016         0.5         1.2
R roof
[DCIA, DNCIA]         [0.2, 0.2]  [0.03, 0.02]  [0.2, 0.2]  [0.8, 1.4]
C roof
[DCIA, DNCIA]         [0.2, 0.2]  [0.03, 0.03]  [0.5, 0.2]  [1.2, 1.2]
[DCIA, DNCIA]         [0.4, 0.3]  [0.08, 0.04]  [0.5, 0.2]  [0.7, 1.2]
Other studies
Temprano et al.,      46           0.3           0.13        1.2
2006, Spain, R
Hood et al., 2007,     0.25        0.0025      500           2.35
Estonia, Mixed
Chow et al., 2012,    [0.3, 0.5,  [0.05, 0.1,   [0.41,      [1.46, 1,
Malaysia, R, C, I      0.3]        0.16]         0.4, 0.8]   1.08]
Koppel et al., 2014,
Estonia, C and road

Land use                              Total suspended solids
                        Build up              Wash-off
                        Coeff.    Exponent    Coeff.       Exponent
This study
R                       3         0.8         0.3          1.5
Mixed R and C           3         0.75        0.3          1.2
I                      10         0.7         0.5          1.5
R roof
[DCIA, DNCIA]          [3, 3]    [0.8, 0.8]  [0.25, 0.2]  [1.5, 1.5]
C roof
[DCIA, DNCIA]          [10, 10]  [0.8, 0.8]  [0.6, 0.2]   [1.5, 1.5]
[DCIA, DNCIA]          [10, 10]  [0.8, 0.8]  [0.6, 0.2]   [1.5, 1.5]
Other studies
Temprano et al.,        46        0.3         0.13         1.2
2006, Spain, R
Hood et al., 2007,      25        1           4.9          1.57
Estonia, Mixed
Chow et al., 2012,      [3, 13,  [0.8, 0.8,  [0.2, 1.4,   [1.4, 0.9,
Malaysia, R, C, I       15]       0.7]        3]           0.6]
Koppel et al., 2014,              0.2                      1
Estonia, C and road

Table 11. Verification of model performance

Event       Rain,  Verification  Observed  Simulated  RE, %
            mm     parameter

                                      Volume, [10.sup.3] L/s
08/06/2014  21.07  Runoff         7443.9    6999.1     6.0
06/11/2014     43  Runoff        12597.7   13372.2    -6.1
21/05/2015   53.4  Runoff        17031.1   17612.1    -3.4
06/08/2015   12.4  Runoff         6118.2    6705.4    -9.6
                                          Load, kg
06/08/2015   12.4  TSS            273.54    246.92     9.7
06/08/2015   12.4  TP              1.540     2.376   -54.3

Event         Observed  Simulated  RE, %

                        Peak flow, L/s
08/06/2014    910.0     839.3        7.8
06/11/2014    736.7     761.3       -3.3
21/05/2015    807.0     888.6      -10.1
06/08/2015    230.1     255.4      -11.0
            Peak flow, mg/L
06/08/2015    111.5     83.4        25.2
06/08/2015      0.5      0.8       -41.5

Table 12. Event-based runoff and concentrations of TSS and TP

Event    Rain,  Runoff     Total flow   Event flow
         mm     duration,  volume,      Mean,       Peak,
                h          [10.sup.6]L  L/s         L/s

Event 1  5.1    19.7       16.3         229.8       1304.9
Event 2  6.2    23.7       16.0         187.6        452.1
Event 3  9.7    13.2       20.0         422.6        915.0

Event    Event TSS loading  Event TP loading
         EMC,  Peak,        EMC,  Peak,
         mg/L  mg/L         mg/L  mg/L

Event 1  33.6  288.9        0.5   2.7
Event 2  50.3  170.9        0.5   1.2
Event 3  69.1  119.2        0.7   1.4

Table 13. Contribution of DCIA to impervious land use

Event                  DCIA/EIA, %
            Max flow,  Flow volume,  TSS,  TP,
            LPS        [10.sup.6] L  kg    kg

Event 1     75.4       50.3          66.3  70.6
Event 2     83.7       89.8          72.8  47.6
Event 3     81.4       85.3          72.4  80.0
All events  80.1       75.1          70.5  66.1
COPYRIGHT 2017 Estonian Academy Publishers
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Maharjan, Bharat; Pachel, Karin; Loigu, Enn
Publication:Proceedings of the Estonian Academy of Sciences
Article Type:Report
Date:Sep 1, 2017
Previous Article:On approximation processes defined by a cosine operator function.
Next Article:Global linearization approach to nonlinear control systems: a brief tutorial overview.

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters