## Abstract

This paper presents a comprehensive framework for the quantitative management of urban runoff. The framework assesses the response of urban catchments to design rainfall events and identifies low-impact development (LID) stormwater control measures (SCMs) for runoff control and flood mitigation. This research's method determines the optimal areas in which to deploy SCMs to control runoff in urban catchments. The optimization method relies on a three-objective simulation-optimization model that (1) minimizes the volume of runoff at the catchment outlet and at flooding nodes, (2) minimizes the implementation and maintenance costs of LID SCMs, and (3) minimizes the service-performance reduction of LID SCMs. The storm water management model (SWMM) is applied for runoff simulation and is coupled with the multi-objective antlion optimization algorithm (MOALOA). The simulation-optimization method is exemplified with an application to District 6 of Tehran's municipality (Iran). The performance of the simulation-optimization method is compared with that of the multi-objective non-dominated sorting genetic algorithm II (NSGAII), and, after confirming the superior capacity of the MOALOA, the latter algorithm is applied to District 6 of Tehran municipality, Iran. The identified optimal LID SCMs are ranked with the technique for order of preference by similarity to ideal solution (TOPSIS) method that reveals the preferences of the runoff managers concerning SCMs choices. The most desirable solution herein found shows the optimal LID SCMs provide a significant reduction in runoff volume at the catchment outlet and flooding nodes.

## INTRODUCTION

Increased urbanization and the replacement of pervious or semi-pervious surfaces with impervious surfaces has been a key factor causing more runoff from precipitation in urban areas (Davis 2005), which increases peak flows and reduces the time to peak flow (Holman-Dodds *et al.* 2003). Urban runoff carries various pollutants from urban areas into the receiving waters (Davis *et al.* 2009; Barkdoll *et al.* 2016). Studies show that in highly urbanized areas, the volume of runoff generation is significant (Xu & Zhao 2016; Yao *et al.* 2018); moreover, as a result of intensifying extreme rainfall events flood risk increases (Hettiarachchi *et al.* 2018). Therefore, appropriate management of urban runoff with low-impact development (LID) stormwater control measures (SCMs) is currently a preferred approach to urban runoff management (Sadeghi *et al.* 2017).

LID SCMs constitute an approach to stormwater management that emphasize the use of small scale, water-retaining features within urban areas to slow, purify, infiltrate and capture urban runoff and precipitation at the point of origin (Davis 2005; Chau 2009; Morsy *et al.* 2016).

There are several types of LID SCMs with varying installation and maintenance costs. Bioretention cells, vegetated swales, rain gardens, porous pavements, dry (infiltration) wells, and infiltration trenches are the most common LID SCMs. They exhibit variable performance in runoff infiltration, detention, and removal of pollutants. Various studies of LID SCMs have been conducted to evaluate their performance in terms of quantitative and qualitative control of urban runoff, such as those dealing with bioretention cells (Davis 2008; Dietz & Clausen 2008; Hunt *et al.* 2008; Carpenter & Hallam 2009; James & Dymond 2011;,Lucke & Nichols 2015; York *et al.* 2015), on porous pavement (Collins *et al.* 2009; Fassman & Blackbourn 2010; Brunetti *et al.* 2016; Huang *et al.* 2016; Alsubih *et al.* 2017; Kamali *et al.* 2017), and infiltration trenches (Emerson *et al.* 2010) and on vegetated swale (Fach *et al.* 2011).

Successful and economical application of LID SCMs requires careful selection, design, and the area over which they are deployed. An effective approach to SCM selection, sizing, and deployment is to combine runoff simulation with optimization of SCMs. In many applications there are many sub-catchments within a catchment, complicating the choice of the most effective layout of candidate LID SCMs without resorting to optimization approaches (Cano & Barkdoll 2016; Eckart *et al.* 2018). Several authors have addressed the problem of SCMs selection, layout, and design with quantitative, qualitative, and economic criteria by applying the simulation-optimization approach. In some studies, minimizing runoff volume or runoff peak flow and LID SCMs costs (quantitative and economic criteria) were considered (Jia *et al.* 2012; Damodaram & Zechman 2012; Tao *et al.* 2014; Eckart 2015; Loáiciga *et al.* 2015; Chui *et al.* 2016; Aminjavaheri & Nazif 2018). Other studies considered the objectives of minimizing runoff pollution and LID SCMs costs (qualitative and economic criteria) (Baek *et al.* 2015; Tobio *et al.* 2015). Several authors combined quantitative, qualitative and economical criteria in the selection, sizing, and location of SCMs for runoff control (Limbrunner *et al.* 2013; Montaseri *et al.* 2015; Sadeghi *et al.* 2017). The cited studies assumed temporal constant performance for each LID SCM neglecting the effect of clogging on their performance over time.

LID SCMs are impaired by clogging due to sedimentation of suspended solids or by the accumulation of oil and grease, blockage by accumulated trash, or by reduction of their infiltration and filtration capacities by biological growth. LID SCMs performance declines over time despite periodic maintenance (Haselbach *et al.* 2016) until replacement takes place. Many experimental, field, and theoretical studies have been conducted to assess clogging effects on LID SCMs and their long term performance (Haselbach *et al.* 2006; Lucke & Beecham 2011). These studies determined the effect of clogging as the ratio of the permeability reduction to the permeability in the unclogged state (Deo *et al.* 2010; Freni & Mannina 2018), as a percentage of the long term hydraulic conductivity reduction (Li & Davis 2008; Pezzaniti *et al.* 2009), or as an analytical probabilistic expression (Zhang & Guo 2014). The reduction of LID SCMs service performance over time varies among LID SCMs depending on site-specific features.

This paper presents a multi-objective simulation-optimization model to determine optimal combination of LID SCMs in urban areas. The model includes economic and quantitative criteria, and considers the reduction of LID SCMs service performance over time. The consideration of service-performance reduction of LID SCMs has been rarely accounted for in previous studies. The simulation-optimization model links the SWMM simulation model to the multi-objective antlion optimization algorithm (MOALOA) by coding them in MATLAB2012b considering three-objective goals: (1) minimizing the runoff volume at the catchment outlet and the overflow volume at nodes to increase urban drainage system reliability that reduces flood potential, (2) minimizing the cost of implementation and maintenance of LID SCMs, and (3) minimizing the service-performance reduction of LID SCMs by choosing efficient SCMs layout. The solutions of the multi-objective optimization of SCMs are expressed in terms of Pareto possibility frontiers, or Pareto fronts, that represent alternative optimal combinations of SCMs with inherent tradeoffs among objectives. The TOPSIS method (Hwang & Yoon 1981) is implemented for selecting the most desirable LID SCMs solution according to weighting criteria.

## MATERIALS AND METHODS

