Next Article in Journal
Promising and Environmentally Friendly Removal of Copper, Zinc, Cadmium, and Lead from Wastewater Using Modified Shrimp-Based Chitosan
Previous Article in Journal
Differentiating Potential Suspended Sediment Sources Using Radionuclide Tracers and Soil Organic Matter Analysis in a Headwater Catchment in Chuncheon, South Korea
Previous Article in Special Issue
Stormwater Treatment Using Natural and Engineered Options in an Urban Growth Area: A Case Study in the West of Melbourne
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Carbon Emission Intensity of Rainwater Bioretention Facilities

1
School of Environmental and Municipal Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China
2
Key Laboratory of Yellow River Water Environment in Gansu Province, Lanzhou 730070, China
3
Tianshui Housing and Urban-Rural Development Bureau, Tianshui 741000, China
4
Beijing General Municipal Engineering Design & Research Institute Co., Ltd., Beijing 100082, China
*
Authors to whom correspondence should be addressed.
Water 2024, 16(1), 183; https://doi.org/10.3390/w16010183
Submission received: 14 December 2023 / Revised: 29 December 2023 / Accepted: 2 January 2024 / Published: 4 January 2024
(This article belongs to the Special Issue Water Sensitive Urban Design and Decentralised Systems)

Abstract

:
To investigate the quantitative relationship between the volume capture of rainfall and carbon emissions from bioretention facilities, this study introduces the concept of the carbon intensity of volume capture of rainfall. The influence of four key factors—climatic conditions, aquifer height, permeability coefficient, and facility area—was investigated using a residential neighborhood in Tianshui, China, as an example. The results reveal that the carbon intensity value is influenced not only by external environmental changes but also by the inherent attributes of bioretention facilities, such as aquifer height, permeability coefficient, and facility area. The maximum carbon intensity value for the volume capture of rainfall was −0.0005 kg CO2/m3, while the minimum was −0.0852 kg CO2/m3, representing a substantial difference of approximately 169 times. Orthogonal experiments identified the facility area as the most significant influencing factor on carbon intensity, with a correlation coefficient of 0.0520. The area of bioretention facilities can be prioritized to meet deployment requirements, taking into account volume capture reduction effects and carbon emissions. For facilities with a high carbon intensity, an emphasis should be placed on enhancing carbon reduction benefits, and various initiatives can be implemented to achieve this goal.

1. Research Background

Along with the rapid development of urbanization in China, a myriad of challenges is increasingly hindering the sustainable progress of cities [1]. These challenges encompass issues such as water scarcity, recurrent waterlogging, severe water pollution, and a continuous decline in groundwater levels [2,3]. The coexistence of these problems in modern cities has produced significant threats to the fundamental aspects of residents’ lives [4]. China proposed in December 2013 to “build a sponge city with natural accumulation, natural infiltration and natural purification” [5]. A sponge city is inspired by the concept of “Low Impact Development (LID)” [6], a sustainable stormwater management concept that minimizes disasters caused by urban rainfall, prevents water loss, enhances groundwater infiltration and recharge, reduces pollutant discharges, improves urban stormwater storage and utilization, and increases urban resilience [7,8,9]. With the rapid advancement of ecological civilization construction, prioritizing a green stormwater infrastructure to address urban stormwater management issues has become a crucial strategy for emphasizing ecological priorities and fostering green development in sustainable urban planning [10].
A sponge city mainly includes the construction of LID facilities, such as green roofs, bioretention facilities (BFs), grass-planted swales, permeable pavements and tiles, and infiltration drains [11,12,13]. In recent years, a substantial number of studies have evaluated the performance of LID facilities individually or in combination with gray infrastructure [14]. Urban flood modeling and flood source tracking methods were used to identify flood source areas and hazard areas, which offered a more nuanced understanding of the optimal scale for LID facility deployment. Also, various other factors such as different costs [15], different climatic conditions [16,17], etc., were considered at this stage of the study. Davis et al. [18] showed that the quantity of treated water and the effectiveness of water quality are highly dependent on the BF’s length, width, slope, type and density of grasses, soil type, quality of runoff, and marsh design. Shao et al. [19] evaluated the potential for green roof deployment by unmanned aerial vehicle, and their results showed that urban-scale LID facilities require different scales of deployment areas to achieve the required scale. Numerous studies have concluded that LID facilities have different stormwater control capabilities depending on their layout (structure, scale) and the environment in which they are located.
Bioretention facilities are stagnant landscaped areas which are designed to attenuate and treat stormwater runoff [20,21]. The bioretention facilities can reduce runoff and peak runoff flows by retaining captured stormwater for a period of time through their filter media [22]. The enhancement of stormwater infiltration by bioretention facilities is primarily achieved by their ability to delay stormwater runoff during peak infiltration periods, and their infiltration capacity is a critical aspect from a hydrological perspective [23]. The soil medium of the bioretention facilities is the most crucial parameter in reducing runoff [24]. Numerous field studies at different sites have confirmed the ability of bioretention systems to reduce stormwater runoff [6,25]. Additionally, the permeable fill media of bioretention basins play a vital role in reducing stormwater pollution [26]. Ou et al. [10] use bioretention facilities as equivalents in an urban sustainability perspective, based on standardized indicators for volume captured per unit area, to simplify and standardize the assessment and planning process for green stormwater infrastructure. Bioretention facilities provide considerable social and environmental advantages through their effective control of stormwater, pollutant reduction, and reasonable investment, establishing them as crucial components of sponge cities [27].
Bioretention facilities, green roofs, and planted swales function as plant carbon sinks, so it is important to continue to conduct in-depth research on the carbon emissions of LID facilities. Her et al. [28] and Cai et al. [29] demonstrated that the carbon reduction pathways of sponge cities, including bioretention facilities, primarily encompass urban runoff volume capture, pollutant emission reduction, air pollution absorption, plant carbon sequestration, urban heat island mitigation, water reuse, etc. Getter et al. [30] observed that green roofs achieve a more effective carbon sink effect. In a related study, Kavehei et al. [31] quantified the carbon emissions and sequestration potential of various LID facilities, including green roofs, rain gardens, bioretention ponds, vegetated buffer strips, and rainwater ponds. Their findings verified that rain gardens possess the highest carbon capacity. Lin et al. [32] developed a comprehensive accounting model based on the Intergovernmental Panel on Climate Change guidelines and life cycle assessment methodology to predict carbon emissions and carbon sinks in sponge cities. The results indicated that sponge cities can achieve a significant amount of carbon sinks during long-term normal operation. It is anticipated that in the second half of the life cycle, the system can attain carbon neutrality. Subsequently, it can be utilized as a carbon emission reduction system to mitigate the greenhouse effect. Su et al. [33] found that sponge cities can reduce carbon emissions from integrated urban drainage systems by an average of 49%, as observed in a comparison with traditional urban drainage systems in Xiamen. Moore et al. [34] conducted a multidimensional comparison of carbon emissions from various LID facilities, encompassing initial emissions, maintenance emissions, and vegetation sequestration over a 30-year period. Their results also showed that rainwater wetlands and bioretention strips exhibited the lowest net carbon emissions per unit area. However, while the calculation methods exhibit some adaptability across different study areas, accurately assessing the magnitude of their carbon emissions and carbon reduction potential remains challenging. Peng et al. [35] calculated the carbon footprint of rain gardens using the life cycle assessment (LCA) methodology, estimated the operational carbon reduction benefits, and concluded that rain gardens have significant potential to mitigate climate change.
The current study focuses on the carbon emission accounting and carbon emission accounting model, exemplified by studies such as Su et al. [33] and Peng et al. [35]. However, due to variations in study area location, size, and climatic conditions, it is challenging to horizontally compare the size of carbon emissions among different study areas. Similarly, within the same study area, the impact of different LID facility layouts on rainfall volume capture varies, and since the volume capture ratio of annual rainfall is the primary goal of sponge city construction, directly comparing the total amount of carbon emissions may lack meaningful significance. Yet, when considering unit quantities, such as intensity, which are often measured in terms of area, there is limited research on the volume capture ratio of annual rainfall (VCRA), a crucial metric for low-impact development facilities. Therefore, we propose the concept of carbon emission intensity based on the volume capture ratio of annual rainfall. This concept can be employed for the comparison of carbon emissions from LID facilities across different study areas and for the evaluation of carbon emissions from various classes of LID facility placement scenarios within the same study area. However, relevant studies on this concept are yet to be seen.
This study innovates by utilizing the storm intensity formula to derive the calculation formula for the total annual runoff control volume. It introduces the concept of carbon emission intensity for the total runoff control volume throughout the entire life cycle and achieves an evaluation of the combined contribution of bioretention facilities to rainfall and carbon emissions.
In this study, considering the effect of sponge city construction, Tianshui City, China, was studied as an example. The variation in carbon emission intensity of volume capture of rainfall from bioretention facilities over the full life cycle under different influencing factors was investigated. The technical route is shown in Figure 1. The results can serve as a reliable foundation for assessing the carbon emissions that come from LID facilities. Additionally, the results can aid in the creation of sponge cities, contributing to their development into more sustainable cities and the implementation of China’s “double carbon” strategy.

