Evaluating the Effects of Low Impact Development Practices on Urban Flooding under Different Rainfall Intensities

Low impact development (LID) is an important control measure against extreme rainfall events and is widely applied to relieve urban flood disasters. To investigate the effects of LID practices on flooding control under different rainfall scenarios, this paper constructs a rainfall–runoff model based on the storm water management model (SWMM) for a typical residential area in Guangzhou, China. The model is calibrated by using observed rainfall and runoff data. A total of 27 rainfall scenarios are constructed to simulate the change characteristics before and after the LID practices. Also, the projection pursuit method based on a particle swarm optimization (PSO) algorithm is used to assess the flooding characteristics. The results show that the constructed rainfall–runoff model can closely reflect the relationship between rainfall and runoff, with all Nash–Sutcliffe coefficients of efficiency (NSE) exceeding 0.7. It was found from the simulation and assessment of the constructed rainfall scenarios that the changes in rainfall characteristics have a considerable impact on the constructed drainage system and that LID practices can properly control the floods. However, with an increase in rainfall peak coefficient, intensity or duration, the control effects of LID tend to reduce. Particularly in the scenario of relatively high rainfall intensity, the impact of rainfall duration and the rainfall peak coefficient on the LID practices is minor.


Introduction
Along with significant changes in global climate and the considerably increasing disturbance caused by human activities, flood disasters in urban areas have become more serious, substantially restricting social and economic development and posing a threat to the safety of human lives and properties [1,2]. In recent years, urban flooding has caused many casualties and property losses in Japan, Singapore, The Netherlands, Britain and other countries [3][4][5][6]. Many major cities in China have also been greatly influenced by urban flooding disasters, including Guangzhou, Wuhan, Shenzhen, Nanjing, Beijing [7]. The key issue is that the drainage systems of these areas cannot adapt to climate change and human activity. Therefore, there is an urgent need to control flood disasters from a design perspective.
In order to alleviate the problem of urban flooding disasters, many strategies of flood mitigation and disaster relief have been put forward [8][9][10]. The traditional efforts include the following: Table 1. Comparison of some Low Impact Development (LID) practices [19]. Compared with traditional flood control measures, LID practices are more flexible in response to climate change, which is a very important aspect for government decision-makers [36]. Also, it is considered to be a pivotal approach to sustainable urban development. Nevertheless, most studies consider only laboratory scenarios or one particular rainfall scenario. Research into the response characteristics of LID under different rainfall scenarios is sparse. In addition, the evaluation of LID practices is mainly based on the evaluation of a single factor (such as increasing peak flow attenuation, increasing lag time, flow attenuation of stormwater runoff, etc.), and it has mainly evaluated at the multidimensional scale, which is likely to cause dimensionality when a multi-factor evaluation is performed.

LID
A large number of studies have shown that the mastering of LID control characteristics and accurately evaluating the control of LID practices under different rainfall scenarios can effectively respond to urban flood disasters and guide the formulation of appropriate measures. Therefore, a rainfall-runoff model is constructed based on the SWMM model [27]. When a large number of scenarios is considered, a better optimization procedure is required. In this paper, 27 rainfall scenarios with different rainfall intensities, rainfall durations and rainfall peak coefficients are considered and a comprehensive assessment of the change characteristics of floods is conducted in various rainfall scenarios by using particle swarm optimization (PSO)-based projection pursuit technology [37] which includes (1) the discovery of the response law of the pipeline network to different rainfall intensity before and after the implementation of LID practices; (2) the qualitative and quantitative evaluation of the variation characteristics of the flood before and after the implementation of the LID practices; and (3) the recognition of the control characteristics of LID under different complex scenarios.