The outline of this work is as follows:

Choosing a simulation model and optimization algorithm.

The simulation-optimization model is developed by coupling the simulation model to the optimization algorithm.

The performance of the optimization algorithm is evaluated by means of an example problem from the SWMM manual.

The simulation model is built for the study area (District 6 of Tehran municipality).

The design rainfall hyetograph is developed and applied to the simulation model.

The catchment's stormwater infrastructure is evaluated.

Suitable LID SCMs are selected for the study area conditions and from results of previous studies.

The simulation-optimization model is implemented and optimal management alternatives are determined.

The most desirable solution is selected among those on the optimal Pareto front using the TOPSIS method.

The next section explains the methodology in detail.

### Rainfall-runoff simulation model for the study area

The SWMM is implemented as the rainfall-runoff model. Afterwards, appropriate LID SCMs are identified for improving the study area's drainage system performance.

#### The SWMM

Among the simulation models, the US Environmental Protection Agency (USEPA) storm water management model (SWMM) has been widely used for qualitative and quantitative simulation of runoff in urban areas and LID SCMs. Jayasooriya & Ng (2014) evaluated and compared 20 models for urban runoff management. They singled out SWMM developed by the United States environmental protection agency (USEPA) as one of the most accurate and practical models.

The site or catchment simulated with SWMM is divided into subcatchments. Subcatchments receive precipitation and generate runoff and pollutant loads. Runoff and pollution loads are routed through channels and pipes to the outlet of the catchment (Rossman 2015).

#### LID SCMs simulation

SWMM 5.1 contains a toolbox to model common types of LID SCMs explicitly and computes the absorption and retention of runoff and pollutants. There are several LID SCMs that are chosen based on site features such as land use, terrain slope, environmental characteristics, and precipitation regime. This study considers several LID SCMs described in the SWMM's user manual. The LID SCMs considered for deployment in the main case study (District 6 of Tehran municipality) are porous pavement and bioretention cells. Advantages of porous pavement are its capacity to reduce the volume of runoff and peak flow (Fassman & Blackbourn 2010; Palla & Gnecco 2015; Alsubih *et al.* 2017). Moreover, it is ideal for sites with limited space for other LID SCMs, making it a good choice for highly urbanized environments. Previous studies have shown considerable reduction of runoff volume and peak flow by bioretention cells (Dietz 2007; Line & Hunt 2009; Trinh & Chui 2013; Eckart *et al.* 2017). Bioretention cells are suitable for implementing alongside streets and on parking lots to control runoff and provide natural green spaces.

The parameters of LID SCMs are determined and defined in the SWMM according to site conditions (District 6 of Tehran municipality) as reported in previous studies (PaDEP 2006; Fallahi Zarandi 2013; Eckart 2015; Chui *et al.* 2016) and the SWMM users' manual (Rossman 2015). The parameter values are listed in Table 1. The soil characteristics are nearly constant in the study area; therefore, similar LID SCMs were considered in all the subcatchments. The maximum percentage of clogging or blockage in porous pavement, bioretention cells, vegetated swales, and infiltration trenches were set equal to 40%, 20%, 20%, and 30%, respectively. These values were assigned considering each LID SCMs performance in previous laboratory studies and field studies in similar study areas (Siriwardene *et al.* 2007; Li & Davis 2008; Pezzaniti *et al.* 2009; Deo *et al.* 2010; Lucke & Beecham 2011; Razzaghmanesh & Becham 2018). The drainage area for each LID SCMs in each subcatchment was estimated based on land use maps. Pedestrian areas, low volume roads, and parking areas for implementing porous pavement, and alongside streets and parking lots were taken into account in relation to bioretention cells.

Parameter | Bioretention cells | Porous pavement | Vegetated swale | Infiltration trench | |
---|---|---|---|---|---|

Surface | Berm height (mm) | 150 | 4 | 200 | 100 |

Vegetation volume fraction (%) | 0.90 | 0 | 0.90 | 0 | |

Surface roughness coefficient | 0.03 | 0.01 | 0.03 | 0.11 | |

Surface slope (%) | 0 | According to subcatchment slope | According to subcatchment slope | ||

Soil | Thickness (mm) | 900 | – | – | – |

Porosity (%) | 0.50 | – | – | – | |

Hydraulic conductivity (mm/h) | 250 | – | – | – | |

Pavement | Thickness (mm) | – | 100 | – | – |

Void ratio | – | 0.15 | – | – | |

Impervious surface fraction (%) | – | According to subcatchment impervious surface fraction | – | – | |

Permeability (mm/h) | – | 500 | – | – | |

Storage | Thickness (mm) | 500 | 350 | – | 700 |

Void ratio | 0.70 | 0.50 | – | 0.7 | |

Seepage rate (mm/h) | 44 | 44 | – | 44 |

Parameter | Bioretention cells | Porous pavement | Vegetated swale | Infiltration trench | |
---|---|---|---|---|---|

Surface | Berm height (mm) | 150 | 4 | 200 | 100 |

Vegetation volume fraction (%) | 0.90 | 0 | 0.90 | 0 | |

Surface roughness coefficient | 0.03 | 0.01 | 0.03 | 0.11 | |

Surface slope (%) | 0 | According to subcatchment slope | According to subcatchment slope | ||

Soil | Thickness (mm) | 900 | – | – | – |

Porosity (%) | 0.50 | – | – | – | |

Hydraulic conductivity (mm/h) | 250 | – | – | – | |

Pavement | Thickness (mm) | – | 100 | – | – |

Void ratio | – | 0.15 | – | – | |

Impervious surface fraction (%) | – | According to subcatchment impervious surface fraction | – | – | |

Permeability (mm/h) | – | 500 | – | – | |

Storage | Thickness (mm) | 500 | 350 | – | 700 |

Void ratio | 0.70 | 0.50 | – | 0.7 | |

Seepage rate (mm/h) | 44 | 44 | – | 44 |

### Optimization

_{ij}= area of LID SCM number j in subcatchment i, n

_{s}= number of subcatchments, j represents the type of LID SCMs 1 = porous pavement, 2 = bioretention cells, 3 = infiltration trench and 4 = vegetated swale (only LID SCMs number 1 and 2 were considered in case study 2).

The decision variables are embedded in the simulation model linked to the optimization algorithm. Equations (2)–(4) express respectively objective functions 1, 2, and 3.

*Z*

_{1}= objective function 1,

*VO*

_{1}= the volume of runoff at the outlet for normal operation when clogging and failure of LID SCMs do not occur,

*VO*= the minimum runoff at the outlet when LID SCMs cover the maximum allowable area,

_{min}*VO*= the maximum volume of runoff at the outlet when ID SCMs are not deployed,

