In this study, scale-based runoff plots of concave grasslands were designed and simulated rainfall experiments were conducted to investigate their retention effectiveness for runoff volume and pollutant loads, and to analyze the influences of concave depths on runoff and pollution retention of grasslands. Results showed that mean time to runoff of concave grasslands was 88.5 minutes, which was 5.3 times than that of flat grassland. Average peak flow rate of concave grasslands was reduced by 36.2% compared with flat grassland. Concaved grasslands averagely retained 58.2% of stormwater runoff. Deeper concave depths significantly increased runoff detention and retention performance of grasslands. Total suspended solids (TSS) load reduction rates of concave grasslands were ranged from 50.8% to 97.3%. Total nitrogen (TN) load reduction rate was 49.8% for concave depth of 10 cm. Total phosphorus (TP) load reduction rates were 45.0% and 93.9% for grasslands with 5 cm and 10 cm concave depths, respectively. Pollution load reduction rates of TSS, TN and TP enhanced along with the increase in concave depths. The estimated minimum area ratios of upslope impervious surface to grasslands of 5 cm and 10 cm concave depths were approximately 1:1 under 20 mm rainfall events, and 38:1 under 5 mm rainfalls, respectively.

  • Concaved grasslands averagely retained 61.1% of stormwater runoff.

  • Mean time to runoff of concave grasslands was 5.3 times than that of flat grassland.

  • Deeper concave depths significantly increased runoff detention and retention performance.

  • TSS, TN and TP load reduction rates enhanced along with increased concave depths.

  • Estimated minimum area ratios of impervious surface to concave grasslands were 0.96–37.6.

Graphical Abstract

Graphical Abstract
Graphical Abstract

Over the past few decades, China's urbanization and impervious surfaces have grown rapidly and led to a worsening situation such as urban pluvial flooding and water environment pollution (Jia et al. 2013; Lu et al. 2016), which is joining force to restrict urban sustainable development. However, traditional improved measures on urban stormwater controls (such as increasing the drainage pipe diameters and extending drainage networks to facilitate the rapid runoff discharge) are prohibitively expensive or not effective at mitigating excessive surface runoff during extreme storm events (Pomeroy 2007; Berland et al. 2017). The challenges which urban planners and landscape architects confronted are retrofitting urban landscape for effectively regulating stormwater runoff (De Greef 2005; Ramyar et al. 2021). However, in practice, rainwater management/treatment facilities are usually designed by civil and environmental engineers and do not focus on the basic elements of urban landscape (Backhaus et al. 2012; Ranzato 2017). To improve this situation and promote sustainable stormwater management, the Chinese government announced a ‘Sponge City’ initiative in building urban rainwater infrastructures in 2013 (Che et al. 2015; Zhang et al. 2019). During 2015 and 2016, Chinese government selected 30 cities as pilot sites, and more than 130 cities have formulated ‘Sponge City’ development plans (Chan et al. 2018).

The concept basis and key technology of ‘Sponge City’ are based on green infrastructure (GI) and low impact development (LID) measures (Jia et al. 2017; Li et al. 2019a). The GI/LID strategies emphasized to utilize on site and small-scale source control techniques to maintain and enhance the pre-development hydrologic function of the urban landscape (Benedict & McMahon 2006; Dietz 2007; Baek et al. 2015; Fletcher et al. 2015; Li et al. 2020). The small-scale on-site GI/LID facilities including a green roof, rain barrel/cistern, bioretention system, permeable pavement, concave grassland and wetland channel, etc. (Jiang et al. 2022). Moreover, the ‘Sponge City’ program also aims to transform the urban landscape whilst promoting the greater use of LID facilities that enhance infiltration and storage of urban stormwater such as infiltration trenches, grassed swales, and storage ponds (Xia et al. 2017; Chan et al. 2018; Bai et al. 2019). Especially, the GI concept attempts to influence urban planning and layouts to maximize the benefits of green space (Fletcher et al. 2015).

Green space has an important role in retaining and detaining stormwater runoff (Yao et al. 2015; Zhang et al. 2015; Syrbe & Chang 2018). However, the capacity of green space in reducing direct runoff has not yet effectively employed because of limitations in urban landscape design (Li et al. 2017; Du et al. 2019). Traditionally, the ground surface of urban green space in China is flat or convex and higher than the surrounding roads (Liu et al. 2014). As a result, these flat or convex grasslands are prone to yield surface runoff under the excess infiltration at heavy rainfall events (Du et al. 2019). Concaved grassland (i.e., low elevation or sunken greenbelt), similar to bioretention, rain garden and grassed swale, which with berms deeper than normal lawns was listed as a key measure of the ‘Sponge City’ initiative; it refers to vegetated land that has a lower elevation than its surroundings and can temporarily retain stormwater (Ministry of Housing and Urban-Rural Development 2014). In contrast to a traditional convex design, concave grassland increases rainwater infiltration capacity of grasslands, as well as reduces peak flow and recharges groundwater aquifer (Hood et al. 2007). Thus, concave grassland is a cost-efficient facility of rainwater infiltration, which plays a significant role in rainwater infiltration, storage and utilization (Zhang et al. 2020).

The hydrologic responses and pollution removal performances of concave grasslands differ with various soil properties, concave depths and rainfall characteristics (Wang & Zhang 2015; Winston et al. 2016; Mai et al. 2018; Du et al. 2019). For example, Du et al. (2019) investigated that concave grassland with a depth of 10–20 cm can mitigate direct runoffs by 23.63–98.35%. Cheng et al. (2007) argued that concave grassland with a concave depth of 10 cm could effectively reduce runoff of pluvial floods with 3-year return period. Ye et al. (2001) found that, at a service impervious area ratio of 1:1, runoff reduction rates of a concave grassland for storms of 10, 50, and 100 years were 87.15%, 58.48% and 50.75%, respectively. Liu et al. (2014) showed that runoff volume and peak discharge reduction were respectively achieved of 23.20% and 29.11% by retrofitting grasslands into a concave morphology with a depth of 5 cm. Winston et al. (2016) showed that concave grassland could reduce peak flow by 24%–96% for a 1-year storm event. Mai et al. (2018) investigated that low-elevation greenbelt with 10 cm depth and 5 cm height gullies achieved 43.54%–94.00% runoff volume reduction and 40.59%–97.48% pollutant removal effectiveness. Hunt et al. (2008) showed bioretention effectively attenuated peak flow rates with a mean reduction percentage of 99%. Shrestha et al. (2018) investigated bioretention cells averagely retained suspended solids loads by 94%. Wang & Zhang (2015) found that concave grasslands significantly reduce phosphate, with average total phosphate (TP) reduction rates of 86.28%. Cheng et al. (2009) investigated that the average removal rate of TP from concave grassland outflow was 47.35%.

