Runoff Effect Evaluation of LID through SWMM in Typical Mountainous, Low-Lying Urban Areas: A Case Study in China

Urban flooding occurs frequently in many regions of China. To reduce the losses caused by urban flooding, sponge city (SPC) and low-impact development (LID) have been carried out in many Chinese cities. However, urban flooding is influenced by various factors, such as climate, land cover characteristics and nearby river networks, so it is necessary to evaluate the effectiveness of LID measures. In this study, the Storm Water Management Model (SWMM) was adopted to simulate historical urban storm processes in the mountainous Fragrance Hills region of Beijing, China. Subsequently, numerical simulations were performed to evaluate how various LID measures (concave greenbelt, permeable pavement, bio-retention, vegetative swales, and comprehensive measures) influenced urban runoff reduction. The results showed that the LID measures are effective in controlling the surface runoff of the storm events with return periods shorter than five years, in particular, for one-year events. Furthermore, the effectiveness on traffic congestion mitigation of several LID measures (concave greenbelt, vegetative swales, and comprehensive measures) was evaluated. However, the effective return periods of storm events are shorter than two years if the effectiveness on traffic congestion relief is considered. In all evaluated aspects, comprehensive measures and concave greenbelts are the most effective, and vegetative swale is the least effective. This indicated that LID measures are less effective for removing ponding from most storm events in a mountainous, low-lying and backward pipeline infrastructure region with pressures from interval flooding and urban waterlogging. The engineering measures including water conservancy projects and pipeline infrastructure construction combined with the non-engineering measures were suggested to effectively control severe urban storms.


Introduction
High-speed urbanization has caused rapid changes to underlying surfaces, resulting in fundamental changes in urban runoff processes.Due to these changes, urban flooding frequently occurs, which is widely influenced by climate change.This has been more distinctly reported as a "sea view in the city" in recent years in China.Urban waterlogging brings a series of socio-economic losses such as traffic paralysis, loss of property, and even human casualties.The effective control of urban waterlogging is both crucial and difficult to manage for urban stormwater runoff.The evaluations of the effectiveness of the various management measures are more important.Various runoff reduction measures have been implemented, including LID, and widely applied throughout the world [1].LID is an innovative urban stormwater management system that was jointly introduced by the stormwater management experts from Programs and Planning Division of Prince George's County Department of Environmental Resources, during the mid-1990s [2].Since the end of the 1990s, with the promotion of LID by the United States Environmental Protection Agency (USEPA), it has been generally recognized and adopted by countries all over the world.A series of related stormwater management regulations have been formulated in Florida, Chicago, and other locales, and remarkable achievements have been made [3,4].The idea behind LID is to depress the negative influence of water quantity as well as the quality of the runoff process caused by urbanization to allow regional runoff processes to return to a natural undeveloped state to the largest degree possible [2].Its meaning has been extended in many countries or regions that have similar ideas regarding stormwater management, such as Water-Sensitive Urban Design (WSUD) in Australia, Sustainable Urban Drainage Systems (SUDS) in the UK, and the natural drainage systems in Seattle [1,[5][6][7].Based on these concepts, SPC for urban stormwater management was proposed in China.The concept of SPC is that the city is similar to a block of sponge that can absorb water in the wet period and release water during a drought, aiming to relieve the waterlogging and water shortage situation in cities.It was first mentioned by the Chinese President in 2013 and proposed in 2014.SPC is a comprehensive concept, focusing on not only the source control (as LID) but also the process control and drainage management.In 2015, SPC construction was proposed, and sixteen cities were selected as the pilot cities in China.The aim of SPC construction is the detention, onsite storage, infiltration, filtration, reuse, and drainage of stormwater [8,9].
The two main purposes of the series of LID measures are controlling the quantity of storm runoff and improving its quality.Regarding the former, LID measures primarily reduce the flood peak and volume by increasing the infiltration and onsite storage.Only the reduction of surface runoff quantity was involved in this study, and the relevant practice measures mainly included permeable pavements, bio-retention facilities, green roofs, and swale systems [10][11][12], measures that have been commonly applied.The permeable pavement is mainly composed of artificial materials (such as porous asphalt and permeable tile) with improvements to the water permeability, and the runoff is reduced by increasing the rainfall infiltration through these measures.Related studies have shown that the runoff reduction has reached 90% in some regions that have applied this measure [13].However, the reduction percentage has been found to vary under different rainfall conditions in different regions [1].Bio-retention facility systems are systems that use plants, soil, and micro-flora for storage infiltration and to purify rainwater runoff.These systems were originally designed for rainwater purification applications [2,3].However, Davis et al. [14] noted that runoff processes can also be reduced by increasing the soil infiltration and evapotranspiration since these measures are close to the natural state.A green roof system is an earlier applied LID measure [15] that mainly cuts runoff by planting vegetation on high-quality roofs.Studies have shown that these measures were able to significantly reduce runoff, and the runoff reduction percentage reached between 29% and 100% in some regions [16].Swale systems, which are different from the three aforementioned systems, are the measures that have demonstrated effectiveness in runoff control, not only at the source but also during runoff processes, especially in terms of reducing the flood velocity [2,12].Cheng et al. [17] examined the impacts of concave greenbelts' design parameters on storm runoffs.The results showed that when the concave greenbelt area ratio was between 10% and 30% and the concave depth ranged from 0.1 m to 0.3 m, the return periods of the designed rainfall would be one year, three years, and five years.
Although the global applications of LID have been relatively successful, the scales of the areas where it is applied are normally units or are experimental, and few have focused on the basin or region scales [18].The reason for this is that the LID was originally designed as a micro-scale measure for small urban areas [2,3].In addition, the areas where LID has been applied have been mainly located in Water 2017, 9, 439 3 of 21 flat, highly urban areas, including parks in downtown areas [19] and city blocks [20]; few have focused on mountain cities [21].The main problem of stormwater management in the urbanization plain areas is local waterlogging, but mountain cities can also suffer floods of the upper reaches or interval floods.China has a varied topography with great differences in climate, and it is also a developing and urbanizing country.These factors make for more complex underlying surface situations, and regional urban waterlogging usually has different causes in China.However, a boom of LID construction has recently occurred in China after advancing SPC construction guidance [22][23][24].It seems that LID measures are universal for urban water problems.In fact, the effects of different LID measures at different rainfall intensities should be evaluated in different regions.The related evaluation should include the most suitable LID measure, its more economic size, its effective rainfall intensity, etc.
Fragrance Hills, located in Beijing, China, was selected as the mountainous test area in this study.The stormwater management in this region includes not only waterlogging caused by unitization and poor drainage but also the threat of interval flash floods.The simulated effectiveness of surface runoff reduction with four different LID measures (concave greenbelt, permeable pavement, bio-retention, vegetative swale), along with their composed measures through SWMM, were evaluated and compared.This study focused on: (1) model construction and validation for regions with no observed discharge; (2) effect evaluation of different LID measures in a semi-natural-semi-urbanized mountainous area in rain absorption and flood peak reduction; and (3) effect evaluation of different LID measures in improving regional traffic jams under different rainfall conditions.This study comprehensively compared the differences of regional water-logging and traffic jams after the implementation of LID measures through SWMM.The results can provide technical support for optimizing and accelerating the construction of regional drainage pipe networks.This study can also provide integrated technical support and decision-making references for optimizing investment in local SPC construction and improvements in the rational planning of SPC implementation schemes.