_{max}*FV*

_{1}= sum of overflow volume at all flooding nodes and for normal situation when clogging and failure of LID SCMs do not occur

*, FV*

_{min}*=*the sum of minimum flood volume at all flooding nodes when LID SCMs cover the maximum allowable area,

*FV*

_{max}*=*the sum of maximum overflow volume at all flooding nodes when LID SCMs are not deployed. The first term on the left-hand side of objective function (1) quantifies reduction of runoff volume at the outlet, while the second term quantifies decreasing flood potential at nodes where drainage system does not have sufficient capacity to pass the design runoff and overflows. Z

_{1}equals zero when all LID SCMs cover their maximum allowable areas, and Z1 becomes 1 when LID SCMs are not deployed. Considering FV and VO is necessary to prevent flooding inside and downstream the catchment respectively. Reducing the amount of runoff at the outlet of the catchment prevents flooding downstream the catchment which is a vulnerable urban area. This consideration is necessary in this study for reducing runoff overflow at flooding nodes inside the catchment.

*Z*

_{2}

*=*objective function 2,

*J*

*=*the number of subcatchments,

*L*

*=*the number of LID SCM types

*, X*(

_{l}*j*) = the area of the

*l*LID SCM in the

^{th}*j*subcatchment, UC

^{th}_{l}= the unit cost of the

*l*LID SCM type,

^{th}*Maximum Cost*= the cost of construction and maintenance of all LID SCMs covering the maximum allowable area.

*Z*

_{3}= objective function 3,

*VO*

_{1}and

*VO*

_{2}are defined identically as in Equation (1),

*VO*

_{2}= the volume of runoff at the outlet when LID SCMs undergo maximum clogging and or blockage, and

*VF*

_{2}= the sum of runoff volume at flooding nodes when LID SCMs undergo maximum expected clogging or blockage.

Concerning the second terms of objective functions 1 and 3 notice that if there are no flooding nodes in the study area then the runoff peak flow in the catchment outlet is considered instead of the sum of overflow volume at flooding nodes (FV) for the purpose of reducing the risk of downstream flooding.

#### The multi-objective antlion optimization algorithm

The antlion optimization algorithm (ALOA) and multi-objective version of the algorithm was introduced by Mirjalili (2015) and Mirjalili *et al.* (2017). It is inspired by the antlion hunting behavior at the larva stage. Mirjalili *et al.* (2017) applied the MOALOA for solving engineering problems and compared the MOALOA performance solving test mathematical functions with that of the particle swarm optimization (PSO) algorithms and the non-dominated sorting genetic algorithm (NSGAII). The MOALOA is implemented in the present study because of its relative superior performance in those comparative studies. The performance of MOALOA in the simulation-optimization model of the present study is first evaluated with the optimization of the example runoff problem found in the SWMM application manual (Gironás *et al.* 2015).

The ALOA and MOALOA apply ants (i.e., the solutions) as search agents; the ants move over the decision space. Antlions (also the solutions because each ant has an antlion counterpart) combat and hunt the ants making the latter become fitter. The position of each ant is updated in each algorithmic iteration. The roulette wheel selection operator selects solutions (ants, antlions) with better fitness-function values (Pencheva *et al.* 2009; Mirjalili 2015). Table 2 lists the characteristics of the ALOA and MOALOA. The MOALOA and ALOA algorithms initialize the position of antlions and ants randomly, and calculates their fitness-function values, thus determining the elite ants and antlions. The roulette wheel operator selects one antlion for each ant in each iteration. Each ant's fitness function value is calculated: if an ant becomes fitter than its corresponding antlion its position becomes the new position of its antlion in the next iteration. The positions of ants are updated according to random walk (Mani *et al.* 2018).

General algorithm | Antlion optimizer |
---|---|

Decision variable | Antlion's and ant's position in each dimension |

Solution | Antlion's position |

Old solution | Old position of antlion |

New solution | New position of antlion |

Best solution | Elite antlion |

Fitness function | Desirability of elite |

Initial solution | Random antlion |

Selection | Roulette wheel |

The process of generating new solution | Random walk over the decision space |

General algorithm | Antlion optimizer |
---|---|

Decision variable | Antlion's and ant's position in each dimension |

Solution | Antlion's position |

Old solution | Old position of antlion |

New solution | New position of antlion |

Best solution | Elite antlion |

Fitness function | Desirability of elite |

Initial solution | Random antlion |

Selection | Roulette wheel |

The process of generating new solution | Random walk over the decision space |

The general framework of population-based evolutionary algorithms is similar (Bozorg-Haddad *et al.* 2017). The algorithms start the optimization process with generated populations of solutions which are compared with each other (Sarzaeim *et al.* 2018). The non-dominated solutions remain in the search algorithm. The evolutionary algorithms attempt to improve remaining solutions in the next iteration. What makes evolutionary algorithms different from one another is the method applied to enhance the non-dominated solutions.

Multi-objective optimization considers diversity and convergence in the improvement of the non-dominated solutions (Mirjalili *et al.* 2017). There are different methods for enhancing convergence such as archive and leaders, the non-dominated sorting mechanism, and niching. The MOALOA applies the archive-and-leaders approach to ensure suitable coverage.

### Selecting an option from a set of alternative optimal solutions

Multi-criteria decision-making (MCDM) methods are implemented to select an option among several suitable alternative choices or solutions. The result of multi-objective optimization algorithms are expressed as a set of non-dominated solutions that form a Pareto possibility frontier or Pareto front. The technique for order of preference by similarity to ideal solution (TOPSIS) is a leading MCDM method used in water resources systems, and it was chosen here for selecting a most desirable solution among several alternative optimal solutions.

#### The TOPSIS method

The TOPSIS method attempts to prioritize management choices on the basis of proximity to the ideal choice and distance from the anti-ideal choice. TOPSIS takes into consideration the weights of criteria and the values of ideal and anti-ideal solutions of an optimization problem and ranks the solutions on the basis of their relative Euclidean distance from the ideal and anti-ideal solutions. TOPSIS identifies the solution closest to the ideal solution and farthest from the anti-ideal solution as the superior solution (Azarnivand & Banihabib 2016). The ideal solution is the most desirable to the decision-maker. Each optimization objective may have a separate ideal solution. The closer a solution is to the ideal value, the more desirable to the decision-maker. The anti-ideal solution is the worst solution for the decision-maker, and the decision-maker invariably prefers the solution obtained to be far from the anti-ideal solution (Banihabib *et al.* 2015; Bozorg-Haddad *et al.* 2016).

The steps of implementation of the TOPSIS method are as follows (Inanlu *et al.* 2012):

Constructing a decision matrix whose rows and columns represent criteria and choices (or solutions) that receive a quantitative value.

