Planning Nature-Based Solutions for Urban Flood Reduction and Thermal Comfort Enhancement

: As a consequence of climate change and urbanization, many cities will have to deal with more ﬂooding and extreme heat stress. This paper presents a framework to maximize the e ﬀ ectiveness of Nature-Based Solutions (NBS) for ﬂood risk reduction and thermal comfort enhancement. The framework involves an assessment of hazards with the use of models and ﬁeld measurements. It also detects suitable implementation sites for NBS and quantiﬁes their e ﬀ ectiveness for thermal comfort enhancement and ﬂood risk reduction. The framework was applied in a densely urbanized study area, for which di ﬀ erent small-scale urban NBS and their potential locations for implementation were assessed. The overall results show that the most e ﬀ ective performance in terms of ﬂood mitigation and thermal comfort enhancement is likely achieved by applying a range of di ﬀ erent measures at di ﬀ erent locations. Therefore, the work presented here shows the potential of the framework to achieve an e ﬀ ective combination of measures and their locations, which was demonstrated on the case of the Sukhumvit area in Bangkok (Thailand). This can be particularly suitable for assessing and planning ﬂood mitigation measures in combination with heat stress reduction.


Introduction
There is an increasing awareness that the interplay between the supposed effects of climate change and global warming combined with rapid and uncontrolled urbanization can lead to serious challenges to urban water managers and city planners. Since the vegetation coverage and green areas are decreasing significantly, the imperviousness rate in different urban areas is increasing [1,2]. As a The results obtained from hazard assessment are used in the selection of small-scale urban NBS ("Which") and the identification of suitable areas for their implementation ("Where"). Best management practices (BMP) Sitting and ArcGIS tools are used to identify the locations suitable for the implementation of the selected NBS measures. The last part ("How much") evaluates the effectiveness of the selected measures for flood reduction and thermal comfort enhancement. Hydrodynamic and microclimatic modeling are used to assess flood and heat stress mitigation effectiveness. The use of micro-climatic models is often used for urban assessments [26,27].
Six parameters are considered for the effectiveness assessment, three related to flood mitigation and three oriented to evaluate heat stress mitigation. Finally, according to the results obtained for these parameters, scores are given to the measures for their comparison. The following sections introduce the case study area framework.

Study Area
The framework was applied in the Sukhumvit area, located in central Bangkok (Figure 2), which is a highly dense urban area of approximately 23 Km 2 . According to The World Bank [28], urban growth in Thailand is mostly situated in the Bangkok urban area, which is among the twenty largest cities in Asia in terms of population, approaching 10 million people. Bangkok is a growing city located in a tropical area, and it is facing many extreme climatic conditions, which will also be intensified in the future as a result of climate change [29].
The annual rainfall in the city is 1651 mm, which mainly takes place in the wet season (from May to October). According to Rehan et al. [30] the number of annual rainy days was increased from 90 to 110 days in the last 30 years. Additionally, the study done by Sheikh [31] shows that the average mean temperature in Bangkok was increased by 0.6 • C between 1985 and 2014. Arifwidodo and Tanaka [32] studied the effect of the Urban Heat Island (UHI) in Bangkok, showing that there is a mean maximum of 5 • C UHI intensity (UHII) between semi urban and urban areas, and a mean maximum UHII of 2 • C between dry and in rainy seasons. In addition, there is a maximum night time UHII of 7 • C during January in Bangkok [29,32]. Flooding is caused by excessive local rainfalls or by the overtopping of embankments due to a high water level in the Chao Phraya River [4,33]. Even though there are numerous pumping stations inside the city to pump the excess of storm water to the river, the city is still highly vulnerable to flooding. The problem is aggravated by over extraction of ground water, which has caused land subsidence of up to 15 cm in many locations [34].

Hazard Assessment
Hazard assessment was undertaken to gain a better understanding of existing conditions, and to identify locations which are more hazardous in terms of urban flooding and heat stress. The choice to undertake a macro scale approach for urban flood modeling and a micro scale approach for the thermal comfort modeling was based on field observations and relevant literature review For assessing urban flood hazards, the existing sewer system was modeled with a hydrodynamic model, considering several scenarios with particular emphasis on areas with a higher frequency of flooding. In terms of hazards due to heat stress, this assessment aimed to identify and evaluate the effects from different urban land uses on heat stress and human thermal comfort. For this purpose, field data was collected from both fixed weather stations and mobile weather measurements (using the instrument Kestrel 5400 Heat Stress Tracker). Five different urban land uses were considered for the heat stress assessment with mobile weather measurements (Table 1), and categorized following the work of Stewart and Oke [35].  Flooding is also a severe problem in Bangkok, causing important economic and health-related problems. This problem has been aggravated in the last years due to urbanization, which has generated land use changes in the city. According to Srivanit et al. [29] the urban/built-up land in Bangkok and its metropolitan area increased by almost three times between 1994 and 2009, growing from 15% in 1994, to 42% in 2009. In contrast, a pronounced decrease in the vegetated area was observed from 1994 (72%) to 2009 (40%). The land use in the Study Area of Sukhumvit is presented in Appendix A.
Flooding is caused by excessive local rainfalls or by the overtopping of embankments due to a high water level in the Chao Phraya River [4,33]. Even though there are numerous pumping stations inside the city to pump the excess of storm water to the river, the city is still highly vulnerable to flooding. The problem is aggravated by over extraction of ground water, which has caused land subsidence of up to 15 cm in many locations [34].

