Analysis of Thermal Pollution Reduction Efﬁciency of Bioretention in Stormwater Runoff under Different Rainfall Conditions

: The thermal pollution of stormwater runoff is an important factor in the degradation of the urban water environment. Bioretention is an effective way to control the thermal pollution of stormwater runoff. To understand the inﬂuence of different rainfall conditions on the reduction of thermal pollution load of rainwater runoff from bioretention facilities and the correlation between the parameters, an experimental study was carried out by using the Beijing P&C rainfall pattern to change the rainfall parameters. By proposing innovative evaluation index parameters and analysis methods, the correlation between different parameters was quantiﬁed. The results showed that the heat pollution load reduction rate (HR) had a strong negative correlation with rainfall inﬂow volume (IV), rainfall duration (RD), and the service area ratio of bioretention facilities (CAR), and the correlation coefﬁcients were − 0.95, − 0.73, and − 0.73, respectively. HR was weakly correlated with rainfall return period (RP), rainfall temperature (RT), and air temperature (At) during the experiment, and the correlation coefﬁcients were − 0.29, 0.20, and 0.20, respectively. The delay between the peak value of heat output from the bioretention and the peak value of heat input from rainfall was about 10 min. Through the research on the change of rainwater heat entering and exiting the facility, we can master some rules of heat reduction in bioretention facilities, which will provide a reference for subsequent research on the difference between the internal heat change and bioretention temperature change.


Introduction
Urbanization construction leads to the increase of urban impervious underlying surface, which changes the thermal performance of the urban surface. The surface temperature of impervious concrete and asphalt is 19~23 • C higher than the air temperature on the premise of sufficient sunlight in summer [1]. According to the data from the National Climatic Data Center of the United States, the average outlet temperature of rainwater pipelines is 8.3 • C higher than that of cold water rivers [2,3]. The rainwater runoff formed by asphalt pavement of urban roads increases the temperature of the receiving water bodies by 5 • C and even 10~12 • C under extreme conditions [4]. The huge heat source formed by these high-temperature objects, under the scouring and convergence of rainfall-runoff, directly causes an instantaneous rise of river temperature [5], destroys the urban water environment, and produces the thermal pollution of rainwater runoff. After the summer rain, the water temperature of rivers and lakes in urban areas often rises sharply, the dissolved oxygen drops, the water quality deteriorates, and a large number of fish die [6]. Under the influence of the thermal pollution of rainwater runoff, the average temperature exchange between media of bioretention facilities in order to understand the difference between heat reduction and temperature reduction and master the variation law of peak value of facility heat discharge curve. It provides theoretical support for the application of bioretention facilities to the removal of thermal pollution from rainwater runoff, and also provides a reference for the subsequent research on the heat change and temperature change between the internal media of bioretention facilities.

Materials and Methods
The artificial rainfall system based on the Beijing P&C rain pattern method was used for rainfall simulation. By changing the rainfall temperature, rainfall duration, rainfall return period, and the service area of bioretention facilities, the thermal pollution control efficiency of facilities under different rainfall conditions was obtained.