Study Area
Guangzhou (22 • 26 -23 • 56 N, 112 • 57 -114 • 03 E) has a tropical and subtropical monsoon climate with extremely significant characteristics of a monsoon climate and subtropical westerly, equatorial westerly and tropical easterly winds. The area is abundant in rainfall, with average annual rainfall of 1675.5 mm, and characterized by a long rainy season and heavy rainfall of great intensity and nonuniform spatial distribution. The rainfall is mainly concentrated in the months from April to September that are affected by warm air. The rainfall from April to September accounts for 80% of annual rainfall. The rainfall in May and June represents 32% of the annual rainfall. Since the reform and opening policies in the 1980s, the forms of land use in Guangzhou have changed substantially.
Land under construction is constantly being expanded, contributing to changes in the characteristics of runoff and to increasingly serious urban flooding.
In order to probe the flood control characteristics of LID practices under different rainfall scenarios, this paper selected a typical highly-developed residential area. It covers a total of 3.13 × 10 5 m 2 and is primarily used for construction with some land set aside for green space and transportation ( Figure 1). The drainage system uses a rainwater and sewage separation system, and its design standard is based on a 2-year rainfall intensity developed by the Guangzhou Water Affairs Bureau in 1993 (with pipe diameter between 900 and 2000 mm and slopes between 0.001 and 0.01).
The study area is divided into 19 sub-catchments, and the drainage system is generalized into 31 drainage pipes, 31 manholes and 1 drainage outlet ( Figure 1) in accordance with land use, the pipe laying situation, the terrain and other data. The rainfall-runoff model is constructed based on the SWMM model to reflect the relationships between runoff and rainfall. Point A (in the downstream area) is selected to record the runoff, and the flow and rainfall are recorded simultaneously in accordance with a flow meter and rain gauge at an interval of 10 min (Figure 1). westerly and tropical easterly winds. The area is abundant in rainfall, with average annual rainfall of 1675.5 mm, and characterized by a long rainy season and heavy rainfall of great intensity and nonuniform spatial distribution. The rainfall is mainly concentrated in the months from April to September that are affected by warm air. The rainfall from April to September accounts for 80% of annual rainfall. The rainfall in May and June represents 32% of the annual rainfall. Since the reform and opening policies in the 1980s, the forms of land use in Guangzhou have changed substantially. Land under construction is constantly being expanded, contributing to changes in the characteristics of runoff and to increasingly serious urban flooding.
In order to probe the flood control characteristics of LID practices under different rainfall scenarios, this paper selected a typical highly-developed residential area. It covers a total of 3.13 × 10 5 m 2 and is primarily used for construction with some land set aside for green space and transportation ( Figure 1). The drainage system uses a rainwater and sewage separation system, and its design standard is based on a 2-year rainfall intensity developed by the Guangzhou Water Affairs Bureau in 1993 (with pipe diameter between 900 and 2000 mm and slopes between 0.001 and 0.01).
The study area is divided into 19 sub-catchments, and the drainage system is generalized into 31 drainage pipes, 31 manholes and 1 drainage outlet ( Figure 1) in accordance with land use, the pipe laying situation, the terrain and other data. The rainfall-runoff model is constructed based on the SWMM model to reflect the relationships between runoff and rainfall. Point A (in the downstream area) is selected to record the runoff, and the flow and rainfall are recorded simultaneously in accordance with a flow meter and rain gauge at an interval of 10 min (Figure 1).

Methodology
In order to investigate the control characteristics of the highly-developed area under different rainfall scenarios after the implementation of LID practices, the observed rainfall-runoff data was employed to calibrate and verify the rainfall-runoff model based on the SWMM model. Also, different rainfall scenarios were constructed based on the rainfall intensity formula in Guangzhou [38] and the Chicago Hydrograph Model [39]. The characteristics of flood change (before and after LID) were analyzed by PSO-based projection pursuit technology.