A series of fundamental and applied questions remain unsolved for the general utilization and development of LID techniques, which actually require further and systematic investigations through extensive practical programs (Vogel et al. 2015; Li et al. 2019b). The main challenge in implementation of concave grasslands is how to efficiently evaluate retention effectiveness, optimize structural configuration and spatial placement (Wang & Banzhaf 2018; Tansar et al. 2022). Although a plethora of studies have been undertaken to assess the runoff reduction effectiveness of LID facilities using pilot experiment and statistical/hydrological model simulation (Li et al. 2019a). However, only few studies were performed through field experiments to investigate the runoff and pollution reduction effectiveness of concave grasslands and explain their runoff and pollution retention mechanism (Elliott & Trowsdale 2007; Ahiablame et al. 2012; Liu et al. 2014). It is difficult to determine the optimal structural configuration and adequate guidelines of concave grasslands without lots of field experiments (Liu et al. 2018). Specifically, the local design criteria and operation conditions have become the critical and urgent questions for the effective design and resilient management of the LID practices (Duan et al. 2016). Therefore, it is necessary to collect practical field data that contribute to finding the most effective strategic solution to overcome barriers for widespread promotion and adoption of LID practices (Shafique & Kim 2015; Eckart et al. 2017).

In the present study, the well-defined scale-based runoff plots of concave grasslands and well-controlled field experimental approaches were conducted in a semi-arid region to investigate the retention performances of runoff quantity and pollution loads of concave grasslands. The impacts of concave depths on stormwater runoff retention and pollution removal rates were quantified, and the minimum area ratios of upslope impervious surface to concave grasslands were estimated. The aim of the present study was: (1) to investigate the effectiveness for runoff volume retention and pollution load reduction of concave grasslands; (2) to analyze the influences of concave depths on runoff and pollution retention of grasslands; (3) to estimate the area ratios of connecting upslope impervious surface to concave grassland. These results are expected to improve our understanding of the hydrological behaviour of concave grasslands, and help urban planners for designing appropriate concave depths and upslope impervious area for retrofitting existing grasslands to be concave.

Runoff plots of concave grasslands

Four scale-based runoff plots were constructed by a custom-made 2 mm thick stain-less steel platform (with internal dimensions of 100 cm wide × 100 cm long × 25–40 cm height) (Figure 1(a)). A marked steel V-flume was created at the outlet of each plot for facilitating the collection of surface runoff. In the inner section of each runoff plot, a soil layer with 15-cm thickness was padded with sandy loam soil (Liu et al. 2020). Physical properties of the soil used in the present study were shown in Table 1. The ground of concave grasslands was designed with lower elevation than the lower edges of the outlet flumes (Figure 1(b)). In addition, the outer edge of runoff plot was also designed 10 cm higher than the grassland ground to prevent rainfall splashing. In this study, three concave grasslands with concave depths of 5 cm (test plot numbered as C-5), 10 cm (C-10) and 15 cm (C-15) were constructed and fully occupied with 1 m2 of the runoff plots (Figure 1(c)). Runoff plots of concave grasslands were individual without connecting to impervious surface inflow. As a control sample, a runoff plot of flat grassland (C-0) was also designed, which following the same structure of the above concave grasslands without concave depth. In order to prevent the soil subsidence caused by the simulated rainfall, the soil was watered several times and filled. Grasslands were planted with tall fescue (Festuca elata Keng ex E. Alexeev) on the soil layers, and the average plant height was about 10 cm during the experimental period. Before the field experiments, runoff plots were placed outdoors about two weeks in summer until the soil moisture returned to the normal condition (Liu et al. 2019a). During the field experiments, the initial substrate moisture condition was controlled by the same rainfall experimental interval (5–6 days) of each rainfall repetitions.
Table 1

Soil physical properties in this study

MeasurementsUnitsValues
Dry bulk density g/cm3 1.35 ± 0.61 
Effective porosity 50.2 ± 6.90 
Saturated water content % (v/v) 51.7 ± 4.24 
Residual water content % (v/v) 26.0 ± 1.30 
Infiltration rate cm/d 23.0 ± 0.48 
MeasurementsUnitsValues
Dry bulk density g/cm3 1.35 ± 0.61 
Effective porosity 50.2 ± 6.90 
Saturated water content % (v/v) 51.7 ± 4.24 
Residual water content % (v/v) 26.0 ± 1.30 
Infiltration rate cm/d 23.0 ± 0.48 
Figure 1

(a) Design of scale-based runoff plots; (b) Structure of concave grasslands; (c) Runoff plots of concave grasslands in the field experiments; (d) Norton rainfall simulator.

Figure 1

(a) Design of scale-based runoff plots; (b) Structure of concave grasslands; (c) Runoff plots of concave grasslands in the field experiments; (d) Norton rainfall simulator.

Close modal

Field rainfall simulation experiments

In this study, a Norton artificial rainfall simulator was set up at 3.5 m above the experimental runoff plots (Figure 1(d)). Veejet 80100 nozzles with 41 kPa water pressure were applied in the spraying systems and spaced 1.1 m apart with a computer that oscillated across the plot to generate a constant rainfall intensity. Specifically, the median volume of raindrop size obtained by this simulator was 2.2 mm, and the uniformity coefficient of rainfall reached more than 0.8. Rainfall depth was monitored by a standard tipping bucket rain gauge (Onset HOBO 0.2 mm Rainfall Smart Sensor, S-RGB), positioned adjacent near and at the same height as that of the experimental runoff plots. In this paper, the runoff outflow from the designated runoff plot was collected in a plastic container below the platform by means of a pipe at its downstream end for continuous monitoring of the weight. In the rainfall experiments, six test repetitions were conducted for each runoff plot. The simulated rainfall events were chosen according to the significant runoff generation process and the capacity of runoff collection by plastic container based on a preliminary experiment (Liu et al. 2019a).

The rainfall intensity was designed based on the local storm formulation (Ji et al. 2002), which was expressed as following:
(1)
where, i represents the rainfall intensity (mm/min), N represents the recurrence period (year), t represents the rainfall duration (min).

Runoff water quality test

