Study on StormWater Management of Grassed Swales and Permeable Pavement Based on SWMM

Grassed swales and permeable pavement that have greater permeable underlying surface relative to hard-pressing surface can cooperate with the city pipe network on participating in urban storm flood regulation. This paper took Nanshan village in Jiangsu Province as an example, the storm-water management model (SWMM) was used to conceptualize the study area reasonably, and the low-impact development (LID) model and the traditional development model were established in the region. Based on the storm-intensity equation, the simulation scene employed the Chicago hydrograph model to synthesize different rainfall scenes with different rainfall repetition periods, and then contrasted the storm-flood-management effect of the two models under the condition of using LID facilities. The results showed that when the rainfall repetition period ranged from 0.33a to 10a (a refers to the rainfall repetition period), the reduction rate of total runoff in the research area that adopted LID ranged from 100% to 27.5%, while the reduction rate of peak flow ranged from 100% to 15.9%, and when the values of unit area were the same, the combined system (permeable pavement + grassed swales) worked more efficiently than the sum of the individuals in the reduction of total runoff and peak flow throughout. This research can provide technical support and theoretical basis for urban LID design.


Introduction
China is in a period of rapid urbanization, with expansion of city scale and concentration of human population.Due to this process, the natural characteristics of the land are changed by various anthropogenic activities.The original forestland, farmland and grassland are gradually transformed into impermeable construction land.These anthropogenic activities are among the most important pollutant sources [1]; numerous pollutants are introduced into the urban environment by vehicular traffic, industrial processes, building construction and commercial activities, and are carried by storm-water.Above all, the impact of urbanization on the water environment includes increased risks in terms of floods, erosion and degradation of stream habitats, and deterioration of water quality [2,3].In addition, the recharge of shallow groundwater resources has been attenuated with the reduction in potential infiltration and recharge.The storm-water runoff originates from residential and commercial areas having large amounts of impervious area connected to the storm-water system known as directly connected impervious areas (DCIAs) [4].Therefore, the control and utilization of rainwater is imperative, and rainwater drainage systems should be changed from the traditional rainwater drainage into systems that mainly rely on permeability and combine with storage and drainage [5,6].The concept of sponge city construction emerges as the times require.
The concept of sponge city refers to the city that can absorb water like a sponge.Sponge cities break through the traditional concept of drainage in urban rainwater management, relying on the sunken lawn, grassed swales, green roof, permeable pavement, bioretention swales and other infrastructure, and take them as carriers, taking full account of the operational safety of urban infrastructure and urban water security in the absence of floods.Reasonable resource utilization of rainwater and maintenance of a good hydrological and ecological environment can be made at the same time as the comprehensive utilization of a variety of ecological technologies including penetration, retention, storage, purification, reuse, efflux and so on, which supplement the groundwater and regulate the water cycle.
Grassed swales are an important part of the sponge city, which is a landscape surface-drainage system planted with vegetation [7,8].Grassed swales allow surface runoff to be detained, filtered and to penetrate at a lower flow rate [9].As a new type of ecological measure, grassed swales can be used to collect road-surface runoff, and replace the gully, ditch or part of the rainwater pipe network.Grassed swales can not only absorb rainwater, but also pretreat rainwater with plant roots, which controls the runoff of contaminants from the source [10].If grassed swales are substituted for a traditional underground drainage system, the problem-that traditional sewage pipes and rainwater pipelines were connected randomly or improperly-can be solved.
Pervious asphalt pavement refers to a kind of asphalt pavement structure that is composed of a highly porous mixture that allows road surface water to enter the subgrade [11].Compared with the present urban roads that mainly have a dense type of pavement, the pervious asphalt pavement can effectively replenish the groundwater; the urban heat island effect would be alleviated, the peak flow during rainstorms and the pressure of the urban drainage system would be reduced, and the vehicle running noise would also be reduced significantly, while the safety and comfort of driving would be improved effectively [12].Since pervious asphalt pavement can significantly improve the ecological and environmental efficiency of the road, it has attracted widespread attention at home and abroad.
Since the design and construction of a sponge project are unlikely to occur simultaneously, and it is hard to adjust after implementation, it comes to be an inevitable trend that software models are used for simulating the effect of the sponge facilities in projects [13,14].Green infrastructure is increasingly being considered for application in urban storm-water management designs.Many municipalities, regulatory agencies and advocacy groups promote the use of low-impact development (LID) to reduce runoff and increase infiltration.To show LID benefits, engineers must quantify the advantages of green versus traditional grey infrastructure [15].Rainfall-runoff models can be classified into three types: physically based models, conceptual models and empirical models.In this latter class of models, the catchment is considered as a black box, without any reference to the internal processes that control the transformation of rainfall to runoff.In recent years, some models derived from studies on artificial intelligence have found increasing use [16].Currently, storm-water management model (SWMM) and Hydro CAD (http://www.hydrocad.net/company.htm)are the most-commonly used simulation software for LID, and the former is the most-widely used.A comprehensive hydrological model, like SWMM, has been widely used for rainfall-runoff simulation.In recent years, simple and effective modern modeling techniques have also brought great attention to the prediction of runoff with rainfall input [17].SWMM is a hydrodynamic rainfall-runoff simulation model developed by the US Environmental Protection Agency that includes a hydrology, hydraulic and water-quality module.It has functions of simulating urban rainfall-runoff processes (including the ground runoff, current in drainage systems, the process of flood regulation and storage), BOD (biochemical oxygen demand), COD (chemical oxygen demand), and the migration and diffusion of total phosphorus, total nitrogen and six other kinds of pollutants.The latest version, SWMM 5.1 (United States Environmental Protection Agency, Washington, DC, USA), set up eight LID technical facilities including grassed swales, roof closure facilities, rain water tanks, permeable pavements, infiltration ditchs, green roofs, rain gardens and bioretention units in the LID module [18].
At present, most of the research on sponge cities in China focus on municipal roads, residential quarters or urban area planning, but are rarely concerned with how to quantitatively analyze the sponge facilities involved in urban rainwater management.It is not difficult to make source reductions for permeable pavement and grassed swales because of the large area and permeability of the underlying surface; the difficulty is how to cooperate with the rainwater pipe network, import exogenous rainwater and perform the functions of halfway transmission, terminal regulation and storage [19].Nanshan village, Tianwang Town, Jurong, Jiangsu Province, is a demonstration point of the National Key Technology Support Program.In this paper, Nanshan village was selected as the research area, where the permeable pavement and grassed swales were arranged reasonably with local conditions in the overall planning of the village [20].The traditional development model and the LID model based on SWMM were established to simulate the runoff regulation effect in the study area under various rainfall repetition periods [21].Finally, comparison of the rainfall regulation effect between permeable pavement and grassed swales was done individually.The results of this study can provide guidance for the selection and arrangement of LID facilities in urban areas, and can also provide some technical support for the construction of sponge cities in China.

Situation of Research Area
Nanshan village lies in the north of Fushan Cherry Garden, adjacent to S340 provincial highway.The research area was named a national demonstration site in rural characteristic tourism and the design area is about four hectares.The village is located in a subtropical zone with subtropical monsoon climate with four distinctive seasons, and the annual average natural precipitation is 1188 mm, which is mainly concentrated from May to September.The plan of the research area is shown in Figure 1.
Water 2017, 9, 840 3 of 12 exogenous rainwater and perform the functions of halfway transmission, terminal regulation and storage [19].Nanshan village, Tianwang Town, Jurong, Jiangsu Province, is a demonstration point of the National Key Technology Support Program.In this paper, Nanshan village was selected as the research area, where the permeable pavement and grassed swales were arranged reasonably with local conditions in the overall planning of the village [20].The traditional development model and the LID model based on SWMM were established to simulate the runoff regulation effect in the study area under various rainfall repetition periods [21].Finally, comparison of the rainfall regulation effect between permeable pavement and grassed swales was done individually.The results of this study can provide guidance for the selection and arrangement of LID facilities in urban areas, and can also provide some technical support for the construction of sponge cities in China.

Situation of Research Area
Nanshan village lies in the north of Fushan Cherry Garden, adjacent to S340 provincial highway.The research area was named a national demonstration site in rural characteristic tourism and the design area is about four hectares.The village is located in a subtropical zone with subtropical monsoon climate with four distinctive seasons, and the annual average natural precipitation is 1188 mm, which is mainly concentrated from May to September.The plan of the research area is shown in Figure 1.

Sub-Catchment Generalization
The topography of the region, low in south and high in north, has a roughly ladder-like distribution, and the region is composed of 25% permeable area and 75% impervious area.The generalized area is divided into nine water catchments and one outfall.A generalized sketch map of sub-basin arrangement can be seen in Figure 2, where c1 to c9 refer to catchment and S1 to S9 mean pipe system.

Sub-Catchment Generalization
The topography of the region, low in south and high in north, has a roughly ladder-like distribution, and the region is composed of 25% permeable area and 75% impervious area.The generalized area is divided into nine water catchments and one outfall.A generalized sketch map of sub-basin arrangement can be seen in Figure 2, where c1 to c9 refer to catchment and S1 to S9 mean pipe system.In the LID model, permeable pavement and grass swales were arranged in the catchment.The area of grassed swales is 800 square meters, accounting for 4% of the total area, while the area of permeable pavement is 2450 square meters, accounting for 12.25% of the total area.After infiltration and treatment of grassed swale and permeable pavement, rainwater overflows into the municipal rainwater pipe network.The drainage path of runoff in the LID model is shown in Figure 3.In the LID model, permeable pavement and grass swales were arranged in the catchment.The area of grassed swales is 800 square meters, accounting for 4% of the total area, while the area of permeable pavement is 2450 square meters, accounting for 12.25% of the total area.After infiltration and treatment of grassed swale and permeable pavement, rainwater overflows into the municipal rainwater pipe network.The drainage path of runoff in the LID model is shown in Figure 3.In the LID model, permeable pavement and grass swales were arranged in the catchment.The area of grassed swales is 800 square meters, accounting for 4% of the total area, while the area of permeable pavement is 2450 square meters, accounting for 12.25% of the total area.After infiltration and treatment of grassed swale and permeable pavement, rainwater overflows into the municipal rainwater pipe network.The drainage path of runoff in the LID model is shown in Figure 3.In the traditional development model, there were no LID facilities.The drainage path of runoff in the traditional development model is shown in Figure 4.In the LID model, permeable pavement and grass swales were arranged in the catchment.The area of grassed swales is 800 square meters, accounting for 4% of the total area, while the area of permeable pavement is 2450 square meters, accounting for 12.25% of the total area.After infiltration and treatment of grassed swale and permeable pavement, rainwater overflows into the municipal rainwater pipe network.The drainage path of runoff in the LID model is shown in Figure 3.

Selection of Model Parameters
The Horton model, with maximum infiltration rate of 76.2 mm/h, minimum infiltration rate of 3.18 mm/h and attenuation constant of 4.12 per hour according to soil characteristics was adopted to simulate the rainfall infiltration [22].The nonlinear reservoir model was chosen to calculate the runoff generation at each grid while the equation of kinetic wave was adopted to calculate the pipeline transmission system [23].In addition, the determination of parameters, which is also determined by the field geotechnical exploration reports and previous settings on the underlying surface, refers to typical values in the user manual and related documents [24].Parameters of grass swale and permeable pavement are shown in Table 1.

Designed Rainfall
Design of rainfall hyetograph refers to the distribution of total rainfall in the designed phase, and it is the boundary condition for delay of simulation based on the hydraulic model; it is also the key module to calculate the designed flow of the drainage pipe network [25].
The Chicago rainfall hydrograph model is a rainfall process based on a rainstorm intensity Equation [26], and Equation (1) has been adopted as a widespread rainstorm intensity equation in China [27].Equation (1) can reflect the regular pattern of rainfall distribution over time: where q refers to the rainfall intensity of rainfall duration (L/s ha); p refers to the repetition period of designing rainfall (a); t d refers to the rainfall duration (min); and A 1 , C, b and n are local parameters, calculated by a statistical method.
Peak coefficient r is the only parameter of the Chicago rainfall hydrograph [29], and the hydrograph can be calculated by Equations ( 3) and (4).Accurate determination of peak coefficient r is the key to fit the curve of designed rainfall: where i before refers to the rainfall intensity before peak flow (mm/min); i after refers to the rainfall intensity after peak flow (mm/min); t b refers to the time before peak flow (min); t a refers to the time after peak flow (min); and r refers to the ratio of peak flow time and rainfall duration (peak coefficient).
In this research, the rainfall duration was set as 3 h and the peak coefficient was set as 0.34 [30].The distribution with time of rainfall under different rainfall repetition periods is shown in Table 2, and the rainfall hydrograph curve of various repeating rainfall periods when the rainfall peak coefficient is 0.34 is shown in Figure 5.In this research, the rainfall duration was set as 3 h and the peak coefficient was set as 0.34 [30].The distribution with time of rainfall under different rainfall repetition periods is shown in Table 2, and the rainfall hydrograph curve of various repeating rainfall periods when the rainfall peak coefficient is 0.34 is shown in Figure 5. Note: * a refers to the rainfall repetition period.

Comparison of LID and Traditional Development
The total runoff of outfall changed with the rainfall intensity under the LID model and the traditional development model, and the simulation results can be seen in Figure 6.The total rainfall here was the product of rainfall and catchment area, so its unit had changed.The rainfall repetition periods were 0.33a, 1a, 3a, 5a and 10a, the rainfall duration was 3 h and the simulation time was 4 h; total runoff, peak flow, runoff-yielding time and peak-flow occurrence time can be seen in Table 3.

Comparison of LID and Traditional Development
The total runoff of outfall changed with the rainfall intensity under the LID model and the traditional development model, and the simulation results can be seen in Figure 6.The total rainfall here was the product of rainfall and catchment area, so its unit had changed.The rainfall repetition periods were 0.33a, 1a, 3a, 5a and 10a, the rainfall duration was 3 h and the simulation time was 4 h; total runoff, peak flow, runoff-yielding time and peak-flow occurrence time can be seen in Table 3.According to Table 3, the total rainfalls over 3 h were 1080, 2052, 2865, 3247 and 3765 m 3 when the rainfall repetition periods were 0.33a, 1a, 3a, 5a and 10a, respectively.After the management effect of LID (i.e., permeable pavement + grassed swales), the outfall runoff was decreased in all groups.Especially, when the rainfall repetition period was 0.33a, the outfall had no flow.According to the comparison from Figure 6, it can be seen that the peak flow was reduced by 0.377, 0.286, 0.335, 0.227 and 0.119 m 3 /s when the rainfall repetition periods were 0.33a, 1a, 3a, 5a and 10a, respectively.With the comparative analysis of sub-catchment runoff control data, it is proposed that the runoff reduction in the traditional development model mainly depends on the infiltration of the surface permeable area, while the runoff reduction in the LID model mainly depends on reduction at the source, infiltration at halfway and regulation on the terminal.
Grassed swale and permeable pavement were conducive to reducing the peak flow under different rainfall repetition periods; however, as the rainfall repetition period and rainfall intensity increased, the runoff reduction efficiency decreased.When the rainfall repetition period was 10a, the difference in values of peak flow, runoff-yielding time and peak-flow occurrence time between the two models was negligible.The result is consistent with the simulation results of the effect of the rainwater garden's regulation and control of the urban rainwater runoff from Jiake Li et al. [31].Furthermore, in Li's design model, the rainwater garden accounts for 2% of the total land area, while the rainfall repetition periods were 2a, 5a, 10a and 20a, the total runoff reduction was more than 25.69% and the peak flow occurrence time was delayed by 5-7 min.
It can be seen from Figure 7 that when the rainfall intensity rises to a certain extent, the difference value of the total runoff between the two development models tends to be downward, indicating that all LID facilities have entered a saturated state, and the regulation and storage of the LID model are not obvious at this point.According to Table 3, the total rainfalls over 3 h were 1080, 2052, 2865, 3247 and 3765 m 3 when the rainfall repetition periods were 0.33a, 1a, 3a, 5a and 10a, respectively.After the management effect of LID (i.e., permeable pavement + grassed swales), the outfall runoff was decreased in all groups.Especially, when the rainfall repetition period was 0.33a, the outfall had no flow.According to the comparison from Figure 6, it can be seen that the peak flow was reduced by 0.377, 0.286, 0.335, 0.227 and 0.119 m 3 /s when the rainfall repetition periods were 0.33a, 1a, 3a, 5a and 10a, respectively.With the comparative analysis of sub-catchment runoff control data, it is proposed that the runoff reduction in the traditional development model mainly depends on the infiltration of the surface permeable area, while the runoff reduction in the LID model mainly depends on reduction at the source, infiltration at halfway and regulation on the terminal.
Grassed swale and permeable pavement were conducive to reducing the peak flow under different rainfall repetition periods; however, as the rainfall repetition period and rainfall intensity increased, the runoff reduction efficiency decreased.When the rainfall repetition period was 10a, the difference in values of peak flow, runoff-yielding time and peak-flow occurrence time between the two models was negligible.The result is consistent with the simulation results of the effect of the rainwater garden's regulation and control of the urban rainwater runoff from Jiake Li et al. [31].Furthermore, in Li's design model, the rainwater garden accounts for 2% of the total land area, while the rainfall repetition periods were 2a, 5a, 10a and 20a, the total runoff reduction was more than 25.69% and the peak flow occurrence time was delayed by 5-7 min.
It can be seen from Figure 7 that when the rainfall intensity rises to a certain extent, the difference value of the total runoff between the two development models tends to be downward, indicating that all LID facilities have entered a saturated state, and the regulation and storage of the LID model are not obvious at this point.

Comparison of Combined System and Single LID Facility
The rainfall repetition individual LID controls, one 5-year 3 h of rainfall event period was set as 5a while the rainfall duration was set as 3 h, the LID facilities (including grassed swales and permeable pavement) were simulated in the study area individually, and the simulation results were compared with the traditional development model and LID model; the final results can be seen in Table 4.It can be seen from Table 4 that the LID facilities have effect in reducing total runoff and peak flow, delaying the runoff-yielding time and peak-flow occurrence time.Particularly, the combined system (i.e., permeable pavement + grassed swales) worked more efficiently than the sum of the individuals in the reduction of total runoff and peak flow, and the delay of peak-flow occurrence time throughout.The reason is that after runoff from sub-catchment areas infiltrates through the grassed swales and permeable pavement surface, the water in storage layers of permeable pavement can infiltrate into the grassed swales through lateral infiltration.The grassed swales reduce the runoff velocity in the transmission of rainwater runoff, which results in the time delay or elimination in runoff-yielding.The connection of an internal infiltration path between the individual LID facilities increases the water-storage space of the sub-catchment area, which leads to a reduction in peak flow and delay in the occurrence time of peak flow.
When the values of unit area are the same, the efficiency of the LID facilities on the total runoff reduction is ranked as: combined system (i.e., permeable pavement + grassed swales) > permeable pavement > grassed swales.The efficiency of the peak-flow reduction is ranked as: combined system (i.e., permeable pavement + grassed swales) > permeable pavement > grassed swales.The efficiency of delay in peak-flow occurrence time is ranked as: combined system (i.e., permeable pavement + grassed swales) > permeable pavement > grassed swales.The reason is that the permeable pavement is composed of surface course, cushion course and porous base course, so that the runoff on the surface can infiltrate freely from top to bottom in the permeable pavement, and finally recharge groundwater.Per unit area, the permeable pavement has a larger reservoir space than grassed swales, the rainwater in the water-storage formation continues to infiltrate after rainfall, and finally the extra water overflows and discharges.The effect of delay in runoff-yielding time is ranked as: combined system (i.e., permeable pavement + grassed swales) > grassed swales > permeable pavement.The reason is that grassed swales have a larger Manning coefficient on the surface due to plantation coverage (the Manning coefficient on the surface is generally 0.3), which leads to the reduction in runoff velocity, and a slower runoff velocity would increase the transmission time of rainwater runoff from runoff-yield to outfall, resulting in the delay or elimination of runoff-yielding time.

Conclusions
This paper presented a research method on storm-water management using SWMM.As a part of this study, the superiority of LID in storm regulation was verified through comparison of the LID model and the traditional development model, and the contrast between the combined system (i.e., permeable pavement + grassed swales) and a single LID facility on the storm-water management effect was further investigated.Based on the results and analysis, the following conclusions were made: (1) The SWMM model can quantify the management effect of LID facilities on storm-water runoff.According to the nonlinear reservoir theory, the model took the extracted information of land surface as the basis to calculate runoff, hydrodynamic theory was used to compute concentration flow unit by unit, and the computational data can reflect the running effect of sponge facilities.In practical engineering, flow meters and level gauges are arranged in the infall and outfall to monitor actual treatment results, and the actual operating results of LID facilities are to be verified in the follow-up field tests for real-world authentication.
(2) Better results of storm-water management could be obtained through LID, including reduction in the total runoff, decrease in the peak flow, elimination or delay in the runoff-yielding time and delay in the peak-flow occurrence time.When the rainfall repetition period ranged from 0.33a to 10a, the reduction rate of total runoff in the research area where LID facilities were arranged ranged from 100% to 27.5%.The reduction rate of peak flow ranged from 100% to 15.9%.The runoff-yielding time was eliminated or delayed by 13 min to 5 min and the peak-flow occurrence time was eliminated or delayed by 17 min to 5 min after LID treatment.
(3) When the values of unit area are the same, the efficiency of the LID facilities on the total runoff reduction is ranked as: combined system (i.e., permeable pavement + grassed swales) > permeable pavement > grassed swales.The efficiency of the peak-flow reduction is ranked as: combined system (i.e., permeable pavement + grassed swales) > permeable pavement > grassed swales.The effect of delay in runoff-yielding time is ranked as: combined system (i.e., permeable pavement + grassed swales) > grassed swales > permeable pavement.The efficiency of delay in peak-flow occurrence time is ranked as: combined system (i.e., permeable pavement + grassed swales) > permeable pavement > grassed swales.The capacity for storm-water management of permeable pavement and grassed swales was limited, and the combined system was well-operated when the rain was light.However, when the area suffered too much rain in a short time, the management was mostly negligible.
(4) The results can provide reference for the selection and arrangement of LID facilities in urban areas; for example, when the regional space is ample, the united arrangement of LID facilities are recommended in the design of the drainage system, whereas when the area is limited, permeable pavement should be a top priority.Grassed swales and permeable pavement emphasize site-specific recommendations, intensive management, improved efficiency and environmentally sound use of inputs, rather than pursuit of the largest area of sinking space.

Figure 1 .
Figure 1.Plan of Nanshan village in Jurong.

Figure 1 .
Figure 1.Plan of Nanshan village in Jurong.

Figure 3 .
Figure 3. Drainage path of runoff in the LID model.

Figure 4 .
Figure 4. Drainage path of runoff in the traditional development model.

Figure 3 .
Figure 3. Drainage path of runoff in the LID model.In the traditional development model, there were no LID facilities.The drainage path of runoff in the traditional development model is shown in Figure 4.

Figure 4 .Figure 3 .
Figure 4. Drainage path of runoff in the traditional development model.

Figure 3 .
Figure 3. Drainage path of runoff in the LID model.In the traditional development model, there were no LID facilities.The drainage path of runoff in the traditional development model is shown in Figure 4.

Figure 4 .Figure 4 .
Figure 4. Drainage path of runoff in the traditional development model.

Figure 5 .
Figure 5. Rainfall hydrograph of various repeating rainfall periods when the rainfall peak coefficient is 0.34.

Figure 5 .
Figure 5. Rainfall hydrograph of various repeating rainfall periods when the rainfall peak coefficient is 0.34.

Figure 6 .
Figure 6.Comparison curve of runoff control effect: (a) when rainfall repetition period is 0.33a; (b) when rainfall repetition period is 1a; (c) when rainfall repetition period is 3a; (d) when rainfall repetition period is 5a; (e) when rainfall repetition period is 10a.

Table 1 .
Table of parameters of LID controls.

Table 2 .
Distribution of 3 h of rainfall under different rainfall repetition periods.

Table 2 .
Distribution of 3 h of rainfall under different rainfall repetition periods.

Table 3 .
Simulated results of the runoff under different rainfall intensities.

Table 3 .
Simulated results of the runoff under different rainfall intensities.

Table 4 .
Simulated results of runoff regulation effect under.