SWMM Model
The SWMM model is a relatively flexible hydrodynamic model developed by USEPA [27], based on the physical mechanisms of mass, energy and momentum conversion, as well as a comprehensive mathematical model constructed for the mechanisms of runoff generation [27]. It can properly simulate the rainfall-runoff relationship, pollutant dispersion, flood control and other hydrological processes in urban areas under natural conditions. Also, the SWMM model is widely used to simulate the hydrologic performance of natural channels, rain and sewage diversion systems or other drainage systems, and it can evaluate the design of LID practices and BMPs. Much research on the SWMM model can be found [40][41][42].
The SWMM model provides different modules and algorithms for simulation as described below.

Surface Runoff Simulation Module
When the SWMM model is included in the calculation of surface production, the sub-basin is regarded as a non-linear reservoir, and the surface runoff refers to the rainfall depth after infiltration and surface water storage. The surface water storage modules are different for different regions. The SWMM model provides the Horton model, the Green-Ampt model and the curve number model to calculate the infiltration. This paper chooses the Horton model for calculating small watersheds. The flow is transformed into the outflow process and the process in which flow changes over time, and the two processes are calculated separately using Manning's equation and the finite difference method.

Transmission Calculation Module
The main method for calculating the flow of the pipe network is the Saint-Venant equation. In the SWMM model, three methods (steady flow routing, kinematic wave routing and dynamic wave routing) are used for pipe network transmission calculation. The kinematic wave routing method is used in this paper.

Flood Calculation Module
If the depth of the node exceeds the maximum available depth, overflow is lost from the system or water is saved at the top of the node, and it re-enters the drainage system under specific conditions. The SWMM model can simulate the approximate location of the overflow, the distribution and the corresponding overflow.

LID Module
The SWMM model can simulate the control effect of LID practices, and seven different types are provided: bio-retention cell, porous pavement, infiltration trench, rain barrel, vegetative swale, rain garden and green roof. In the SWMM model, different types of LID practices are divided into several vertical layers and represented by five to 23 parameters (e.g., thickness, surface roughness, offset height). More details about this can be found in [27].
The SWMM model's user manual [27] provides a range of values for sensitive parameters based on a large number of simulation studies; this provides references for calibrating the model's parameters. Thus, the recommended range of the parameters and the trial and error approach are used.

Evaluation of Model Accuracy
The Nash-Sutcliffe Efficiency (NSE) coefficient, which is often considered to be one of the most important indices for measuring the model's simulation accuracy [43], was used to evaluate the accuracy of the constructed model simulation results. Its expression is as follows: where n means the total number of time-steps; S i represents the simulated value at time-step i; O i represents the observed flow at time-step i; and O is the average value of the observed flow. When NSE = 1, the observed flow is consistent with the simulated flow. If NSE is within the range of 0 to 1, the simulated flow is acceptable. If NSE is greater than 0.5, the simulation results are satisfactory [44,45].

