Performance Evaluation of Enhanced Bioretention Systems in Removing Dissolved Nutrients in Stormwater Runo ﬀ

: Bioretention has great potential in managing and purifying urban stormwater runo ﬀ . However, information regarding the removal of nutrients in bioretention systems with the use of media, plants, and saturated areas is still limited. In this study, three devices of control, conventional bioretention (DS), and strengthened bioretention (SZ) were investigated to enhance the simultaneous removal of nitrogen and phosphorus. The experimental column SZ showed the best performance for total phosphorus (TP), ammonia (NH 4 + -N) and total nitrogen (TN) removal (85.6–92.4%, 83.1–92.7%, 57.1–74.1%, respectively), whereas DS columns performed poorly for NH 4 + -N removal (43.6–81.2%) under di ﬀ erent conditions. For the removal of nitrate, the columns of Control and DS exhibited negative performance ( − 14.3% and − 8.2%) in a typical event. Further evaluation of water quality revealed that in the early stages of rainfall, the e ﬄ uent of the SZ column was able to reach quality standards of Grade IV for surface water in China. Moreover, although the ion-exchange and phosphate precipitation occurred on the surface of the media, which were placed in the saturation zone, it did not change the surface crystal structure.


Introduction
With the growth of urbanization, natural green spaces have been gradually replaced by impervious pavements, resulting in reduced nitrogen retention capacity in many urban areas [1][2][3]. Urban rainfall-runoff often carries high concentrations of nitrogen, phosphorus, nutrients, organic pollutants, and heavy metals [4][5][6]. Studies have shown that half of the nitrogen pollutants in urban rivers come from surface runoff formed during urban rainfall and the concentrations of metal elements in some runoff can exceed more than ten times the limitations for surface water quality [7,8]. Nitrogen, heavy metals, and phosphorus present in the surface runoff lead to the functional degradation of about two-thirds of urban rivers, and the pollutants discharged from this surface runoff present a serious threat to the urban water system [9][10][11][12]. In order to control the amount of nitrogen chemical oxygen demand (COD) and phosphorus in stormwater runoff and prevent the further deterioration of the urban water ecosystem, decentralized management of rainwater becomes critically important [13]. Many countries have adopted different rainwater management systems based on the existing urban intensities, and concentrations; (c) fuzzy mathematics was used to evaluate the effluent quality of bioretention systems in regular rainfall events; and, (d) static adsorption, scanning electron microscope (SEM), energy-dispersive X-ray spectroscopy (EDS), and X-ray diffraction (XRD) were combined to determine the pollutant removal mechanisms of different combinations of the bioretention system. On this basis, an optimization strategy is proposed for the simultaneous removal of nitrogen and phosphorus in the bioretention system. The results provide a reliable and scientific reference for improving the performance of the bioretention system and removing organic matter, nitrogen, and phosphorus from the initial road runoff.

Set-Up of the Bioretention Columns
The simulated bioretention system includes a peristaltic pump, a water intake bucket, and a bioretention system. The bioretention system is composed of a plexiglass cylinder and a polyethylene pipe with an inner diameter of 10 mm and a height of 80 cm. On 24 May 2019, three setups of simulated bioretention systems were constructed, including the control group, the conventional bioretention (DS) group, and the strengthened bioretention (SZ) group ( Figure 1). The control group and the DS group drained water from the bottom. The SZ group drained water 40 cm above the base, and the 40 cm column was in a long-term saturated state. Since phosphorus is easily captured by the physical and chemical adsorption of media during infiltration, only a portion of the nitrogen can be absorbed by the medium, and the rest is removed by nitrification and denitrification. Based on the above theory, the media layer was arranged from top to bottom as the water-retaining layer, adsorption layer, and catchment layer. Landscape plant Ophiopogon japonicus was selected as the model plant for bioretention. The fibrous root of Ophiopogon japonicus is well developed and has a fast growth rate. It can adapt to periodic changes on the surface runoff. In the control group, no plants were grown, while in the DS and SZ groups, Ophiopogon japonicus was planted. The specific configuration is shown in Table 1.

Synthetic Stormwater Runoff
The concentration of phosphorus and nitrogen in the simulated urban surface runoff refers to the monitoring results of the surface runoff in Shanghai, Nanjing, Wuhan, and Zhenjiang in China [34][35][36]. The rates of removal of DP (PO 4 3− ) and dissolved nitrogen (NH 4 + -N and NO 3 − -N) were analyzed. In order to determine the performance of bioretention in the case of excessive pollutants in contaminated Appl. Sci. 2020, 10, 3148 4 of 20 stormwater, two different concentrations of synthetic water runoff were used in the study, as presented in Table 2. The concentration of pollutants in the synthetic runoff for the regular event was designed to be similar to stormwater runoff, as suggested in the literature [31,37]. The experimental conditions for the extreme event were used to evaluate the performance of the bioretention in the case of sudden high concentrations in stormwater. Filter layer A non-woven permeable geotextile membrane / Catchment layer Crushed limestone gravel (grain size of 5-10 mm, void space of 34%) 5 cm Table 2. Two types of concentrations for synthetic runoff used in this study.