Hazard Assessment
Hazard assessment was undertaken to gain a better understanding of existing conditions, and to identify locations which are more hazardous in terms of urban flooding and heat stress. The choice to undertake a macro scale approach for urban flood modeling and a micro scale approach for the thermal comfort modeling was based on field observations and relevant literature review.
For assessing urban flood hazards, the existing sewer system was modeled with a hydrodynamic model, considering several scenarios with particular emphasis on areas with a higher frequency of flooding. In terms of hazards due to heat stress, this assessment aimed to identify and evaluate the effects from different urban land uses on heat stress and human thermal comfort. For this purpose, field data was collected from both fixed weather stations and mobile weather measurements (using the instrument Kestrel 5400 Heat Stress Tracker). Five different urban land uses were considered for the heat stress assessment with mobile weather measurements (Table 1), and categorized following the work of Stewart and Oke [35]. The assessment of heat stress variation was based on measurements of different weather parameters taken at five different locations. These weather parameters were: Air temperature (Ta), mean radiant temperature (Tmrt), predicted mean vote (PMV) and physiological equivalent temperature (PET). The rationale behind this assessment was to analyze the effect of different urban land uses on heat stress and thermal comfort.
In terms of the selection of parameters, according to Coccolo et al. [36], "the Tmrt is considered to be as an artificial measure to express the degree of exposure to the environmental radiation. The radiant temperature is related to the amount of radiant heat transferred from a surface, and it depends on the material's ability to absorb or emit heat, or its emissivity". As Höppe [37] stated, "PET is defined as the air temperature at which, in a typical indoor setting (without wind and solar radiation), the heat budget of the human body is balanced with the same core and skin temperature as under the complex outdoor conditions to be assessed". Predicted mean vote (PMV) is one of the most used thermal indices by researchers. It is initially based on Fanger's heat balance model, and is an index ranging from −3 for cold weather and +3 for hot weather). This index is an outcome of A result of the perceived sensation of the thermal environment of a group of people, and was initially developed for indoor environments [38].

Selection of NBS Types and Their Suitable Sites: Which and Where
Four measures were selected based on several factors, including: Recommendations in the literature [39][40][41], feasibility of implementation in the case study area and the capability to be modeled in Mike Urban (for urban flood analysis) and the ENVI-met microclimatic simulation model (for thermal comfort analysis). The selected measures are: green roofs, pervious pavements, bio retentions and rain gardens. Details about each measure and the parameters used in the models are provided in Appendix B.
For suitability analysis, these measures were categorized in two groups: Green roof (GR) and pervious pavement (PP) in the first category, and bio-retention (BR) and rain garden (RG) in the second one. The analysis of possible locations and the maximum application of these measures was done using satellite images and the geographic information system (GIS) data. For the second category, the BMP sitting tool Sustain [42] was also used (see Appendix C).
The boundary of the case study area (macro scale) is depicted with a green line in Figure 3 (left). Two distinct locations inside the case study area of Sukhumvit were chosen for microclimatic simulations, and they are shown in Figure 3 (within the red rectangles). These locations were categorized as open low-rise (A) and compact high-rise (B), and are examples of two most common urban configurations in the case study area: A site with low-rise buildings and a site with dense high-rise buildings. This selection was made to achieve a comprehensive analysis of how the measures will be effective in each of these conditions. The two selected micro scale sites represent two different urban configurations for urban climate zones based on the research done by Stewart and Oke [35] (see Appendix D).

Model Development and Data Analysis
In order to quantify the effectiveness of NBS for urban flood reduction, a flood hazard map was produced by applying a 1D/2D modeling approach within Mike Urban software [43] (see Appendix C), see also Vojinovic and Tutulic [44]. The model was run for different rainfall return periods and different cases of NBS measures applications in the study area. The input data to create the model was collected from the Bangkok Metropolitan Administration office, Hydro and Agro Informatics Institute and the Department of Drainage and Sewerage office in Bangkok. The flooding reduction in the study area was assessed considering storm water runoff reduction, peak flow reduction and time to peak delay.
In terms of the heat stress, the ENVI-met v4.1.3 model [45] was used for assessing the effectiveness of measures in relation to thermal comfort. ENVI-met is a three-dimensional computational fluid dynamics non-hydrostatic S.V.A.T. (soil, vegetation, atmosphere, and transfer) model (Appendix C). This software is commonly used for modeling surface-plant-air interactions in urban environments, and it can also simulate flows around buildings, heat and vapor transfer at urban surfaces, turbulence and exchanges of energy and mass between the vegetation and its surroundings, and simple chemical reactions [46][47][48][49].
The input data and parameters used in each of the models are presented in Table 2. The changes made in each model in order to represent the application of NBS are presented as inputs, while the outputs explain the type of results obtained from each model. The results are compared with the case that no such measures are applied (also referred to as a 'baseline scenario' or 'business as usual') in order to evaluate the change of conditions obtained from each alternative.