Study Area
The area in this study is located at the foot of the eastern side of Fragrance Hills, which is a closed mountain catchment (hereafter, FHC).The total area is 11.75 km 2 , the highest elevation is 667 m, and the lowest elevation is 70 m.It is mainly a mountainous area that belongs to Fragrant Hill Town, Haidian District of Beijing, China.Its administrative divisions and geographical location are shown in Figure 1. have focused on mountain cities [21].The main problem of stormwater management in the urbanization plain areas is local waterlogging, but mountain cities can also suffer floods of the upper reaches or interval floods.China has a varied topography with great differences in climate, and it is also a developing and urbanizing country.These factors make for more complex underlying surface situations, and regional urban waterlogging usually has different causes in China.However, a boom of LID construction has recently occurred in China after advancing SPC construction guidance [22][23][24].It seems that LID measures are universal for urban water problems.In fact, the effects of different LID measures at different rainfall intensities should be evaluated in different regions.The related evaluation should include the most suitable LID measure, its more economic size, its effective rainfall intensity, etc. Fragrance Hills, located in Beijing, China, was selected as the mountainous test area in this study.The stormwater management in this region includes not only waterlogging caused by unitization and poor drainage but also the threat of interval flash floods.The simulated effectiveness of surface runoff reduction with four different LID measures (concave greenbelt, permeable pavement, bioretention, vegetative swale), along with their composed measures through SWMM, were evaluated and compared.This study focused on: (1) model construction and validation for regions with no observed discharge; (2) effect evaluation of different LID measures in a semi-natural-semi-urbanized mountainous area in rain absorption and flood peak reduction; and (3) effect evaluation of different LID measures in improving regional traffic jams under different rainfall conditions.This study comprehensively compared the differences of regional water-logging and traffic jams after the implementation of LID measures through SWMM.The results can provide technical support for optimizing and accelerating the construction of regional drainage pipe networks.This study can also provide integrated technical support and decision-making references for optimizing investment in local SPC construction and improvements in the rational planning of SPC implementation schemes.

Study Area
The area in this study is located at the foot of the eastern side of Fragrance Hills, which is a closed mountain catchment (hereafter, FHC).The total area is 11.75 km 2 , the highest elevation is 667 m, and the lowest elevation is 70 m.It is mainly a mountainous area that belongs to Fragrant Hill Town, Haidian District of Beijing, China.Its administrative divisions and geographical location are shown in Figure 1.The FHC is in a warm, semi-humid zone of the East Asia monsoon region, with a significant continental monsoon climate and a large annual temperature range.The average annual temperature is 11.7 • C. According to statistics, the annual mean precipitation of the Haidian District from 1956 to 2010 was 525.4 mm, and it was 614.3 mm during normal years and 336.4 mm during dry years.The annual precipitation distribution has been extremely uneven, with a maximum of 1115.7 mm (1956), and a minimum of only 281.4 mm (1965).The precipitation during the wet years can be almost four times more than that during the dry years, with an extremely uneven distribution within the year.In other words, the precipitation during the flood season accounts for over 80% of the total year's precipitation.The groundwater recharge mainly depends on precipitation and river infiltration, and the ground runoff flows from the northwest to southeast in the ground.

Observed Precipitation Data
The observed precipitation processes were based on the gauge monitoring in the FHC and supplied by Beijing Water Authority (BWA).Three typical precipitation processes were selected for model calibration and validation.The storm of 21 July 2012 (hereafter, "7.21" storm) was a typical storm event in Beijing beginning at 12:20 on 21 July and ending at approximately 5:00 on 22 July.Its total precipitation was approximately 264.5 mm in less than 17 h.This storm had significant effects and 1.9 million people suffered from the storm, of whom at least 79 died in Beijing.The storm of 20 July 2016 (hereafter "7.20" storm) was a typical moderate rain event and the storm of 7 September 2016 (hereafter, "9.07" storm) was a light event, compared with "7.21".The characteristics of these rain events are listed in Table 1.In this study, typical precipitation processes of one-year, two-year, five-year, and ten-year return periods (five-minute precipitation processes, which lasted 24 h in total) were selected as the storm inputs for LID effectiveness evaluation.They were calculated according to the local hydrologic handbook, technical code and standard [25][26][27], and the related processes are shown in Supplementary Material S1.The precipitation processes with different return periods are shown in Figure 2.