Experimental Design
After planting the Ophiopogon japonicus, the plants were watered every three days. The concentrations of TP and TN were 0.062 ± 0.021 mg/L and 0.87 ± 0.26 mg/L, respectively. After two months of cultivation as the plants gradually grew, the bioretention system was tested and monitored in July 2018. The removal efficiency of pollutants was investigated based on rainfall intensity, sunny weather, and concentration.
Based on the empirical model of rainfall intensity for Nanjing (Jiangsu Province, China) in Equation (1), we simulated the rainfall patterns using the Chicago Storm Theory, as provided by Equations (2) and (3) [30]. The rainfall duration (t) was assumed to be 360 min with rain peak using a rain peak coefficient of 0.4. Three different recurrence intervals in years (1a, 5a, 10a) were used, as presented in Table 3 and Figure 2.  Next, we calculated the rainwater runoffs associated with different recurrence intervals, based on the drainage criterion (GB50014-2006) in China. The rainwater design flow rate was calculated using Equation (4). (4) where Qs (mL/min) is the design flow rate of the rainwater, φ is the runoff coefficient (assumed to be 0.3 for soil road such as in a suburb), and F (cm 2 ) is the percentage of the catchment area (assumed to be 8% as the ratio of the bioretention area to the catchment area).
In Nanjing, 60% of rainfall is not more than 25.0 mm, daily rainfall is not less than 5.0 mm, and the average interval time between rainfall is seven days. The simulation interval days were set at three levels (i.e., 3, 8, and 14 days), which were used in analyzing the effect of sunny weather on the rate of removal of nitrogen and phosphorus.

Evaluation of Removal Efficiency by Bioretention
Water sample testing was based on the Standard of Monitoring and Analysis Method of Water and Wastewater (Fourth Edition). The water sample was filtered using a 0.45 µm membrane, and subsequently, ultraviolet spectrophotometry was performed to determine the TN. Phenol disulfonic Before Peak: After Peak: where i (mm/min) is the rainfall intensity, t (min) is the rainfall duration, p is the recurrence interval, and A, b, and n are the constants dependent on the units employed and the recurrence interval of the storm, r is the coefficient of rain peak. Next, we calculated the rainwater runoffs associated with different recurrence intervals, based on the drainage criterion (GB50014-2006) in China. The rainwater design flow rate was calculated using Equation (4).
where Q s (mL/min) is the design flow rate of the rainwater, ϕ is the runoff coefficient (assumed to be 0.3 for soil road such as in a suburb), and F (cm 2 ) is the percentage of the catchment area (assumed to be 8% as the ratio of the bioretention area to the catchment area). In Nanjing, 60% of rainfall is not more than 25.0 mm, daily rainfall is not less than 5.0 mm, and the average interval time between rainfall is seven days. The simulation interval days were set at three levels (i.e., 3, 8, and 14 days), which were used in analyzing the effect of sunny weather on the rate of removal of nitrogen and phosphorus.

Evaluation of Removal Efficiency by Bioretention
Water sample testing was based on the Standard of Monitoring and Analysis Method of Water and Wastewater (Fourth Edition). The water sample was filtered using a 0.45 µm membrane, and subsequently, ultraviolet spectrophotometry was performed to determine the TN. Phenol disulfonic acid spectrophotometry was employed to determine the NO 3 − -N content. Nessler's reagent colorimetric method was used to analyze the NH 4 + -N content. TP was determined using potassium persulfate oxidation and the molybdenum antimony anti-spectrophotometric method. During the test, the cumulative load of pollutants in the bioretention system was calculated as follows: Appl. Sci. 2020, 10, 3148 6 of 20 The removal ratio for different forms of pollutants was calculated by Equation (5): where η (%) is the removal ratio, C 0 (mg/L) is the concentrations in the runoff before entering the bioretention system, and C e (mg/L) is the concentration in the outflow.

Surface Characterization and Chemical Compositions
To evaluate the changes in morphology, crystal structure, and the chemical compositions of the surface, the materials in the bioretention before and after the simulated rainfall tests were characterized using the SEM, EDS, and XRD. The compositions and distribution of elements on the surface of the substrate were determined by EDS analysis. Before all analytical tests, the materials were dried in an oven at 105 • C for 24 h.

Water Quality Assessment Based on Fuzzy Water Pollution Index Method
Comprehensive evaluation theory is a method evaluated by fuzzy mathematics. Using the fuzzy theory, the method from a number of factors undertakes a comprehensive evaluation according to subordinate degrees of arbitrated aspects. The fuzzy evaluation method can rank and evaluate objects by the size of the overall evaluation scores and also divide them into different grades [16]. This study uses an analytical decision-making method based on the fuzzy congruity theory. By using the fuzzy theory, some qualitative factors are transformed into quantitative parameters, and the fuzziness and uncertainty of determination are solved [38,39].

Establishing the Model of the Fuzzy Comprehensive Index (FCI) Method
The procedure for constructing the comprehensive evaluation model is as follows [40,41]: (1) Identify the factor set U = {u 1 , u 2 , ···, u n }. Assume that the number of factors related to the evaluated object is n. (2) Identify the remark set V = {v 1 , v 2 , ···, v n }. Assume that the number of all the possible remarks is m.
(3) Process the single factor evaluation, which means evaluating every factor u i {i = 1, 2, ···, n}, and getting the single factor evaluation set r i = {r i1 , r i2 , ···, r im }. Determine the evaluation matrix by the mapping relationship.
(4) Since every factor's status is not necessarily equal, there is a need for weighting factors. Assume that the fuzzy set A = {a 1 , a 2 , ···, a n } represents the weight distribution of every factor and n. The comprehensive evaluation set is: 2.6.2. Calculation of the FCI of Water Quality 1. Fuzzy Relation Matrix R The fuzzy relation matrix represents the subjection degree of each evaluation factor to each water quality grade. In this study, the subjection degree matrix R is calculated by using the simplified half-trapezoidal distribution function and the linear triangle function [42]. The surface water is categorized into six grades (I, II, III, IV, V, worse than V), which is based on the "Environmental quality standard for surface water". The standard values for each grade are Q(1), Q(2), Q(3), Q(4), Q(5), Q(6). The membership pollution factor r ij can be defined as follows in Equation (6): where x i represents the measured value of the i pollution factor; Q (j) represents the j level quality standard of the i pollution factor; Q (j+1) represents the j+1 level quality standard of the i pollution factor.