2. Overview of the Study Area and Research Methods

2.1. Overview of the Study Area

The study area is in a newly built residential area in Tianshui, Gansu Province. It is a catchment area within the residential area, and the study simulates the performance of a bioretention facility in this catchment area. The study area is 2930 m2, consisting of built-up areas, green spaces, and roads (including concrete and paved roads). Stormwater pipes and outlets were set up based on the elevation and pipe layout of the study area. The natural drainage direction is from northwest to southeast, as shown in Figure 2.

2.2. Model Construction and Research Methods

2.2.1. Design Rainfall Intensity

The rainfall intensity in Tianshui was calculated by Equations (1) and (2):
q = 712.900 ( 1 + 1.90 × lg ( p ) ) ( t + 8.711 ) 0.742 ( p < 20   a )
q = 1336.703 ( 1 + 1.96 × lg ( p ) ) ( t + 12.94 0 ) 0.840 ( p 20   a )
where q represents the rainfall intensity in L/(s·ha), p represents the design rainfall return period in a, and t represents the rainfall duration in min.
Using the Chicago Rainfall Type Generation Software (1.0), a rainfall intensity curve with a peak coefficient of 0.398 was generated based on Equation (1) for design rainfall return periods of 0.5 a, 1 a, 2 a, 3 a, 5 a, and 10 a, with a rainfall duration of 120 min (Figure 3). Similarly, a long-duration rainfall curve with a rainfall duration of 1440 min was generated based on Equation (2) for design rainfall return periods of 20 a, 30 a, and 50 a, as shown in Figure 4.

2.2.2. InfoWorks ICM Model Construction and Parameter Calibration

The InfoWorks ICM (2023.1) was selected for simulation. This software is a comprehensive urban watershed drainage model that can simulate the entire urban rainwater cycle system. It can simulate and evaluate the performance of various rainwater storage and utilization facilities, providing technical support for preventing urban waterlogging and controlling and utilizing rainwater resources [36]. The rainfall–runoff model was constructed based on the topography and design drawings of the residential area, and it was generalized into 8 sub-catchment areas, 16 nodes, and 16 pipes. The study area is one of the sub-catchment areas, as shown in Figure 1.
The runoff coefficient method is the parameter rate determination method of the urban stormwater model selected in this paper [37]. According to the “Standard for Design of Outdoor Wastewater Engineering (GB 50014-2021)” and considering the characteristics of the site’s underlying surface and soil, the runoff coefficients for roof surfaces, concrete roads, paved roads, and green spaces were set as 0.90, 0.85, 0.60, and 0.15, respectively. The comprehensive runoff coefficient for the site was calculated through area weighting, resulting in a value of 0.66. After establishing the model for the site, a simulation verification was conducted using a rainfall of 26.58 mm with a two-year return period under normal rainfall conditions. The resulting runoff coefficient for the site was also 0.66, confirming the rationality of the model parameters that were set in this study.

2.2.3. Experimental Design

The basic parameters of the BF were determined based on the relevant specifications for sponge city construction [38] and the existing literature [28,39]. The selected bioretention facility is an anti-seepage-shaped BF, referred to as the BF. The influencing factors and scope of bioretention facilities are defined by domestic regulations and relevant international norms. These guidelines emphasize the significance of considering aquifer height, permeability coefficient, and facility area as crucial factors in the construction process [40,41]. Additionally, within the context of climate change, the study of the impact of climate conditions on bioretention facilities has been introduced. The climate conditions varied within ±20% of the rainfall amount, calculated using the local rainfall intensity formula with a step size of 10%. The ranges for the aquifer height and permeability coefficient were set at 100 mm to 300 mm (with a step size of 50 mm) and 50 mm/h to 150 mm/h (with a step size of 25 mm/h), respectively. The facility areas are 45.5 m2, 72.8 m2, 91 m2, 136.5 m2, 182 m2, and 273 m2, which corresponds to 5%, 8%, 10%, 15%, 20%, and 30% of the proportion of the facility area to the green area. The corresponding values for each influencing factor are shown in Table 1.
In order to determine the value of the carbon intensity and carbon reduction effect under different levels of single influencing factors, 19 simulation scenarios under different influencing factors were established through the method of single-factor experiment (scenario numbers are 1–18 and 36). Afterwards, 24 scenarios were established through orthogonal experiments, which were scenarios 19–42, to determine the strength of the correlation between the individual influencing factors on the carbon intensity and carbon reduction benefits. The details of each scenario are shown in Table 2.

2.2.4. Calculation Method for Rainwater Control Performance

The rainwater control performance is represented by the volume capture ratio of annual rainfall (VCRA) and the volume capture of annual rainfall (VCA).
The VCRA is calculated as the weighted average of the difference between the rainfall and runoff volumes, as shown in Equation (3). The VCA is calculated as the VCRA divided by the annual rainfall volume, as shown in Equation (4).
V C R A = j = 1 9 V r a i n , j V r u n o f f , j P j A j = 1 9 H j P j × 100 %
V C A = V C R A × H a × A
where V C R A represents the volume capture ratio of annual rainfall in %, P j represents the return period (j = 1–9 represents return periods of 0.5, 1, 2, 3, 5, 10, 20, 30, and 50 years), V r a i n , j represents the rainfall volume for a return period of P j in m3, V r u n o f f , j represents the runoff volume generated when the LID facility is in place for a return period of P j (i.e., the overflow from the BF) in m3, H a is the annual average rainfall volume in the study area in mm (set as 505 mm), H j indicates the rainfall intensity for a return period of P j in mm, V C A represents the volume capture of annual rainfall in m3, and A indicates the size of the study area in m2 (set as 2930 m2). The value of V r u n o f f , j was calculated using the InfoWorks ICM model.

2.2.5. Full Life Cycle Costing Calculation Method

This study only analyzed the full life cycle costing of LID facilities, including construction, operation and maintenance, and removal costs, without considering other required costs (Figure 5). The life cycle of the BF is set to 30 years in this paper. The calculation formula refers to existing research results [42,43,44,45,46], such as Equations (5)–(8):
F L C C = C C a p i t a l + t = 1 n k P V O & M k , t + C r e m o v e
P V O & M t = F V O & M t 1 + d t ,   t
F V O & M t = C C a p i t a l × p × 1 + r t ,   t
C d i s m a n t l e = y × C C a p i t a l S
where F L C C represents the full life cycle costing of the BF, in CNY. n k refers to the life cycle of the LID facility in years, taking n k = 30. C C a p i t a l is the construction cost of the facility in CNY, and the unit price is referenced from relevant research [47] and engineering examples in the northwest region of China. P V O & M t represents the present value of the maintenance cost of the LID facility in year t in CNY. F V O & M t represents the future value of the maintenance cost of the LID facility in year t in CNY. d is the discount rate in %, taking 8% as the current social discount rate. r is the average inflation rate in %, taking 3%. p is the proportion of annual operation and maintenance costs to construction costs in %, taking p = 8. C d i s m a n t l e is the cost of facility dismantle in CNY. y is the proportion of removal cost to construction cost in %, taking 10%. S is the salvage value, taking S = 0.