Device Design
As shown in Figure 1, build 1 m × 1 m × 1.6 m (long × wide × High) stainless steel bioretention test box, the internal cylinder was a 500 mm diameter cylinder, and the filler was quartz sand with a particle size of 2~4 mm. The exterior was wrapped with aluminum foil high-density rubber plastic heat-absorbing cotton, and the insulation layer was made between the box and the internal soil column and filled with insulation asbestos. After the verification of the experimental device in the early stage, the duration of each group of experiments was 2~3 h, the internal temperature loss degree of the device was 0.11 °C/h, and the temperature loss rate accounted for about 1.04% of the inlet water temperature during the experiment. From top-to-bottom inside the soil column were a 150 mm water-storage layer, 250 mm planting-soil layer, 1000 mm internal-filler layer, and 200 mm bottom-gravel layer [18,19]. The soil particle size of the planting soil layer was 0.5~1 mm. Ophiopogon japonicas, which is common in northern China, was used as a plant. The particle size of the gravel layer at the bottom was 8~15 mm. DN32 seepage discharge pipes were respectively arranged in the aquifer and the gravel layer at the bottom, and thermometers were arranged at the bottom seepage discharge pipe orifice to record the outflow water temperature. From top-to-bottom inside the soil column were a 150 mm water-storage layer, 250 mm planting-soil layer, 1000 mm internal-filler layer, and 200 mm bottom-gravel layer [18,19]. The soil particle size of the planting soil layer was 0.5~1 mm. Ophiopogon japonicas, which is common in northern China, was used as a plant. The particle size of the gravel layer at the bottom was 8~15 mm. DN32 seepage discharge pipes were respectively arranged in the aquifer and the gravel layer at the bottom, and thermometers were arranged at the bottom seepage discharge pipe orifice to record the outflow water temperature. An artificial rainfall system was set above the bioretention column, including four square sprinkler heads with a side length of 150 mm and their supporting pipes, flow meters, centrifugal pumps, and water tanks. A thermostatic heater was used inside the water tank to control the artificial rainfall temperature. The rainfall interval between each experiment was 20~22 h, and 10 groups of experiments were conducted. The inflow flow was recorded by the flowmeter, the rainfall time interval was recorded by the computer, and the outflow flow of the bioretention facility was measured manually by the volumetric method.

Conceptual Design
On the premise that the rainwater runoff coefficient Ψ is certain, the rainfall volume was calculated by combining the rainfall intensity Formulas (1) and (2) of zone II in Beijing: where q is the design storm intensity [L/(s·hm 2 )]; t is the rainfall duration (min); V is the total design rainfall volume (L); Ψ is the runoff coefficient; F is the catchment area (hm 2 ); and P is the design return period (a). Since the shortest time for the Code for Design of Stormwater Management and Harvest Engineering to allocate rainfall patterns is 180 min, to obtain the distribution ratio of rainfall patterns under 60 min and 120 min rainfall duration, the Beijing P&C rain pattern method was selected for experimental research [20]. During the experiment, the artificial rainfall system took 5 min as a rainfall cycle, and the rainfall was calculated according to the design rainfall parameters in Formulas (1) and (2). As the biological detention facility is usually used as the source treatment facility for small and medium-sized rainfall, the rainfall return period (RP) was 1 year/2 years/5 years, the rainfall duration (RD) was 60 min and 120 min, the catchment area ratio (CAR) was 1:5 and 1:10, and the total influent volume (IV) during the experiment was recorded. As the temperature of summer rainwater runoff in Beijing is 22.8~26.6 • C [16], the experimental rainfall temperatures (RT) were selected as 22 • C, 25 • C, and 28 • C, respectively, the experiment was carried out from low to high, and the air temperature (AT) during the experiment was recorded. After the rainfall stopped, we continued to observe and record the outflow temperature and volume within one hour.
Due to a large number of experimental variables, the orthogonal experiment method was used to select representative combinations for the experiment. We used SPSS software to carry out orthogonal experimental design, added the four abovementioned parameter variables (RT/RP/CAR/RD), and set their values to obtain the orthogonal scheme design table. See Table 1 for the specific experiment contents. The reduction of runoff thermal pollution by green rainwater infrastructure is mainly reflected in two aspects: one is the reduction of runoff discharge volume, and the other is the reduction of outlet water temperature of the facility through heat exchange between media. The reduction rate of thermal pollution of each rainfall of the facility can be calculated according to Formulas (3) and (4): where EMT is the average temperature of field rainfall-runoff ( • C); Q t is the runoff generated by this rainfall (m 3 /s); T t is the runoff temperature generated by this rainfall ( • C); δt is the time increment (min); and T d is the total duration of runoff (min).
where r is the reduction rate of heat pollution load of rainwater runoff of a single facility (%); EMT out is the average temperature of the discharged runoff after passing through the bioretention facilities ( • C); EMT in is the average temperature of surface runoff on the underlying surface of the catchment served by the bioretention facilities ( • C); V out is the volume of rainwater runoff discharged by the bioretention facilities of this rainfall (m 3 ); V in is the volume of rainwater runoff from this rainfall into the bioretention facilities (m 3 ); C is the specific heat capacity of water (4.2 kJ/(kg· • C)); and ρ is the density of water (1 × 10 3 kg/m 3 ).