Weight Coefficient
The weight coefficient of pollution factors is used to measure the influence of pollution factors on environmental water quality. In this study, the weight value is determined by calculating the excess ratio. The higher the weight coefficient, the stronger the influence of the factor on water quality. The formula is as follows in Equation (7): where C i is the measured value of the pollution factor i; S i is the arithmetic mean value of environmental quality standards for pollution factors at all levels in i, such that in Equation (8): 3. Fuzzy Composite Index B* In this study, we use the weighted average principle to calculate the Fuzzy composite index B* in Equation (9): where B* is the relative position of the object in each quality grade; b j is the membership of grade j; and, k is an undetermined coefficient which is used to control the value of b j .

Effect of Rainfall Intensity on the Retention of Pollutants by Bioretention
All three groups with bioretention systems had high rates of removal of TP for the whole rainfall duration (Figure 3), indicating that the bioretention technology is able to effectively reduce the content of phosphorus in the rainfall-runoff. The conventional bioretention system (control group) had the lowest removal efficiency of TP, while the DS and SZ columns remained less than 0.4 mg/L of TP in the initial effluent and showed a higher purification effect. This could be attributed to the plant uptake, media adsorption, and chemical precipitation that occurred within the system [30]. With the increase in rainfall time, the concentration of phosphorus in the rainwater outflow increased gradually. One possible explanation is that the initial fast-reversible adsorption sites on the media surface were gradually occupied, and then the phosphate in the rainwater moved slowly towards the irreversible adsorption sites during infiltration [43]. At the same time, the released adsorption reversible sites could regain the adsorptive capacity of the media for the bioretention system before the beginning of the next rainfall.  The removal rate of TP decreased with increasing rainfall time. The removal rate of TP in the SZ column was still highest under these three different rainfall intensities. When the rainfall recurrence interval was 1a, the bioretention columns with aquifers (SZ) was able to achieve a removal rate of more than 92.4% in the initial effluent. However, with an increase in rainfall intensity, the removal rate of TP in the effluent had decreased gradually. The removal of TP was mainly due to the physical and chemical adsorption of the media. At the same time, higher rainfall intensity would have led to more considerable runoff and higher total pollution load. The hydraulic retention time of rainwater in the bioretention system was shortened, and the phosphorus adsorption by the fillers with limited adsorption capacity was easily saturated. The traditional bioretention system (no saturated area) was affected by the rainfall intensity. When the recurrence interval was increased to 10a, the removal rate of TP decreased to 72.6%. Therefore, even under extreme weather conditions, the bioretention system with a saturated zone proved to be more reliable and stable.
Under different rainfall intensities, the removal efficiency of NH4 + -N was relatively different ( Figure 4). With continued rainfall, the concentration of NH4 + -N in the effluent had gradually increased, with the SZ column exhibiting the highest removal efficiency. When the rainfall recurrence interval was 1a, the removal rate of NH4 + -N reached 90.1%, and NH4 + -N concentration in the effluent was lower than 0.85 mg/L by the end of the rainfall. This is because higher rainfall intensities lead to greater runoffs and higher hydraulic loads. Rainwater runoff could quickly pass through the planting layer, which decreases the hydraulic retention time of rainwater runoff in the planting layer and reduces the fast adsorption and filtration of NH4 + -N by the plant roots and soil layers [30,35]. At the same time, the reaction time between nitrifying bacteria and NH4 + -N in the rainwater runoff is reduced, which subsequently diminishes the treatment efficiency of the bioretention system for the rainwater runoff.   The removal rate of TP decreased with increasing rainfall time. The removal rate of TP in the SZ column was still highest under these three different rainfall intensities. When the rainfall recurrence interval was 1a, the bioretention columns with aquifers (SZ) was able to achieve a removal rate of more than 92.4% in the initial effluent. However, with an increase in rainfall intensity, the removal rate of TP in the effluent had decreased gradually. The removal of TP was mainly due to the physical and chemical adsorption of the media. At the same time, higher rainfall intensity would have led to more considerable runoff and higher total pollution load. The hydraulic retention time of rainwater in the bioretention system was shortened, and the phosphorus adsorption by the fillers with limited adsorption capacity was easily saturated. The traditional bioretention system (no saturated area) was affected by the rainfall intensity. When the recurrence interval was increased to 10a, the removal rate of TP decreased to 72.6%. Therefore, even under extreme weather conditions, the bioretention system with a saturated zone proved to be more reliable and stable.
Under different rainfall intensities, the removal efficiency of NH 4 + -N was relatively different ( Figure 4). With continued rainfall, the concentration of NH 4 + -N in the effluent had gradually increased, with the SZ column exhibiting the highest removal efficiency. When the rainfall recurrence interval was 1a, the removal rate of NH 4 + -N reached 90.1%, and NH 4 + -N concentration in the effluent was lower than 0.85 mg/L by the end of the rainfall. This is because higher rainfall intensities lead to greater runoffs and higher hydraulic loads. Rainwater runoff could quickly pass through the planting layer, which decreases the hydraulic retention time of rainwater runoff in the planting layer and reduces the fast adsorption and filtration of NH 4 + -N by the plant roots and soil layers [30,35]. At the same time, the reaction time between nitrifying bacteria and NH 4 + -N in the rainwater runoff is reduced, which subsequently diminishes the treatment efficiency of the bioretention system for the rainwater runoff. The removal rate of TP decreased with increasing rainfall time. The removal rate of TP in the SZ column was still highest under these three different rainfall intensities. When the rainfall recurrence interval was 1a, the bioretention columns with aquifers (SZ) was able to achieve a removal rate of more than 92.4% in the initial effluent. However, with an increase in rainfall intensity, the removal rate of TP in the effluent had decreased gradually. The removal of TP was mainly due to the physical and chemical adsorption of the media. At the same time, higher rainfall intensity would have led to more considerable runoff and higher total pollution load. The hydraulic retention time of rainwater in the bioretention system was shortened, and the phosphorus adsorption by the fillers with limited adsorption capacity was easily saturated. The traditional bioretention system (no saturated area) was affected by the rainfall intensity. When the recurrence interval was increased to 10a, the removal rate of TP decreased to 72.6%. Therefore, even under extreme weather conditions, the bioretention system with a saturated zone proved to be more reliable and stable.
Under different rainfall intensities, the removal efficiency of NH4 + -N was relatively different ( Figure 4). With continued rainfall, the concentration of NH4 + -N in the effluent had gradually increased, with the SZ column exhibiting the highest removal efficiency. When the rainfall recurrence interval was 1a, the removal rate of NH4 + -N reached 90.1%, and NH4 + -N concentration in the effluent was lower than 0.85 mg/L by the end of the rainfall. This is because higher rainfall intensities lead to greater runoffs and higher hydraulic loads. Rainwater runoff could quickly pass through the planting layer, which decreases the hydraulic retention time of rainwater runoff in the planting layer and reduces the fast adsorption and filtration of NH4 + -N by the plant roots and soil layers [30,35]. At the same time, the reaction time between nitrifying bacteria and NH4 + -N in the rainwater runoff is reduced, which subsequently diminishes the treatment efficiency of the bioretention system for the rainwater runoff.   rate of removal of NO 3 − -N increased from −19.3% to 18.9%, and the bioretention system (SZ) with a saturated zone showed a higher rate of removal of TN. This is because forming a submerged area at the bottom of the retention system would be difficult with small runoffs. There were many aerobic zones in the system. The nitrification process converts NH 4 + -N into NO 3 − -N. However, the system cannot provide suitable anaerobic conditions, thereby limiting the process of denitrification [29,32]. As a result, the rate of removal of NO 3 − -N was extremely low, and even the concentration of the effluent was higher than that of the influent. Therefore, further studies are required to improve the removal efficiency of NO 3 − -N in the bioretention system.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 20 Figure 5 shows the operational effect of the bioretention system on NO3 − -N under different rainfall intensities. The concentration of NO3 − -N in the outflow declined with increasing sampling time. The rate of removal of NO3 − -N increased from −19.3% to 18.9%, and the bioretention system (SZ) with a saturated zone showed a higher rate of removal of TN. This is because forming a submerged area at the bottom of the retention system would be difficult with small runoffs. There were many aerobic zones in the system. The nitrification process converts NH4 + -N into NO3 − -N. However, the system cannot provide suitable anaerobic conditions, thereby limiting the process of denitrification [29,32]. As a result, the rate of removal of NO3 − -N was extremely low, and even the concentration of the effluent was higher than that of the influent. Therefore, further studies are required to improve the removal efficiency of NO3 − -N in the bioretention system.  Figure 6 shows TN removal in the different types of bioretention systems with the same rainfall duration. The effect of TN treatment was similar to that of NO3 -N, with the change in rainfall, showing a positive correlation trend with a 15% difference in the rate of removal. Regardless of the change in rainfall intensity, the SZ column with a saturation function had the highest rate of TN removal (more than 40%). When the rainfall recurrence interval was 10a, the TN removal reached 60%. This enhanced function may be attributed to the denitrification process in the aquifer under anaerobic conditions. NO3 − -N was denitrified into nitrogen in the deep anaerobic zone and removed by the bioretention system [24].