Normalizing the metrical elements' values to construct a normalized data matrix.

- Determining the weight of each criteria and computing the weighted values for each matrix element through Equation (5): where
*F*= the weighted value of the_{ij}*i*row (criterion) and^{th}*j*column (choices or solutions) of the decision matrix, = the normalized value of the^{th}*i*row and^{th}*j*column of the decision matrix, and^{th}*w*= the_{i}*i*criterion weight.^{th} Ranking the solutions on the basis of the value of

*CL*such that the closer this amount is to zero, the lower the rank; and the closer it is to 1, the higher the rank._{i}

The decision-making must choose from the set of optimal solutions of the multi-objective optimization problem, each of which is superior to others at least in terms of one objective. This study applies the TOPSIS method with equal weights to the three-objective functions to choose a most desirable solution among those present in the MOALOA-calculated Pareto front of solutions.

## CASE STUDIES

Two case studies are solved in this work. The first is the example appearing in the application manual of the SWMM. The solutions calculated with the MOALOA are compared with those of the NSGAII to serve as a test of the MOALOA relative to the well-established NSGAII.

The second case study deals with District 6 of Tehran municipality. This district is predominantly residential with predominant impervious surface, which causes flooding during heavy rains. MahabGhods Consulting Engineering Company (2011a, 2011b) recommended the use of LID SCMs in all areas of the city of Tehran.

### Sample case of the SWMM application manual

The SWMM's example is a catchment of 117,359 m^{2} divided into seven subcatchments. The areas of the subcatchments range between 8012 and 27,478 m^{2}, and their slopes range between 2% and 3%. The fraction of impervious areas in the subcatchments varies between 0% and 40%. The input data including the design rainfalls with return periods of 2, 5, and 10 years and subcatchments' features are found in the SWMM manual (Gironás *et al.* 2015). Here the SWMM is implemented with a 2-hour design rainfall with a 10-year return period, in which 43.45 mm of rainfall occurred within 2 hours.

### District 6 of Tehran municipality

District 6 of Tehran municipality has an area equal to 2,137.18 hectares (1 ha = 10^{4} m^{2}), and is in central Tehran. This is a highly urbanized site including residential areas, streets, parks, and the University of Tehran. It is shown in Figure 3 that District 6 is divided into 45 subcatchments (Tehran University-Graduate Faculty of Environment 2015). The slope of the subcatchments varies between 1% and 4%, and the impervious fraction near 100% in most of the subcatchments and declining to 30–40% in a few subcatchments where there are parks with pervious surfaces. The drainage system of this area includes 30 circular buried channels with 1.6-meter depth, and one rectangular open channel with 0.7-meter depth. The Manning roughness coefficient of channels are in the range of 0.012–0.035.

*i*= rainfall intensity (mm/h) and

*D*= rainfall duration (minutes).

### Calculating costs of implementation and maintenance of LID SCMs

The expense for LID SCMs includes all costs of their construction and maintenance during their 20-year service life. It is determined based on the list of values issued by the Management and Planning Organization of Iran (2017) and the study by Fallahi Zarandi (2013). The costs from the study by Fallahi Zarandi are scaled to 2017 costs by applying the cost-appreciation rate applicable in Iran (Central Bank of the Islamic Republic of Iran 2017).

*et al.*2017). Equation (12) converts these costs to current value: where

*C*= current equivalent cost,

*A*= annual cost,

*N*represents the service life of LID SCMs, in years.

## RESULTS

### Evaluating the performance of the SWMM-MOALOA simulation-optimization model with a SWMM sample problem

The relative superior performance of the MOALOA compared with the NSGAII and the multi-objective particle swarm optimization (MOPSO) algorithm has been demonstrated with several optimization problems (Mirjalili 2015; Mirjalili *et al.* 2017). The SWMM was linked with the MOALOA and with the NSGAII in MATLAB Version 2012. The SWWM manual's sample problem was solved with the two simulation-optimization models. The simulation-optimization model was run 10 times with the MOALOA-SWWM and the NSGAII-SWWM under identical conditions using a Core i7 computer with a processing speed of 3.5 GHz and a memory of 8 GB. The sample catchment is divided into seven subcatchments (Gironás *et al.* 2015), one of which is 100% pervious and does not require runoff control. Four LID SCMs were considered for deployment, namely, porous pavement, bioretention cells, infiltration trench, and vegetated swales, in the six remaining subcatchments. This means there are 24 decision variables to be optimized. These decision variables are the areas over which each LID SCM would be constructed within each subcatchment. The unit costs of implementation and maintenance for porous pavement, bioretention cells, vegetated swales, and infiltration trench were set equal to $10, $19.75, $7.64, and $16.15, respectively. The objective functions defined for this sample problem are those expressed by Equations (2)–(4), except that the parameter *FV* (the sum of overflow volume at flooding nodes) is replaced by the peak flow at the outlet in Equations (2) and (4). There are no flooding nodes in this example; therefore, peak flow was deemed the appropriate variable to simulate in this example to reduce the risk of downstream flooding. This problem was solved using MOALOA-SWMM and NSGAII-SWMM with 100 iterations and 50 prospective solutions in the initial population of solutions. Figure 1 illustrates the Pareto front resulting from the combination of the 10 runs of each algorithm.

As seen in Figure 1(a), the solutions of the Pareto front from 10 runs of the MOALOA-SWMM (whose calculated solution are markers shaped as stars) have a closer distance to the origin of coordinate than the NSGAII-SWMM solutions, thus the former solutions are closer to the ideal solution. Table 3 lists various features of the simulation-optimization runs. As seen in Table 3 the average Euclidean distances of the MOALOA-SWMM Pareto fronts (recall there were 10 runs) from the origin of coordinate vary between 0.4 and 0.62, while those of the NSGAII-SWMM vary between 0.53 and 0.69. Evidently the Euclidean distances of the best and the worst runs of MOALOA-SWMM are 0.13 and 0.07 less (better) than the corresponding values of NSGAII-SWMM, respectively, and the Pareto front resulting from the best run of the MOALOA-SWMM is approximately 1.72 times nearer the coordinate's origin than the best Pareto front calculated with the NSGAII-SWMM. The hyper volume index for objective functions 1 and 2, and objective functions 2 and 3 was calculated considering (0,0) as the base point. The average hyper volume values of objective functions 1 and 2 for 10 runs equal 0.251 and 0.276 respectively for MOALOA and NSGA II. The average hyper volume values of objective functions 2 and 3 equal 0.115 and 0.114 respectively for MOALOA and NSGA II. The results show overall better performance of the MOALOA-SWMM.