Volume Control Ratio of Rainwater Runoff Heat Load
Bioretention facilities are widely used in various source volume-control measures. For the current status of heat pollution, the concept of the ratio of heat load-control rate to volume-control rate (HV control ratio) is proposed by referring to the cumulative load cumulative volume fraction curve (M (V) curve) to intuitively compare the heat pollution volume control effects of different bioretention facilities. The HV is the ratio of the heat load reduction rate η and the runoff volume reduction rate µ of this rainfall. This ratio is a dimensionless number, and its meaning is such that whenever 1% of the runoff volume can be controlled to control 1% of the heat load, the HV control ratio is 1.
Only one HV control ratio can be obtained from a single rainfall. Multiple HV control ratios of the bioretention facility can be obtained by changing the rainfall conditions to simulate multiple rainfalls, forming an HV control ratio curve. By comparing the HV control ratio data of different bioretention facilities under the same rainfall conditions, the optimal bioretention facilities for controlling the unit volume of thermal pollution of rainwater runoff under certain rainfall conditions can be obtained. Multiple groups of HV control ratio data of the same bioretention facility in a certain range of rainfall conditions are calculated, which can be used to analyze the impact of different rainfall parameters on the HV control ratio. Therefore, the HV of bioretention facilities can be expressed by Formula (5): where η is the ratio of the heat load reduction rate of this rainfall (%); µ is the runoff volume reduction rate of this rainfall (%); C is the specific heat capacity of water (4.2 × 10 3 J/(kg· • C)); ρ is the influent density (1 × 10 3 kg/m 3 ); V in is the total volume of influent water (m 3 ); V out is the total volume of effluent (m 3 ); V i in is the inflow volume per minute (m 3 ); V i out is the effluent volume per minute (m 3 );T i in is the inlet water temperature per minute ( • C); and T i out is the outlet water temperature per minute ( • C).

Temperature Control Ratio of Rainwater Runoff Heat Load
Similar to the concept of the HV control ratio, the ratio of heat load control ratio-totemperature control ratio (HT) is proposed for the relationship between runoff temperature controlled by bioretention facilities and heat load-reduction rate. HT is the ratio of the heat load-reduction rate η and the runoff volume reduction rate ϕ of this rainfall. This ratio is a dimensionless number. Its meaning is such that whenever 1% runoff temperature can be controlled to control 1% heat load, the control ratio of HT is 1.
The HT of bioretention facilities can be expressed by Formula (6): where η is the ratio of the heat load reduction rate of this rainfall (%); ϕ is the runoff volume reduction rate of this rainfall; C is the specific heat capacity of water (4.2 × 10 3 J/(kg· • C)); ρ is the influent density (1 × 10 3 kg/m 3 ); V i in is the inflow volume per minute (m 3 ); V i out is the effluent volume per minute (m 3 ); T i in is the inlet water temperature per minute ( • C); T i out is the outlet water temperature per minute ( • C); EMT in is the average temperature of surface runoff on the underlying surface of the catchment served by the bioretention facilities ( • C); and EMT out is the average temperature of the discharged runoff after passing through the bioretention facilities ( • C).