Effect of Antecedent Dry Periods on Pollutants Retention in Bioretention Systems
The results suggest that dry periods before rainfall have an important influence on the accumulation, migration, and transformation of nitrogen in the bioretention system. Overall, during the rainfall event, the effluent nitrogen concentration in the different bioretention systems showed a  a positive correlation trend with a 15% difference in the rate of removal. Regardless of the change in rainfall intensity, the SZ column with a saturation function had the highest rate of TN removal (more than 40%). When the rainfall recurrence interval was 10a, the TN removal reached 60%. This enhanced function may be attributed to the denitrification process in the aquifer under anaerobic conditions. NO 3 − -N was denitrified into nitrogen in the deep anaerobic zone and removed by the bioretention system [24].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 20 Figure 5 shows the operational effect of the bioretention system on NO3 − -N under different rainfall intensities. The concentration of NO3 − -N in the outflow declined with increasing sampling time. The rate of removal of NO3 − -N increased from −19.3% to 18.9%, and the bioretention system (SZ) with a saturated zone showed a higher rate of removal of TN. This is because forming a submerged area at the bottom of the retention system would be difficult with small runoffs. There were many aerobic zones in the system. The nitrification process converts NH4 + -N into NO3 − -N. However, the system cannot provide suitable anaerobic conditions, thereby limiting the process of denitrification [29,32]. As a result, the rate of removal of NO3 − -N was extremely low, and even the concentration of the effluent was higher than that of the influent. Therefore, further studies are required to improve the removal efficiency of NO3 − -N in the bioretention system.  Figure 6 shows TN removal in the different types of bioretention systems with the same rainfall duration. The effect of TN treatment was similar to that of NO3 -N, with the change in rainfall, showing a positive correlation trend with a 15% difference in the rate of removal. Regardless of the change in rainfall intensity, the SZ column with a saturation function had the highest rate of TN removal (more than 40%). When the rainfall recurrence interval was 10a, the TN removal reached 60%. This enhanced function may be attributed to the denitrification process in the aquifer under anaerobic conditions. NO3 − -N was denitrified into nitrogen in the deep anaerobic zone and removed by the bioretention system [24].