Algorithm | Feature | Run 1 | Run 2 | Run 3 | Run 4 | Run 5 | Run 6 | Run 7 | Run 8 | Run 9 | Run 10 |
---|---|---|---|---|---|---|---|---|---|---|---|

MOALOA | Mean Euclidean distance between Pareto front and the origin of the coordinate system | 0.40 | 0.61 | 0.58 | 0.62 | 0.53 | 0.57 | 0.50 | 0.47 | 0.45 | 0.51 |

Number of points on the Pareto front | 100 | 93 | 94 | 92 | 98 | 89 | 92 | 100 | 100 | 86 | |

Run Time (s) | 2,950 | 2,743 | 2,695 | 2,849 | 2,615 | 2,754 | 2,917 | 2,825 | 2,990 | 2,890 | |

Hyper Volume indicator for objective functions (1) and (2) | 0.2 | 0.28 | 0.25 | 0.29 | 0.26 | 0.28 | 0.24 | 0.23 | 0.22 | 0.26 | |

Hyper Volume indicator for objective functions (2) and (3) | 0.1 | 0.13 | 0.13 | 0.12 | 0.11 | 0.11 | 0.12 | 0.1 | 0.11 | 0.12 | |

NSGAII | Mean Euclidean distance between Pareto front and the origin of the coordinate system | 0.57 | 0.53 | 0.58 | 0.69 | 0.58 | 0.61 | 0.59 | 0.60 | 0.60 | 0.60 |

Number of points on the Pareto front | 20 | 18 | 18 | 18 | 18 | 16 | 20 | 18 | 15 | 15 | |

Run Time (s) | 2,985 | 2,864 | 2,689 | 2,762 | 2,889 | 2,850 | 2,828 | 2,735 | 2,963 | 2,912 | |

Hyper Volume indicator for objective functions (1) and (2) | 0.28 | 0.25 | 0.26 | 0.3 | 0.28 | 0.28 | 0.28 | 0.27 | 0.28 | 0.28 | |

Hyper Volume indicator for objective functions (2) and (3) | 0.1 | 0.1 | 0.11 | 0.13 | 0.1 | 0.13 | 0.11 | 0.13 | 0.11 | 0.12 |

Algorithm | Feature | Run 1 | Run 2 | Run 3 | Run 4 | Run 5 | Run 6 | Run 7 | Run 8 | Run 9 | Run 10 |
---|---|---|---|---|---|---|---|---|---|---|---|

MOALOA | Mean Euclidean distance between Pareto front and the origin of the coordinate system | 0.40 | 0.61 | 0.58 | 0.62 | 0.53 | 0.57 | 0.50 | 0.47 | 0.45 | 0.51 |

Number of points on the Pareto front | 100 | 93 | 94 | 92 | 98 | 89 | 92 | 100 | 100 | 86 | |

Run Time (s) | 2,950 | 2,743 | 2,695 | 2,849 | 2,615 | 2,754 | 2,917 | 2,825 | 2,990 | 2,890 | |

Hyper Volume indicator for objective functions (1) and (2) | 0.2 | 0.28 | 0.25 | 0.29 | 0.26 | 0.28 | 0.24 | 0.23 | 0.22 | 0.26 | |

Hyper Volume indicator for objective functions (2) and (3) | 0.1 | 0.13 | 0.13 | 0.12 | 0.11 | 0.11 | 0.12 | 0.1 | 0.11 | 0.12 | |

NSGAII | Mean Euclidean distance between Pareto front and the origin of the coordinate system | 0.57 | 0.53 | 0.58 | 0.69 | 0.58 | 0.61 | 0.59 | 0.60 | 0.60 | 0.60 |

Number of points on the Pareto front | 20 | 18 | 18 | 18 | 18 | 16 | 20 | 18 | 15 | 15 | |

Run Time (s) | 2,985 | 2,864 | 2,689 | 2,762 | 2,889 | 2,850 | 2,828 | 2,735 | 2,963 | 2,912 | |

Hyper Volume indicator for objective functions (1) and (2) | 0.28 | 0.25 | 0.26 | 0.3 | 0.28 | 0.28 | 0.28 | 0.27 | 0.28 | 0.28 | |

Hyper Volume indicator for objective functions (2) and (3) | 0.1 | 0.1 | 0.11 | 0.13 | 0.1 | 0.13 | 0.11 | 0.13 | 0.11 | 0.12 |

Figure 1(b) depicts a two-dimensional view of the Pareto fronts calculated with the MOALOA-SWMM and NSGAII-SWMM models. The value of the third objective function is represented in Figure 1(b) by the diameter of the circles chosen as markers to represent the calculated solutions. It is seen in Figure 1(b) that there are a few points on the MOALOA Pareto that are located above those of the NSGAII Pareto in relation to the *x* and *y* axes, yet the value of the third objective function is smaller (and thus preferable) for the MOALOA solutions.

The results from this SWWM sample problem demonstrate the capacity of the MOALOA-SWMM model when tested against the well-known NSGAII model coupled with the SWMM. The next sections describe the implementation of the MOALOA-SWMM model to the LID SCMs optimization in District 6 of Tehran.

### The design hyetograph, and runoff generation in District 6 of Tehran municipality

The application of Equation (10) and the alternating block method produce the design hyetograph shown in Figure 2. The 25-year return period design hyetograph was developed based on the Tehran Stormwater Management Master Plan (MahabGhods Consulting Engineering Company 2011a, 2011b). The rainfall duration is longer than the concentration time because the largest peak runoff rate occurs when the entire drainage area is contributing (USDA SCS 1986). The SWWM simulated runoff associated with the design hyetograph in the absence of LID SCMs. The input parameters for SWMM such as curve number, urban percent of imperviousness, depression storage, and Manning's roughness coefficient were calibrated by Tehran University, Faculty of Environment (Tehran University-Graduate Faculty of environment 2015). The results show that runoff volume produced by the design rainfall exceeds the conduit capacity in the four nodes. Table 4 lists the duration, discharge, and volume of overflow at four flooding nodes. Figure 3 depicts the locations of flooding nodes within the study catchment. The runoff volume and peak flow at the site outlet equal 581.05 × 10^{3} m^{3} and 71.66 × 10^{3} m^{3}/s, respectively. These runoff volumes and peak flows at the catchment outlet and flooding nodes inside the catchment may induce flooding in the downstream urban area and within the catchment. Therefore, LID SCMs are deployed to improve the performance of the drainage system and decrease the risk of downstream flooding.

Node no. | Duration (h) | Maximum overflow discharge (m^{3}/s) | Total overflow volume (m^{3}10^{3}) |
---|---|---|---|

10 | 0.69 | 5.01 | 6.55 |

16 | 3.41 | 7.57 | 30.64 |

31 | 2.74 | 11.70 | 33.40 |