Projection Pursuit Method
The projection pursuit method is mainly used to project high-dimensional data to low dimensions, to avoid "dimension disaster". The characteristics of the original data are described by the optimal projection direction and projection values [46]. It has been widely applied to water resource assessment, engineering stability assessments and many other applications [47][48][49]. Compared with traditional methods (e.g., neural networks and principal component analysis), this method can eliminate the influence of subjective factors.
Since the flooding of an urban drainage system involves a large number of indices, it is difficult to comprehensively assess all of the flood characteristics through high-dimensional indices and data. According to the mechanisms and objectives of urban flooding analysis, the volume of flooding, the maximum flooding rate and the duration of flooding are selected as indices used to assess the characteristics of the flooding. A high-dimensional flood evaluation system may be reduced to a one-dimensional system in order to comprehensively evaluate the regional flood characteristics using the projection pursuit method. The calculation steps are as follows: Step (1): Normalize the evaluation index. As the evaluation indices have different dimensions and lack comparability, they must be normalized. Assume an evaluation index set x * = x ij i = 1, 2, . . . , n; j = 1, 2, 3 , where x ij is the j th index value of the i th sample, and n is the sample size. A large number of studies have shown that the selected indices are positively related to flood evaluation. Therefore, all indices are used in Equation (2) to normalize x ij to X ij .
where x jmax and x jmin are the maximum and the minimum of evaluation index j th of all samples respectively.
Step (2): Find the function of the construct projection index. The flooding evaluation index is expressed as a one-dimensional projection value Therefore, the construction of the projection index function Q(a) is as follows: where S Z and D Z indicate the standard deviation and local density of Z i , respectively, i.e., where R is radius of the window, and it is usually 0.1 ; P is the number of indices; r(i, j) is the distance between samples; and u(R − r(i, j)) is the unit step function.
The change in the projection index function Q(a) is determined by the projection direction. When Q(a) is maximized, it is considered to be the most likely function to represent the structural features of high dimensional data, and the corresponding vector is the best projection direction α * , where α * is the projection direction vector involving the maximum number of high-dimensional data characteristics. The maximization of the objective function and its corresponding constraints are denoted by (6) and (7), respectively.
Step (3): Solving the objective function. The objective function (6) is a complex nonlinear system which needs to be solved by advanced optimization algorithms. A large number of studies have shown that the PSO algorithm has significant versatility in high-dimensional global optimization with simple and effective operation [50][51][52][53]. It is therefore used to solve the objective function in this study. The PSO algorithm is a heuristic search method that treats a solution to a problem as a "particle" of the search space. The basic idea is to initialize a group of random particles (i.e., random solutions), and then find the optimal solution by iterative calculations. The velocity, position and fitness of the particle N), N is the pre-defined number of particles. During each iteration, the optimal solution (individual optimal value p d (t)) will be found by the particle through the particles themselves, and the optimal solution (global optimal value g d (t)) found among the entire population of particles during each calculation updates the position and velocity of the particle. Its update rules are as follows: where, t is the current iteration as inertia weight; c 1 and c 2 represent individual learning factor and group learning factor, respectively; and b 1 (t) and b 2 (t) represent the individuals and groups of particles, respectively. In this paper, 150 random particles are initialized, and the maximum number of iterations is set to 1500. The maximum and minimum weighting factors are 0.9 and 0.4, respectively. c 1 = c 2 = 2 gave the best overall performance [54,55]. Therefore, c 1 and c 2 are kept as constant of 2.
In the iterative calculation process, the update rules of the individual optimal value and the global optimal value are shown in (10) and (11).
where f d (t + 1) is the fitness value of the iteration of particle d in t + 1.
To solve the coupling between PSO and the projection pursuit method, we use the following procedure: Step a j x ij (i = 1, 2, . . . , n) to calculate the value of Z i Step (2) Use Equation (4) and (5) for calculating S z and D z Step (3)  When f d (t + 1) = f d (t), or it achieves the maximum number of iterations, the resulting g d (t) is equal to α * .
Step (4) Analyze the flood characteristics. α * will be substituted into the formula Z i = p ∑ j=1 a j x ij (i = 1, 2, . . . , n) to obtain the best projection value for each manhole Z i * ; the greater the projection value is, the larger the flood.

Setting of Rainfall Scenarios
Rainfall design is the basis and premise for planning regional drainage systems and flood control measures. To probe the impact of climate change on rainfall, the Guangzhou Water Affairs Bureau analyzed serial information about rainfall from the most recent 60 years in 2011 [38], compared the information with the Guangzhou rainstorm formula enacted in 1993 [56] and introduced a rainfall intensity formula more in line with the rainstorm characteristics of Guangzhou (Table 2). A total of 27 designed rainfall events, which are composed of different return periods (one year, five years and ten years), rainfall duration (1 h, 1.5 h and 2 h) and rainfall peak coefficients (0.375, 0.5 and 0.8), were selected to investigate the control effects of LID practices. Here, the rainfall peak coefficient means the time that the rainfall peak occurred divided by the rainfall duration. Note: i is design rainfall intensity, t (hour) is rainfall duration and P is rainfall return period.