Effect of Antecedent Dry Periods on Pollutants Retention in Bioretention Systems
The results suggest that dry periods before rainfall have an important influence on the accumulation, migration, and transformation of nitrogen in the bioretention system. Overall, during the rainfall event, the effluent nitrogen concentration in the different bioretention systems showed a

Effect of Antecedent Dry Periods on Pollutants Retention in Bioretention Systems
The results suggest that dry periods before rainfall have an important influence on the accumulation, migration, and transformation of nitrogen in the bioretention system. Overall, during the rainfall event, the effluent nitrogen concentration in the different bioretention systems showed a downward trend, followed by an upward drift. The difference in the effluent nitrogen concentrations of the different structures was found insignificant (Figure 7). The number of adsorption sites vacated by the media was not significantly different. The initial effluent concentration of NH 4 + -N in the SZ group was higher than that in the later stage. The main explanation is that the submerged layer is usually anoxic, which is unconducive to aerobic nitrification [24,29,32]. The effluent NH 4 + -N concentration in the bioretention system was lower than 0.2 mg/L in the experiment, which showed that the bioretention system had higher removal efficiency of NH 4 + -N in the rainwater runoff.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 20 downward trend, followed by an upward drift. The difference in the effluent nitrogen concentrations of the different structures was found insignificant (Figure 7). The number of adsorption sites vacated by the media was not significantly different. The initial effluent concentration of NH4 + -N in the SZ group was higher than that in the later stage. The main explanation is that the submerged layer is usually anoxic, which is unconducive to aerobic nitrification [24,29,32]. The effluent NH4 + -N concentration in the bioretention system was lower than 0.2 mg/L in the experiment, which showed that the bioretention system had higher removal efficiency of NH4 + -N in the rainwater runoff.   The effluent TN concentration also showed an upward trend. The difference in the effluent TN concentration in the different bioretention systems was not significant. The overall change in the trend was similar to that of the NO3 − -N. The removal capacity of NH4 + -N by the filter layer in the different systems exhibited minor differences. The longer the antecedent dry period, the more NO3 − -N was removed by plants or by abiotic fixation. Given the combined effect of these two  Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 20 downward trend, followed by an upward drift. The difference in the effluent nitrogen concentrations of the different structures was found insignificant (Figure 7). The number of adsorption sites vacated by the media was not significantly different. The initial effluent concentration of NH4 + -N in the SZ group was higher than that in the later stage. The main explanation is that the submerged layer is usually anoxic, which is unconducive to aerobic nitrification [24,29,32]. The effluent NH4 + -N concentration in the bioretention system was lower than 0.2 mg/L in the experiment, which showed that the bioretention system had higher removal efficiency of NH4 + -N in the rainwater runoff.   The effluent TN concentration also showed an upward trend. The difference in the effluent TN concentration in the different bioretention systems was not significant. The overall change in the trend was similar to that of the NO3 − -N. The removal capacity of NH4 + -N by the filter layer in the different systems exhibited minor differences. The longer the antecedent dry period, the more NO3 − -N was removed by plants or by abiotic fixation. Given the combined effect of these two The effluent TN concentration also showed an upward trend. The difference in the effluent TN concentration in the different bioretention systems was not significant. The overall change in the trend was similar to that of the NO 3 − -N. The removal capacity of NH 4 + -N by the filter layer in the different systems exhibited minor differences. The longer the antecedent dry period, the more NO 3 − -N was removed by plants or by abiotic fixation. Given the combined effect of these two measures, the difference in TN removal by the different systems was found to be minimal. Although a longer dry period means that more denitrification could be carried out and would result in a lower NO 3 − -N content, the results show that the difference in the amount of NO 3 − -N between the different systems was not significant (Figure 9). This means that the comprehensive effect of the filter layer and the aquifer is able to reduce the difference in TN removal rates [30].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 11 of 20 measures, the difference in TN removal by the different systems was found to be minimal. Although a longer dry period means that more denitrification could be carried out and would result in a lower NO3 − -N content, the results show that the difference in the amount of NO3 − -N between the different systems was not significant (Figure 9). This means that the comprehensive effect of the filter layer and the aquifer is able to reduce the difference in TN removal rates [30].  Figure 10 shows the relationship between the dry period duration and TP removal. The TP concentration increased with a longer dry period. When the dry period was three days, the TP removal efficiency rates of the three systems was between 84.7% to 91.4%, and the effluent TP concentration was lower than 0.2 mg/L. From the 8 th to the 14 th day, the removal efficiency of TP decreased from 86.2% to 71.7%. The main reason is that the media adsorption serves as the primary mode of phosphate removal. Since high phosphorus content is intercepted by the infiltration layer and the media layer, the phosphorus in the system is slowly released outward with a prolonged dry period. At the same time, the anoxic state in the submerged area steers the phosphorus-accumulating bacteria to again release the fixed phosphorus [27,30]. Thus, in the next rainfall, the water in the submerged area would flow out.   Figure 10 shows the relationship between the dry period duration and TP removal. The TP concentration increased with a longer dry period. When the dry period was three days, the TP removal efficiency rates of the three systems was between 84.7% to 91.4%, and the effluent TP concentration was lower than 0.2 mg/L. From the 8th to the 14th day, the removal efficiency of TP decreased from 86.2% to 71.7%. The main reason is that the media adsorption serves as the primary mode of phosphate removal. Since high phosphorus content is intercepted by the infiltration layer and the media layer, the phosphorus in the system is slowly released outward with a prolonged dry period. At the same time, the anoxic state in the submerged area steers the phosphorus-accumulating bacteria to again release the fixed phosphorus [27,30]. Thus, in the next rainfall, the water in the submerged area would flow out.