Drainage Pipeline Data
The drainage pipe networks have not been constructed in the FHC.Only one drainage pipeline is laid beneath Fragrant Hill South Road; the pipeline length is 438 m and the diameter is between 1 m and 2 m.Due to the large topographic slope, overland flow frequently occurs in the mountainous areas, moves along the roads or pavements toward the low-lying areas, and flows into canals or rivers.Thus, the roads as the drainage channels are the main part of the drainage system in this study (Figure 3a).sub-catchments were further divided according to the drainage directions and boundaries (rivers, mountain gullies and roads in Figure 3a), which were both identified by the RSI.Third, the region was divided into small sub-catchments based on the types of underlying surfaces identified by the RSI, historical investigations data and the related research [28].Finally, based on the balance between simulation precision and elapsed time of computer, the entire modeling region was divided into 72 sub-catchments, and the confluence path and confluence nodes were also built (Figure 3b).

Investigation, Survey and Interview
The supplementary investigation, survey and interview (ISI for short) were implemented.Because the observed discharge data are absent in the FHC, five positions (Figure 4) were selected sub-catchments were further divided according to the drainage directions and boundaries (rivers, mountain gullies and roads in Figure 3a), which were both identified by the RSI.Third, the region was divided into small sub-catchments based on the types of underlying surfaces identified by the RSI, historical investigations data and the related research [28].Finally, based on the balance between simulation precision and elapsed time of computer, the entire modeling region was divided into 72 sub-catchments, and the confluence path and confluence nodes were also built (Figure 3b).

Investigation, Survey and Interview
The supplementary investigation, survey and interview (ISI for short) were implemented.Because the observed discharge data are absent in the FHC, five positions (Figure 4) were selected

Elevation Data and Sub-Catchment Division
The original elevation data include Digital Topographic Map (DTM) 1:2000 and Remote Sensing Images (RSI) 1:2000, both supplied by Beijing Institute of Surveying and Mapping (BISM).The digital elevation model (DEM) was compiled based on DTM.RSI, DEM and the underlying surface properties were applied for the sub-catchment division.
The principles of the sub-catchment division are the elevation, the uniformity of underlying surface, and the nearest drainage in this study.First, the entire region was roughly divided into three large sub-catchments (areas with red boundaries in Figure 3a) based on the DEM.Second, the large sub-catchments were further divided according to the drainage directions and boundaries (rivers, mountain gullies and roads in Figure 3a), which were both identified by the RSI.Third, the region was divided into small sub-catchments based on the types of underlying surfaces identified by the Water 2017, 9, 439 6 of 21 RSI, historical investigations data and the related research [28].Finally, based on the balance between simulation precision and elapsed time of computer, the entire modeling region was divided into 72 sub-catchments, and the confluence path and confluence nodes were also built (Figure 3b).