Analytical Method
Because the rainfall parameters and their experimental results are relatively independent, Spearman's Rank Correlation Coefficient Method was selected as the correlation analysis method among the parameters in this study.
Spearman's Rank Correlation Coefficient Method is used to assign a value between −1 and 1 by measuring the degree of correlation between the two waveforms. When the assignment is positive, there is a positive correlation, and when the assignment is negative, there is a negative correlation, and the higher the absolute value assigned, the stronger the correlation.
The specific method of Spearman's Rank Correlation Coefficient is to arrange the waveform data X in ascending or descending order to obtain the sequence A and record the position of each element in the sequence X in A as r i to obtain the rank-order column r of X. Another set of waveform data Y is arranged in the same way to obtain the rankorder column s of Y. The rank difference sequence d is obtained by successively making a difference between the elements in the sequence r and s, and the elements in the sequence d are successively brought into Spearman's Rank Correlation Coefficient Formula (7) [21]: where n is the sample size of waveform data; σ Is Spearman's rank correlation coefficient; and i is the ith sample.
To better distinguish the correlation, this paper classified the data results of Spearman's Rank Correlation Coefficient [22], as shown in Table 2.

Influence of Rainfall Pattern Parameters on the Reduction of Heat Pollution Load of Rainwater Runoff
By changing the rain pattern parameters, the experimental study on heat-pollution reduction of rainwater runoff from bioretention facilities was carried out. During the experiment, the air temperature was 16~27 • C, and the heat-reduction rate (HR), volumereduction rate (VR), the contribution of volume reduction (CVR), and contribution of heat exchange (CHE) were obtained. See Table 3 below. It can be seen from Table 3 that, based on the abovementioned rain pattern parameters, the VR of the bioretention facilities was 13~61%, which is close to the 27~76% rainfall effluent results obtained by Jones M [23] in the study on the thermal pollution of the bioretention facilities in the riparian biological reserve, so the experimental results are reliable. The reduction range of thermal pollution load of rainwater runoff was 21~63%. In addition to volume reduction, there was also heat exchange between the heat in runoff rainwater and the internal medium of bioretention facilities, so the overall reduction of thermal pollution load was slightly higher than the volume reduction rate. We created an HR and VR relation diagram, and obtained the results presented in Figure 2. With the gradual increase of VR, the rising speed of HR began to accelerate gradually. When VR was greater than 30%, the rising speed slowed down gradually. According to the data in Table 2, under the same rainfall temperature, with the increase of total rainfall volume, HR/VR also gradually increased, indicating that with the increase of rainfall volume, the contribution rate of volume reduction in heat exchange gradually decreases. The rainfall volume is related to the rainfall return period, rainfall duration, and the catchment area of facility service. Draw Figure 3 according to the sequence of temperature from low to high and volume from small to large. It can be found that with the increase  The rainfall volume is related to the rainfall return period, rainfall duration, and the catchment area of facility service. Draw Figure 3 according to the sequence of temperature from low to high and volume from small to large. It can be found that with the increase of the total volume of rainfall, the reduction rate of heat pollution decreased on the whole. Under the condition of similar rainfall volume, combined with the parameters of rainfall temperature in Table 3, it can be known that the higher the rainfall temperature, the higher the HR of facilities. The preliminary analysis of the cause of this phenomenon is that under the premise that the internal temperature of the bioretention facility is close, the RT increases, increasing the temperature difference between the rainwater and the internal medium of the facility, and enhancing heat exchange. The rainfall volume is related to the rainfall return period, rainfall dur catchment area of facility service. Draw Figure 3 according to the sequence o from low to high and volume from small to large. It can be found that with of the total volume of rainfall, the reduction rate of heat pollution decreased Under the condition of similar rainfall volume, combined with the paramet temperature in Table 3, it can be known that the higher the rainfall temperatu the HR of facilities. The preliminary analysis of the cause of this phenomenon the premise that the internal temperature of the bioretention facility is clo creases, increasing the temperature difference between the rainwater and th dium of the facility, and enhancing heat exchange. The three figures, representing the reduction rate of heat pollution, th duction rate, and the volume reduction percentage, respectively, trend mo gether in Figure 3. On this basis, taking the volume reduction rate of VR as linear fitting relationship between VR and CVR can be obtained, as show below. The three figures, representing the reduction rate of heat pollution, the volume reduction rate, and the volume reduction percentage, respectively, trend more closely together in Figure 3. On this basis, taking the volume reduction rate of VR as the x-axis, the linear fitting relationship between VR and CVR can be obtained, as shown in Figure 4 below.  CVR and VR showed a linear correlation with a slope of 0.68. Therefore, CVR and VR showed a linear correlation with a slope of 0.68. Therefore, with the gradual increase of VR, the growth rate of CVR slowed down, slightly lower than that of HR. As the reduction of thermal pollution is mainly caused by the volume reduction and the heat exchange inside the facility, with the increase of VR, the proportion of heat exchange (CHE) between the medium inside the bioretention facility and rainwater gradually decreased.