Model Development and Data Analysis
In order to quantify the effectiveness of NBS for urban flood reduction, a flood hazard map was produced by applying a 1D/2D modeling approach within Mike Urban software [43] (see Appendix C), see also Vojinovic and Tutulic [44]. The model was run for different rainfall return periods and different cases of NBS measures applications in the study area. The input data to create the model was collected from the Bangkok Metropolitan Administration office, Hydro and Agro Informatics Institute and the Department of Drainage and Sewerage office in Bangkok. The flooding reduction in the study area was assessed considering storm water runoff reduction, peak flow reduction and time to peak delay.
In terms of the heat stress, the ENVI-met v4.1.3 model [45] was used for assessing the effectiveness of measures in relation to thermal comfort. ENVI-met is a three-dimensional computational fluid dynamics non-hydrostatic S.V.A.T. (soil, vegetation, atmosphere, and transfer) model (Appendix C). This software is commonly used for modeling surface-plant-air interactions in urban environments, and it can also simulate flows around buildings, heat and vapor transfer at urban surfaces, turbulence and exchanges of energy and mass between the vegetation and its surroundings, and simple chemical reactions [46][47][48][49].
The input data and parameters used in each of the models are presented in Table 2. The changes made in each model in order to represent the application of NBS are presented as inputs, while the outputs explain the type of results obtained from each model. The results are compared with the case that no such measures are applied (also referred to as a 'baseline scenario' or 'business as usual') in order to evaluate the change of conditions obtained from each alternative.

Surface and Soil parameters
The Green area percentage is increase by 5%. The inputs are Number of trees, LAD, RAD, plant height, Albedo and Leaf type.

Scenarios Development
Further to the above, ten scenarios were considered in each case, to evaluate urban flooding and thermal comfort (macro and micro scales respectively). Table 3 shows the scenarios for urban flood simulations, which include two different precipitation return periods for each of the four selected NBS, in addition to the baseline scenario. Table 3. Nature-Based Solutions (NBS) implementation scenarios for urban flood assessment.  Table 4 shows scenarios for microclimatic thermal comfort simulations, which include two different site characteristics and the four selected NBS, in addition to the baseline scenario. Figure 4 shows the overall framework for assessment of effectiveness from NBSs in relation to urban flood and thermal comfort. For each case, the models used, scenarios, variables and outputs, are presented. Table 4. NBS implementation scenarios for thermal comfort assessment.

Comparative Effectiveness of NBSs
The main purpose of the comparative performance scoring performed here is to identify the most effective measure, taking into consideration all parameters and scenarios. The effectiveness of measures for flooding reduction and thermal comfort enhancement are evaluated using the six parameters shown in Figure 4. According to the results obtained, the measures are scored from 1 to 4 according to the relation to their performance compared with other measures for the same scenario. This implies that the measure with the highest effectiveness will be scored as 4, and the one with the least effectiveness (comparably) will be scored as 1.

Comparative Effectiveness of NBSs
The main purpose of the comparative performance scoring performed here is to identify the most effective measure, taking into consideration all parameters and scenarios. The effectiveness of measures for flooding reduction and thermal comfort enhancement are evaluated using the six parameters shown in Figure 4. According to the results obtained, the measures are scored from 1 to 4 according to the relation to their performance compared with other measures for the same scenario. This implies that the measure with the highest effectiveness will be scored as 4, and the one with the least effectiveness (comparably) will be scored as 1.

Hazard Assessment
The initial heat stress assessment was performed from mobile weather measurements. These measurements were performed in five different urban land uses within the case study area (M sites in Figure 5). Values of air temperature, wind speed, glob temperature and humidity were collected; and some thermal comfort indices including PMV and PET were calculated. The Kestrel 5400 Heat Stress Tracker was used for this purpose. This procedure showed that urban parks (M4) and open low-rise sites (M5) can be respectively 1.4 • C and 1.1 • C cooler than a highly dense urban site in the study area (M2 and M3). Furthermore, the results showed that the PMV and PET can be also lower in these two sites compared to highly dense urban sites. Notice that M3 and M5 are the sites chosen to perform the thermal modeling, presented as areas B and A, respectively, in Figure 3.
Additionally, by overlaying the flood hazard map with the buildings data ( Figure 5), it was observed that less floods and more heat stress are likely to occur in the highly dense upstream area (A1 in Figure 5), when compared other parts of the study area. However, less urbanized areas (A2) located downstream of this high-rise area present more of the flood-related issues. Since in a highly dense urban area there are more impervious surfaces, it can be expected that this can have a great impact upon flood-related issues in downstream areas.