41 | 0.26 | 1.70 | 1.07 |

Node no. | Duration (h) | Maximum overflow discharge (m^{3}/s) | Total overflow volume (m^{3}10^{3}) |
---|---|---|---|

10 | 0.69 | 5.01 | 6.55 |

16 | 3.41 | 7.57 | 30.64 |

31 | 2.74 | 11.70 | 33.40 |

41 | 0.26 | 1.70 | 1.07 |

### Results from the NSGAII-SWMM and MOALOA-SWMM for LID SCMs optimization in District 6 of Tehran

The three-objective problem of the present study was solved with MOALOA and NSGAII applying 100 iterations and an initial population size of 150 solutions. The number of decision variables equals 90 (=45 subcatchments × 2 possible SCMs per subcatchment). The decision variables represent the areas over which a type of LID SCM is deployed within each subcatchment.

Figure 4(a) illustrates the Pareto fronts calculated with the MOALOA-SWMM and NSGAII-SWMM models in three dimensions. Figure 4(b) depicts the same results in a two-dimensional visualization. The diameters of the circles represent the value of the third objective function in Figure 4(b). The smaller diameters signify better solutions under minimization. It is seen in Figure 4(a) and 4(b) that the MOALOA-SWMM solutions have wider distribution (that is, greater diversity of choices or solutions) and greater density than the NSGAII-SWMM solutions. The MOALOA-SWMM model produces 135 non-dominated Pareto solutions compared with 50 by the NSGAII-SWMM model. Figure 4(b) shows that the NSGAII-SWMM model produced no solutions for *Z*_{1} < 0.1, while MOALOA produced non-dominated solutions for objectives *Z*_{1} and *Z*_{2} for the entire defined feasible range of solutions. This demonstrates the capacity of the MOALOA-SWMM model to more precisely and broadly (diversity) search the multi-objective space of the optimization problem and calculate non-dominated Pareto solutions. It is worth noting that the average Euclidean distance from the origin of the coordinate system for the final Pareto front calculated with the MOALOA-SWMM equals 0.637, and that of the Pareto calculated with the NSGAII-SWMM model equaled 0.764. This indicates the Pareto front produced by MOALOA-SWMM is closer to the ideal solution (located at the origin of the coordinate system under minimization).

Figure 5 depicts the non-dominated Pareto solutions of the three-objective optimization. The numbers on the x-axis and the y-axis represent the values of the first objective function (minimization of the volume of runoff at the outlet of the catchment and at flooding nodes) and the second objective function (minimization of the construction and maintenance cost of LID SCMs), respectively. The points or solutions on the Pareto front are represented by circles with different diameters, whose sizes denote the value of the third objective function (reduction of system performance due to failure of LID SCMs by clogging and blockage). The smaller diameters signify better solutions under minimization.

It is seen in Figure 5 that as the second objective function (which represents cost) increases (worsens) from 0 to 0.5, the first objective function decreases (improves) from 0.9 to about 0.15, and, thenceforth, the decreasing trend of the first objective function continues with a lower slope. In fact, as the second objective function increases from 0.5 to 0.9, the first objective function decreases by only 0.15. This indicates that marginal improvement of the first objective drastically decreases for Z2 > 0.5. Figure 5 shows that there are solutions with different diameters that signify the value of the third objective function (or reduction of LID SCMs performance) on the Pareto front.

The third objective function plays an important role in the selection of LID SCMs. This is so because for close values of the first and second objective function the magnitude of the difference between the values of the third objective function signifies the difference in reduction performance of the LID SCMs due to clogging and blockage. Table 5 lists the values of two solutions on the Pareto fronts with similar first and second objective functions. According to Table 5 for the two solutions numbered 74 and 75, signaled by arrows in Figure 6, the first and second objective functions have approximately equal values, but the third objective function has different values for these solutions, which indicates a difference in the reduction of service performance of the LID SCMs associated with each solution. After blockage and reduction of the performance of LID SCMs the volume of runoff at the outlet of the subcatchments and the volume of runoff at flooding nodes increased by 12.8 and 1.41 thousand cubic meters, respectively, for Solution 74, while they increased by 24.2 and 5.47 thousand cubic meters, respectively, for Solution 75. The runoff at the outlet is twice larger for Solution 75 than that of Solution 74, and the runoff at the flooding nodes is five times larger for Solution 75 than that of Solution 74. This considerable difference in runoff underlines the importance of including the third objective function in selecting SCMs with better service performance, which is called service reliability. Otherwise Solutions 74 and 75 would be nearly indistinguishable from each other, yet, by including objective function 3 (reduction of performance) Solution 74 is clearly superior in the long run.

Solution point no. | Objective function 1 | Objective function 2 | Objective function 3 | Difference between VO_{2} and VO_{1} (10^{3} m^{3}) | Difference between VF_{2} and VF_{1} (10^{3} m^{3}) |
---|---|---|---|---|---|

74 | 0.331 | 0.347 | 0.104 | 12.803 | 1.412 |

75 | 0.332 | 0.340 | 0.163 | 24.289 | 5.470 |

Solution point no. | Objective function 1 | Objective function 2 | Objective function 3 | Difference between VO_{2} and VO_{1} (10^{3} m^{3}) | Difference between VF_{2} and VF_{1} (10^{3} m^{3}) |
---|---|---|---|---|---|

74 | 0.331 | 0.347 | 0.104 | 12.803 | 1.412 |

75 | 0.332 | 0.340 | 0.163 | 24.289 | 5.470 |

The MOALOA-SWMM model was run for District 6 of Tehran employing the first objective function (Equation (2)), reduction of the runoff volume at the catchment outlet and at flooding nodes, and the second objective function (Equation (3), minimization of the SCMs construction and maintenance cost). Figure 6 shows that the Pareto front generated by this two-objective optimization with MOALOA-SWMM has suitable density, coverage, and incorporates all the solutions from the highest to the lowest values of both objective functions, demonstrating the accurate and effective performance of the simulation-optimization model. These Pareto solutions reduce the volume of outlet runoff and the volume of runoff at flooding nodes between 0 and 215.8 × 10^{3} m^{3} and between 0 and 51.243 × 10^{3} m^{3} for a cost of 0 to $47,932. The next section applies the TOPSIS method to select the most desirable solutions employing two-objective and three-objective optimization and compares the impact of considering objective 3 (reduction of LID SCMs performance over time) in the deployment of LID SCMs.

### Determining the most desirable solutions with the TOPSIS method