Calibration and Verification of Model Parameters
According to the observed rainfall-runoff data, the principles of the SWMM model and the recommended parameters, five rainfalls were selected to calibrate the parameters. Based on the sensitivity analysis of the SWMM model parameters, the trial and error approach was used to calibrate the parameter values repeatedly until the simulation results agreed with the observation results. The determined model's main sensitivity parameters are shown in Table 3. The simulation runoff process was obtained by repeatedly calibrating the model's parameters, and the change process of the observed runoff was consistent with the simulated runoff ( Figure 2 and Table 4). The correlation coefficient (R 2 ) between the simulated flow and the observed flow and NSE coefficient exceed 0.91, and the peak relative error remained within 5%. From the rainfall-runoff process, the constructed regional rainfall-runoff model accurately reflects the relationship between runoff and rainfall.
To verify the rationality and reliability of the calibrated parameters, five different rainfall scenarios were selected for parameter verification ( Table 4). The results show that the NSE coefficients between the simulated runoff and observed runoff exceed 0.7, that R 2 is greater than 0.86, and that the discharge relative error and the flood peak relative error are less than 5%. Therefore, the constructed rainfall-runoff model is believed to be highly accurate. Water 2017, 9, 548 9 of 17 and the peak relative error remained within 5%. From the rainfall-runoff process, the constructed regional rainfall-runoff model accurately reflects the relationship between runoff and rainfall. To verify the rationality and reliability of the calibrated parameters, five different rainfall scenarios were selected for parameter verification ( Table 4). The results show that the NSE coefficients between the simulated runoff and observed runoff exceed 0.7, that 2 R is greater than 0.86, and that the discharge relative error and the flood peak relative error are less than 5%. Therefore, the constructed rainfall-runoff model is believed to be highly accurate.

Analyze the Existing Situation of the Drainage System
To adapt to the changes of rainfall characteristics and land use, the constructed rainfall-runoff model simulates the functioning of the pipe network in different scenarios [38]. As can be seen from Table 5, there are seven to ten overload pipes in the two-year rainfall scenario, and the overload time is more than one hour ( Table 5). The drainage system cannot satisfy the designed standard. With the increased rainfall duration or the rainfall peak coefficient, the drainage system is under increasing pressure. To ensure safety, appropriate flood controls must be implemented. Note: In the Rainfall Scenarios row, y indicates rainfall return period, h indicates rainfall duration; and the number below y and h indicates rainfall peak coefficient.

Control Characteristics of LID Practices under Complex Scenarios
In order to ensure safety and resource conservation, it was evaluated how well LID can contribute to flood control and it's characteristics under complex scenarios were analyzed. Since rain gardens and bioretention swales have better control ability and visual effects than other LID practices, the control characteristics of these two LID practices are analyzed under complex

Analyze the Existing Situation of the Drainage System
To adapt to the changes of rainfall characteristics and land use, the constructed rainfall-runoff model simulates the functioning of the pipe network in different scenarios [38]. As can be seen from Table 5, there are seven to ten overload pipes in the two-year rainfall scenario, and the overload time is more than one hour ( Table 5). The drainage system cannot satisfy the designed standard. With the increased rainfall duration or the rainfall peak coefficient, the drainage system is under increasing pressure. To ensure safety, appropriate flood controls must be implemented. Note: In the Rainfall Scenarios row, y indicates rainfall return period, h indicates rainfall duration; and the number below y and h indicates rainfall peak coefficient.

Control Characteristics of LID Practices under Complex Scenarios
In order to ensure safety and resource conservation, it was evaluated how well LID can contribute to flood control and it's characteristics under complex scenarios were analyzed. Since rain gardens and bioretention swales have better control ability and visual effects than other LID practices, the control characteristics of these two LID practices are analyzed under complex scenarios. The calculation is based on the volume of runoff produced by the rainfall as measured once in a two-year period [38]. The layout position is as shown in Figure 1.