Hazard Assessment
The initial heat stress assessment was performed from mobile weather measurements. These measurements were performed in five different urban land uses within the case study area (M sites in Figure 5). Values of air temperature, wind speed, glob temperature and humidity were collected; and some thermal comfort indices including PMV and PET were calculated. The Kestrel 5400 Heat Stress Tracker was used for this purpose. This procedure showed that urban parks (M4) and open low-rise sites (M5) can be respectively 1.4 °C and 1.1 °C cooler than a highly dense urban site in the study area (M2 and M3). Furthermore, the results showed that the PMV and PET can be also lower in these two sites compared to highly dense urban sites. Notice that M3 and M5 are the sites chosen to perform the thermal modeling, presented as areas B and A, respectively, in Figure 3.
Additionally, by overlaying the flood hazard map with the buildings data ( Figure 5), it was observed that less floods and more heat stress are likely to occur in the highly dense upstream area (A1 in Figure 5), when compared other parts of the study area. However, less urbanized areas (A2) located downstream of this high-rise area present more of the flood-related issues. Since in a highly dense urban area there are more impervious surfaces, it can be expected that this can have a great impact upon flood-related issues in downstream areas.

NBS Types Selection and Suitability Analysis (Which and Where)
Four different small-scale NBS measures were selected from the analysis of local characteristics, namely the availability of possible locations for green roofs implementation, or low slope and low traffic pavements, which could be changed into pervious pavements:
Bio-retention (BR), with shrubs as topping vegetation at a height of 1.2 m 4.
Rain garden (RG), with street trees and lawn as topping vegetation at the height of 6 m From the suitability analysis, maximum rates of measures application where obtained. The results of these analysis showed that for the macro scale simulations (urban flooding), green roofs and pervious pavements can be implemented as an average on 27% and 15% of the whole study area, respectively. In addition, bio-retention and rain gardens can both be implemented as an average on 4% of the whole study area. For the micro scale simulations (thermal comfort), green roofs and pervious pavements can be implemented within a maximum of 36% and 27% of the selected micro scale study areas, respectively. While bio-retentions and rain gardens can both be implemented as a maximum of 5% of the selected microscale study area.
The two different results from suitability analysis are due to two different scales used in this study. Urban flood assessment and analysis requires a macro scale simulation study. However, the microclimatic thermal comfort assessment requires a microscale simulation study. The maximum possible application of the measures will be different when studying a macro scale site as a whole, and when studying the selected micro scale sites within the whole. Figure 6 shows the effectiveness of each NBS for flood mitigation according to the results obtained from the hydrodynamic model for two rainfall scenarios. The parameters presented are runoff volume reduction ( Figure 6a) and peak discharge reduction (Figure 6b). From the analysis of results we can observe that the effectiveness of the measures is reduced when the rainfall return period increases.

Effectiveness of NBS's on Urban Flooding (How Much Impact on Flood Reduction)
Additionally, it appears that 'green roofs' is the most efficient NBS type for this case study area, having effectiveness of up to 39% and 40% in reduction for total runoff volume and peak discharges, respectively, for a two-year return period rainfall. On the other hand, 'pervious pavements' was found to be the least effective NBS type. The main reason for green roofs for being the most effective is the relatively large suitable area for its implementation. According to the suitability analysis, around 27% of the area was considered suitable for implementing green roofs.

The Effectiveness of NBS on Thermal Comfort Enhancement (How much Impact on Thermal Comfort)
Regarding the microclimatic situation, model results for the base case and for the case of implementing rain gardens in open low-rise sites are shown as an example in Figure 7. The effectiveness of NBS was measured by comparing the results of variation in Air temperature (Ta) and Mean radiant temperature (Tmrt) with the base case scenario. In addition, variations of Predicted mean vote (PMV) in relation with the base case scenario were considered. The obtained results show that the cooling effect of trees, which were used in rain gardens, was widely dispersed if we analyze the reduction in air temperature. Whereas, Tmrt and PMV were significantly reduced, but only in the The results also show that there was no significant change in the delay for time to peak in most of the cases; thus this result was not plotted in Figure 6. One reason for this could be the existence of numerous small catchments in the case study area.