Figures 5 and 6 show the most desirable solutions selected with the TOPSIS method corresponding to the three-objective and two-objective optimizations, respectively. The two-objective solution is such that the reductions of the volumes of runoff at the outlet and at the flooding nodes equal 163.4 × 10^{3} and 43.24 × 10^{3} m^{3}, respectively, for a cost of $18.54 × 10^{6} of LID SCMs. By considering maximum blockage and the failure percentage of LID practices the amount of volume of runoff at outfall and the amount of runoff at flooding nodes are increased by 27.3 × 10^{3} and 7.2 × 10^{3} m^{3}, respectively, in comparison with the assumed condition in which clogging does not occur. The three-objective optimization solution produced reductions in the volumes of runoff at the outlet and at the flooding nodes equal to 121 × 10^{3} and 32.69 × 10^{3} m^{3}, respectively, with a cost of $14 × 10^{6}. By considering the maximum percentage of blockage and failure of LID SCMs practices runoff at the outlet and at flooding nodes increases by 12.8 × 10^{3} and 1.41 × 10^{3} m^{3} respectively in comparison with the assumed condition in which clogging does not occur. This considerable reduction of runoff between the cases in which clogging of LID SCMs is assumed not to occur and that in which it does occur demonstrates the necessity of taking account of the third objective function, which improves the selection of LID SCMs by leading to greater service reliability.

TOPSIS method's most desirable solution associated with the three-objective optimization selected deployment of porous pavement in 34% of the maximum area that could be allocated to this LID SCM, which equals 808,135 m^{2}, and allocated bioretention cells over 35% of the maximum area that could be allocated to this LID SCM, which equals 300,878 m^{2}. Bioretention cells are costlier than porous pavement, yet, previous studies have demonstrated that they are less vulnerable to blockage and clogging (Li *et al.* 2009).

Table 6 lists the maximum allowable areas' percentages for implementing LID SCMs, and their optimal percentages of the total areas of subcatchments. The maximum possible areas for implementing LID SCMs is determined based on land use maps of the study area. The optimal areas for LID SCMs in each subcatchment are the result of three-objective optimization, which are shown as the percentage of the total subcatchment areas as listed in Table 6 and Figure 7.

Subcatchment | ||||
---|---|---|---|---|

Porous pavement | Bioretention cells | Porous pavement | Bioretention cells | |

1 | 11 | 5 | 5.7 | 1.6 |

2 | 11 | 6 | 4.9 | 1.3 |

3 | 10 | 5 | 3.4 | 1.7 |

4 | 11 | 5 | 0.5 | 1.9 |

5 | 12 | 4 | 4.3 | 1.8 |

6 | 9 | 5 | 2.7 | 1.9 |

7 | 11 | 5 | 2.6 | 2.1 |

8 | 12 | 4 | 6.1 | 2.3 |

9 | 11 | 5 | 1.6 | 2.0 |

10 | 8 | 6 | 2.9 | 1.9 |

11 | 12 | 4 | 4.7 | 1.9 |

12 | 10 | 6 | 3.0 | 1.7 |

13 | 11 | 5 | 3.0 | 0.4 |

14 | 14 | 5 | 6.4 | 2.5 |

15 | 12 | 5 | 4.5 | 1.0 |

16 | 10 | 5 | 6.0 | 1.4 |

17 | 11 | 5 | 3.7 | 1.7 |

18 | 10 | 7 | 7.6 | 0.5 |

19 | 10 | 7 | 5.0 | 1.6 |

20 | 12 | 6 | 3.7 | 0.3 |

21 | 10 | 5 | 3.2 | 1.2 |

22 | 15 | 7 | 3.0 | 0.7 |

23 | 11 | 6 | 2.0 | 1.4 |

24 | 11 | 5 | 2.3 | 2.1 |

25 | 11 | 5 | 1.0 | 2.2 |

26 | 11 | 5 | 4.9 | 0.9 |

27 | 9 | 5 | 3.8 | 1.8 |

28 | 8 | 4 | 1.2 | 2.0 |

29 | 12 | 5 | 4.6 | 0.3 |

30 | 11 | 5 | 5.5 | 0.1 |

31 | 11 | 6 | 0.1 | 1.0 |

32 | 12 | 3 | 1.6 | 0.3 |

33 | 13 | 7 | 3.2 | 0.7 |

34 | 10 | 6 | 2.9 | 1.8 |

35 | 11 | 7 | 4.6 | 1.6 |

36 | 11 | 6 | 7.3 | 2.3 |

37 | 11 | 6 | 0.3 | 1.8 |

38 | 12 | 5 | 6.0 | 1.7 |

39 | 11 | 5 | 5.9 | 2.1 |

40 | 10 | 6 | 5.1 | 1.5 |

41 | 10 | 8 | 1.7 | 1.6 |

42 | 10 | 6 | 5.4 | 0.3 |

43 | 11 | 5 | 4.2 | 3.0 |

44 | 10 | 5 | 3.4 | 1.4 |

45 | 11 | 5 | 4.0 | 1.6 |

Subcatchment | ||||
---|---|---|---|---|

Porous pavement | Bioretention cells | Porous pavement | Bioretention cells | |

1 | 11 | 5 | 5.7 | 1.6 |

2 | 11 | 6 | 4.9 | 1.3 |

3 | 10 | 5 | 3.4 | 1.7 |

4 | 11 | 5 | 0.5 | 1.9 |

5 | 12 | 4 | 4.3 | 1.8 |

6 | 9 | 5 | 2.7 | 1.9 |

7 | 11 | 5 | 2.6 | 2.1 |

8 | 12 | 4 | 6.1 | 2.3 |

9 | 11 | 5 | 1.6 | 2.0 |

10 | 8 | 6 | 2.9 | 1.9 |

11 | 12 | 4 | 4.7 | 1.9 |

12 | 10 | 6 | 3.0 | 1.7 |

13 | 11 | 5 | 3.0 | 0.4 |

14 | 14 | 5 | 6.4 | 2.5 |

15 | 12 | 5 | 4.5 | 1.0 |

16 | 10 | 5 | 6.0 | 1.4 |

17 | 11 | 5 | 3.7 | 1.7 |

18 | 10 | 7 | 7.6 | 0.5 |

19 | 10 | 7 | 5.0 | 1.6 |

20 | 12 | 6 | 3.7 | 0.3 |

21 | 10 | 5 | 3.2 | 1.2 |

22 | 15 | 7 | 3.0 | 0.7 |

23 | 11 | 6 | 2.0 | 1.4 |

24 | 11 | 5 | 2.3 | 2.1 |

25 | 11 | 5 | 1.0 | 2.2 |

26 | 11 | 5 | 4.9 | 0.9 |

27 | 9 | 5 | 3.8 | 1.8 |

28 | 8 | 4 | 1.2 | 2.0 |

29 | 12 | 5 | 4.6 | 0.3 |

30 | 11 | 5 | 5.5 | 0.1 |