In this study, we used the collected rainwater from a cement ground at the experiment site, and replenished with tap water as the source of supply water for the rainfall simulator. Water quality parameters of concave grassland runoff in this study were focused on total suspended solids (TSS), total nitrogen (TN) and TP concentrations. The outflow from the grassland was collected in a plastic container below the runoff plot through a drainage pipe at its downstream end. Before runoff sampling, the plastic containers and sampling bottles were rinsed three times with distilled water. In each rainfall event, the composite water sample was manually collected after the runoff flow is over. In addition, in order to determine the water quality of input rainwater, six rainwater samples were also collected. The water samples were collected with 0.5 L pre-cleaned polyethylene bottles and immediately stored in a refrigerator at the experimental station. After field experiments, the collected water samples were carried to the laboratory and prepared for testing water quality. The tested methods of water quality were according to the Standard Methods for the Examination of Water and Wastewater, which was published by the American Public Health Association (APHA et al. 2005). TSS concentration was measured using the filtering, drying, and weighing method (Little et al. 2005). TN concentration was measured using the alkaline potassium persulfate digestion and UV spectrophotometric method. TP concentration was measured using the persulfate digestion spectrophotometric method (Liu et al. 2019b).

Data analysis methods

The relationship describing the runoff retention from the concave grassland was computed as follows (Liu et al. 2019a):
(2)
where, Rr was the runoff retention percentage (%), P referred to the rainfall volume actually received by concave grassland (mm), R was the runoff depth of concave grassland (mm).
The formula for calculating the pollutant load reduction rate is shown in following (Liu et al. 2020):
(3)
where, Rp is the runoff pollutant load reduction rate (%), Cin is the average pollutant concentration of rainwater input (mg/L), Vin is the total rainfall input volume (L), Cout is the average pollutant concentration of runoff (mg/L), and Vout is the total runoff volume (L).

On the basis of the monitored rainfall and runoff data collected in this study, selected hydrological indicators including time to runoff, runoff discharge depth, mean flow rate, peak flow rate, runoff retention percentage, runoff coefficient and pollutant load reduction rate were calculated and utilized for further analysis. The mean and standard deviation of each runoff indicator and pollutant load reduction rate were quantified for detecting the general features of the runoff retention and related pollution removal characteristics (Liu et al. 2020). The one-way analysis of variance (ANOVA) followed by the Tukey post-hoc test was used to compare the differences in means of time to runoff, runoff retention percentage and peak flow rate for different concave depths and analyzed in the SPSS 17.0. Histograms and line charts were all prepared in the SigmaPlot 14.0.

Calculating the area ratio of impervious area to concave grassland

To adequately infiltrate and retain stormwater runoff under limited land resources, concave grasslands should be positioned where they can capture inflow runoff from impervious surfaces (Yang & Chui 2018). In smaller rainfall events, concave grassland usually has a void storage volume that can accommodate and store inflow runoff from the connecting upslope impervious surface. According to maximum runoff retention capacity of flat and concave grasslands, the minimum area ratio of upslope impervious surface to concave grassland can be calculated as following:
(4)
(5)
where, is the inflow from impervious surface (mm), P is the rainfall depth (mm), is the initial loss of impervious surface (mm); it assumed as 3 mm in this study, is the area of imperious surface (m2), is the maximum retention capacity of concave grassland (mm), is the maximum retention capacity of flat grassland (mm), is the area of concave grassland (m2), r is the minimum area ratio of upslope imperious surface to concave grassland.

General hydrologic performance of concave grasslands

During the field rainfall experiments, the simulated rainfall intensity was recorded with a mean value of 0.62 mm/minute, and with a standard deviation of 0.04. The rainfall amount had a mean value of 84.1 mm, and the average rainfall duration was 136.7 minutes. These simulated rainfall events were all regarded as heavy storms. A hydrograph of the monitored typical rainfall-runoff process of the flat and concave grasslands is shown in Figure 2, which intuitively illustrates the performance on runoff retention, runoff onset delay and peak flow reduction of the concave grasslands compared to the flat grassland. Table 2 showed the statistical characteristics of overall hydrological performances for the flat and concave grasslands. Mean value of time to runoff of concave grasslands was 88.5 minutes, which was 5.3 times that of flat grassland (16.8 minutes). Mean flow rate of flat grassland was 0.6 L/minute, and mean flow rate of the concave grasslands achieved was 0.4 L/minute. The average peak flow rate of concave grasslands was reduced by 38.1% compared with flat grassland. Under these heavy storms, runoff retention of flat grassland was only 5.0%. The concave grasslands dramatically retained runoff, mean runoff retention was achieved of 58.2%. Average runoff coefficient of the concave grasslands was 0.89, which was half of that of the flat grassland (0.42). In this study, average runoff retention of concave grasslands was occupied a lower range of the reported values. The retention percentages of concave grasslands in this study were much lower than Du et al. (2019), Mai et al. (2018), and Ye et al. (2001) who have reported runoff retention ranged from about 23%-94%. Higher efficiency was also found in DeBusk & Wynn (2011), who observed bioretention could reduce 97% to 99% of surface runoff. And de Macedo et al. (2019) found bioretention presents an average runoff retention efficiency of 70%. A plausible reason for this is the fact that we have used large storms in the rainfall simulation experiments. The runoff detention effects were aligned with Liu et al. (2018), it concluded that the concave grasslands can significantly delay time to runoff 7.3 minutes and 5.2 minutes compared with the impervious surface. As well as, Winston et al. (2016) showed three bioretention cells reduce runoff by 59, 42, and 36% over the monitoring period. Therefore, if the flat or convex grasslands are changed to be concave to increase infiltration and retain partial impervious runoff, their capacity of stormwater runoff detention and retention as well as peak flow reduction would be greatly enhanced. These various findings about hydrological performance of concave grasslands and bioretention cells are mainly due to the large variations in soil media composition and thickness, concave depth, and climate conditions (Davis et al. 2009; Skorobogatov et al. 2020).
Table 2

Characteristics of runoff responses of the flat and concave grasslands