The Effectiveness of NBS on Thermal Comfort Enhancement (How Much Impact on Thermal Comfort)
Regarding the microclimatic situation, model results for the base case and for the case of implementing rain gardens in open low-rise sites are shown as an example in Figure 7. The effectiveness of NBS was measured by comparing the results of variation in Air temperature (Ta) and Mean radiant temperature (Tmrt) with the base case scenario. In addition, variations of Predicted mean vote (PMV) in relation with the base case scenario were considered. The obtained results show that the cooling effect of trees, which were used in rain gardens, was widely dispersed if we analyze the reduction in air temperature. Whereas, Tmrt and PMV were significantly reduced, but only in the shaded areas of the trees, as a result of the prevention of direct sun radiation. Figure 7 shows the distribution of the air temperature at 4:00 p.m. and at the height of 1.0 m from the ground. According to the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) [50], the thermal comfort measurements have to be done at the center of the human body, which internationally is established at 1.10 m of height. However, a vertical grid size configuration was used in this model, and the software only provided the vertical elevation data at the height of 1.0 m and 1.4 m. Therefore, the results in this work are presented for an elevation of 1.0 m, which is considered the pedestrian level.
The results from the microclimatic modeling for this case established that 'rain gardens' was the most effective measure in the open low-rise buildings (site A), with a maximum reduction of 0.66 • C in air temperature (Ta) (see Table 5). The reason why rain gardens performed best is that trees were considered as vegetation for this measure. Similar studies on the effect of different vegetation on thermal comfort also revealed that 'trees' was the most effective measure in providing outdoor thermal comfort, and this effectiveness can be enhanced by increasing the LAD (Leaf Area Density) and height of the trees [51][52][53]. In contrary, in a compact high-rise building (site B) rain gardens did not have the same effect. In site B, even though this measure was still the most effective in air temperature reduction, it had an impact of only 0.25 • C maximum reduction in Ta. This difference on the impact of the same measure in sites A and B is because more paved surfaces were exposed to sunlight in the low-rise site compared to the high-rise site. As a result, the application of trees prevented the sunlight to reach the paved surfaces in site A, which consequently lead to a higher decrease in the sensible heat fluxes. In other words, in site B the shadow provided by the tall buildings did not give the chance for the trees to further decrease the temperature by providing their own shadow.
Regarding bio-retention, this measure was the least effective in the open low-rise building or site A, with a maximum reduction capacity of 0.16 • C. This result confirmed the role of the tree's height and leaf area density (LAD) upon thermal comfort, since by reducing these parameters for the bio-retention case compared to the rain gardens case, the air temperature reduction was significantly less. On the other hand, pervious pavements had a good impact of 0.41 • C on air temperature reduction in the site A. The implementation of pervious pavements in the model was represented by changing the albedo of pavements from 0.4 to 0.8, in order to reduce the absorption of sun short wave radiation. It was observed that this measure had a better performance during the time that the sun was shining almost vertically, around 1:00 p.m. However, the other measures had their highest performance during the heat stress peak daytime, at 4 p.m. Sustainability 2019, 11, x FOR PEER REVIEW 12 of 28  Table 3); (b) Effectiveness of RG implementation on Ta reduction for scenario A-RG; (c) Effectiveness of RG implementation on PMV reduction for scenario A-RG; (d) Effectiveness of RG implementation on Tmrt reduction for scenario A-RG.  Table 3); (b) Effectiveness of RG implementation on Ta reduction for scenario A-RG; (c) Effectiveness of RG implementation on PMV reduction for scenario A-RG; (d) Effectiveness of RG implementation on Tmrt reduction for scenario A-RG.  Table 5 also shows that the effectiveness of the measures on Ta reduction was quite different when the site characteristic changed from open low-rise (A) to compact high-rise buildings (B). For instance, 'green roofs' had almost no effect on the reduction of Ta in scenario B, at 4:00 p.m. and at the pedestrian level. However, in site A, 'green roofs' had an effectiveness of 0.17 • C on reduction of air temperature at the pedestrian level. In fact, as the buildings in site B are very high (around 70 to 104 m), the effect of vegetation on top of such buildings did not reach the pedestrian level. There is previous research that revealed that the green roofs on top of high-rise buildings do not have any significant effect on air temperature for pedestrians [18,[54][55][56].
An interesting finding of this research was that in a site of high-rise buildings, tree plantation and vegetation can be more effective during the night time. The results clearly show that both measures, bio-retention and rain gardens at site B, had better effectiveness at 1:00 a.m., with 0.43 • C and 0.15 • C reduction in air temperature, respectively. Therefore, it can be derived that vegetation planting (especially in the high-rise building site) is the best practice for the urban heat island control, which is at its peak during night time. High-rise buildings absorb sun energy during the day and release it during the night. As a result, the temperature does not get reduced enough during the night time in these sites, and it causes an urban heat island (see also [55]), which can be controlled relatively by implementing NBS measures.
Regarding other parameters, the effectiveness of bio-retention and rain gardens on the reduction of Tmrt were almost the same in both site characteristics. The results show that these two measures had a maximum reduction in Tmrt of around 16 • C and 19 • C, respectively (at 4:00 p.m.) in both sites.
Whereas, 'green roofs' had also a reasonable contribution to Tmrt reduction in site A (Max 17.81 • C from 51.18 • C), it had a very low effect on Tmrt reduction at pedestrian level at site B (0.10 • C from 52.31 • C), which is logical, according to the previous discussion on the height of the buildings. Similarly, its effect on PMV max reduction was also 0.8 at site A and 0.01 at site B, which again does not show any significant change at site B.
The effect of pervious pavements on reduction of Tmrt is quite controversial. The results show that pervious pavements not only had no contribution in reduction of the Tmrt at any of the study areas, but it even increased slightly the value of Tmrt by 0.6 • C and 0.69 • C at sites A and B, respectively. This result shows that by changing the albedo of the streets from 0.4 to 0.8, the sun radiation was being more reflected to the atmosphere. Therefore, the pedestrians experienced more radiation, which led to a growth in Tmrt. Similarly, the impact on PMV for pervious pavements application at site B was slightly negative. However, there was a positive impact in site A, showing a change of 0.68 from 4.67 on PVM.