022, 14, x FOR PEER REVIEW
To further analyze the correlation between the parameters, based on the experimental results in Table 3 and the design parameters in Table 1, the Spearman's Rank Correlation Coefficient Method was used to obtain Figure 5.  CVR and VR showed a linear correlation with a slope of 0.68. Therefore, with the gradual increase of VR, the growth rate of CVR slowed down, slightly lower than that of HR. As the reduction of thermal pollution is mainly caused by the volume reduction and the heat exchange inside the facility, with the increase of VR, the proportion of heat exchange (CHE) between the medium inside the bioretention facility and rainwater gradually decreased.
To further analyze the correlation between the parameters, based on the experimental results in Table 3 and the design parameters in Table 1, the Spearman's Rank Correlation Coefficient Method was used to obtain Figure 5. It can be seen from Figure 5 that HR was affected by IV, RD, and CAR, showing a strong negative correlation. Among them, the correlation coefficient between IV and HR was as high as −0.95, and RD and CAR had the same impact on HR, both of which were It can be seen from Figure 5 that HR was affected by IV, RD, and CAR, showing a strong negative correlation. Among them, the correlation coefficient between IV and HR was as high as −0.95, and RD and CAR had the same impact on HR, both of which were −0.73. This result is consistent with the existing research contents [24]. RP had a weak negative correlation with HR, which was −0.29. RT and AT had a weak positive correlation with HR, with a coefficient of 0.20. It can be found that during the experiment, the shorter the RD and the smaller the CAR, the more obvious the thermal pollution reduction effect of the bioretention facilities. However, the smaller the RP and the higher the RT were, the higher the thermal pollution reduction rate. Additionally, due to their weak positive correlation, the HR performance was still slightly improved. In the actual scene, the rainfall temperature is affected by the air temperature, and the positive correlation coefficient between AT and VR is 0.46. Therefore, in the summer when the AT is high, the thermal pollution of surface rainwater runoff will be better removed by the bioretention facilities due to the increase of VR. According to the analysis, the cause of this phenomenon is that as the air temperature rises, the surface evaporation rate increases and the water content inside the bioretention facilities decrease before rainfall. Therefore, more rainwater can be absorbed during rainfall, and the volume reduction rate is increased.
There are two ways to reduce the thermal pollution of bioretention facilities: (1) reducing the volume of rainfall; and (2) heat exchange between rainwater and internal media of bioretention facilities. It is obvious from Figure 5 that VR and HR have a strong positive correlation, and the correlation coefficient was as high as 0.95. According to the data in Table 3, the smaller the IV, the higher the VR, the higher the corresponding HR, and the larger the proportion of volume reduction in the removal of heat pollution load; when HR is low, IV is large, the heat exchange between rainwater and internal medium of bioretention facilities is more obvious, and the proportion of CHE gradually increases. Since RT also had a weak positive correlation with CHE with a coefficient of 0.15, it can be interpreted that the contribution of CHE will increase slightly with the increase in rainfall temperature. Therefore, the previous analysis suggesting that higher RT leads to higher HR is tenable.