Investigation, Survey and Interview
The supplementary investigation, survey and interview (ISI for short) were implemented.Because the observed discharge data are absent in the FHC, five positions (Figure 4) were selected according to the flood drainage routes and the waterlogging situation in the low-lying sub-catchments: ( 1  The flood marks of typical runoff processes ("7.21", "7.20" and "9.07") in each position were selected as the key factor of ISI and the corresponding accurate values were collected through detailed and repeated ISI.Because the roads functioned as drainage channels, the height of flood marks was equal to the maximum runoff depth (MRD) in this study.The ISI MRD (hereafter, MRDI) of different runoff processes in each position were collected and selected as the calibration and validation variable.
The peak time (PT) is an important characteristic of the runoff processes and the approximate values of the typical processes in each position were collected through interview (hereafter, PTI).PTI was selected as the reference variable for model calibration and validation.
In addition, the types of underlying surface and the geometric characteristics of sub-catchments were also investigated, surveyed and applied for setting the initial values of the model parameters and the correct sub-catchment division.The flood marks of typical runoff processes ("7.21", "7.20" and "9.07") in each position were selected as the key factor of ISI and the corresponding accurate values were collected through detailed and repeated ISI.Because the roads functioned as drainage channels, the height of flood marks was equal to the maximum runoff depth (MRD) in this study.The ISI MRD (hereafter, MRD I ) of different runoff processes in each position were collected and selected as the calibration and validation variable.
The peak time (PT) is an important characteristic of the runoff processes and the approximate values of the typical processes in each position were collected through interview (hereafter, PT I ).PT I was selected as the reference variable for model calibration and validation.
In addition, the types of underlying surface and the geometric characteristics of sub-catchments were also investigated, surveyed and applied for setting the initial values of the model parameters and the correct sub-catchment division.

Hydrological Model Selection
There are many types of methods for evaluating the effects of LID measures; the most simple and accurate method is using models.The relevant models include the Long-Term Hydrologic Impact Assessment-Low Impact Development (L-THIA-LID) model [29]; the improved SCS-CN model [30,31]; the System for Urban Storm-water Treatment and Analysis Integration (SUSTAIN) model [32]; and the SWMM [21,[33][34][35].In 2016, DHI-China developed the Sponge City Aided Design (SCAD) module and added it to the Mike Urban model directed against the requirements of the Chinese market [36].
Among these models, the SWMM is widely applied in surface runoff simulations.It is a physically based simulation model in which the storm runoff quantity can be simulated based on the physical processes of surface runoff, infiltration, surface ponding and flow routing.It has a flexible structure and an open code.SWMM has been updated several times, and its functions have been improved significantly [37].Compared with other models, SWMM has a better simulation effect [38], and it has been widely used in the evaluation and forecasting of the surface runoff process of storm floods by relevant scholars.Marsalek et al. [33] conducted in-depth studies on 12 storm events in three typical small watersheds in the United States and found that the SWMM simulation results fit well with the measured runoff results.Daeryong et al. [34] established a SWMM model for Ulsan, Korea, and analyzed the required scales and storage ranges of three types of storage ponds (with return periods of 2 years, 10 years, and 100 years).Furthermore, they assessed and compared the construction and land costs, as well as the benefits of these ponds.It was found that the storage pond with a return period of 2 years was the most profitable.Villarreal et al. [35] simulated the effectiveness after stormwater and flood management measures were applied in a central urban area, and the rainfall return periods were 0.5, 2, 5, and 10 years.
Within LID simulation, SWMM is one of the earliest hydrological models to be supplemented with an LID module [2].The LID module of SWMM provides five single measures for storm runoff control.The hydrological process and its key factors, such as regional surface runoff and peak discharge, can be simulated by applying the LID module combined with the hydraulic module.Hence, the runoff reduction effect of LID can be evaluated through these key factors.The SWMM for different LID measures has also been extensively used in different regions and countries, especially in China.Using SWMM, Qin et al. [39] evaluated the effect of LID measures (including the LID compositions) under different rainfall types in the Guangming New District of Shenzhen, and Hu et al. [40] analyzed the implementation effect of LIDs for municipal roads in this region.Jia et al. [19] evaluated the implementation effect of LID-BMPs for Beijing Olympic Park.Li et al. [41] performed a simulation analysis on a rainwater garden of a residential area in Xi'an.Wang et al. [42] assessed the drainage pipeline reconstruction and LID planning in a residential zone of Beijing as well.In addition, relevant research has shown that SWMM could be useful in areas including plains [19], mountains [12,21] and karsts [43].
Because the SWMM was well applied in both storm runoff simulations and LID evaluations and was appropriate for different types of areas, it was selected as the simulation model.SWMM 5.0.022released in 2011 [12,44] was selected in this study.

Initial Value of Model Parameter
The initial values of the sensitive parameters are crucial for model simulation.The related studies [45,46] indicated that the characteristic width (W) of the region, impervious coefficient, Manning W is determined by the shape and size of the sub-catchment.Because the study area is large and its topography and underlying surface are complex, the flow length is difficult to determine, i.e., the calculation method based on flow length is inappropriate in this study.Therefore, W was calculated according to Equation (1) [46], where K is the shape correction coefficient and A is the area of the sub-catchment.K is crucial for calculating W and was calibrated based on the measured geometric characteristics.According to the shape of the sub-catchments in the FHC (Figure 3b), the initial value of K was assigned as 0.9 in this study.A and W of each sub-catchment were listed in Table S1.
The impervious coefficient can directly influence the surface runoff quantity.In this study, the impervious areas of each sub-catchment were calculated based on the RSI and investigation results, and then the corresponding impervious coefficients were obtained (Table 2).n and D of the impervious and pervious areas are influenced by the properties of the underlying surface, and the setting of initial values are all empirical.Because the roads were generalized as channels and the FHC is a shrub area [28], the initial n of impervious areas (channels) were set as 0.017, and the initial n of pervious areas were set as 0.4 (Table 2), according to the manual of SWMM 5.0 [46].In this study, D of the impervious and pervious areas were set as 2 and 7 mm, respectively (Table 2), based on the SWMM manual, previous studies [12] and characteristics of the FHC.These sensitive parameters were chosen as the optimization parameters for the model calibration and validation.
The model supplies three methods for infiltration processes, and the Horton model was selected in this study.The infiltration (f ) was calculated as Equation (2) [44], where f c is the maximum infiltration rate, f 0 is the minimum infiltration rate, and k is the attenuation coefficient.
Although f c and f 0 are less sensitive compared to the above parameters, they are quite critical in infiltration calculation.The soil in the FHC is mainly composed of loamy soil and cinnamon soil, and related research supplied the plant species of study area [28].Due to these properties, the initial values of f c in each sub-catchment, the initial k and the initial f 0 were set according to the SWMM manual (Table 3).The other parameters of SWMM such as cross-section size, the characteristics of the drainage pipeline (channel), and slope, were calculated based on the DEM, RSI and investigation data (Table S1).

Model Calibration and Validation
Because the discharge data are absent in the FHC, the catchment model was calibrated and validated based on the agreement between the simulated MRD (hereafter, MRD s ) and MRD I in each ISI position.The absolute value of relative error (|RE|) was chosen as the measure for evaluating the agreement, described as Equation (3).The agreement between the simulated PT (hereafter, PT s ) and PT I in each position was chosen as the reference for model calibration and validation.
The "7.21" storm was chosen for the model calibration and its simulated runoff processes in different positons are shown in Figure 5.The optimal parameters values in each sub-catchment were finally selected (Table 4) according to the minimum |RE| (Table 5) in each ISI position.It was shown that MRD s was similar to MRD I in most of the positions and the related average |RE| was only 7.90% (Table 5) and each PT s was close to PT I (Table 6).
Water 2017, 9, 439 9 of 20 ( ) The "7.21" storm was chosen for the model calibration and its simulated runoff processes in different positons are shown in Figure 5.The optimal parameters values in each sub-catchment were finally selected (Table 4) according to the minimum |RE| (Table 5) in each ISI position.It was shown that MRDs was similar to MRDI in most of the positions and the related average |RE| was only 7.90% (Table 5) and each PTs was close to PTI (Table 6).Table 4. Calibrated results of the parameters.The "7.20" and "9.07" storms were chosen for model validation and their simulated runoff processes are shown in Figures 6 and 7, respectively.It was illustrated that none of the |RE| were larger than 10% and the average value was 7.64% and 2.88%, respectively (Table 5).The values of PT s were all closed to the corresponding PT I (Table 6).The calibration and validation results indicated that the model simulation results were consistent with the investigation data, according to the national standard [47].The model can be feasible for flood simulation under rainfall with different return periods.

LID Selection and Design Scenarios
There are different LID components in SWMM 5.0.022, of which concave greenbelts (rain gardens), permeable pavement, bio-retention cells, and vegetative swales were chosen for the simulation according to the suggestion of some relevant authority research [48] and the local technical code [27].These single LIDs and their composition measures were designed as different scenarios.The parameters of each LID measure were assigned according to the underlying surface types of each  The calibration and validation results indicated that the model simulation results were consistent with the investigation data, according to the national standard [47].The model can be feasible for flood simulation under rainfall with different return periods.

LID Selection and Design Scenarios
There are different LID components in SWMM 5.0.022, of which concave greenbelts (rain gardens), permeable pavement, bio-retention cells, and vegetative swales were chosen for the simulation according to the suggestion of some relevant authority research [48] and the local technical code [27].These single LIDs and their composition measures were designed as different scenarios.The parameters of each LID measure were assigned according to the underlying surface types of each sub-catchment in the FHC and the local technical regulation [27].The absence of LID was designed as another scenario for comparison.

Concave Greenbelt Scenario
Concave greenbelts were set in 20 different typical residential communities (locating in different sub-catchments).The designed concave depth was 10 cm and the concave ratios were 50%, 70% and 90% as three typical sub-scenarios.The other related parameters are listed in Table S2.

Permeable Pavement Scenario
Permeable pavements were set in the roads of 20 residential communities (the same with the concave greenbelts scenario) and in the municipal sidewalks of 30 sub-catchments.The permeable pavement ratio in municipal sidewalks was set as 100% and that in the roads of residential community were set as 50%, 60% and 70%, respectively, also as three typical sub-scenarios.The thickness of the surface layer was 60 mm with 0.2 porosity; the thickness of the entire layer was 200 mm with 0.75 porosity; and the infiltration rate was assigned as 0.1 mm/s according to the local technical regulation [27].The other related parameters were listed in Table S3 (residential community roads) and Table S4 (municipal sidewalks).

Bio-Retention Cell Scenario
Bio-retention cells were set in two typical residential communities (located in S21 and S69 sub-catchment).The depth of the surface storage was 150 mm, with 90% vegetation coverage; the soil thickness was 500 mm, with 0.5 porosity.The field capacity was 0.35 and the wilting point was 0.187.The conductivity gradient of the bio-retention cell was 10, with 36 mm/h hydraulic conductivity and the height of the gravel layer was 400 mm with 210 mm suction of the soil moisture.The other related parameters were listed in Table S5.

Vegetative Swale Scenario
Five vegetative swales were set on the sides of Fragrant Hill Road (Figure 4), because the slope of this road is smaller in the FHC and its ponding situation was severe, of which the length was 213.43 m, 313.86 m, 63.41 m, 130.63 m and 275.92 m, respectively.The total swale length was more than 997 m and the total storage volume was approximately 157 m 3 .All the sections were U shaped with 0.1 m storage depth.The average suction of soil moisture along the capillary was 169.9 mm, the hydraulic conductivity of saturated soil moisture was 0.6138 mm/h, and the initial loss of soil moisture loss (difference value between the soil porosity and initial soil moisture) was 39.25 mm.The storage areas of different water depth in each swale were listed in Table S6.

Comprehensive Measures Scenario
The comprehensive measures were set as the combination of a concave greenbelt scenario with 50% concave ratio, permeable pavement scenario (of which the ratio was 70% in residential community roads and 100% in municipal sidewalks), a bio-retention cell scenario and a vegetative swale scenario.

Basic Scenario
There are three outfalls in the FHC (marked as O1, O2 and O3 in Figure 4), because "O1" is the outfall nearest Xiangquan Roundabout and in the low-lying area, the runoff processes of no LID at O1 outfall were selected as the base case for comparison.
The Permeable pavements and bio-retention cells were both applied in the part of the FHC and evaluating their impacts on the overall study area was inappropriate.Therefore, the runoff processes of no LID in two typical sub-catchments were designed as the base case for evaluating and comparing the effectiveness of these two measures.S38 with a larger impervious area was selected for permeable pavement scenario and S69 close to Xiangquan Roundabout was selected for the other one (Figure 4).

Results and Discussion
The different runoff processes in different scenarios were simulated through the validation model, while the above designed precipitation processes (Figure 2) were input.It is evident that the reduction of runoff volume and peak discharge is crucial to the effectiveness evaluation of LID.However, traffic congestion is also an apparent problem influenced by urban flooding.The FHC is located on the top of the Xiangquan Roundabout, which is an important western transportation hub of the Haidian District.This hub is in the low-lying zone and has frequently encountered waterlogging and flooding from the FHC in recent years.Thus, the effectiveness of relieving traffic congestion through the LID measure application was also involved in this study.

Results in Basic Scenario
The runoff processes under different return periods of O1 outfall (Figure 8), S38 and S69 were simulated in the model.Their corresponding runoff depth and the MRDs of O1 are listed in Table 7. Table 7 shows that the runoff increased sharply when the return period increased: the runoff depth is only approximately 3 mm in the one-year event but more than 80.0 mm in the ten-year event.The "O1" runoff processes of no LID (Figure 8) shows that the runoff volume initially raised slowly, then after 16:30 raised sharply, and then declined quickly.These trends of runoff processes are consistent with those of the corresponding rainfall processes (Figure 2).In addition, the PTs of runoff processes are all later than those of the rainfall processes, which agrees with the law of surface runoff and drainage.ratio is only 0.1% in the five-year and ten-year rainfall events.
There are several reasons behind these results: the first and foremost is the limit of the size and numbers of vegetative swales that influenced the water storage capacity, resulting in the ineffectiveness on peak discharge reduction in the study area; the second reason is that the coupling reaction of local waterlogging and the interval flood in low-lying areas would be more prominent when the storm intensity increased.This explanation is concordant with the analysis of the runoff volume reduction.

Traffic Congestion Relief
The ponding situation of "O1" is paramount to relieve the traffic congestion on the Xiangquan Roundabout during storms, due to the station of O1 outfall (Figure 4).Therefore, the duration and depth of ponding were both involved in this section.The critical value of the ponding depth was assigned as 15 cm for the drainage design of municipal roads, considering the limit of the ponding height for the vehicles.Thus, the durations of ponding depth over 15 cm at O1 in different LID scenarios and the basic scenario were calculated (Table 10) and can reflect different traffic situations on the Xiangquan Roundabout.Table 10 and Figure 8 illustrate that the frequency of traffic  Based on these results in basic scenarios, the evaluation and comparison of different scenarios (LID measures) in the above aspects were studied and illustrated as follows.

Runoff Volume Reduction
The percentages of runoff reduction in different LID scenarios were calculated and analyzed for evaluating their runoff reduction effectiveness (Table 8).Among all of the single measures, concave greenbelt is the most effective in runoff reduction, especially when the return period of precipitation is shorter than two years.The runoff reduction percentages of concave greenbelt ranged from approximately 5.2% to 57.3% in all types of rainfall events and the reduction percentages are much higher (20.5-57.3%) in the one-year and two-year rainfall events.It is notable that when the concave ratio is 50%, the total runoff can be reduced to 52.5% for the one-year rainfall event, proving the high effectiveness of concave greenbelts in flood mitigation.By contrast, when the concave ratio was 70% or 90%, the reduction percentage increased slightly in comparison with that of the 50% ratio, within a range of 1.5-5.0% in different types of rainfall.That is, the concave ratio of 50% is better for stormwater detention.Therefore, there is no need to expand the facility scale.Vegetative swale is the least significant among all the single measures, and its entire runoff reduction percentages at the "O1" outlet are only 0.3-3.0%.There are two reasons for this: the first reason is the limits of the application area and frequency in mountain areas (only five lengths of one road applied this measure in the whole region); the other reason is the fast flooding of mountain floods.These results indicated that the vegetative swale is not suitable in mountains areas, although it is effective in controlling runoff in plain areas [49][50][51].
In the permeable pavement scenario, the runoff reduction percentage of S38 ranged from 7.3% to 12.2% in the one-year rainfall event and from 1.9% to 7.0% in the ten-year rainfall event when the permeable pavement ratio of the roads ranged from 50% to 70% (Table 8).This result indicated 70% permeable pavement is much more effective than 50% permeable pavement when it was applied in the roads of a typical residential community.
In the bio-retention measure scenario, Table 8 illustrated that the total runoff of S69 was reduced significantly to 10.2-12.1%.It is remarkable that when the return period of rainfall increased, the decrease trend of runoff reduction percentage was smoother than that in the other scenarios, which implies that this measure has a better reduction effectiveness in various rainfall intensities.
In the comprehensive measures scenario, the percentages of runoff reduction ranged from 15.4% to 57.3%The results compared with the other scenarios demonstrated that the reduction effectiveness (57.3%, the maximum one) is equal to that of 90% concave greenbelt in the one-year rainfall event and the comprehensive utilization of various measures can further improve the effectiveness in the other rainfall events.It is notable that the effectiveness in the comprehensive measures scenario is not equivalent to the summation of each single measure.A possible reason for this phenomenon is that applying various measures changes the state and process of the runoff in part of the study area, which likely offsets some of the overall effectiveness.As a result, the comprehensive utilization of various measures is less effective than the total effectiveness of various measures.
Overall, all the above LID measures can reduce the runoff volume to a certain magnitude in the rainfall events with the return periods shorter than five years and they are less effective in the precipitation events in the return periods longer than five years.The concave greenbelt and comprehensive measures were more effective, while the vegetative swale was the least effective measure.Permeable pavement can have a certain effect in typical sub-catchment and the bio-retention has better reduction effectiveness for a variety of rainfall intensities in the typical catchment.The effectiveness of runoff reduction in all of the scenarios decreased as the rainfall intensity increased.

Peak Discharge Reduction
The different runoff processes of "O1" in the concave greenbelt scenario (with 50% ratio), vegetative swale scenario and comprehensive measures scenario were simulated (Figure 8) and the percentage of MRD reduction was chosen as the evaluating variable.Then, the effectiveness of these three LID measures on peak discharge reduction was calculated (Table 9).Bio-retention and permeable pavement were not involved because they were only applied in several typical sub-catchments.Figure 8 and Table 9 demonstrate that the comprehensive measures are the most significant in the effectiveness of peak discharge in all rainfall events.By contrast, vegetative swale is the least significant.The concave greenbelt is more significant in the one-year and two-year rainfall events.The comprehensive measures and concave greenbelt are particularly effective on the peak discharge reduction in the one-year rainfall event, and the reduction percentage is 55.7% and 53.9%, respectively.However, the vegetative swale plays a poor role in this aspect, and the related reduction ratio is only 0.1% in the five-year and ten-year rainfall events.
There are several reasons behind these results: the first and foremost is the limit of the size and numbers of vegetative swales that influenced the water storage capacity, resulting in the ineffectiveness on peak discharge reduction in the study area; the second reason is that the coupling reaction of local waterlogging and the interval flood in low-lying areas would be more prominent when the storm intensity increased.This explanation is concordant with the analysis of the runoff volume reduction.

Traffic Congestion Relief
The ponding situation of "O1" is paramount to relieve the traffic congestion on the Xiangquan Roundabout during storms, due to the station of O1 outfall (Figure 4).Therefore, the duration and depth of ponding were both involved in this section.The critical value of the ponding depth was assigned as 15 cm for the drainage design of municipal roads, considering the limit of the ponding height for the vehicles.Thus, the durations of ponding depth over 15 cm at O1 in different LID scenarios and the basic scenario were calculated (Table 10) and can reflect different traffic situations on the Xiangquan Roundabout.Table 10 and Figure 8 illustrate that the frequency of traffic congestion resulting from flooding is seldom in the one-year event in all the scenarios; the situation of traffic congestion becomes more severe with the growing of the rainfall return periods.The shortening percentage in different LID scenarios (Table 11) were calculated based on the duration in basic scenario, for the convenience of evaluation and comparison.Tables 10 and 11 indicate that the concave greenbelt and comprehensive measures are both effective in relieving the traffic congestion in all of the rainfall events, and the least percentage is 29.1%.In the two-year rainfall event, the corresponding effectiveness for both scenarios was significant.These two durations of ponding depth over 15 cm are 48 min and 30 min, respectively (both less than one hour), and the corresponding shortening percentage is 51.0% and 69.4%, respectively.The result indicated only by these percentages is concordant with evaluation results in the other aspects.
However, the vegetative swale played less of a role in relieving traffic congestion in all of the rainfall events.Its largest shortening duration is only 2 min in the two-year event, and it is ineffective in the ten-year event.The ineffectiveness of applying vegetative swale in traffic congestion relief is inevitable when this measure showed its ineffectiveness both in runoff volume reduction and perk discharge reduction.
Tables 10 and 11 also indicate that the effectiveness in traffic congestion relief worsened when the return periods of rainfall increased, and this decline trend is similar to those in the other aspects (Tables 8 and 9).In addition, the comprehensive measures, as the most effective LID measures in this aspect, decreased most sharply when the return periods increased.Its shortening percentage decreased from 69.4% to 34.0%.
However, all of the durations of over 15 cm are more than one hour in the five-year and ten-year rainfall events.This means that traffic congestion would also exist in these rainfall events in the Xiangquan Roundabout if these LID measures are constructed.Table 10 indicates that the duration over 15 cm is approximately 50 min in two-year event, if the concave greenbelt is applied.This means that the traffic congestion would also be serious in the two-year event, although the related reduction percentages is high, which is different from the evaluation of the other methods.
Because Xiangquan Roundabout is located in the low-lying region of the FHC and its pipeline infrastructure is deficient, LID measures have less of an effect when the precipitation intensity is large; and even if LIDs could have a greater effect, they would only be effective in the storage of the local storm, while the flood from the upper reaches and interval region can also remarkably result in the waterlogging of this transportation hub.However, the substantial controlling measures for the mountain flood should combine water conservancy projects and drainage pipeline networks, and only LID measures are less effective.This analysis would supply the different decision-making for the local government.
In summary, only the implementation of LID measures was far from ideal for solving the interval problems of both flood and waterlogging on the mountainous areas.The analysis agrees with the results of the other relevant research [52].For the mountainous areas, rational analysis and decision-making in accordance with the regional characteristics and the above-described results are required before the implementation of stormwater management, flood management and SPC.In this study, the critical goal of regional waterlogging control is to guarantee the normal operation of transportation hubs of the Xiangquan Roundabout and relieving the traffic congestion.Accordingly, the construction of drainage pipe networks in the downstream low-lying areas should be first considered in making SPC schemes, i.e., the condition of the flood pouring down the street should be avoided.If well-funded, some retaining facilities such as flood-cutting channels should be built in upstream areas.Meanwhile, the drainage network should be expanded and the drainage canals should be dredged on both sides of the roads in middle-reach areas.Moreover, the LID measures coupled with underground tunnels should be built in downstream areas.Only when the runoff during the entire process was controlled, the dilemma of Xiangquan Roundabout being frequently flooded after storms would be fundamentally improved.
In addition, effective non-engineering measures should be employed in the FHC.A regional water-logging risk map should be developed and announced according to the simulation results of the ponding depth under different storm scenarios.The awareness of risk evasion should be acknowledged and build up to the public, especially the inhabitants in the low-lying areas (also the high-risk areas).Meanwhile, the corresponding countermeasures for different levels of ponding risk should be formulated.Although the risk map is static, it should be the crucial reference to the local regulation and planning for the future urban development and construction.Based on this static map, the regional real-time water-logging warning would be achieved combining the dynamic rainfall forecasting through satellite and radar.Furthermore, due to the importance of Xiangquan Roundabout for local transportation, a regional risk map for traffic jam and the corresponding countermeasures for various traffic jam risks should also be formulated according to the suggested traffic alleviation measures in this paper.In this way, the real-time warning of traffic congestion would be achieved, based on the risk map for traffic congestion and the dynamic warning of water-logging risk.If the congestion risk was high, the public would be notified in advance.Thus, the residents could optimize the traffic routes and avoid driving in the high-risk area.Vehicles and people that have entered the high risk area should be guided to the safe spots in accordance with the established evacuation routes, in order to minimize the property losses and lower the risk of traffic congestion.
In the future, field experiments and effect analysis on different LID measures should be carried out based on the related research results and can be convenient for comparing the effectiveness of the model simulation compared with practical engineering.The related results would supply more significant reference for decision-making, to further improve future stormwater and flood management.Additionally, the expected economic benefits of employing LID measures would be analyzed and discussed.These two aspects should be investigated in the future to propose reasonable construction scales of LID measures and then optimize regional the SPC scheme.
This study applied the data and information of detailed investigations, surveys and interviews for model calibration and validation due to the absence of observed runoff data in the study region.In the future, the pipeline and canal flow should be observed and monitored through a triangular weir, hydrometric propeller, and Acoustic Doppler Current Profilers (ADCP) to improve the modeling accuracy.
This study evaluated the effectiveness of the LID measures in reducing the quantity of runoff.In the future, the effectiveness on controlling non-point source pollution from storms should be analyzed and discussed.The above conclusions were applicable to mountainous urban units with semi-arid climate zones.In the future, the feasibility of the LID measures in urban units of humid zones and arid zones should be evaluated.

Conclusions
This study applied the SWMM for numerical simulations in the Fragrance Hills catchment (FHC) of Beijing, China, which is an urban mountainous low-lying region with no drainage pipe networks.This study extended the SWMM model to the data-sparse mountainous area, which was different from the previous applications of SWMM, which were mainly applied in highly urbanized regions.
Based on SWMM, the roads were generalized as drainage pipelines since the surface confluences in the region were flooded along the roads; then, according to DEM, RSI and field investigation data, the region was sub-divided into 72 sub-catchments, and thus the digital input of the drainage pipe networks for the model was finished.Because the observed discharge was absent, MRD, which is easily measured by a supplementary survey, was selected as the calibration variable.The parameters in the model were then calibrated using the typical storm process on 21 July 2012, and validated using a storm runoff process and a small rain runoff process on 20 July 2016, and 7 September 2016,

Figure 2 .Figure 3 .
Figure 2. Typical design precipitation processes of different return periods of: (a) one year; (b) two years; (c) five years; and (d) ten years.

Figure 2 .
Figure 2. Typical design precipitation processes of different return periods of: (a) one year; (b) two years; (c) five years; and (d) ten years.

Figure 2 .Figure 3 .
Figure 2. Typical design precipitation processes of different return periods of: (a) one year; (b) two years; (c) five years; and (d) ten years.

Figure 3 .
Figure 3. Drainage system (a); and sub-catchments (b) of study area.In this drainage system map (a), rea lines are mountain ridges, blue lines are rivers, black lines are mountain gullies and pink lines are roads.
) at the southeast gate of the Fragrant Hill Park and in the upper reaches of the entire valley, and the sub-catchment area mainly covered the Fragrant Hill Hotel and the upper sub-catchment areas; (2) at the intersection between the upper sector of Fragrant Hill and Zhenghuangqicun Roads, below the southeast sections of Fragrant Hill Park; the sub-catchment area mainly covered Fragrant Hill Park, as well as both sides of Fragrant Hill Road; (3) at the upper end of Fragrant Hill Valley, and the basin covered the entire southern sub-catchment area; (4) at the portal of the Fragrant Hill Road Tunnel; and (5) in the northern sub-catchment area at the exit of the Beihan River.Water 2017, 9, 439 6 of 20 according to the flood drainage routes and the waterlogging situation in the low-lying subcatchments: (1) at the southeast gate of the Fragrant Hill Park and in the upper reaches of the entire valley, and the sub-catchment area mainly covered the Fragrant Hill Hotel and the upper subcatchment areas; (2) at the intersection between the upper sector of Fragrant Hill and Zhenghuangqicun Roads, below the southeast sections of Fragrant Hill Park; the sub-catchment area mainly covered Fragrant Hill Park, as well as both sides of Fragrant Hill Road; (3) at the upper end of Fragrant Hill Valley, and the basin covered the entire southern sub-catchment area; (4) at the portal of the Fragrant Hill Road Tunnel; and (5) in the northern sub-catchment area at the exit of the Beihan River.

Figure 4 .
Figure 4. Distribution of the ISI positions.

Figure 4 .
Figure 4. Distribution of the ISI positions.
) and depression storage (D) of the impervious and pervious areas are sensitive parameters in SWMM.