Effect of Nutrient Concentration on Pollutant Retention by Bioretention Systems
Appl. Sci. 2020, 10, x FOR PEER REVIEW 11 of 20 measures, the difference in TN removal by the different systems was found to be minimal. Although a longer dry period means that more denitrification could be carried out and would result in a lower NO3 − -N content, the results show that the difference in the amount of NO3 − -N between the different systems was not significant (Figure 9). This means that the comprehensive effect of the filter layer and the aquifer is able to reduce the difference in TN removal rates [30].  Figure 10 shows the relationship between the dry period duration and TP removal. The TP concentration increased with a longer dry period. When the dry period was three days, the TP removal efficiency rates of the three systems was between 84.7% to 91.4%, and the effluent TP concentration was lower than 0.2 mg/L. From the 8 th to the 14 th day, the removal efficiency of TP decreased from 86.2% to 71.7%. The main reason is that the media adsorption serves as the primary mode of phosphate removal. Since high phosphorus content is intercepted by the infiltration layer and the media layer, the phosphorus in the system is slowly released outward with a prolonged dry period. At the same time, the anoxic state in the submerged area steers the phosphorus-accumulating bacteria to again release the fixed phosphorus [27,30]. Thus, in the next rainfall, the water in the submerged area would flow out.  Figure 10. The efficiency of TP removal using bioretention systems with the Control column, DS column, and SZ column during different antecedent dry periods. Figure 10. The efficiency of TP removal using bioretention systems with the Control column, DS column, and SZ column during different antecedent dry periods. Figure 11a,b shows the TP and NH 4 + -N removal efficiency rates of bioretention systems for normal and extreme events with the recurrence interval of the rainfall being 1a. In general, influent stormwater with higher nutrient concentrations would consequently result in higher concentrations in the bioretention system's effluent [44]. In this study, the TP removal efficiency was found to increase in higher concentrations (extreme event) for both the Control and SZ columns. A higher phosphorus load means a larger concentration gradient, which in turn can promote the adsorption process between pollutants and media in the solution given the greater availability of chelate sites on the media.

Effect of Nutrient Concentration on Pollutant Retention by Bioretention Systems
concentration. This suggests that the bioretention system with a saturated area could still possess sufficient treatment capacity to handle the excessive concentration of nutrients, at least within a short period of time. Figure 11c,d compares the NH4 + -N and TN removal efficiency rates of different bioretention systems for standard and extreme concentrations. For NH4 + -N, the removal efficiency in the high nutrient concentration (extreme event) was higher than that of a normal event, especially for the column with the saturated zone (>85%). Similar observations were found for the TN removal in extreme events, where the columns with the saturated zone exhibited a higher removal efficiency of over 55%. According to the study by Zhao et al. [46], the removal ratio of TN in inflows with high nutrient concentrations was about 40%, which was higher compared to low nutrient inflows. Note that the improved bioretention in the SZ column (combines the saturated zone) could enable the simultaneous removal of phosphorus and nitrogen from stormwater runoff.  The extent of phosphorus removal is markedly influenced by the ratio of phosphorus in the solution (i.e., DP) and the adsorbed phosphorus [45]. Similarly, for the NO 3 − -N removal, all columns exhibited better removal efficiency for high nutrient concentration compared with the normal concentration. This suggests that the bioretention system with a saturated area could still possess sufficient treatment capacity to handle the excessive concentration of nutrients, at least within a short period of time. Figure 11c,d compares the NH 4 + -N and TN removal efficiency rates of different bioretention systems for standard and extreme concentrations. For NH 4 + -N, the removal efficiency in the high nutrient concentration (extreme event) was higher than that of a normal event, especially for the column with the saturated zone (>85%). Similar observations were found for the TN removal in extreme events, where the columns with the saturated zone exhibited a higher removal efficiency of over 55%. According to the study by Zhao et al. [46], the removal ratio of TN in inflows with high nutrient concentrations was about 40%, which was higher compared to low nutrient inflows. Note that the improved bioretention in the SZ column (combines the saturated zone) could enable the simultaneous removal of phosphorus and nitrogen from stormwater runoff.