Discussion on NBS's Performance
In general, the simulations of the base case scenario for both sites (A and B) on thermal comfort show that during the day, and specially in the peak temperature time of the day (4:00 p.m.), site B (compact high-rise) provided a better thermal comfort at the pedestrian level, compared to site A (open low-rise). The main reason for this was related to the shadow provided during the day by the high-rise buildings. On the other hand, during the nights this issue reversed, and the site A had a better thermal comfort, as there was more air ventilation and less structures to release the heat accumulated during the day, to the environment during the night time. A research done by Hedquist and Brazel [57] on a case study of Arizona, U.S.A., presented a similar result.
An overall comparison of the measure's effectiveness in both scenarios (A and B) interestingly shows that the NBS measures had a better performance in the low-rise buildings (A) site during the daytime and a better performance in the high-rise buildings (B) site during the night time. This is because the more direct sunshine in site A during the day provided a better opportunity for NBS measures to improve the existing thermal comfort of the site. This was achieved by either providing more shadow or reradiating the sun's shortwaves. However, in site B, the shadow provided by the buildings during the day gave less opportunity for the NBS to show their effectiveness. On the other hand, as during the night, site B was hotter compared to site A, which is a result of the urban heat island [58], the implementation of NBS was more effective at site B during the night, which could contribute to the reduction of urban heat island intensity at these locations.
In conclusion, the effectiveness assessment of NBS measures using a microclimatic modeling clearly showed that the effectiveness of the measures on thermal comfort enhancement depends on several factors. These factors are: The characteristics of the site, the type of NBS, the coverage and location of the measures and the time of the day.

Comparative Effectiveness for Urban Flood Reduction
Analyzing the effectiveness scoring of the measures in relation to urban flood reduction (Table 6), it can be observed that 'green roofs' represent the most effective measure for this purpose while 'pervious pavements' appears as the least effective measure. This is in line with other studies, e.g., Carpenter and Kaluvakolanu [59] and Berardi et al. [60]. However, this effectiveness scoring does not consider the differences on percentages of implementation area for each measure. It is to be noted that one of the main reasons why 'green roofs' is the most effective NBS for control of urban flooding, is the fact that this measure is the most feasible to be implemented at such an urbanized part of the city. Therefore, again it is that site characteristics determine the performance of the measures.

Comparative Effectiveness for Thermal Comfort Enhancement
In relation to thermal control, 'rain gardens' is found to be the most effective measure in order to provide more thermal comfort (Table 7). Besides, 'pervious pavements' followed by 'green roofs' are the least effective measures for providing thermal comfort at the pedestrian level. Similar results were obtained in previous studies [61,62]. Total comparative performance score 4 16 9 12