Influence of Rain Pattern Parameters on HV and HT
During the experiment, EMT of effluent from bioretention facilities in each group was obtained, and HV and HT were obtained based on HR, as shown in Table 4 below. Based on the data in Table 4, Spearman's Rank Correlation Coefficient Formula was applied to obtain the experimental results in Figure 6. The average EMT of the effluent from the bioretention facilities was closely related to RD and AT during the experiment, showing a strong positive correlation, and the Spearman's Coefficients were 0.73 and 0.60 respectively. Among them, the longer the RD, the higher the EMT of the outflow. During the experiment, the temperature was related to the experimental results. Through analysis, it can be understood that the higher the temperature, the higher the surface plant and soil temperature of the bioretention facili- The average EMT of the effluent from the bioretention facilities was closely related to RD and AT during the experiment, showing a strong positive correlation, and the Spearman's Coefficients were 0.73 and 0.60 respectively. Among them, the longer the RD, the higher the EMT of the outflow. During the experiment, the temperature was related to the experimental results. Through analysis, it can be understood that the higher the temperature, the higher the surface plant and soil temperature of the bioretention facilities, and the reduction effect of the bioretention on the rainwater temperature will be reduced. In the previous section, it was found that the higher the temperature is, the higher the VR of the bioretention facilities will be, and the better the effect of heat pollution reduction will be. Therefore, the reduction capacity of rainwater runoff heat pollution load and outlet temperature can not be generalized. HR is mainly affected by volume reduction. When the temperature is high, the contribution of volume reduction to the reduction of heat pollution load is enhanced, the heat exchange between the media is reduced, and the bottom outlet temperature of the bioretention facility will be higher than the bottom outlet temperature when the temperature is low, but the outlet volume will be reduced.
Based on the data in Table 4 and Figure 6, the following conclusions can be drawn based on removing the data with large deviation. Within the range of rainfall temperature of 22~28 • C and rainfall volume of 41.75~145.05 L, 1.235~1.696% of heat load could be controlled whenever 1% of runoff volume was controlled, and 1.517~2.778% of heat load could be controlled whenever 1% of runoff temperature was controlled. There was a strong negative correlation between HV and HT, and the coefficient was −0.87. With the increase of HV, the value of HT decreased, and both of them were affected by the catchment area. HV was also negatively correlated with AT; HT was greatly affected by RT, showing a moderate negative correlation and a correlation coefficient of −0.49. The outlet EMT of the facility also had a moderate positive impact on HT, with a coefficient of 0.39.