2.2.6. Full Life Cycle Carbon Emission Accounting Method

The full life cycle carbon emission accounting of LID facilities includes carbon emissions during the construction, operation and maintenance, and dismantling stages of the LID facility (Figure 6), calculated using the carbon emission factor method. The life cycle carbon emissions are calculated using Equations (9)–(12).
C E = 30 n k C E C S + C E O M + C E D S
C E C S = A D × E F C S
C E O M = V C A × E F O + A D × E F M A D × E F H
C E D S = A D × E F D S
where C E represents the full life cycle carbon emissions of the BF in kg CO2. C E C S and C E D S represent the carbon emissions during the construction stage and dismantling stage of the LID facility in kg CO2. C E O M represents the carbon emission during the operation and maintenance stage of the LID facility’s full life cycle in kg CO2. A D represents the area of the BF in m2. E F C S and E F D S represent the carbon emission factors during the construction and dismantling stages of the LID facility in kg CO2/m2. E F O represents the operational carbon emission factor of the LID facility during the full life cycle in kg CO2/m3. E F M represents the maintained carbon emission factor of the LID facility during the full life cycle in kg CO2/m2. And E F H represents the carbon sink factor of the carbon sequestration for the green space of the LID facility full life cycle in kg CO2/m2. The values of some constants in the formulas are shown in Table 3.

2.2.7. Concept and Calculation Method of Carbon Intensity of VCA Based on Full Life Cycle

Based on the LCA theory, the carbon emission intensity of VCA, defined as the ratio of carbon emissions of the facility over the runoff control volume throughout its life cycle, is proposed. It is used to quantitatively evaluate the relationship between the rainwater control capacity and carbon emissions of LID facilities. It is calculated according to Equations (13) and (14). A smaller value of this indicator indicates fewer carbon dioxide emissions per unit of controlled runoff, leading to better environmental benefits.
E = C E C V 30
C V 30 = 30 × V C A
where E represents the carbon emission intensity of the runoff control volume in kg CO2/m3 of volume capture of annual rainfall. C V 30 is the VCA of full life cycle in m3. Other symbols have the same meanings as before.

2.2.8. Concept and Calculation Method of Carbon Reduction Benefit Based on Full Life Cycle

The concept of a carbon reduction benefit is proposed based on the cost-effectiveness [51] calculation method and the LCA theory. It refers to the ratio of carbon reduction of the LID facility over its cost throughout its life cycle. It is used to evaluate the relationship between the cost and carbon emission reduction of the LID facility and is calculated according to Equations (15) and (16). A higher value of carbon reduction benefit indicates higher economic effectiveness and value of carbon emission reduction in the construction of LID facilities.
B = C E R L C C
C E R = A D × E F H + V C A × E F H L
where B represents the carbon reduction benefit in kg CO2/CNY, C E R represents the carbon emission reduction over the full life cycle of the LID facility, and the value is equal to the carbon emission reduction generated by activities such as carbon sequestration in green spaces, rainwater utilization, runoff reduction, and rainwater purification in kg CO2, and E F H L indicates the combined carbon emission factor for rainwater utilization and rainwater purification over the full life cycle of the LID facility in kg CO2/m3. Other symbols have the same meanings as before. The values of some constants in the formulas are shown in Table 3.

3. Results and Discussion

3.1. Compositional Analysis of Carbon Emission from BF in Full Life Cycle

The proportion of carbon emissions from stages or activities of the full life cycle of the BF in scenario 32 is shown in Figure 7. According to Equations (9)–(12), the carbon emission in the full life cycle can be calculated as follows: 5895.33 kg CO2 for the construction stage, 921.56 kg CO2 and 700.25 kg CO2 for the operation activity and maintenance activity in the operation and maintenance stage, and 80.54 kg CO2 for the dismantling phase. In addition, we calculated −9131.85 kg CO2 for the carbon sink of the full life cycle. The pure carbon emission (sum of carbon emissions and sinks) during the operation and maintenance stage is −7524.16 kg CO2, and the pure carbon emission during the full life cycle is −1458.10 kg CO2. The carbon emission accounting results of the full life cycle of bioretention facilities in this study differ from some existing research results [32,52] because of the inconsistency in the definition of carbon sinks in the boundary of full life cycle accounting. The related studies’ results include carbon sequestration in green spaces, rainwater purification, and rainwater utilization of the operation and maintenance phase of the bioretention facilities as carbon sinks in the accounting system. In this paper, referring to the relevant literature [53], only the activities that produce resources or energy in the production process and can be delivered to the outside are defined as carbon sinks, and the actions that meet this condition are the greenfield carbon sequestration in the operation and maintenance stage of bioretention facilities. Other activities in the operation and maintenance phase of bioretention facilities, such as rainwater purification and utilization, do not satisfy this study’s definition of carbon sinks. Together with the sink activities, they are defined as carbon emission reductions, which is outside the boundary of the full life cycle accounting and is discussed separately in Section 3.3.

3.2. Analysis of Carbon Emission Intensity under Different Influencing Factors

3.2.1. Climate Condition

As shown in Figure 8a, climate change affects rainwater control performance: as the intensity of storms increases, the VCA of the facility increases, which shows that the BF has a stormwater control capacity for increasing rainfall. The life cycle carbon emissions under different climatic conditions ranged from −874.1175 kg CO2 to −657.3433 kg CO2. The carbon intensity ranged from −0.0306 kg CO2/m3 to −0.0176 kg CO2/m3, with a changing trend that is consistent with the life cycle carbon emissions and a variation range of 74%. The full life cycle carbon emissions gradually increased with the increase in rainfall because of the constant parameters of the BF. Carbon emissions and carbon sinks stay unaltered during the construction and demolition stages. However, the change in climatic conditions led to a shift in rainfall, and the organic substance in the corresponding runoff was treated by the BF, which increased the direct carbon emissions during the operation and maintenance stage, increasing life cycle carbon emissions.

3.2.2. Aquifer Height

As shown in Figure 8b, the larger the aquifer height is, the stronger the ability of rainwater capture is, which is mainly reflected in the increase in VCA. According to the increasing height of the aquifer, VCA showed an increasing trend, and the higher the height of the aquifer was, the more significant the benefit of the facility for runoff control was, which is similar to the results of the study by Tu et al. [54]. As rainwater runoff flows into the BF, one part infiltrates through the soil layer, and the other part is stored in the aquifer. The higher the height of the aquifer is, the more stormwater runoff can be held in the facility, and the greater the storage capacity for stormwater runoff is. The life cycle carbon emissions of the BF increased from −827.2060 kg CO2 at an aquifer height of 100 mm to −695.6845 kg CO2 at an aquifer height of 300 mm. The carbon emission intensity increased monotonically with an increase in aquifer height, ranging from −0.0271 kg CO2/m3 to −0.0194 kg CO2/m3, with a variation range of 40%. The increased aquifer height allows the BF to retain more rainwater during the operation and maintenance stage. Therefore, the amount of rainwater purification is high, and the direct carbon emissions during the operation phase are correspondingly high.

3.2.3. Permeability Coefficient

As shown in Figure 8c, the PC of the BF follows the same trend as the stormwater control performance. The PC resulting from 10–100 mm/h changes in the stormwater control effect have the greatest impact. With the increase in the infiltration coefficient, the impact of bioretention facilities on the stormwater volume capture effect gradually decreases. This is because during rainfall, as runoff rainwater flows into the soil void, the soil infiltration coefficient gradually decreases and eventually reaches a stabilized value. Higher initial PCs can slow this process, and the runoff control effect improves. This is similar to the findings of Pan et al. [55]. The life cycle carbon emissions of the BF ranged from −886.6840 kg CO2 to −676.5319 kg CO2 under soil PCs of 50 mm/h to 150 mm/h. The carbon emission intensity values increased monotonically with an increase in soil PC, ranging from −0.0316 kg CO2/m3 to −0.0185 kg CO2/m3. The carbon emission intensity at a 200 mm/h PC increased by 71% compared to that at a 50 mm/h permeability coefficient because changes in soil permeability affected the stormwater management capacity of the BF, as demonstrated by Haaland et al. [56]. This, in turn, affected the amount of stormwater that was purified by the BF. And it could also affect carbon emissions during the operation and maintenance stage. In the life cycle assessment boundary of this study, changes in soil permeability had no effect on the variance in carbon emission intensity during the construction and demolition stages. Therefore, an increase in soil PC led to a rise in life cycle carbon intensity.