Overall Analysis of Effectiveness and Recommendation for NBS Application
From the comparative rankings of the measures for flood reduction and thermal comfort enhancement in the Sukhumvit area in Bangkok (Thailand), it can be observed that 'rain gardens' are likely to be the most effective NBS type with respect to both criteria, flood reduction and thermal comfort enhancement (Table 8). It is interesting to observe that trees (with average height of six meters), which were part of the design of rain gardens in this case, represent the most influencing factor for urban microclimate conditions. As such, for this particular case study area this should be maintained in order to achieve better thermal comfort. Furthermore, design of rain gardens also included greater depth of soil zone for infiltration (e.g., 80 cm) purposes and the inclusion of a storage zone (e.g., 18 cm depth), which is likely to contribute towards higher effectiveness in flood volume reduction. For the Sukhumvit area, this measure should play an important role in urban planning activities. Green roofs and rain gardens are found to be as the second and third most effective NBS types, respectively.
Regardless of the overall effectiveness results, it is important to consider that different measures performed best on the two different criteria. As a result, to achieve the best performance in both objectives for the Sukhumvit case study area, it is recommendable to implement a combination of different NBS types. Thermal comfort enhancement (T) Low rise (A) 5 7 6 12 High rise (B) 5 5 8 12 Overall comparative score 14 28 23 36 This is possible, since the NBS types analyzed here do not compete for free spaces. For example, green roofs and rain gardens can be applied in a site simultaneously, since one uses roofs and the other one uses spaces alongside the streets.
Additionally, the location of the site as well as its land use type play a role on the performance of each NBS type, and as such, they should be considered for implementing the right measure. In this case, measures with higher effectiveness in flood reduction should be applied in the high-rise site, since it is at the upstream area, and has higher imperviousness. Reducing runoff through providing more disconnections in this case study area is likely to have a positive impact downstream, where flood impacts are higher. For Sukhumvit, the NBS type which showed best flood reduction performance is 'green roofs'. This can be explained due to its enhanced suitability that results on increased disconnection area; that leads to a better effectiveness for urban flood reduction when compared to other NBS types. However, this NBS type has low effectiveness on heat stress reduction at this compact high-rise site. Therefore, the best alternative for the site B is the combination of green roofs and rain gardens, which together would have a good performance on both criteria. Moreover, even though green roofs are less effective in relation to outdoor thermal comfort enhancement, their application may have many other benefits, such as energy consumption reduction, air pollution mitigation, storm water management, sound absorption and ecological preservation [63].

Conclusions
In this paper, a novel framework that can be used for the selection of small-scale urban nature-based solutions to reduce flooding and enhance human thermal comfort has been presented. The framework has been applied in the Sukhumvit area in Bangkok (Thailand). The obtained results show that the combined implementation of different NBS types is likely to have a good potential to make this area more resilient and sustainable to cope with future challenges related to climate change and the high rate of urbanization. By applying this novel framework, it was possible to identify the most promising NBS types that can be applied in different parts of the area to effectively achieve both objectives at the same time.
Several interesting findings were obtained from the present work. For instance, green roofs are likely to achieve better performance in the reduction of urban flooding when compared to the other NBS types studied. However, this particular NBS type is not effective in thermal comfort enhancement in sites with compact, high-rise buildings. Regarding the effectiveness of the NBS measures for thermal comfort, the results showed that this is mainly related to the provision of shadows from trees. Therefore, rain gardens with street trees as covering vegetation would have the best performance for the open low-rise scenario. Although, the results also showed that the effectiveness of different NBS types changes according to the site characteristics and time of the day. Therefore, a combined application of green roofs and rain gardens is recommended in compact high-rise building areas. We conclude then that this method is very helpful to identify adequate measures according to local characteristics, and to choose a combination of measures that is best for each particular site.
The importance of combing micro scale and macro scale effectiveness assessment of NBS was demonstrated through the present work. The use of microclimatic modeling in this framework showed that the effectiveness of NBS for thermal comfort enhancement depends on several factors. These factors are: The characteristics of the implementation site, the type of NBS, the coverage and location of the measures, and the time of day. Moreover, the results illustrate the usefulness of macro scale urban flood modeling for an assessment of the effectiveness of different NBS types for the reduction of flood impacts. Consequently, the present work proves the effectiveness of different NBS types taken from a micro and macro scale perspectives.
The outcome of the present research aims to provide some additional knowledge to city planners and decision makers in gaining better understanding of the effectiveness of different NBS measures for different sites and local conditions.

Appendix B
Pervious Pavements Figure A2. Schematic view of a pervious pavement system   Green Roofs Figure A3. Schematic view of a green roof system.

Appendix C
Site selection assessment for implementing the measures A GIS extension tool called the Sustain-BMP siting tool, developed by the United States Environmental Protection Agency (US EPA), was used to analyze possible locations of NBS measures in the area. This tool was used only in the case of BR and RG, because the implementation of these two measures needs more detailed site feasibility assessment than in the case of GR and PP. This tool has been used for similar studies in several previous cases [67][68][69]. Figure A6 shows the overall methodology for the suitability analysis of the four selected NBS measures.
The input data for this tool, such as land use and a two meters' resolution DEM, were provided for the Bangkok Metropolitan Administration office. Data about soil types, imperviousness and ground water level were obtained from the Department of Drainage and Sewage in Bangkok. The input data for this tool, such as land use and a two meters' resolution DEM, were provided for the Bangkok Metropolitan Administration office. Data about soil types, imperviousness and ground water level were obtained from the Department of Drainage and Sewage in Bangkok.