Figure 8 .
Figure 8. Simulated runoff processes in different scenarios with different return-period events: (a) one year; (b) two years; (c) five years; and (d) ten years.

Figure 8 .
Figure 8. Simulated runoff processes in different scenarios with different return-period events: (a) one year; (b) two years; (c) five years; and (d) ten years.

Table 1 .
Characteristics of the three chosen rainfall events.

Table 2 .
Initial values of sensitive parameters.

Table 3 .
Initial values of the infiltration parameters.

Table 4 .
Calibrated results of the parameters.

Table 5 .
Comparison between MRDS and MRDI of typical rainfall processes.
Note:1Calibration spot is the corresponding ISI position in this study.The "7.20" and "9.07" storms were chosen for model validation and their simulated runoff

Table 5 .
Comparison between MRD S and MRD I of typical rainfall processes.

Table 6 .
Comparison between PT S and PT I of typical runoff processes.

Table 6 .
Comparison between PTS and PTI of typical runoff processes.

Table 9 .
Percentages of MRD reduction in different scenarios.

Table 7 .
The simulated runoff and MRD in basic scenario.

Table 8 .
Percentage of runoff reduction under different scenarios.

Table 9 .
Percentages of MRD reduction in different scenarios.

Table 10 .
Durations of ponding depth over 15 cm in different scenarios.

Table 11 .
Shortening percentages of durations of over 15 cm in different scenarios.