3.2.4. The Facility Area

The FA will directly affect the facility’s water intake, which is one of the leading design parameters of the BF. As seen in Figure 8d, the greater its value is the better the stormwater control performance eventually becomes. This is because the change in the layout area does not affect the rate of soil infiltration, and overall, with the increase in the layout area ratio, the inlet flow decreases, and the facility’s ability to regulate runoff tends to increase. The full life cycle carbon emissions increased monotonically with an increase in the ratio of the BF area to the green space area. When the BF accounted for 30% of the green space area, the life cycle carbon emission was −0.0852 kg CO2/m3. The carbon emission intensity changes with a runoff reduction volume from 5% to 30% follow a diverse trend in terms of the life cycle carbon emissions, which decrease from −0.0040 kg CO2/m3 to −0.0852 kg CO2/m3 with a variation range of 95%. This phenomenon may be explained via the life cycle theory. The increase in the facility area lead to an increase in the carbon emission intensity during the construction, operation and maintenance, and dismantling stages. At the same time, the expansion of the bioretention facility area results in an increase in the amount of greenery and impermeable surfaces, which strengthens the carbon sink effect [57]. The increased carbon sink throughout the life cycle can offset part of the increased carbon emission intensity in the construction, operation and maintenance, and demolition phases, slowing down the trend of increasing full life carbon emissions. Therefore, the larger the area of bioretention facilities is, the more significant the carbon sequestration effect is and the lower the carbon emissions are.

3.3. Analysis of Carbon Reduction Benefits of Different Influencing Factors

3.3.1. Climate Condition

Figure 9a shows that during the change from climate condition C−20 to C+20, the carbon emission reduction throughout the life cycle of the BF ranged from 4.2899 × 104 kg CO2 to 5.4251 × 104 kg CO2. The carbon reduction benefit increased from 1.5363 kg CO2/CNY in the C−20 period to 1.9429 kg CO2/CNY in the C+20 period, with an increase of 30.6%. The influence of changes in the climate condition on different carbon emission reduction measures is observed, as the carbon emission reduction from green space carbon sequestration remains constant due to the constant area of deployment. The carbon emission reduction from rainwater utilization is the highest, the rainwater purification carbon emission reduction is the second highest, and both increase with the increase in rainfall. These findings suggested that variations in external CSs provided an increase in the BF’s brilliant carbon reduction capability. This is because the facility’s structure, including storage, infiltration, and drainage layers, can preserve rainwater and indirectly contribute to carbon reduction through factors such as rainwater utilization. This observation is consistent with the findings of Cai [29]. Taking C+10 as an example, under this climatic condition, the BF can achieve a carbon reduction benefit of 1.8422 kg CO2 per unit cost, indicating good economic value.

3.3.2. Aquifer Height

Figure 9b shows that the carbon reduction throughout the full life cycle of the BF gradually increased with the increase in aquifer height, ranging from 4.5356 × 104 kg CO2 to 5.2244 × 104 kg CO2. The carbon reduction benefit increased gradually as the aquifer height increased from 1.6243 kg CO2/CNY for an aquifer height of 100 mm to 1.8710 kg CO2/CNY for an aquifer height of 300 mm. When compared to a storage layer height of 100 mm, each increase of 50 mm can result in an increase of approximately 5.0% to 14.4% in carbon reduction benefits, because the BF can retain more rainwater and contribute to increased carbon reduction through rainwater purification and utilization, resulting in a higher carbon reduction potential. Increasing the height of the aquifer will significantly reduce the use of tap water and conserve rainwater resources [29]. This result was also demonstrated by Moore et al. [58].

3.3.3. Permeability Coefficient

Figure 9c shows that as the PC increased, the carbon emission reduction throughout the full life cycle of the BF increased, and the carbon reduction benefit increased from 1.5128 kg CO2/CNY to 1.9069 kg CO2/CNY, with an increase of 21%. With the increase in the PC of the BF, permeability pairs affected the carbon reduction capacity. The infiltration rate of rainwater increased, and the rainwater infiltrated into the drainage layer and was discharged through blind pipes, which can reduce the total runoff and peak flow, resulting in increased carbon sequestration. Plants and the composition of the soil medium are significant factors influencing infiltration coefficients and are critical to the ability of bioretention areas to remove pollutants from stormwater and retain stormwater. Therefore, it can be inferred that carbon reduction activities are equally important for both stormwater purification and stormwater utilization in these areas [59].

3.3.4. The Facility Area

Figure 9d shows that as the FA increased, the carbon emission reduction throughout the life cycle of the BF gradually increased from 3.8515 × 104 kg CO2 to 7.3734 × 104 kg CO2. This is due to an increase in the area of bioretention facilities, an increase in the carbon sink capacity of the greenfield carbon sink, and an increase in the carbon emission reduction capacity of rainwater purification and rainwater utilization, all of which contribute to a gradual increase in full-life carbon emission reductions. The carbon reduction benefits for a layout area ratio of 5% to 30% are 2.7586 kg CO2/CNY, 2.0103 kg CO2/CNY, 1.7499 kg CO2/CNY 1.3560 kg CO2/CNY, 1.1460 kg CO2/CNY, and 0.8802 kg CO2/CNY, and the carbon reduction benefit gradually decreases with the increase in the deployment area, with a decrease of 213%. Compared to the growth in carbon reduction across the life cycle, the increase in the FA resulted in a more noticeable increase in the overall cost of the life cycle. Therefore, the layout area’s ratio has a substantial impact on the overall life cycle cost, which is consistent with the conclusions drawn by Wang et al. [60].

3.4. Orthogonal Experiment

3.4.1. Carbon Intensity

Figure 10 illustrates the outcomes in terms of carbon emission of the volume capture and carbon intensity throughout the full life cycle for various scenarios related to the bioretention facility. These scenarios correspond to numbers 19–42, as outlined in Table 2.
In the 25 sets of scenarios, the reduction in the full life cycle runoff control values varied between 16,582.2064 m3 and 50,697.7611 m3. Additionally, the full life cycle carbon emissions exhibited a range from −2357.9866 kg CO2 to −14.2365 kg CO2, while the carbon intensity of the volume capture of rainfall ranged from −0.0728 kgCO2/m3 to −0.0005 kgCO2/m3. To validate the impacts of various influencing factors on the bioretention facility’s performance, orthogonal experiments were conducted using a level difference analysis, and the results are shown in Table 4.
According to the results of the extreme difference analysis for carbon intensity in Table 4, the analysis reveals that the coefficients (R) for different factors—climatic condition, aquifer height, permeability coefficient, and the facility area—are 0.0150, 0.0090, 0.0150, and 0.0520, respectively. A higher R value in the extreme difference indicates a greater impact of the corresponding factor.
The influence of the four factors on the volume capture and carbon emission intensity is analyzed in the following order: the facility area, permeability coefficient, climate condition, and aquifer height. Therefore, in the design of bioretention facilities, priority can be given to the area ratio of bioretention facilities, considering its significant impact on the volume capture and carbon emissions. This study also found that, in comparison to other design characteristics, the area occupied by bioretention facilities has a greater impact on the intensity of carbon emissions.

3.4.2. Carbon Reduction Benefit