Mike Urban Simulation
A 1D/2D modeling simulation for producing the flood hazard map had been run by Mike Urban in different scenarios of implementing the multifunctional measures in the case study area. This is a well-recognized software for hydrodynamic modeling [71,72]. Additionally, the existing sewer network of Sukhumvit was also initially modeled using Mike Urban, which facilitated the building of the model for this study. Two modifications to the previously built model of Sukhumvit were introduced: Change in runoff routing and change in the simulation engine. Figure A7 shows the overall methodology followed in Mike Urban simulations to study the effectiveness of NBS measures on urban flooding.

Mike Urban simulation
A 1D/2D modeling simulation for producing the flood hazard map had been run by Mike Urban in different scenarios of implementing the multifunctional measures in the case study area. This is a well-recognized software for hydrodynamic modeling [70,71]. Additionally, the existing sewer network of Sukhumvit was also initially modeled using Mike Urban, which facilitated the building of the model for this study. Two modifications to the previously built model of Sukhumvit were introduced: Change in runoff routing and change in the simulation engine. Figure A7 shows the overall methodology followed in Mike Urban simulations to study the effectiveness of NBS measures on urban flooding. The input data for this tool, such as land use and a two meters' resolution DEM, were provided for the Bangkok Metropolitan Administration office. Data about soil types, imperviousness and ground water level were obtained from the Department of Drainage and Sewage in Bangkok.

Mike Urban Simulation
A 1D/2D modeling simulation for producing the flood hazard map had been run by Mike Urban in different scenarios of implementing the multifunctional measures in the case study area. This is a well-recognized software for hydrodynamic modeling [71,72]. Additionally, the existing sewer network of Sukhumvit was also initially modeled using Mike Urban, which facilitated the building of the model for this study. Two modifications to the previously built model of Sukhumvit were introduced: Change in runoff routing and change in the simulation engine. Figure A7 shows the overall methodology followed in Mike Urban simulations to study the effectiveness of NBS measures on urban flooding.  ENVI-met microclimatic simulation ENVI-met is a three-dimensional non-hydrostatic climate model for the simulation of surface-plant-air interactions, especially for conditions inside urban environments. Since the model is designed for the microscale, the resolution output is high, ranging from 0.5 to 10 m horizontally, with a temporal resolution of 10 sec and the ability to simulate timeframes from 24 to 48 h. The model requires the user to input certain parameters, such as defining the model area (area input file) and configuring the initial atmospheric conditions, surfaces (including soils), vegetation, and time intervals [72]. Figure A8 shows the conceptual framework describing the required procedure and steps for running the microclimatic simulations in ENVI-met and visualizing and analyzing the results using LEONARDO 2014 and Biomet, respectively.

ENVI-met microclimatic Simulation
ENVI-met is a three-dimensional non-hydrostatic climate model for the simulation of surfaceplant-air interactions, especially for conditions inside urban environments. Since the model is designed for the microscale, the resolution output is high, ranging from 0.5 to 10 m horizontally, with a temporal resolution of 10 sec and the ability to simulate timeframes from 24 to 48 h. The model requires the user to input certain parameters, such as defining the model area (area input file) and configuring the initial atmospheric conditions, surfaces (including soils), vegetation, and time intervals [73]. Figure A8 shows the conceptual framework describing the required procedure and steps for running the microclimatic simulations in ENVI-met and visualizing and analyzing the results using LEONARDO 2014 and Biomet, respectively.

Appendix D
The selection of urban climate zones follows the urban climate zone categorization proposed by Stewart and Oke [35]. The zone numbers 1 & 6 of this categorization (see Figure A9) were used as sites B & A in microclimatic simulation. Figure A9. Characteristics of zones chosen for microclimatic simulation (modified from [35]).

Appendix D
The selection of urban climate zones follows the urban climate zone categorization proposed by Stewart and Oke [35]. The zone numbers 1 & 6 of this categorization (see Figure A9) were used as sites B & A in microclimatic simulation.

ENVI-met microclimatic Simulation
ENVI-met is a three-dimensional non-hydrostatic climate model for the simulation of surfaceplant-air interactions, especially for conditions inside urban environments. Since the model is designed for the microscale, the resolution output is high, ranging from 0.5 to 10 m horizontally, with a temporal resolution of 10 sec and the ability to simulate timeframes from 24 to 48 h. The model requires the user to input certain parameters, such as defining the model area (area input file) and configuring the initial atmospheric conditions, surfaces (including soils), vegetation, and time intervals [73]. Figure A8 shows the conceptual framework describing the required procedure and steps for running the microclimatic simulations in ENVI-met and visualizing and analyzing the results using LEONARDO 2014 and Biomet, respectively.

Appendix D
The selection of urban climate zones follows the urban climate zone categorization proposed by Stewart and Oke [35]. The zone numbers 1 & 6 of this categorization (see Figure A9) were used as sites B & A in microclimatic simulation. Figure A9. Characteristics of zones chosen for microclimatic simulation (modified from [35]). Figure A9. Characteristics of zones chosen for microclimatic simulation (modified from [35]).