31 | 11 | 6 | 0.1 | 1.0 |

32 | 12 | 3 | 1.6 | 0.3 |

33 | 13 | 7 | 3.2 | 0.7 |

34 | 10 | 6 | 2.9 | 1.8 |

35 | 11 | 7 | 4.6 | 1.6 |

36 | 11 | 6 | 7.3 | 2.3 |

37 | 11 | 6 | 0.3 | 1.8 |

38 | 12 | 5 | 6.0 | 1.7 |

39 | 11 | 5 | 5.9 | 2.1 |

40 | 10 | 6 | 5.1 | 1.5 |

41 | 10 | 8 | 1.7 | 1.6 |

42 | 10 | 6 | 5.4 | 0.3 |

43 | 11 | 5 | 4.2 | 3.0 |

44 | 10 | 5 | 3.4 | 1.4 |

45 | 11 | 5 | 4.0 | 1.6 |

Figure 7 depicts the spatial distribution of the LID SCMs associated with the most desirable solution obtained with three-objective optimization (calculated percentage area of LID SCMs to their maximum allowable percentage area). It can be seen in Figure 7 the upstream subcatchment encompassing node 31 has a low-percentage area of LID SCMs because the cause of flooding at node 31 is runoff entering as inflow from an upstream catchment. Therefore, the deployment of LID SCMs in subcatchment of node 31 would not be effective. In contrast, the subcatchments with nodes 10 and 16 have large areas covered by LID SCMs. Evidently, the selected solution provides optimal LID SCMs based on their capacity to reduce runoff, their costs, and their service reliability, features that have a strong geographical context given the heterogeneity of land-use conditions. Generally, subcatchments in which runoff passes through the flooding nodes are more suitable for allocating LID SCMs and produce more effective designs.

Table 7 reports the runoff volume at flooding nodes before and after applying the LID SCMs. These values demonstrate the runoff volume at flooding nodes is reduced considerably. Yet, at node 31 the runoff volume is only slightly reduced because it enters as inflow from upstream areas. Table 8 lists the peak discharge, the runoff volume, and the time to peak discharge at the catchment outlet before and after LID SCMs deployment. It is evident from Table 8 that LID SCMs reduce the runoff volume by 121,000 m^{3} and delays the time to peak discharge by 20 minutes. Figure 8 depicts the outflow hydrograph before and after applying LID SCM. It is evident from Figure 8 the reduction of the peak discharge and the delay in the time to peak discharge.

Node no. | Overflow volume before applying LID SCMs (10^{3} m^{3}) | Overflow volume after applying LID SCMs (10^{3} m^{3}) | Reduced runoff volume (10^{3} m^{3}) | Fraction of reduced runoff volume (%) |
---|---|---|---|---|

10 | 6.55 | 1.38 | 5.17 | 78 |

16 | 30.64 | 16.445 | 14.20 | 46 |

31 | 33.40 | 21.15 | 12.25 | 37 |

41 | 1.07 | 0 | 1.07 | 100 |

Node no. | Overflow volume before applying LID SCMs (10^{3} m^{3}) | Overflow volume after applying LID SCMs (10^{3} m^{3}) | Reduced runoff volume (10^{3} m^{3}) | Fraction of reduced runoff volume (%) |
---|---|---|---|---|

10 | 6.55 | 1.38 | 5.17 | 78 |

16 | 30.64 | 16.445 | 14.20 | 46 |

31 | 33.40 | 21.15 | 12.25 | 37 |

41 | 1.07 | 0 | 1.07 | 100 |

Status | Peak discharge at outlet (10^{3} m^{3}/s) | Runoff volume at outfall (10^{3} m^{3}/s) | Time to peak flow (min) |
---|---|---|---|

Before applying LID SCMs | 61.21 | 581.05 | 200 |

After applying LID SCMs | 47.70 | 460.06 | 220 |

Status | Peak discharge at outlet (10^{3} m^{3}/s) | Runoff volume at outfall (10^{3} m^{3}/s) | Time to peak flow (min) |
---|---|---|---|

Before applying LID SCMs | 61.21 | 581.05 | 200 |

After applying LID SCMs | 47.70 | 460.06 | 220 |

## CONCLUSION AND DISCUSSION

The transition of native landscapes to urban areas increase the impervious surfaces. These changes reduced infiltration, which leads to increased runoff and flood risks in urban areas. LID SCMs are deployed to decrease the stress on urban stormwater infrastructure, improve their performance, purifying urban runoff, and for urban embellishment. In this study, a comprehensive approach has been introduced to determining the optimal areas of deployment of LID SCMs which relies on runoff simulation using SWMM coupled to multi-objective optimization to produce Pareto fronts of alternative solutions by considering three objectives and ranking of solutions with the TOPSIS method.

Objective 1: Reduction of stormwater overflows volume at flooding nodes inside the catchment and runoff volume at the catchment outlet to prevent flooding inside and downstream of the catchment respectively. Taken into account the reduction of runoff volume at the catchment outlet for cases where vulnerable urban areas exist downstream the catchment is of particular significance.

Objective 2: Cost minimization, which results in the most affordable alternative for the best performance of LID SCMs and the largest reduction in runoff volume, which is a key objective in most stormwater optimization problems due to financial limitations.

Objective 3: Minimize the service-performance reduction of LID SCMs considering this objective function improves their selection, assuring with longer-term performance. The third objective function improves the long term performance of designed alternatives is herein demonstrated by (1) comparing the results of the two-objective optimization without considering this objective function with the three-objective optimization results and (2) by comparing the results of the three-objective optimizations of two case studies in which the values of their first and second objective functions were approximately equal, while their third objective function value was significantly different, thus demonstrating the marginal value added by addition of the third objective.

The combination of porous pavement and bioretention cell as herein implemented to improve the performance of Tehran's District 6 stormwater infrastructure whose results indicate the ability of these two LID SCMs to control runoff volume.

The MOALOA solved the three-objective optimization, and its results were compared with those of well-established NSGAII using 10 runs. The results show the superiority of MOALOA. Thus, the MOALOA was shown to be an effective algorithm in conjunction with SWMM for solving complex multi-objective optimization problems involving LID SCMs deployment.

This study's approach was tested and successfully implemented with sample problems from the SWMM manual and with a catchment encompassing District 6 of Tehran municipality (Iran). The proposed simulation-optimization approach can assist stormwater experts the world over in choosing optimal combinations of LID SCMs with confidence to reduce urban flooding and improve water quality.

## ACKNOWLEDGEMENTS

The authors thank Iran's National Science Foundation (INSF) for its financial support of this research. Special thanks to Dr Abdolreza Karbasi for sharing data and information for District 6 of Tehran municipality.

## CONFLICTS OF INTEREST

None.