Evaluation of Water Quality in Bioretention Systems
The calculation results of membership and weight coefficient are shown in Tables 4 and 5. In the three different bioretention systems, the TN had the highest weighted value in terms of the pollution factor. Focusing on the TN values, the weighted value of bioretention with the saturated zone was found to be higher than the other two, with the weight of the initial outflow above 0.6. This means that dealing with TN in bioretention systems is more complicated. As the rainfall extends, the weights of Ammonia and Nitrogen in the control and DS groups are increased, while the TN weight decreases. However, in the SZ group, the weights of Ammonia and Nitrogen are stable, which suggests that in the precipitation process, bioretention is relatively stable given the saturated zone's ability to purify rainwater.   Figure 12 shows the effluent's water quality evaluation grade in the different bioretention systems. With the increase in rainfall time, the bioretention systems' effluent water quality was found to have a declining trend. The order of effluent water quality for the complete rainfall process is SZ > DS > Control. In the first three hours, the B* of SZ was 3.57, and the water quality reached the Grade IV standard for surface water, which is higher than the effluent water quality in the Control (below Grade V) and in the DS (Grade V). By the 4th hour of rainfall, the effluent water quality of all the three detention ponds was below Grade V. Although the water quality evaluation result of SZ was below the Grade V level, the FCI B* was measured at 5.401, which is close to the water quality standard of Grade V. This suggests that, in general, the SZ system can effectively control the amount of ammonia-nitrogen, TN, and TP in stormwater runoff. The application of the FCI method in evaluating water quality in bioretention systems is able to show the changes in effluent water quality during a rainfall event and determine the primary components of water pollution at specific time periods.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 14 of 20  Figure 12 shows the effluent's water quality evaluation grade in the different bioretention systems. With the increase in rainfall time, the bioretention systems' effluent water quality was found to have a declining trend. The order of effluent water quality for the complete rainfall process is SZ > DS > Control. In the first three hours, the B* of SZ was 3.57, and the water quality reached the Grade IV standard for surface water, which is higher than the effluent water quality in the Control (below Grade V) and in the DS (Grade V). By the 4 th hour of rainfall, the effluent water quality of all the three detention ponds was below Grade V. Although the water quality evaluation result of SZ was below the Grade V level, the FCI B* was measured at 5.401, which is close to the water quality standard of Grade V. This suggests that, in general, the SZ system can effectively control the amount of ammonia-nitrogen, TN, and TP in stormwater runoff. The application of the FCI method in evaluating water quality in bioretention systems is able to show the changes in effluent water quality during a rainfall event and determine the primary components of water pollution at specific time periods.

Study on the Adsorption Mechanism
SEM was performed to further explore the micro-morphology of the media and its changes in the operations of the bioretention systems (Figure 13a-d). Zeolite has a rough surface, folded, and a compact micropore structure. The micropore structure of 1-5 µm is distributed all over the media surface. It can provide an ideal environment for the growth of the biofilm and can improve the bio-load and ion exchange performance. Figure 12b-d are the surfaces of post-experiment zeolite, sponge iron, and limestone, respectively. Compared to the original zeolite (Figure 13a), the zeolite used after the experiment exhibited a rough surface and more crystals. Polyhedral crystal aggregates and phosphate precipitates were formed on the surface of the media due to the adsorption of phosphorus. Previous studies have also reported the formation of precipitates on the surface of