Analysis of the Heat-Emission Curve of Bioretention
Due to the different rainfall durations, the proportion of rainfall distribution every 5 min changed. Therefore, the occurrence time of rainfall peak in this experiment is only related to RD and has nothing to do with other rainfall parameters. On this basis, the experimental data were sorted out, and the peak time of rainfall heat (PTRH) and the peak time of outflow heat (PTOH) of 10 groups of experiments were obtained, as shown in Table 5 below. It can be seen from Table 5 that under the conditions of 60 min and 120 min rainfall duration, the peak value of the rainfall heat of the P&C rainfall pattern in Beijing occurred every at 25 min at the beginning of rainfall. It can be seen from the analysis that under the above rainfall duration, the peak rainfall happened to occur at the same time, which is when the rainfall started at 25 min, while the inflow temperature of the bioretention facility remained unchanged during the experiment, so PTRH also occurred when the experiment started for 25 min.
According to the preliminary analysis, most of the PTOH during the experiment occurred at 35 min after the rainfall. When RT and IV gradually increased, the total heat of rainfall increased, and the PTOH of bioretention facilities moved back, which occurred 40 min after the rainfall began. The occurrence time of PTOH was related to RT and IV and had no obvious relationship with rainfall duration.
Based on the abovementioned results, to understand the relationship between the heat of rainwater entering and flowing out of the facility during the process from rainfall to outflow of the bioretention facility, two typical rainfall experiments of HR were selected for research under the rainfall duration of 60 min and 120 min respectively. The selected rainfall pattern parameters are shown in Table 6 below. We calculated and counted the rainfall inflow heat and bottom outflow heat every 5 min during the experiment, and drew a cumulative curve of their percentage in the total rainfall heat and total emission heat to obtain the results shown in Figure 7. It could be found from the experimental results that under the condition of 60 mi rainfall duration, the peak value of rainfall heat was 25 min at the beginning of rainfal and the peak value of outflow heat was 35 min at the beginning of rainfall, with a dela It could be found from the experimental results that under the condition of 60 min rainfall duration, the peak value of rainfall heat was 25 min at the beginning of rainfall, and the peak value of outflow heat was 35 min at the beginning of rainfall, with a delay of about 10 min. When the rainfall ended, the accumulated outflow heat at the bottom outlet accounted for 91.28% of the total outflow heat of the experiment; under the condition of 120 min rainfall duration, the peak value of rainfall heat and outflow heat was the same as that at 60 min, and the peak delay was 10 min. When the rainfall ended, the accumulated outflow heat at the bottom accounted for 96.84% of the total outflow heat of the experiment.
The experiment was ended by taking the accumulated outflow heat accounting for more than 95% of the total outflow heat as the standard. The 60 min rainfall duration experiment should be continued to observe and record for at least 10 min after the rainfall stops. At this time, the accumulated outflow heat accounted for 95.79% of the total outflow heat. In the 120 min rainfall duration experiment, when the rainfall stopped, the accumulated outflow heat accounted for more than 95% of the total outflow heat, which meets the standard, and the experimental observation could be ended.

Summary and Conclusions
The bioretention facilities can reduce the heat-pollution load by 21~63% and the rainfall volume by 13~61% when the RP is within 5 years; the RD was 60 min, 120 min, and the CAR was 1:5 and 1:10. The reduction of the thermal pollution load of rainwater runoff by bioretention facilities was mainly in the form of volume reduction.
HR had a strong negative correlation with IV, RD, and CAR in bioretention facilities, with correlation coefficients of −0.95, −0.73, and −0.73, respectively, and a weak negative correlation with RP of −0.29. Both RT and AT had a weak positive correlation with HR with a correlation coefficient of 0.20. The increase of AT would increase VR, and the removal effect of bioretention facilities on thermal pollution would be better. The rise of RT would increase the temperature difference between rainwater and the medium inside the bioretention, and the rise of CHE would slightly increase HR.
The heat pollution load reduction capacity of rainwater runoff was not equivalent to the outlet temperature reduction capacity, and HR was mainly affected by volume reduction. With the increase of temperature, the contribution of volume reduction to HR could be enhanced, which would reduce the heat exchange between the media, resulting in the outlet temperature at the bottom of the bioretention facility being higher than that when the temperature is low under the same conditions, but the outlet volume would be reduced more, and the total outlet heat would still be reduced.
During the experiment, when the RD was 60 min and 120 min, the peak value of rainfall heat was 25 min from the beginning of rainfall, and the peak value of outflow heat was 35 ± 5 min from the beginning of rainfall. Under the condition that the cumulative rainfall was close, RD had little effect on the arrival time of the peak rainfall and the peak heat emission of bioretention facilities, and the peak delay of both was about 10 min. The end standard of the experiment was that the accumulated outflow heat at the bottom of the bioretention accounted for more than 95% of the total outflow heat. When RD was 60 min, we continued to observe and record for at least 10 min after the rainfall stopped. For the experiment with an RD of 120 min, the end standard of the experiment was met when the rainfall stopped.
It is believed that the development of the research on thermal pollution of rainwater runoff by bioretention facilities will be gradually deepened. At present, the removal effect of bioretention facilities on thermal pollution of rainwater runoff is recognized by many experts and scholars, but the research on temperature change between different levels within bioretention facilities is still blank. How to build a temperature change model within the facilities based on experimental research will be the next research direction of bioretention facilities in the study of thermal pollution of rainwater runoff.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.