IndicatorsFlat grasslandConcaved grasslands
Time to runoff (minutes) 16.83 ± 1.05 88.50 ± 48.23 
Runoff discharge depth (mm) 77.36 ± 1.93 32.75 ± 33.18 
Mean flow rate (L/minute) 0.61 ± 0.19 0.44 ± 0.28 
Peak flow rate (L/minute) 0.97 ± 0.05 0.60 ± 0.42 
Runoff retention percentage (%) 5.34 ± 3.85 58.20 ± 35.12 
Runoff coefficient 0.89 ± 0.16 0.42 ± 0.35 
IndicatorsFlat grasslandConcaved grasslands
Time to runoff (minutes) 16.83 ± 1.05 88.50 ± 48.23 
Runoff discharge depth (mm) 77.36 ± 1.93 32.75 ± 33.18 
Mean flow rate (L/minute) 0.61 ± 0.19 0.44 ± 0.28 
Peak flow rate (L/minute) 0.97 ± 0.05 0.60 ± 0.42 
Runoff retention percentage (%) 5.34 ± 3.85 58.20 ± 35.12 
Runoff coefficient 0.89 ± 0.16 0.42 ± 0.35 
Figure 2

Typical rainfall and runoff processes of flat and concave grasslands.

Figure 2

Typical rainfall and runoff processes of flat and concave grasslands.

Close modal

Effects of concave depths on runoff retention

Figure 3 showed the effects of concave depths on time to runoff, runoff retention and peak flow rate of concave grasslands. In the simulated rainfall events, no runoff was discharged from the grassland with concave depth of 15 cm (C-15). As the concave depth increased from 5 to 10 cm, mean time to runoff of concave grassland C-10 was 2.9 times that of the concave grassland C-5; however, their difference was not statistically significant. A significant difference of time to runoff was investigated between grasslands C-0 and C-10 for the analyzed rainfall events registered p < 0.05. The runoff retention was significantly increased as concave depths increased (p < 0.05). A little runoff volume was discharged from the concave grassland C-10, thus the runoff retention achieved up to 97.1%, which was 3.9 times than that of the concave grassland C-5. Peak flow rate of the concave grassland C-10 was one third of the grassland C-0 (p< 0.05). Moreover, no statistically significant difference in the mean peak flow rate was found between the grasslands C-0 and C-5 in this study period. Therefore, the deeper concave depths led to a delay in runoff generation, and increased the capacity of runoff retention for concave grasslands. This is owing to runoff volume retention of the concave grasslands was mainly depended on their rainwater storage capacity, which was directly determined by their concave depths. Similarly, to maximize runoff retention effectiveness, proportionally large bioretention surface area, deeper storage depths, deeper soil thickness and high soil infiltration rate are desired (Davis et al. 2009; Li et al. 2019c).
Figure 3

Influences of the concave depth on: (a) time to runoff, (b) runoff retention percentage, (c) peak flow rate. Mean values followed by the different letters are significantly different (p < 0.05) as determined by t-test followed by Tukey post-hoc test.

Figure 3

Influences of the concave depth on: (a) time to runoff, (b) runoff retention percentage, (c) peak flow rate. Mean values followed by the different letters are significantly different (p < 0.05) as determined by t-test followed by Tukey post-hoc test.

Close modal
Results of the regression analysis showed that the concave depths and time to runoff had a significant positive correlation (p < 0.001). There was also a positive correlation between the concave depths and runoff retention (see Figure 4), and the non-linear regression equation expressed as: y = 0.98x2 − 0.83x + 5.02 (p < 0.001). It indicated that a polynomial equation is applicable to predict runoff retention capacity of concave grasslands with different concave depths under similar structural details and rainfall conditions. The determination coefficient (R2) in the regression analysis of concave depths with time to runoff and retention percentages of concave grasslands were 0.60 and 0.99, respectively, which indicated that the concave depths accounting for 60% of the variance in time to runoff and about 99% of the variance in runoff retention percentages of concave grasslands.
Figure 4

Nonlinear regression of the concave depths with: (a) time to runoff, and (b) runoff retention percentage.

Figure 4

Nonlinear regression of the concave depths with: (a) time to runoff, and (b) runoff retention percentage.

Close modal

Runoff water quality of concave grasslands

Mean values and standard deviations of the tested water quality parameters of rainwater input and concave grasslands outflow were shown in Figure 5. The averaged water quality concentration in rainwater input was 41.9 mg/L of TSS, 0.6 mg/L of TN and 0.03 mg/L of TP. Mean TSS levels of three concave grasslands ranged from 11.2 mg/L to 21.7 mg/L. Mean TSS levels in outflow from the concave grasslands were lower than their levels of input rainwater, indicated that the concave grasslands can effectively filter the small soil particles through sedimentation and filtration. However, except TP in grassland C-5, mean TN and TP concentrations in concave grasslands outflow were higher than that of input rainwater. Average TN concentration of concave grasslands was elevated with concave depths increased, and varied from 0.9 mg/L to 5.7 mg/L. All concave grasslands had low TP concentration, the mean value ranged from 0.02 mg/L to 0.04 mg/L. Moreover, TN and TP concentrations of concave grassland C-10 were higher than that of grassland C-5, indicating that the pollutants were readily leached out from the deeper concave depths. In this study, except the higher TSS concentration, the collected rainwater had a relatively lower nutrient concentration. Thus, nitrogen and phosphorus will be leaching from the soil layer of concave grasslands into discharged runoff. An explanation is that deeper concave depth may increase the duration of runoff retention, which increased the nitrogen leaching from soil layer, resulting in an increase in nitrogen concentration in the outflow runoff from concave grasslands.
Figure 5

Averaged concentrations of water quality parameters of rainwater and concave grasslands.

Figure 5

Averaged concentrations of water quality parameters of rainwater and concave grasslands.

Close modal

Impacts of concave depths on pollutant load removal

The reduction effects on TN and TP loads by concave grasslands were attributed to their abilities of runoff volume retention. As demonstrated in Table 3, TSS load reduction rates of concave grasslands were ranged from 50.8% to 97.3%. TN load reduction rate was 49.8% for concave depth of 10 cm. TP load reduction rates were 45.0% and 93.9% for C-5 and C-10 grasslands, respectively. However, TN load reduction rates were ineffective for C-0 and C-5 grasslands, and TP load reduction was ineffective for C-0 grassland. The results were inconsistent with the research presented by Hou et al. (2014), who found the average TP and TN concentrations were reduced by 75.8% and 66.6% by a sunken lawn infiltration system, respectively. As well, Davis et al. (2006) reported good reductions of bioretention in phosphorus (65%–87%) and moderate reductions in nitrogen (49%–59%) concentrations. Fahui et al. (2013) examined the bioretention media performed well with exceptional removal of over 95% of TSS and 82%–96% of phosphorus. Zhang et al. (2020) simulated the water quality effects of rain gardens exhibited the maximum contaminant reduction rates with TSS of 15.5%, TN of 17.3%, and TP of 19.1%. Overall, runoff pollutants load reduction rates of TSS, TN and TP were enhanced along with the increase in concave depths for grasslands. This trend was consistent with the high runoff retention of deeper concave grassland. Consequently, pollutant retention by the concave grasslands was strongly driven by the retention of runoff water volume. The water quality improvement of concave grasslands related with soil media and vegetation species by complex physicochemical reaction, thus their performance will be varied among pollutants types and different structural combination. Therefore, specifically optimizing the growth media, addition of amendments, screening of suitable vegetation, and design alterations should be further researched (Vijayaraghavan et al. 2021).