Influence of Rainfall Intensity on Flooding and LID Practices
To investigate the comprehensive characteristics of floods with changes in rainfall intensity before and after the implementation of LID practices, the comprehensive characteristics of flood distribution under different rainfall scenarios (2y-2h-0.375, 5y-2h-0.375, 10y-2h-0.375) were analyzed. As can be seen from Figure 3, before LID practices were taken, the drainage system failed to drain the stormwater runoff produced by the rainfall intensity once in a two-year period and above. As the rainfall intensity increases, the projection values also rise. Figure 3 also illustrates that the flood-affected area increases and the flood become more serious with increasing rainfall intensity. After LID practices are implemented, the drainage system was able to promptly drain the runoff produced by the 2y-2h-0.375 rainfall scenario, and it was rendered free from flooding. However, with the increase in the intensity of the rainfall, the stormwater runoff generated by the 5y-2h-0.375 rainfall scenario caused relatively serious flooding in J9 and flooding in J15. As rainfall intensity increased to once in ten years, the flood-affected area and flood manholes increased. However, LID practices can alleviate the drainage pressure on the pipe network. It is effective in controlling the flood-affected area and the number of flood manholes, especially under the low intensity rainfall scenario.
Water 2017, 9, 548 11 of 17 scenarios. The calculation is based on the volume of runoff produced by the rainfall as measured once in a two-year period [38]. The layout position is as shown in Figure 1.

Influence of Rainfall Intensity on Flooding and LID Practices
To investigate the comprehensive characteristics of floods with changes in rainfall intensity before and after the implementation of LID practices, the comprehensive characteristics of flood distribution under different rainfall scenarios (2y-2h-0.375, 5y-2h-0.375, 10y-2h-0.375) were analyzed. As can be seen from Figure 3, before LID practices were taken, the drainage system failed to drain the stormwater runoff produced by the rainfall intensity once in a two-year period and above. As the rainfall intensity increases, the projection values also rise. Figure 3 also illustrates that the flood-affected area increases and the flood become more serious with increasing rainfall intensity. After LID practices are implemented, the drainage system was able to promptly drain the runoff produced by the 2y-2h-0.375 rainfall scenario, and it was rendered free from flooding. However, with the increase in the intensity of the rainfall, the stormwater runoff generated by the 5y-2h-0.375 rainfall scenario caused relatively serious flooding in J9 and flooding in J15. As rainfall intensity increased to once in ten years, the flood-affected area and flood manholes increased. However, LID practices can alleviate the drainage pressure on the pipe network. It is effective in controlling the flood-affected area and the number of flood manholes, especially under the low intensity rainfall scenario.

Influence of Rainfall Duration on Flooding and LID Practices
In order to analyze the comprehensive characteristics of the influence of rainfall duration on flooding and LID practices, three rainfall scenarios (5y-1h-0.375, 5y-1.5h-0.375, 5y-2h-0.375) are selected. Before the implementation of LID practices, the flood-affected area and manholes with high projection values change slightly with the increase in rainfall duration (Figure 4). However, flood-affected areas tend to increase, indicating that long-duration rainfall increases flooding. After the implementation of LID practices, the flooding is significantly controlled, and the flood-affected area is drastically reduced. The shorter the rainfall duration, the smaller the flood-affected area. Yet, as rainfall duration increases, the flood-affected range gradually increases.

Influence of Rainfall Duration on Flooding and LID Practices
In order to analyze the comprehensive characteristics of the influence of rainfall duration on flooding and LID practices, three rainfall scenarios (5y-1h-0.375, 5y-1.5h-0.375, 5y-2h-0.375) are selected. Before the implementation of LID practices, the flood-affected area and manholes with high projection values change slightly with the increase in rainfall duration (Figure 4). However, flood-affected areas tend to increase, indicating that long-duration rainfall increases flooding. After the implementation of LID practices, the flooding is significantly controlled, and the flood-affected area is drastically reduced. The shorter the rainfall duration, the smaller the flood-affected area. Yet, as rainfall duration increases, the flood-affected range gradually increases.