The full life cycle carbon reduction and carbon reduction benefit of different experimental schemes of bioretention facilities are shown in Figure 11. The full life cycle carbon reduction is 2.4385 × 104–7.7423 × 104 m3, and the full life cycle carbon reduction benefits are 0.8802–3.1223 kg CO2/CNY. To validate the impact of different influencing factors on the bioretention facilities, an analysis of the orthogonal experiment was carried out to perform the level difference analysis, and the results are shown in Table 5.
According to the results of the extreme difference analysis for carbon emission intensity in Table 5, it is observed that the R of the different factors climatic condition, aquifer height, permeability coefficient, and the facility area are 0.4690, 0.2700, 0.4160, and 1.5020, respectively. Specifically, the order is the facility area, permeability coefficient, climate condition, and aquifer height. Consequently, in the design of bioretention facilities, expanding the installation area complements efforts to enhance both the total runoff control and carbon emission reduction.

3.5. Analysis of the Relationship between VCRA and Carbon Emission Intensity

The VCRA and carbon intensity for scenarios 1–42 were analyzed by categorizing the test results into 11 groups based on preopening and closing intervals. The intervals were determined by starting from 45% and progressing in 5% increments. The data were then analyzed in Figure 12 using box-and-whisker plots.
As depicted in Figure 12, the distribution of VCRA varied widely, spanning from 45% to 100%. Data points within different intervals of VCRA exhibited a non-uniform distribution, with the intervals of 65% to 70% and 80% to 85% having the largest number of distributions, each containing eight test scenarios. Within the 65% to 70% interval, scenario 17 exhibited the minimum carbon intensity at −0.0852 kgCO2/m3, while scenario 39 reached the maximum carbon intensity of −0.0005 kgCO2/m3. In the 80% to 85% interval, carbon intensity increased from −0.0553 kgCO2/m3 for scenario 33 to −0.0185 kgCO2/m3 for scenario 12. The intervals with the fewest scenarios were 45~50% and 50~55%, each with only one test scenario. The outlier point (in Figure 12, point a) is attributed to the constancy of carbon emissions during the two major phases of construction and demolition when the area of the bioretention facility remains constant. The impact of other factors influencing stormwater retention capacity on the carbon intensity gradually becomes more apparent. This is particularly notable in scenarios where the infiltration coefficient undergoes orders of magnitude-sized changes, leading to a reduction in the amount of stormwater treatment. The dual factor changes result in a noticeable effect on carbon intensity. Therefore, this point is disregarded in the analysis.
These findings suggest two key points. Firstly, certain bioretention facility scenarios can effectively capture rainwater volume within a specified range. Some scenarios meet local requirements, with the VCRA exceeding 80% and eliminating the need for additional LID facilities in the study area. Secondly, within the same VCRA interval, carbon intensity values varied significantly among schemes, indicating that different scenarios with the same range yield different carbon emission effects. Therefore, studying how to minimize carbon intensity while ensuring the VCRA is a direction of research that is worthy of exploration. The carbon intensity exhibits a decreasing trend with the increase in the VCRA. As the deployment area significantly influences carbon intensity throughout the entire life cycle, the schemes in the 75~80% interval have the same deployment area, resulting in a stable carbon emission intensity. The interval of 60~65% experiences the largest fluctuation in carbon intensity, with scenario 25 demonstrating the maximum value at 182 units. Kaykhosravi et al. [61] demonstrated that BFs were ranked highest for their social, economic, and environmental benefits, with environmental benefits centered around factors such as the amount of pollutants treated. This study also revealed that bioretention facilities exhibit a greater potential in terms of carbon intensity. In this interval, the infiltration coefficient and aquifer height are at lower values, highlighting the significant impact of the deployment area on the carbon emission intensity [33].

3.6. Relationship between FLCC and Bioretention Facility Performance in Terms of Carbon Emissions

A comparative analysis of the size and cost of bioretention facilities and all of their capacities in different scenarios was carried out. The results of the analysis are shown in Figure 13. Figure 13a illustrates that the water control capacity of the BF gradually increases with the rise in cost, leading to a simultaneous decrease in the full life carbon emission and carbon emission intensity. The primary driver for the increased cost is the expansion of the facility area. The figure illustrates a direct proportional relationship between the size of the BF and the cost. The relationship between the cost and the full life cycle carbon emissions is inverse, and it can be learned from Figure 13b,c that if the area of the bioretention facility deployment increases, the volume of water that can be captured by the facility during rainfall increases, the area of green space can lead to carbon sinks, and the amount of carbon sequestered increases with the increase in the area of deployment. These carbon sinks effectively offset the heightened carbon emissions incurred during the construction, operation, and maintenance of bioretention facilities. Han et al. [62] indicate that the larger the size of the BF, the higher the social benefits. However, their findings did not include results related to carbon emissions within the context of social benefits. Furthermore, the change in cost and carbon intensity displayed an inverse proportionality, indicating that the carbon emissions per unit of water controlled by bioretention facilities decreased as the deployment area increased.

3.7. Relationship between Carbon Emission Intensity and Carbon Emission Reduction

In order to further explore the relationship between carbon emission intensity and carbon reduction benefits, this study conducted a linear fitting of the two quantities for all scenarios, and the results are shown in Figure 14.
As illustrated in Figure 14, there is a well-established fitting relationship between carbon intensity and carbon reduction benefit. This relationship is attributed to the higher proportion of costs that are associated with increased volume capture of rainfall per unit area, indicating a positive correlation between water volume control and carbon intensity. As carbon emissions increased, the carbon abatement benefits of bioretention facilities improved. These findings align with the results reported by Su et al. [33] in the literature. Therefore, stronger carbon abatement methods and emphasizing an intensified approach are recommended for bioretention facilities with greater carbon emission intensity. Attention should be directed towards optimizing carbon abatement effectiveness in these scenarios.

4. Conclusions and Prospects

4.1. Conclusions

In this study, the carbon emissions of bioretention facilities in different scenarios were calculated using the full life cycle model, and the carbon intensity of bioretention facilities in the full life cycle was evaluated using the carbon intensity of volume capture of rainfall. The influencing factors were analyzed to refine scenarios and optimize facilities. The results are shown in the following:
(1)
The carbon intensity of the volume capture of rainfall effectively assesses the carbon emission levels of bioretention facilities, providing a theoretical foundation for the study of carbon emissions in sponge cities.
(2)
The carbon intensity value ranges from a maximum of −0.0005 kg CO2/m3 to a minimum of −0.0852 kg CO2/m3, exhibiting a significant difference of approximately 169 times. This value is not only affected by the external environmental changes, but also by the bioretention facility’s own attributes such as the aquifer height, permeability coefficient, and facility area.
(3)
The results of orthogonal experiments show that the strongest influence on the carbon intensity of the volume capture of rainfall is the facility area, with a correlation coefficient of 0.0520. Under the consideration of the total runoff reduction effect and the carbon emission situation, the bioretention facilities can be prioritized by adjusting the deployment area to satisfy the requirements of the deployment.
(4)
The maximum carbon reduction benefit of bioretention facilities is 3.1223 kg CO2/CNY, differing approximately 2.55 times from the minimum value of 0.8802 kg CO2/CNY. For bioretention facilities with a higher carbon emission intensity, emphasis should be placed on carbon emission reduction efforts, and various initiatives can be implemented to enhance their carbon reduction benefits.

4.2. Prospects

This study has focused on a singular bioretention facility, necessitating additional research to substantiate the carbon intensity of LID facilities. The research directions encompass:
(1)
Investigating the varied impact of different LID facilities on rainwater control, prompting further exploration of carbon intensity for individual LID facilities.
(2)
Conducting a study on the carbon intensity of combined LID arrangements at the parcel level, capitalizing on their synergistic effect in enhancing rainfall and flood control.
(3)
Climate conditions exert a significant influence on stormwater runoff capture at LID facilities. Employing more accurate climate prediction methods can facilitate research on the carbon emission intensity across various climate conditions.
(4)
Constructing a carbon emission model for LID facilities based on data from prior studies. The model will consider varying climatic conditions, utilizing the total runoff control rate as a target. This exploration aims to unveil the potential for carbon emission reduction and strategies for sponge city construction.

Author Contributions