Table 3

Pollutants load reduction rates of the concave grasslands

Concaved grasslandsAverage pollutants reduction rates (%)
TSSTNTP
C-0 50.85 − 34.68 − 36.19 
C-5 80.02 − 31.89 45.05 
C-10 97.27 49.76 93.91 
Concaved grasslandsAverage pollutants reduction rates (%)
TSSTNTP
C-0 50.85 − 34.68 − 36.19 
C-5 80.02 − 31.89 45.05 
C-10 97.27 49.76 93.91 

The area ratios of upslope impervious surface to concave grasslands

The present study designed the individual runoff plots of concave grasslands without connecting to surrounding impervious surface inflow. In general, concave grasslands allow stormwater to flow into grasslands from the surrounding roads, pavements and squares, etc. Ideally, the minimum area ratio of upslope impervious surface to concave grassland can be estimated according to maximum runoff retention capacity of flat and concave grasslands. Based on the experimental rainfall-runoff data of grasslands, the calculated maximum runoff retention capacity of C-0, C-5, and C-10 grasslands were 4.2 mm, 20.5 mm, and 79.3 mm, respectively. Table 4 showed the estimated minimum area ratios of upslope impervious surface to grasslands varied from 0.96 to 37.6 under 5, 10, 20 mm rainfall events. The minimum area ratio of upslope impervious surface to grassland C-5 was approximately 1:1 under 20 mm rainfall events. The grassland C-10 could fully control the inflow from a 37.6 times area of upslope impervious surface under rainfall depth of 5 mm. Hou et al. (2020) simulated that when the ponding depth was 15 cm and the flow area ratio was 15:1, runoff control rate of rain garden was 31.9%–100% for 0.5 year to 50 years return periods. Richards et al. (2015) reported that concave grassland with a ratio of 7.5% to catchment area could reduce the flood volume by more than 90%. In this study, the concave grasslands can fully accommodate the assumed inflow runoff from minimum upslope impervious areas without outflow discharge. In application, the estimated minimum impervious area ratios should be enlarged with partial retention goals. These results can help urban planners to design the appropriate criteria of area ratio for upslope impervious surface to promote more efficient concave grassland retrofitting.

Table 4

Estimated area ratios of upslope impervious surface to concave grasslands

Concaved grasslandsMaximum runoff retention capacity (mm)Area ratios of upslope impervious surface under different rainfall depths
5 mm10 mm20 mm
C − 5 20.5 8.2 2.3 0.96 
C − 10 79.3 37.6 10.7 4.4 
Concaved grasslandsMaximum runoff retention capacity (mm)Area ratios of upslope impervious surface under different rainfall depths
5 mm10 mm20 mm
C − 5 20.5 8.2 2.3 0.96 
C − 10 79.3 37.6 10.7 4.4 

The present study experimentally verified runoff volume and pollution load reduction effectiveness of concave grasslands through scale-based runoff plots and simulated storm events, and quantified the impacts of concave depths on stormwater runoff retention and pollution removal rates. Compared with flat grassland, the concave grasslands dramatically delayed runoff generation, retained runoff volume and reduced peak flow rate. Due to runoff volume reduction of concave grasslands depended on their water storage and infiltration capacities, the deeper concave depths led to a delay in runoff generation, and increase the capacity of runoff retention. Results of the regression analysis showed that the concave depths significantly increased time to runoff, and there was also a positive correlation between the concave depths and runoff retention. The developed polynomial equation is applicable to predict runoff retention performance for a concave grassland under similar structural details and rainfall conditions, which help to understand the contribution of concave depth on runoff retention performance. Except TP concentration in grassland with concave depth of 5 cm, mean TSS levels in outflow from concave grasslands were lower as well as TN and TP concentrations were higher than that of input rainwater. Pollution load reduction rates of TSS, TN and TP enhanced along with the increase in concave depths. It indicated that pollutant retention by the concave grasslands was strongly driven by the retention of runoff water quantity. The estimated minimum area ratios of upslope impervious surface to concave grasslands were varied from 0.96 to 37.6 under different rainfall events. In application, the estimated area ratios should be enlarged with partial runoff retention goals.

These experimental results are potentially helpful in better understanding of the mechanism of runoff and pollution retention of concave grasslands, and it would be beneficial for urban planners to raise awareness of retrofitting grasslands to be concave. In stormwater management practices, the concave design is highly proposed in newly built urban grasslands and retrofitting of existing grasslands. Moreover, spatial distribution planning of grasslands should consider draining portion of adjacent impervious surface runoff into the grasslands. These results also assist urban managers and planners for designing an appropriate criteria of concave depth and upslope impervious area ratio to promote more efficient concave grassland retrofitting. As limitations of the results application, results in the present study may largely depend on design structures of concave grasslands, rainfall characteristics, and climatic regions. Future study is needed to experimentally investigate the runoff retention capacity of concave grasslands in different climate regions, continuously simulate the long-term retention performance of concave grasslands, and assess the reduction effectiveness of large-scale implementation of concave grasslands in urban areas for flooding risk control and water quality management.

This study was supported by the National Natural Science Foundation of China [42071051], and the Major Science and Technology Projects of Gansu Province [21ZD4FA008].

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Ahiablame
L.
,
Engel
B. A.
&
Chaubey
I.
2012
Effectiveness of low impact development practices: literature review and suggestions for future research
.
Water Air and Soil Pollution
223
(
7
),
4253
4273
.
https://doi.org/10.1007/s11270-012-1189-2
.
APHA, AWWA. & WEF
2005
Standard Methods for the Examination of Water and Wastewater
, 21st edn.
American Public Health Association
,
Washington, DC
,
USA
.
Backhaus
A.
,
Dam
T.
&
Jensen
M. B.
2012
Stormwater management challenges as revealed through a design experiment with professional landscape architects
.
Urban Water Journal
9
(
1
),
29
43
.
https://doi.org/10.1080/1573062X.2011.633613
.
Baek
S.
,
Choi
D.
,
Jung
J.
,
Lee
H.
,
Lee
H.
,
Yoon
K.
&
Cho
K. H.
2015
Optimizing low impact development (LID) for stormwater runoff treatment in urban area, Korea: experimental and modeling approach
.
Water Research
86
,
122
131
.
https://doi.org/10.1016/j.watres.2015.08.038
.
Benedict
M. A.
&
McMahon
E.
2006
Green Infrastructure: Linking Landscapes and Communities
.
Island Press
,
Washington, DC
,
USA
.
Berland
A.
,
Shiflett
S. A.
,
Shuster
W. D.
,
Garmestani
A. S.
,
Goddard
H. C.
,
Herrmann
D. L.
&
Hopton
M. E.
2017
The role of trees in urban stormwater management
.
Landscape and Urban Planning
162
,
167
177
.
https://doi.org/10.1016/j.landurbplan.2017.02.017
.
Chan
F. K.
,
Griffiths
J.
,
Higgitt
D.
,
Xu
S.
,
Zhu
F.
,
Tang
Y.
,
Xu
Y. Y.
&
Thorne
C. R.
2018
‘Sponge City’ in China-A breakthrough of planning and flood risk management in the urban context
.
Land Use Policy
76
,
772
778
.
https://doi.org/10.1016/j.landusepol.2018.03.005
.
Che
W.
,
Zhao
Y.
,
Li
J. Q.
,
Wang
W. L.
,
Wang
J. L.
,
Wang
S. S.
&
Gong
Y. W.
2015
Explanation of sponge city development technical guide: basic concept and comprehensive goals
.
Journal of China Water & Wastewater
31
(
8
),
1
5
.
(in Chinese)
.
Cheng
J.
,
Xu
Q.
,
Yang
K.
,
Liu
L.
&
Li
B.
2007
Rainfall penetration and storage effect and influencing factors of sunken green space
.
Water & Wastewater Engineering
33
(
5
),
45
49
.
(in Chinese)
.
Cheng
J.
,
Yang
K.
,
Huang
M.
,
Xie
B.
&
Li
X.
2009
Reduction effect of sunken green space on urban rainfall-runoff pollution
.
China Environmental Science
29
(
6
),
611
616
.
(in Chinese)
.
Davis
A. P.
,
Shokouhian
M.
,
Sharma
H.
&
Minami
C.
2006
Water quality improvement through bioretention media: nitrogen and phosphorus removal
.
Water Environment Research
78
(
3
),
284
293
.
https://doi.org/10.2175/106143005X94376
.
Davis
A. P.
,
Hunt
W. F.
,
Traver
R. G.
&
Clar
M.
2009
Bioretention technology: overview of current practice and future needs
.
Journal of Environmental Engineering
135
(
3
),
109
117
.
https://doi.org/10.1061/(ASCE)0733-9372(2009)135:3(109)
.
DeBusk
K. M.
&
Wynn
T. M.
2011
Storm-water bioretention for runoff quality and quantity mitigation
.
Journal of Environmental Engineering
137
(
9
),
800
808
.
https://doi.org/10.1061/(ASCE)EE.1943-7870.0000388
.
De Greef
P.
2005
Rotterdam Waterstad 2035
.
Episode Publishers
,
Rotterdam
,
Netherlands
.
de Macedo
M. B.
,
do Lago
C. A. F.
&
Mendiondo
E. M.
2019
Stormwater volume reduction and water quality improvement by bioretention: potentials and challenges for water security in a subtropical catchment
.
Science of the Total Environment
647
,
923
931
.
https://doi.org/10.1016/j.scitotenv.2018.08.002
.
Dietz
M. E.
2007
Low impact development practices: a review of current research and recommendations for future directions
.
Water, Air, & Soil Pollution
186
(
1
),
351
363
.
https://doi.org/10.1007/S11270-007-9484-Z
.
Du
S.
,
Wang
C.
,
Shen
J.
,
Wen
J.
,
Gao
J.
,
Wu
J.
,
Lin
W. P.
&
Xu
H.
2019
Mapping the capacity of concave green land in mitigating urban pluvial floods and its beneficiaries
.
Sustainable Cities and Society
44
,
774
782
.
https://doi.org/10.1016/j.scs.2018.11.003
.
Duan
H. F.
,
Li
F.
&
Yan
H.
2016
Multi-objective optimal design of detention tanks in the urban stormwater drainage system: LID implementation and analysis
.
Water Resources Management
30
(
13
),
4635
4648
.
https://doi.org/10.1007/s11269-016-1444-1
.
Eckart
K.
,
McPhee
Z.
&
Bolisetti
T.
2017
Performance and implementation of low impact development–A review
.
Science of the Total Environment
607
,
413
432
.
https://doi.org/10.1016/j.scitotenv.2017.06.254
.
Elliott
A. H.
&
Trowsdale
S.
2007
A review of models for low impact urban stormwater drainage
.
Environmental Modelling & Software
22
(
3
),
394
405
.
https://doi.org/10.1016/j.envsoft.2005.12.005
.
Fahui
N.
,
Zhanmeng
L.
,
Luyan
M.
&
Jie
L.
2013
Role of media: implication for application of urban storm water runoff bioretention systems
. In
2013 Third International Conference on Intelligent System Design and Engineering Applications
. pp.
1229
1232
.
Fletcher, T. D., Shuster, W. D., Hunt, W. F., Ashley, R., Butler, D., Arthur, S., Trowsdale, S., Barraud, S., Semadeni-Davies, A., Bertrand-Krajewski, J., Mikkelsen, P. S., Rivard, G., Uhl, M., Dagenais, D. & Viklander, M.
2015
SUDS, LID, BMPs, WSUD and more – The evolution and application of terminology surrounding urban drainage
.
Urban Water Journal
12
(
7
),
525
542
.
https://doi.org/10.1080/1573062X.2014.916314
.
Hood
M.
,
Clausen
J. C.
&
Warner
G. S.
2007
Comparison of stormwater lag times for low impact and traditional residential development
.
Journal of the American Water Resources Association
43
(
4
),
1036
1046
.
https://doi.org/10.1111/j.1752-1688.2007.00085.x
.
Hou
L.
,
Yang
H.
&
Li
M.
2014
Removal of chemical oxygen demand and dissolved nutrients by a sunken lawn infiltration system during intermittent storm events
.
Water Science and Technology
69
(
2
),
398
406
.
https://doi.org/10.2166/wst.2013.728
.
Hou
J.
,
Liu
F.
,
Tong
Y.
,
Guo
K.
,
Ma
L.
&
Li
D.
2020
Numerical simulation for runoff regulation in rain garden using 2D hydrodynamic model
.
Ecological Engineering
153
,
105794
.
https://doi.org/10.1016/j.ecoleng.2020.105794
.
Hunt
W. F.
,
Smith
J. T.
,
Jadlocki
S. J.
,
Hathaway
J. M.
&
Eubanks
P. R.
2008
Pollutant removal and peak flow mitigation by a bioretention cell in urban Charlotte, NC
.
Journal of Environmental Engineering
134
(
5
),
403
408
.
https://doi.org/10.1061/(ASCE)0733-9372(2008)134:5(403)
.
Ji
R. C.
,
Guo
X. D.
&
Liu
Y. L.
2002
A study of the rain storm intensity relation in Lanzhou City
.
Journal of Lanzhou Railway University (Natural Sciences)
21
(
6
),
65
68
.
(in Chinese)
.
Jia
H.
,
Yao
H.
&
Yu
S. L.
2013
Advances in LID BMPs research and practice for urban runoff control in China
.
Frontiers of Environmental Science & Engineering
7
(
5
),
709
720
.
https://doi.org/10.1007/s11783-013-0557-5
.
Jia
H.
,
Wang
Z.
,
Zhen
X.
,
Clar
M.
&
Yu
S. L.
2017
China's sponge city construction: a discussion on technical approaches
.
Frontiers of Environmental Science & Engineering
11
(
4
),
18
45
.
https://doi.org/10.1007/s11783-017-0984-9
.
Jiang
C.
,
Li
J.
,
Hu
Y.
,
Yao
Y.
&
Li
H.
2022
Construction of water-soil-plant system for rainfall vertical connection in the concept of sponge city: a review
.
Journal of Hydrology
605
,
127327
.
https://doi.org/10.1016/j.jhydrol.2021.127327
.
Li
F.
,
Zhang
F.
,
Li
X.
,
Wang
P.
,
Liang
J.
,
Mei
Y.
,
Chen
W. W.
&
Qian
Y.
2017
Spatiotemporal patterns of the use of urban green spaces and external factors contributing to their use in central Beijing
.
International Journal of Environmental Research and Public Health
14
(
3
),
237
.
https://doi.org/10.3390/ijerph14030237
.
Li
C.
,
Peng
C.
,
Chiang
P.
,
Cai
Y.
,
Wang
X.
&
Yang
Z.
2019a
Mechanisms and applications of green infrastructure practices for stormwater control: a review
.
Journal of Hydrology
568
,
626
637
.
https://doi.org/10.1016/j.jhydrol.2018.10.074
.
Li
F.
,
Yan
X. F.
&
Duan
H. F.
2019b
Sustainable design of urban stormwater drainage systems by implementing detention tank and LID measures for flooding risk control and water quality management
.
Water Resources Management
33
(
9
),
3271
3288
.
https://doi.org/10.1007/s11269-019-02300-0
.
Li
Q.
,
Wang
F.
,
Yu
Y.
,
Huang
Z.
,
Li
M.
&
Guan
Y.
2019c
Comprehensive performance evaluation of LID practices for the sponge city construction: a case study in Guangxi, China
.
Journal of Environmental Management
231
,
10
20
.
https://doi.org/10.1016/j.jenvman.2018.10.024
.
Li
Y.
,
Li
H. X.
,
Huang
J.
&
Liu
C.
2020
An approximation method for evaluating flash flooding mitigation of sponge city strategies – a case study of Central Geelong
.
Journal of Cleaner Production
257
,
120525
.
https://doi.org/10.1016/j.jclepro.2020.120525
.
Little
J. L.
,
Bennett
D. R.
&
Miller
J. J.
2005
Nutrient and sediment losses under simulated rainfall following manure incorporation by different methods
.
Journal of Environmental Quality
34
(
5
),
1883
1895
.
https://doi.org/10.2134/jeq2005.0056
.
Liu
W.
,
Chen
W.
&
Peng
C.
2014
Assessing the effectiveness of green infrastructures on urban flooding reduction: a community scale study
.
Ecological Modelling
291
(
1
),
6
14
.
https://doi.org/10.1016/j.ecolmodel.2014.07.012
.
Liu
W.
,
Chen
W.
&
Feng
Q.
2018
Field simulation of urban surfaces runoff and estimation of runoff with experimental curve numbers
.
Urban Water Journal
15
(
5
),
418
426
.
https://doi.org/10.1080/1573062X.2018.1508597
.
Liu
W.
,
Feng
Q.
,
Chen
W.
&
Deo
R. C.
2020
Stormwater runoff and pollution retention performances of permeable pavements and the effects of structural factors
.
Environmental Science Pollution Research
27
(
24
),
30831
30843
.
https://doi.org/10.1007/s11356-020-09220-2
.
Liu
W.
,
Feng
Q.
,
Chen
W.
,
Wei
W.
&
Deo
R. C.
2019a
The influence of structural factors on stormwater runoff retention of extensive green roofs: new evidence from scale-based models and real experiments
.
Journal of Hydrology
569
,
230
238
.
https://doi.org/10.1016/j.jhydrol.2018.11.066
.
Liu
W.
,
Wei
W.
,
Chen
W.
,
Deo
R. C.
,
Si
J.
,
Xi
H.
,
Li
B.
&
Feng
Q.
2019b
The impacts of substrate and vegetation on stormwater runoff quality from extensive green roofs
.
Journal of Hydrology
576
,
575
582
.
https://doi.org/10.1016/j.jhydrol.2019.06.061
.
Lu
S.
,
Zhang
X.
,
Bao
H.
&
Skitmore
M.
2016
Review of social water cycle research in a changing environment
.
Renewable & Sustainable Energy Reviews
63
,
132
140
.
https://doi.org/10.1016/j.rser.2016.04.071
.
Mai
Y.
,
Zhang
M.
,
Chen
W.
,
Chen
X.
,
Huang
G.
&
Li
D.
2018
Experimental study on the effects of LID measures on the control of rainfall runoff
.
Urban Water Journal
15
(
9
),
827
836
.
https://doi.org/10.1080/1573062X.2018.1561912
.
Ministry of Housing and Urban-Rural Development
2014
Technical Guidance for Sponge City: Constructing Low Impact Development Rainwater System
.
MHURD
,
Beijing
,
China
.
(in Chinese)
.
Pomeroy
C. A.
2007
Evaluating the Impacts of Urbanization and Stormwater Management Practices on Stream Response
.
PhD diss.
,
Colorado State University
,
USA
.
Ramyar
R.
,
Ackerman
A.
&
Johnston
D. M.
2021
Adapting cities for climate change through urban green infrastructure planning
.
Cities
117
,
103316
.
https://doi.org/10.1016/j.cities.2021.103316
.
Ranzato
M.
2017
Landscape elements as a basis for integrated water management
.
Urban Water Journal
14
(
7
),
694
703
.
https://doi.org/10.1080/1573062X.2016.1240807
.
Richards
P. J.
,
Farrell
C.
,
Tom
M.
,
Williams
N. S.
&
Fletcher
T. D.
2015
Vegetable raingardens can produce food and reduce stormwater runoff
.
Urban Forestry & Urban Greening
14
(
3
),
646
654
.
https://doi.org/10.1016/j.ufug.2015.06.007
.
Shafique
M.
&
Kim
R.
2015
Low impact development practices: a review of current research and recommendations for future directions
.
Ecological Chemistry and Engineering
22
(
4
),
543
.
https://doi.org/10.1515/eces-2015-0032
.
Shrestha
P.
,
Hurley
S. E.
&
Wemple
B. C.
2018
Effects of different soil media, vegetation, and hydrologic treatments on nutrient and sediment removal in roadside bioretention systems
.
Ecological Engineering
112
,
116
131
.
https://doi.org/10.1016/j.ecoleng.2017.12.004
.
Skorobogatov
A.
,
He
J.
,
Chu
A.
,
Valeo
C.
&
van Duin
B.
2020
The impact of media, plants and their interactions on bioretention performance: a review
.
Science of The Total Environment
715
,
136918
.
https://doi.org/10.1016/j.scitotenv.2020.136918
.
Syrbe
R. U.
,
Chang
J.
,
2018
Options and challenges for implementing green spaces in urban development
. In:
Towards Green Cities
(
Grunewald
K.
,
Li
J.
,
Xie
G.
&
Kümper-Schlake
L.
, eds).
Springer
,
Heidelberg
,
Germany
.
Tansar
H.
,
Duan
H. F.
&
Mark
O.
2022
Catchment-scale and local-scale based evaluation of LID effectiveness on urban drainage system performance
.
Water Resources Management
36
(
2
),
507
526
.
https://doi.org/10.1007/s11269-021-03036-6
.
Vijayaraghavan, K., Biswal, B. K., Adam, M. G., Soh, S. H., Tsen-Tieng, D. L., Davis, A. P., Chew, S. H., Tan, P. Y., Babovic, V. & Balasubramanian, R.
2021
Bioretention systems for stormwater management: recent advances and future prospects
.
Journal of Environmental Management
292
,
112766
.
https://doi.org/10.1016/j.jenvman.2021.112766
.
Vogel, J. R., Moore, T. L., Coffman, R. R., Rodie, S. N., Hutchinson, S. L., McDonough, K. R., McLemore, A. J. & McMaine, J. T.
2015
Critical review of technical questions facing low impact development and green infrastructure: a perspective from the Great Plains
.
Water Environment Research
87
(
9
),
849
862
.
https://doi.org/10.2175/106143015X14362865226392
.
Wang
J.
&
Banzhaf
E.
2018
Towards a better understanding of green infrastructure: a critical review
.
Ecological Indicators
85
,
758
772
.
https://doi.org/10.1016/j.ecolind.2017.09.018
.
Wang
S.
&
Zhang
L.
2015
The cutting effect of phosphorus in simulated road runoff by concave herbaceous field
.
Environmental Science & Technology
38
(
9
),
119
122
.
(in Chinese)
.
Winston
R. J.
,
Dorsey
J. D.
&
Hunt
W. F.
2016
Quantifying volume reduction and peak flow mitigation for three bioretention cells in clay soils in northeast Ohio
.
Science of the Total Environment
553
,
83
95
.
https://doi.org/10.1016/j.scitotenv.2016.02.081
.
Xia
J.
,
Zhang
Y.
,
Xiong
L.
,
He
S.
,
Wang
L.
&
Yu
Z.
2017
Opportunities and challenges of the Sponge City construction related to urban water issues in China
.
Science China Earth Sciences
60
(
4
),
652
658
.
https://doi.org/10.1007/s11430-016-0111-8
.
Yang
Y.
&
Chui
T. F. M.
2018
Optimizing surface and contributing areas of bioretention cells for stormwater runoff quality and quantity management
.
Journal of Environmental Management
206
,
1090
1103
.
https://doi.org/10.1016/j.jenvman.2017.11.064
.
Yao
L.
,
Chen
L.
,
Wei
W.
&
Sun
R.
2015
Potential reduction in urban runoff by green spaces in Beijing: a scenario analysis
.
Urban Forestry & Urban Greening
14
(
2
),
300
308
.
https://doi.org/10.1016/j.ufug.2015.02.014
.
Ye
S.
,
Liu
H.
&
Meng
G.
2001
Analysis of rain storage and infiltration in Sunken-Lawn under the condition of designed storm
.
Journal of China Agricultural University
6
(
6
),
53
58
.
(in Chinese)
.
Zhang
B.
,
Xie
G.
,
Li
N.
&
Wang
S.
2015
Effect of urban green space changes on the role of rainwater runoff reduction in Beijing, China
.
Landscape and Urban Planning
140
,
8
16
.
https://doi.org/10.1016/j.landurbplan.2015.03.014
.
Zhang
L.
,
Sun
X.
&
Xue
H.
2019
Identifying critical risks in Sponge City PPP projects using DEMATEL method: a case study of China
.
Journal of Cleaner Production
226
,
949
958
.
https://doi.org/10.1016/j.jclepro.2019.04.067
.
Zhang
L.
,
Ye
Z.
&
Shibata
S.
2020
Assessment of rain garden effects for the management of urban storm runoff in Japan
.
Sustainability
12
(
23
),
9982
.
https://doi.org/10.3390/su12239982
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).