Study on the Adsorption Mechanism
SEM was performed to further explore the micro-morphology of the media and its changes in the operations of the bioretention systems (Figure 13a-d). Zeolite has a rough surface, folded, and a compact micropore structure. The micropore structure of 1-5 µm is distributed all over the media surface. It can provide an ideal environment for the growth of the biofilm and can improve the bio-load and ion exchange performance. Figure 12 are the surfaces of post-experiment zeolite, sponge iron, and limestone, respectively. Compared to the original zeolite (Figure 13a), the zeolite used after the experiment exhibited a rough surface and more crystals. Polyhedral crystal aggregates and phosphate precipitates were formed on the surface of the media due to the adsorption of phosphorus. Previous studies have also reported the formation of precipitates on the surface of similar materials [30,47]. EDS results showed that the main elements on the surface of limestone and sponge iron were calcium, iron, aluminum, and silicon ( Figure 12). The EDS peak intensity of Ca, Al, and Fe in the media decreased after the rainfall test ( Figure 12). Combined with the results of SEM, the formation of phosphorus peaks in the media after rainfall indicates that phosphate precipitation has occurred on the surface of the media, and the possible reactions are as follows in Equations (10)- (12): Appl. Sci. 2020, 10, x FOR PEER REVIEW 15 of 20 similar materials [30,47]. EDS results showed that the main elements on the surface of limestone and sponge iron were calcium, iron, aluminum, and silicon (Figure 12g,h). The EDS peak intensity of Ca, Al, and Fe in the media decreased after the rainfall test (Figure 12e,f). Combined with the results of SEM, the formation of phosphorus peaks in the media after rainfall indicates that phosphate precipitation has occurred on the surface of the media, and the possible reactions are as follows in Equations (10)-(12):  The XRD patterns, as shown in Figure 14, indicate that both zeolite and limestone have regular crystal structures, before and after the simulated rainfall. The characteristic diffraction peaks did not change significantly, which suggests that the crystal structure of the zeolite and the limestone did not transform during rainwater purification. Although the results of the SEM and EDS suggest that phosphate precipitates are present on the surface of zeolite and limestone, the characteristic peaks for phosphate were not found by XRD. This indicates that the chemical interaction between the cations and the zeolites may have occurred mainly on the surface of the media [12]. It is also possible that the characteristic diffraction peaks of phosphate coincide with those of zeolite and limestone [30,48].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 16 of 20 Figure 13. The figure shows the SEM images for (a) the zeolite before tests; (b) the zeolite after tests; (c) the post-experiment sponge iron; and, (d) the post-experiment limestone. The EDS spectrum is also illustrated for the following samples: (e) pre-experiment zeolite; (f) post-experiment zeolite; (g) post-experiment sponge iron; and, (h) post-experiment limestone.
The XRD patterns, as shown in Figure 14, indicate that both zeolite and limestone have regular crystal structures, before and after the simulated rainfall. The characteristic diffraction peaks did not change significantly, which suggests that the crystal structure of the zeolite and the limestone did not transform during rainwater purification. Although the results of the SEM and EDS suggest that phosphate precipitates are present on the surface of zeolite and limestone, the characteristic peaks for phosphate were not found by XRD. This indicates that the chemical interaction between the cations and the zeolites may have occurred mainly on the surface of the media [12]. It is also possible that the characteristic diffraction peaks of phosphate coincide with those of zeolite and limestone [30,48]. The media in the bioretention system serves as the carrier of microorganism and facilitator for plant growth and is also where physical, chemical, and biological reactions take place. Plants absorb pollutants directly from the water, while the presence of microorganisms in the root zone is beneficial for plant growth, as illustrated in Figure 15. The removal of pollutants by the bioretention system uses the media as the attachment point and the reaction interface. Through the interaction of microorganisms and plants, the pollutants are removed through microbial decomposition and utilization, media adsorption and precipitation, and plant harvesting. The media in the bioretention system serves as the carrier of microorganism and facilitator for plant growth and is also where physical, chemical, and biological reactions take place. Plants absorb pollutants directly from the water, while the presence of microorganisms in the root zone is beneficial for plant growth, as illustrated in Figure 15. The removal of pollutants by the bioretention system uses the media as the attachment point and the reaction interface. Through the interaction of microorganisms and plants, the pollutants are removed through microbial decomposition and utilization, media adsorption and precipitation, and plant harvesting. Figure 13. The figure shows the SEM images for (a) the zeolite before tests; (b) the zeolite after tests; (c) the post-experiment sponge iron; and, (d) the post-experiment limestone. The EDS spectrum is also illustrated for the following samples: (e) pre-experiment zeolite; (f) post-experiment zeolite; (g) post-experiment sponge iron; and, (h) post-experiment limestone.
The XRD patterns, as shown in Figure 14, indicate that both zeolite and limestone have regular crystal structures, before and after the simulated rainfall. The characteristic diffraction peaks did not change significantly, which suggests that the crystal structure of the zeolite and the limestone did not transform during rainwater purification. Although the results of the SEM and EDS suggest that phosphate precipitates are present on the surface of zeolite and limestone, the characteristic peaks for phosphate were not found by XRD. This indicates that the chemical interaction between the cations and the zeolites may have occurred mainly on the surface of the media [12]. It is also possible that the characteristic diffraction peaks of phosphate coincide with those of zeolite and limestone [30,48]. The media in the bioretention system serves as the carrier of microorganism and facilitator for plant growth and is also where physical, chemical, and biological reactions take place. Plants absorb pollutants directly from the water, while the presence of microorganisms in the root zone is beneficial for plant growth, as illustrated in Figure 15. The removal of pollutants by the bioretention system uses the media as the attachment point and the reaction interface. Through the interaction of microorganisms and plants, the pollutants are removed through microbial decomposition and utilization, media adsorption and precipitation, and plant harvesting. The removal efficiency of NH 4 + -N in the present study was studied by the single factor experiment method. Generally, there are multiple factors simultaneously influencing the removal performance of bioretention system with plants and saturated areas in practical conditions. For example, adsorption, nitrification, denitrification, and the removal efficiency are influenced by temperature, and the bioretention system can reach very low temperatures on cold winter days. Therefore, for practical considerations, further research focusing on the influence of local weather conditions, stormwater composition, and load levels, on the removal of nutrients and heavy metals from stormwater runoff are recommended.

Conclusions
Understanding the dynamics and effectivity of bioretention systems is an essential step in creating effective measures to control stormwater runoff pollution. In this study, saturated zone and plants were introduced into the conventional bioretention system to determine the influence of rainfall intensity, concentration, and antecedent dry periods on stormwater runoff treatment performance, in terms of effluent quality, adsorption, and complexation. The results show that nitrogen and phosphorus removal could significantly be improved by introducing saturated zones and plants for different rainfall intensities, concentrations, and wet-dry cycles. In the early stages of rainfall, the quality of effluent water from bioretention systems can reach water quality standards of Grade III and IV for surface water in China, but the effluent quality deteriorated over time. The results indicate that the combined use of saturated areas, media, plants, and microorganisms in the bioretention system is an effective measure for ecological purification. According to the SEM-EDS analysis, phosphate precipitates were found to have formed on the media surface after the adsorption. The XRD analysis shows that the original crystal structure of the media does not change during rainwater infiltration. Rainfall intensity can significantly influence the removal of pollutants in the bioretention systems. The results suggest longer rainfall recurrence intervals would result in reduced treatment capacity of bioretention systems in normal events. In terms of removal efficiency, the dry period was found to have a negative correlation for TP and a positive correlation for TN; however, it showed little effect on the removal efficiency for NH 4 + -N. Further research on road surface runoff management can focus on the microbial mechanisms towards long-term purification and the transformation mechanism of nitrogen.