Conceptualization, B.W.; Methodology, B.W. and D.W.; Formal Analysis, D.W. and H.C.; Investigation, H.C., W.L. and L.C.; Data Curation, X.L. and H.L.; Writing—Original Draft Preparation, D.W.; Writing—Review and Editing, T.Z.; Supervision, X.W. and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

The National Natural Science Foundation of China (No. 52060014); National Key Research and Development Program of China (No. 2018YFE0206200); Gansu Province Construction Science and Technology Project (No. JK2023-18); Gansu Science and Technology Plan (No. 21JR7RA340).

Data Availability Statement

Data are available upon reasonable request to the corresponding author.

Conflicts of Interest

Author Chen Hai was employed by the company Tianshui Housing and Urban-Rural Development Bureau. Authors Li Wei, Cao Lianbao, and Liu Jianlin were employed by the company Beijing General Municipal Engineering Design & Research Institute Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Chen, D.; Hu, W.; Li, Y.; Zhang, C.; Lu, X.; Cheng, H. Exploring the temporal and spatial effects of city size on regional economic integration: Evidence from the Yangtze River Economic Belt in China. Land Use Policy 2023, 132, 106770. [Google Scholar] [CrossRef]
  2. Cui, T.; Long, Y.; Wang, Y. Choosing the LID for Urban Storm Management in the South of Taiyuan Basin by Comparing the Storm Water Reduction Efficiency. Water 2019, 11, 2583. [Google Scholar] [CrossRef]
  3. Chen, D.; Lu, X.; Hu, W.; Zhang, C.; Lin, Y. How urban sprawl influences eco-environmental quality: Empirical research in China by using the Spatial Durbin model. Ecol. Indic. 2021, 131, 108113. [Google Scholar] [CrossRef]
  4. Azari, B.; Tabesh, M. Urban storm water drainage system optimization using a sustainability index and LID/BMPs. Sustain. Cities Soc. 2022, 76, 103500. [Google Scholar] [CrossRef]
  5. Tang, S.; Jiang, J.; Zheng, Y.; Hong, Y.; Chung, E.-S.; Shamseldin, A.Y.; Wei, Y.; Wang, X. Robustness analysis of storm water quality modelling with LID infrastructures from natural event-based field monitoring. Sci. Total Environ. 2021, 753, 142007. [Google Scholar] [CrossRef] [PubMed]
  6. Tang, S.; Luo, W.; Jia, Z.; Liu, W.; Li, S.; Wu, Y. Evaluating Retention Capacity of Infiltration Rain Gardens and Their Potential Effect on Urban Stormwater Management in the Sub-Humid Loess Region of China. Water Resour. Manag. 2016, 30, 983–1000. [Google Scholar] [CrossRef]
  7. Yang, D.; Zhao, X.; Anderson, B.C. Integrating Sponge City Requirements into the Management of Urban Development Land: An Improved Methodology for Sponge City Implementation. Water 2022, 14, 1156. [Google Scholar] [CrossRef]
  8. Zhang, S.; Li, Y.; Ma, M.; Song, T.; Song, R. Storm Water Management and Flood Control in Sponge City Construction of Beijing. Water 2018, 10, 1040. [Google Scholar] [CrossRef]
  9. Yang, Y.; Li, J.; Huang, Q.; Xia, J.; Li, J.; Liu, D.; Tan, Q. Performance assessment of sponge city infrastructure on stormwater outflows using isochrone and SWMM models. J. Hydrol. 2021, 597, 126151. [Google Scholar] [CrossRef]
  10. Ou, J.; Li, J.; Li, X.; Zhang, J. Planning and Design Strategies for Green Stormwater Infrastructure from an Urban Design Perspective. Water 2024, 16, 29. [Google Scholar] [CrossRef]
  11. Dutta, A.; Torres, A.S.; Vojinovic, Z. Evaluation of Pollutant Removal Efficiency by Small-Scale Nature-Based Solutions Focusing on Bio-Retention Cells, Vegetative Swale and Porous Pavement. Water 2021, 13, 2361. [Google Scholar] [CrossRef]
  12. Liu, Z.; Li, W.; Wang, L.; Li, L.; Xu, B. The scenario simulations and several problems of the Sponge City construction in semi-arid loess region, Northwest China. Landsc. Ecol. Eng. 2022, 18, 95–108. [Google Scholar] [CrossRef]
  13. Raimondi, A.; Becciu, G. Performance of Green Roofs for Rainwater Control. Water Resour. Manag. 2021, 35, 99–111. [Google Scholar] [CrossRef]
  14. Qi, W.; Ma, C.; Xu, H.; Chen, Z.; Zhao, K.; Han, H. Low Impact Development Measures Spatial Arrangement for Urban Flood Mitigation: An Exploratory Optimal Framework based on Source Tracking. Water Resour. Manag. 2021, 35, 3755–3770. [Google Scholar] [CrossRef]
  15. Lee, H.; Woo, W.; Park, Y.S. A User-Friendly Software Package to Develop Storm Water Management Model (SWMM) Inputs and Suggest Low Impact Development Scenarios. Water 2020, 12, 2344. [Google Scholar] [CrossRef]
  16. Sun, Y.-w.; Pomeroy, C.; Li, Q.-y.; Xu, C.-d. Impacts of rainfall and catchment characteristics on bioretention cell performance. Water Sci. Eng. 2019, 12, 98–107. [Google Scholar] [CrossRef]
  17. Yang, Y.; Chui, T.F.M. Optimizing surface and contributing areas of bioretention cells for stormwater runoff quality and quantity management. J. Environ. Manag. 2018, 206, 1090–1103. [Google Scholar] [CrossRef] [PubMed]
  18. Davis Allen, P.; Hunt William, F.; Traver Robert, G.; Clar, M. Bioretention Technology: Overview of Current Practice and Future Needs. J. Environ. Eng. 2009, 135, 109–117. [Google Scholar] [CrossRef]
  19. Shao, H.; Song, P.; Mu, B.; Tian, G.; Chen, Q.; He, R.; Kim, G. Assessing city-scale green roof development potential using Unmanned Aerial Vehicle (UAV) imagery. Urban For. Urban Green. 2021, 57, 126954. [Google Scholar] [CrossRef]
  20. Mahmoud, A.; Alam, T.; Yeasir, A.; Rahman, M.; Sanchez, A.; Guerrero, J.; Jones, K.D. Evaluation of field-scale stormwater bioretention structure flow and pollutant load reductions in a semi-arid coastal climate. Ecol. Eng. 2019, 142, 100007. [Google Scholar] [CrossRef]
  21. Rodrigues, A.L.M.; da Silva, D.D.; de Menezes Filho, F.C.M. Methodology for Allocation of Best Management Practices Integrated with the Urban Landscape. Water Resour. Manag. 2021, 35, 1353–1371. [Google Scholar] [CrossRef]
  22. Zhang, X.; Chen, L.; Guo, C.; Jia, H.; Shen, Z. Two-scale optimal management of urban runoff by linking LIDs and landscape configuration. J. Hydrol. 2023, 620, 129332. [Google Scholar] [CrossRef]
  23. Johnson Jeffrey, P.; Hunt William, F. Field Assessment of the Hydrologic Mitigation Performance of Three Aging Bioretention Cells. J. Sustain. Water Built Environ. 2020, 6, 04020017. [Google Scholar] [CrossRef]
  24. Liu, J.; Sample, D.J.; Bell, C.; Guan, Y. Review and Research Needs of Bioretention Used for the Treatment of Urban Stormwater. Water 2014, 6, 1069–1099. [Google Scholar] [CrossRef]
  25. Alikhani, J.; Nietch, C.; Jacobs, S.; Shuster, B.; Massoudieh, A. Modeling and Design Scenario Analysis of Long-Term Monitored Bioretention System for Rainfall-Runoff Reduction to Combined Sewer in Cincinnati, OH. J. Sustain. Water Built Environ. 2020, 6, 04019016. [Google Scholar] [CrossRef] [PubMed]
  26. Smyth, K.; Drake, J.; Li, Y.; Rochman, C.; Van Seters, T.; Passeport, E. Bioretention cells remove microplastics from urban stormwater. Water Res. 2021, 191, 116785. [Google Scholar] [CrossRef] [PubMed]
  27. Peng, Y.; Wang, Y.; Chen, H.; Wang, L.; Luo, B.; Tong, H.; Zou, Y.; Lei, Z.; Chen, S. Carbon reduction potential of a rain garden: A cradle-to-grave life cycle carbon footprint assessment. J. Clean. Prod. 2024, 434, 139806. [Google Scholar] [CrossRef]
  28. She, L.; Wei, M.; You, X.-y. Multi-objective layout optimization for sponge city by annealing algorithm and its environmental benefits analysis. Sustain. Cities Soc. 2021, 66, 102706. [Google Scholar] [CrossRef]
  29. Cai, Y.; Zhao, Y.; Wei, T.; Fu, W.; Tang, C.; Yuan, Y.; Yin, Q.; Wang, C. Utilization of constructed wetland technology in China’s sponge city scheme under carbon neutral vision. J. Water Process Eng. 2023, 53, 103828. [Google Scholar] [CrossRef]
  30. Getter, K.L.; Rowe, D.B.; Robertson, G.P.; Cregg, B.M.; Andresen, J.A. Carbon sequestration potential of extensive green roofs. Environ. Sci. Technol. 2009, 43, 7564–7570. [Google Scholar] [CrossRef]
  31. Kavehei, E.; Jenkins, G.; Adame, M.; Lemckert, C. Carbon sequestration potential for mitigating the carbon footprint of green stormwater infrastructure. Renew. Sustain. Energy Rev. 2018, 94, 1179–1191. [Google Scholar] [CrossRef]
  32. Lin, X.; Ren, J.; Xu, J.; Zheng, T.; Cheng, W.; Qiao, J.; Huang, J.; Li, G. Prediction of life cycle carbon emissions of sponge city projects: A case study in Shanghai, China. Sustainability 2018, 10, 3978. [Google Scholar] [CrossRef]
  33. Su, X.; Shao, W.; Liu, J.; Jiang, Y.; Wang, J.; Yang, Z.; Wang, N. How does sponge city construction affect carbon emission from integrated urban drainage system? J. Clean. Prod. 2022, 363, 132595. [Google Scholar] [CrossRef]
  34. Moore, T.L.C.; Hunt, W.F. Predicting the carbon footprint of urban stormwater infrastructure. Ecol. Eng. J. Ecotechnol. 2013, 58, 44–51. [Google Scholar] [CrossRef]
  35. Peng, J.; Cao, Y.; Rippy, M.A.; Afrooz, A.R.M.N.; Grant, S.B. Indicator and Pathogen Removal by Low Impact Development Best Management Practices. Water 2016, 8, 600. [Google Scholar] [CrossRef]
  36. Sidek, L.M.; Jaafar, A.S.; Majid, W.H.A.W.A.; Basri, H.; Marufuzzaman, M.; Fared, M.M.; Moon, W.C. High-resolution hydrological-hydraulic modeling of urban floods using InfoWorks ICM. Sustainability 2021, 13, 10259. [Google Scholar] [CrossRef]
  37. Liu, X. Parameter calibration method for urban rainfall-runoff model based on runoff coefficient. Water Wastewater Eng. 2009, 45, 213–217. [Google Scholar] [CrossRef]
  38. Ministry of Housing Urban-Rural Construction of the People’s Republic of China. Technical Guide for Sponge City Construction-Construction of Rain Water System for Low Impact Development; MOHURD: Beijing, China, 2014; p. 88.
  39. Zhang, J. InfoWorks ICM Stormwater Models in the Application and Practice of Sponge City Construction in Mountainous Cities. Master’s Thesis, Chongqing University, Chongqing, China, 2018. [Google Scholar]
  40. T/CUWA 40052-2022; Technical Regulations for Stormwater Bioretention. China Urban Water Association: Beijing, China, 2022; p. 24.
  41. Cunningham, A.; Colibaba, A.; Hellberg, B.; Roberts, G.S.; Simcock, R.; Speed, S.R.; Vigar, N.; Woortman, W. Stormwater Management Devices in the Auckland Region; Auckland Council: Auckland, New Zealand, 2017. [Google Scholar]
  42. Chui, T.F.M.; Liu, X.; Zhan, W. Assessing cost-effectiveness of specific LID practice designs in response to large storm events. J. Hydrol. 2016, 533, 353–364. [Google Scholar] [CrossRef]
  43. Cao, S.; Dong, C. Life cycle cost benefit assessments for gerrn buildings. J. Tsinghua Univ. 2012, 52, 843–847. [Google Scholar] [CrossRef]
  44. Eckart, K.; McPhee, Z.; Bolisetti, T. Multiobjective optimization of low impact development stormwater controls. J. Hydrol. 2018, 562, 564–576. [Google Scholar] [CrossRef]
  45. Zhou, G.; Mei, C.; Liu, J.; Wang, H.; Shao, W.; Li, Z.; Wang, D.; Fu, I. Analysis on hydrological response and cost-benefit of sponge city construction in Ximen District of Pingxiang City. Water Resour. Hydropower Eng. 2019, 50, 10–17. [Google Scholar] [CrossRef]
  46. Xiao, X. Study on Life Cycle Carbon Emission and Life Cycle Cost of Green Buildings. Master’s Thesis, Beijing Jiaotong University, Beijing, China, 2021. [Google Scholar]
  47. Chongqing Municipal Design and Research Institute Co., Ltd. Indicators for Estimating Investment in Sponge City Construction Projects: ZYA1-02(01)-2018; China Planning Press: Beijing, China, 2018. [Google Scholar]
  48. Li, J.; Zhang, X.; Huimin, L. Study on carbon emission accounting in construction and operation of a sponge city in Beijing. Water Resour. Prot. 2023, 39, 86–93. [Google Scholar]
  49. Ma, J. Carbon Source Analysis and Carbon Emission Study of Typical Measures for Sponge City Construction. Master’s Thesis, Shanxi Agricultural University, Taiyuan, China, 2018. [Google Scholar]
  50. Li, C.; Zheng, T.; Peng, K.; Cheng, W.; Xu, J.; Qiao, J.; Huang, J. Study on carbon emission of sponge city stormwater system based on life cycle assessment. Environ. Sustain. Dev. 2019, 44, 132–137. [Google Scholar] [CrossRef]
  51. Banting, D.; Doshi, H.; Li, J.; Missios, P.; Au, A.; Currie, B.A.; Verrati, M. Report on the Environmental Benefits and Costs of Green Roof Technology for the City of Toronto; City of Toronto and Ontario Centres of Excellence—Earth and Environmental Technologies; Ryerson University: Toronto, ON, Canada, 2005. [Google Scholar]
  52. Zheng, T. Estimation of carbon emission during sponge city recoSnstruction of residential community. China Water Wastewater 2021, 37, 112–119. [Google Scholar] [CrossRef]
  53. China Urban Water Association. Guidelines for Carbon Accounting and EmissionReduction in the Urban Water Sector; China Architecture & Building Press: Beijing, China, 2022; p. 202. [Google Scholar]
  54. Tu, A.; Li, Y.; Mo, M.; Nie, X. Hydrological effects of design parameters optimization of bioretention facility based on RECARGA model. J. Soil Water Conserv. 2020, 34, 149–153. [Google Scholar] [CrossRef]
  55. Pan, J.; Ni, R.; Zheng, L. Influence of In-situ Soil and Groundwater Level on Hydrological Effect of Bioretention. Pol. J. Environ. Stud. 2022, 31, 3745–3753. [Google Scholar] [CrossRef]
  56. Haaland; Christinevan Den Bosch; Konijnendijk, C. Challenges and strategies for urban green-space planning in cities undergoing densification: A review. Urban For. Urban Green. 2015, 14, 760–771. [Google Scholar] [CrossRef]
  57. Liu, J.; Wang, J.; Ding, X.; Shao, W.; Mei, C.; Li, Z.; Wang, K. Assessing the mitigation of greenhouse gas emissions from a green infrastructure-based urban drainage system. Appl. Energy 2020, 278, 115686. [Google Scholar] [CrossRef]
  58. Moore, T.L.; Hunt, W.F. Ecosystem service provision by stormwater wetlands and ponds–a means for evaluation? Water Res. 2012, 46, 6811–6823. [Google Scholar] [CrossRef]
  59. Li, G.; Xiong, J.; Zhu, J.; Liu, Y.; Dzakpasu, M. Design influence and evaluation model of bioretention in rainwater treatment: A review. Sci. Total Environ. 2021, 787, 147592. [Google Scholar] [CrossRef]
  60. Wang, M.; Zhang, D.; Cheng, Y.; Tan, S.K. Assessing performance of porous pavements and bioretention cells for stormwater management in response to probable climatic changes. J. Environ. Manag. 2019, 243, 157–167. [Google Scholar] [CrossRef] [PubMed]
  61. Kaykhosravi, S.; Khan, U.T.; Jadidi, M.A. A simplified geospatial model to rank LID solutions for urban runoff management. Sci. Total Environ. 2022, 831, 154937. [Google Scholar] [CrossRef] [PubMed]
  62. Han, R.; Li, J.; Li, Y.; Xia, J.; Gao, X. Comprehensive benefits of different application scales of sponge facilities in urban built areas of northwest China. Ecohydrol. Hydrobiol. 2021, 21, 516–528. [Google Scholar] [CrossRef]