Influence of Rainfall Peak Coefficient on Flooding and LID Practices
The rainfall peak coefficient is one of the most important factors which influences flooding. Three rainfall scenarios with the same rainfall intensity and duration but different rainfall peak coefficients (5y-2h-0.375, 5y-2h-0.5, 5y-2h-0.8) were selected to analyze the influence of the rainfall peak coefficient on flooding and LID practices. As can be seen from Figure 5, when no LID control practices were taken, the impact of flooding is small with the changes in rainfall peak coefficient. The flood-affected area and the projection values are virtually unchanged. However, the rainfall peak coefficient has some influence on LID practices. The control effects of LID practices are diminished with the increase in the rainfall coefficient ( Figure 5). The flood-affected area and the projection

Influence of Rainfall Peak Coefficient on Flooding and LID Practices
The rainfall peak coefficient is one of the most important factors which influences flooding. Three rainfall scenarios with the same rainfall intensity and duration but different rainfall peak coefficients (5y-2h-0.375, 5y-2h-0.5, 5y-2h-0.8) were selected to analyze the influence of the rainfall peak coefficient on flooding and LID practices. As can be seen from Figure 5, when no LID control practices were taken, the impact of flooding is small with the changes in rainfall peak coefficient. The flood-affected area and the projection values are virtually unchanged. However, the rainfall peak coefficient has some influence on LID practices. The control effects of LID practices are diminished with the increase in the rainfall coefficient ( Figure 5). The flood-affected area and the projection values of manholes J9 and J15 tended to increase ( Figure 5). They are the largest under the 5y-2h-0.8 rainfall scenario.

Changes in the Comprehensive Flood and LID Control Characteristics under Different Rainfall Scenarios
According to the comprehensive distribution characteristics of flooding, the projection values of the manholes J9 and J15 and the downstream area are high. The reasons for this are as follows: (i) for J15, it is located in the joint point of two pipes, and therefore its convergence area is large; (ii) J9 covers a large proportion of impervious surface in this area, and therefore rapid convergence speed and the severe overflow of the downstream pipe network contribute to grave flooding. The downstream flooding of the drainage system is caused by the failure of the timely drainage of stormwater, the (overly) full drainage network and the rise of the water level in the manholes.
In order to more comprehensively analyze the change of projection values and the LID control characteristics under complex scenarios, areas in J9 and J15 that have the greatest projection values were selected for analyzing the change characteristics of flooding under different complex rainfall scenarios (Figures 6 and 7). After the implementation of LID practices, the projection values of these two manholes were significantly reduced by at least 28%. However, the projection values increase with the increase in rainfall intensity, rainfall duration or rainfall peak coefficient. With the increase in rainfall intensity, there is a more obvious reduction in the impact of rainfall duration and rainfall peak coefficient on projection in the area (J9) with the greater amount of impervious areas. Some similar conclusions are obtained by Qin et al. [29].

Changes in the Comprehensive Flood and LID Control Characteristics under Different Rainfall Scenarios
According to the comprehensive distribution characteristics of flooding, the projection values of the manholes J9 and J15 and the downstream area are high. The reasons for this are as follows: (i) for J15, it is located in the joint point of two pipes, and therefore its convergence area is large; (ii) J9 covers a large proportion of impervious surface in this area, and therefore rapid convergence speed and the severe overflow of the downstream pipe network contribute to grave flooding. The downstream flooding of the drainage system is caused by the failure of the timely drainage of stormwater, the (overly) full drainage network and the rise of the water level in the manholes.
In order to more comprehensively analyze the change of projection values and the LID control characteristics under complex scenarios, areas in J9 and J15 that have the greatest projection values were selected for analyzing the change characteristics of flooding under different complex rainfall scenarios (Figures 6 and 7). After the implementation of LID practices, the projection values of these two manholes were significantly reduced by at least 28%. However, the projection values increase with the increase in rainfall intensity, rainfall duration or rainfall peak coefficient. With the increase in rainfall intensity, there is a more obvious reduction in the impact of rainfall duration and rainfall peak coefficient on projection in the area (J9) with the greater amount of impervious areas. Some similar conclusions are obtained by Qin et al. [29].
LID implementation in J9 accounts for only a small proportion, and its effects on mitigating flooding are limited, making the control characteristics of LID practice not significant. When the rainfall duration is greater than 1.5 h or the rainfall peak coefficient is greater than 0.5, the projection value of J9 increases slightly. However, the upstream and downstream LID practices can alleviate the drainage pressure on the drainage network, thereby reducing the projection values of the areas where no LID practices are taken. On the contrary, the percentage of LID is relatively high upstream (above manhole J15), and the flood control characteristics of the LID practices are more obvious. The projection value of J15 is greater with changes in the rainfall duration and rainfall peak coefficient, suggesting that the area where many LID practices were implemented is significantly affected by rainfall duration and rainfall peak coefficient.   LID implementation in J9 accounts for only a small proportion, and its effects on mitigating flooding are limited, making the control characteristics of LID practice not significant. When the rainfall duration is greater than 1.5 h or the rainfall peak coefficient is greater than 0.5, the projection value of J9 increases slightly. However, the upstream and downstream LID practices can alleviate the drainage pressure on the drainage network, thereby reducing the projection values of the areas where no LID practices are taken. On the contrary, the percentage of LID is relatively high upstream (above manhole J15), and the flood control characteristics of the LID practices are more obvious. The projection value of J15 is greater with changes in the rainfall duration and rainfall peak coefficient, suggesting that the area where many LID practices were implemented is significantly affected by rainfall duration and rainfall peak coefficient.   It can be seen from the comparison between J9 and J15 that the influence of rainfall intensity, rainfall duration and rainfall peak coefficient on J15 are more obvious than that on J9. This is because LID implementation in J9 accounts for only a small proportion, and its effects on mitigating flooding are limited, making the control characteristics of LID practice not significant. When the rainfall duration is greater than 1.5 h or the rainfall peak coefficient is greater than 0.5, the projection value of J9 increases slightly. However, the upstream and downstream LID practices can alleviate the drainage pressure on the drainage network, thereby reducing the projection values of the areas where no LID practices are taken. On the contrary, the percentage of LID is relatively high upstream (above manhole J15), and the flood control characteristics of the LID practices are more obvious. The projection value of J15 is greater with changes in the rainfall duration and rainfall peak coefficient, suggesting that the area where many LID practices were implemented is significantly affected by rainfall duration and rainfall peak coefficient.

Conclusions
Rainfall has a great impact on urban flooding in the drainage system. In this study, the effect of LID practices on flood reduction in a typical residential area of Guangzhou were investigated under different rainfall scenarios. The flooding hot spots and their comprehensive characteristics were analyzed by the SWMM model with the PSO-based projection pursuit technique. The main conclusions were that: (i) the projection values of different manholes have different responses to the changes in rainfall scenarios, but the values of manholes J9 and J15 were the largest before and after the implementation of LID practices. Therefore, manholes J9 and J15 were flooding hotspots and should be paid more attention; (ii) LID practices are effective in flood reduction. They can control the flooding under the 2y-2h-0.375 rainfall scenario. However, the control effects of LID practices tend to decrease as the rainfall intensity, rainfall duration or rainfall peak coefficient increase; (iii) compared with different rainfall scenarios, the control ability of LID practices are more effective in flood reduction for shorter duration, lower intensity and smaller peak coefficient rainfall events; (iv) as the rainfall intensity increases, the impact of rainfall duration and rainfall peak coefficient on projection values is reduced, particularly in the area with the greater amount of impervious areas. Among rainfall intensity, rainfall duration and rainfall peak coefficient, the control effects of LID practices are most affected by rainfall intensity.