Figure 1. Technology roadmap.
Figure 1. Technology roadmap.
Water 16 00183 g001
Figure 2. Schematic map of the study area.
Figure 2. Schematic map of the study area.
Water 16 00183 g002
Figure 3. 120 min short-duration rainfall intensity.
Figure 3. 120 min short-duration rainfall intensity.
Water 16 00183 g003
Figure 4. 1440 min long-duration rainfall intensity.
Figure 4. 1440 min long-duration rainfall intensity.
Water 16 00183 g004
Figure 5. Full life cycle costing accounting boundary and activities.
Figure 5. Full life cycle costing accounting boundary and activities.
Water 16 00183 g005
Figure 6. Full life cycle carbon emission accounting boundary and activities.
Figure 6. Full life cycle carbon emission accounting boundary and activities.
Water 16 00183 g006
Figure 7. Schematic diagram of the composition of carbon emissions of BF over the full life cycle.
Figure 7. Schematic diagram of the composition of carbon emissions of BF over the full life cycle.
Water 16 00183 g007
Figure 8. Impact of (a) climatic conditions, (b) aquifer height, (c) permeability coefficient, and (d) the facility area on carbon emissions and carbon intensity.
Figure 8. Impact of (a) climatic conditions, (b) aquifer height, (c) permeability coefficient, and (d) the facility area on carbon emissions and carbon intensity.
Water 16 00183 g008aWater 16 00183 g008b
Figure 9. Effects of (a) climatic conditions, (b) aquifer height, (c) permeability coefficient, and (d) the facility area on the carbon reduction benefit.
Figure 9. Effects of (a) climatic conditions, (b) aquifer height, (c) permeability coefficient, and (d) the facility area on the carbon reduction benefit.
Water 16 00183 g009
Figure 10. Results of scenarios 19–42 on carbon intensity of volume capture of rainfall.
Figure 10. Results of scenarios 19–42 on carbon intensity of volume capture of rainfall.
Water 16 00183 g010
Figure 11. Results of scenarios 19–42 for carbon reduction benefit.
Figure 11. Results of scenarios 19–42 for carbon reduction benefit.
Water 16 00183 g011
Figure 12. Carbon intensity under different levels of VCRA.
Figure 12. Carbon intensity under different levels of VCRA.
Water 16 00183 g012
Figure 13. FLCC and (a) water control capacity, (b) carbon emission, and (c) carbon intensity.
Figure 13. FLCC and (a) water control capacity, (b) carbon emission, and (c) carbon intensity.
Water 16 00183 g013
Figure 14. The relationship between carbon emission intensity and carbon emission reduction benefits.
Figure 14. The relationship between carbon emission intensity and carbon emission reduction benefits.
Water 16 00183 g014
Table 1. Values of influencing factors.
Table 1. Values of influencing factors.
Influencing FactorValue
Climatic condition 1 (CC)C−20C−10CC+10C+20
Aquifer height (AH) (mm)100150200250300
Permeability coefficient (PC) (mm/h)102550100200
The facility area (FA) (m2)45.572.891136.5182273
Note: 1 C−20, C−10, C, C+10, and C+20 represent rainfall amounts corresponding to 80%, 90%, 100%, 110%, and 120% of the rainfall calculated using the local rainfall intensity formula, respectively.
Table 2. Simulation scenarios.
Table 2. Simulation scenarios.
Scenario NumberFactorsScenario NumberFactors
CCAH (mm)PC (mm/h)FA
(m2)
CCAH (mm)PC (mm/h)FA
(m2)
10.82005091220.830010091
20.92005091230.9100200136.5
31.12005091240.91502572.8
41.22005091250.9200100182
511005091260.92501091
611505091270.93005045.5
71250509128110010072.8
81300509129115010182
9120010913012005091
101200259131125020045.5
1112001009132130025136.5
12120020091331.110050182
1312005045.5341.115020091
1412005072.8351.12002545.5
15120050136.5361.1250100136.5
16120050182371.13001072.8
17120050273381.21002591
180.81001045.5391.215010045.5
190.815050136.5401.220010136.5
200.820020072.8411.22505072.8
210.825025182421.2300200182
Table 3. Constant values.
Table 3. Constant values.
Constant E F C S
(kg CO2/m2)
E F D S
(kg CO2/m2)
E F O
(kg CO2/m3)
E F M
(kg CO2/m2)
E F H
(kg CO2/m2)
E F H L
(kg CO2/m3)
Value 144.25235.67103.78905.130066.900044.6160
Note: 1 The constant values concerning the relevant carbon emission factors are taken from studies such as Li [48], Ma [49], and Li [50] by accounting.
Table 4. Orthogonal experiment results for carbon intensity.
Table 4. Orthogonal experiment results for carbon intensity.
CCAHPCFA
K1−0.1513−0.1891−0.1813−0.0379
K2−0.1637−0.1488−0.1565−0.0730
K3−0.1603−0.1885−0.1476−0.1101
K4−0.1489−0.1240−0.1611−0.2351
K5−0.1639−0.1377−0.1415−0.3320
R0.01500.00900.01500.0520
Table 5. Orthogonal experiment results for carbon reduction benefits.
Table 5. Orthogonal experiment results for carbon reduction benefits.
CCAHPCFA
K18.33367.50717.961612.8031
K29.70389.25468.798610.5100
K39.05678.25759.04179.0959
K48.49909.63968.91386.3631
K58.34369.27799.22105.1646
R0.46900.27000.41601.5020
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, D.; Liu, X.; Li, H.; Chen, H.; Wang, X.; Li, W.; Cao, L.; Liu, J.; Zhang, T.; Wei, B. The Carbon Emission Intensity of Rainwater Bioretention Facilities. Water 2024, 16, 183. https://doi.org/10.3390/w16010183

AMA Style

Wang D, Liu X, Li H, Chen H, Wang X, Li W, Cao L, Liu J, Zhang T, Wei B. The Carbon Emission Intensity of Rainwater Bioretention Facilities. Water. 2024; 16(1):183. https://doi.org/10.3390/w16010183

Chicago/Turabian Style

Wang, Deqi, Xuefeng Liu, Huan Li, Hai Chen, Xiaojuan Wang, Wei Li, Lianbao Cao, Jianlin Liu, Tingting Zhang, and Bigui Wei. 2024. "The Carbon Emission Intensity of Rainwater Bioretention Facilities" Water 16, no. 1: 183. https://doi.org/10.3390/w16010183

APA Style

Wang, D., Liu, X., Li, H., Chen, H., Wang, X., Li, W., Cao, L., Liu, J., Zhang, T., & Wei, B. (2024). The Carbon Emission Intensity of Rainwater Bioretention Facilities. Water, 16(1), 183. https://doi.org/10.3390/w16010183

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop