Multicriteria GIS-Based Approach in Priority Areas Analysis for Sustainable Urban Drainage Practices: A Case Study of Pato Branco, Brazil

Geographic information systems (GIS) and multicriteria decision methods are robust techniques for supporting the urban planning process, including urban drainage. New low-impact management approaches (LID) for rainwater have been investigated and have become increasingly used. Considering the central region of Pato Branco city, which suffers recurrent problems related to flooding, this work presents a method to identify potential areas for the application of LIDs, such as rainwater collection tanks, permeable pavements, green roofs, and rain gardens. The identification of these areas is based on the analysis in a GIS environment considering criteria related to both the land slope, the characteristics of land use and technical parameters. Thus, we observe that rainwater collection tanks are indicated for all habitations, permeable pavements are recommended for 6.30% of the study area, while green roofs can be implemented in 3.97% of the area. Finally, 3.03% can receive rain gardens. In total, 13.30% of the central region of Pato Branco can receive LIDs. The results obtained reveal that the use of the GIS tool associated with multicriteria analysis is efficient in choosing locations for the implementation of LIDs as alternatives for the management of urban drainage.


Introduction
The urbanization process impacts on the urban environment and has caused recurrent damages and losses to the population. Among the main problems are the changes in land use [1], waterproofing of the soil surface [2], floods [3], and changes in the urban hydrological cycle [4]. Consequently, these impacts increase the volume of runoff [5], changing the quality of the water that drains superficially [6] and the decrease in underground recharge [7].
Floods are worrisome because approximately 80% of the world population [8] is subject to financial damage, destruction of infrastructure, and deaths due to this phenomenon [9,10]. In addition, the socioeconomic risks associated with floods are more evident when the phenomenon occurs in an urban area compared to rural areas [11], these risks have led to a diversity of studies on the causes of floods [10,12,13].
Climate change can increase both frequency and magnitude of floods on a global scale [14]. The increase in emissions of gases that contribute to the greenhouse effect can contribute to reducing the period of return of extreme precipitation, making flash floods more frequent [15,16].
Besides, the physical factors that are generally linked to the occurrence of floods are the gradient of the hydrographic basin, the use and cover of the soil, the permeability of the soil, and the physical Eng 2020, 1 98 aspects of the implementation of sustainable drainage systems in consolidated urban areas. The authors performed a multiscale analysis, that is, at the city, local and microscale levels.
Thus, given the problems related to recurrent flooding in urbanized areas, this work aims to identify potential areas for the application of LIDs, using the central region of the Pato Branco city as a study area. We use a methodology using GIS-based multicriteria analysis for the integration of geospatial data and technical parameters. To enable data integration, the data were normalized by fuzzy logic. The exclusion of inappropriate areas was applied using Boolean logic in map algebra. The sustainable urban drainage practices chosen (permeable pavement, green roof, rain garden, and RHW) were chosen due to easy implementation in areas with consolidated urbanization. The identification was performed using GIS combined with the analysis of the slope and the characteristics of land use.

Materials and Methods
In this work, we investigate a method for selection of areas with the potential for an application of four types of LID: 1. rain gardens, 2. permeable pavements, 3. rainwater harvesting, and 4. green roofs. We have opted for these LIDs because of the ease of implementation in areas with consolidated urbanization, with retrofit capacity in a wide range of land use configurations [24,44,45]. In addition, they are sustainable practices, common for rainwater management, with benefits in reducing flow volumes and peak hydrograph [46][47][48][49]. The method workflow is represented in Figure 1.
Eng 2020, 1, FOR PEER REVIEW 3 environmental, social, and economic aspects of the implementation of sustainable drainage systems in consolidated urban areas. The authors performed a multiscale analysis, that is, at the city, local and microscale levels. Thus, given the problems related to recurrent flooding in urbanized areas, this work aims to identify potential areas for the application of LIDs, using the central region of the Pato Branco city as a study area. We use a methodology using GIS-based multicriteria analysis for the integration of geospatial data and technical parameters. To enable data integration, the data were normalized by fuzzy logic. The exclusion of inappropriate areas was applied using Boolean logic in map algebra. The sustainable urban drainage practices chosen (permeable pavement, green roof, rain garden, and RHW) were chosen due to easy implementation in areas with consolidated urbanization. The identification was performed using GIS combined with the analysis of the slope and the characteristics of land use.

Materials and Methods
In this work, we investigate a method for selection of areas with the potential for an application of four types of LID: 1. rain gardens, 2. permeable pavements, 3. rainwater harvesting, and 4. green roofs. We have opted for these LIDs because of the ease of implementation in areas with consolidated urbanization, with retrofit capacity in a wide range of land use configurations [24,44,45]. In addition, they are sustainable practices, common for rainwater management, with benefits in reducing flow volumes and peak hydrograph [46][47][48][49]. The method workflow is represented in Figure 1.

Study Area
The study area comprises the central urban area of Pato Branco city ( Figure 2). Pato Branco is located at coordinates 52°43' S and 52°37' W in Paraná State, Brazil's south region. The total geographical area is 539,029 km 2 . The central urban area was chosen because that region is susceptible to recurrent flooding and it is a place with high traffic of people and economic, commercial, and residential interest.

Study Area
The study area comprises the central urban area of Pato Branco city ( Figure 2). Pato Branco is located at coordinates 52 • 43 S and 52 • 37 W in Paraná State, Brazil's south region. The total geographical area is 539,029 km 2 . The central urban area was chosen because that region is susceptible to recurrent flooding and it is a place with high traffic of people and economic, commercial, and residential interest. According to the Brazilian Institute of Geography and Statistics (IBGE), Pato Branco has an estimated population of 82,881 inhabitants in 2019 [50], with 94.09% residing in the urban area and 5.91% in the rural area. Therefore, the municipality has an approximate demographic density of 134.25 inhabitants/km 2 . This region has a temperate climate of the Cfb type, with well-distributed rains throughout the year. According to the database of the Instituto de Águas do Paraná (Paraná Water Institute), Pato Branco has an average annual rainfall of 2089 mm (considering 30 years of historical data).

Identification of Potential Areas
The methodology proposed here is based on the multicriteria analysis in a Geographic Information System (GIS) environment. For this, the ArcGIS 10 software was used. In the GIS environment, all the necessary information was organized in layers, sharing a single database. To ensure geospatial compatibility, all information layers were spatially referenced according to the Brazilian Geodetic System, which adopts the Reference System for the Americas-2000 (SIRGAS 2000).
The slope map (Figure 3b) of the study area was generated from the elevation contour every 5 m. In this step, the Digital Elevation Model (DEM) (Figure 3a) was generated from the interpolation by the topo to raster method (ANUDEM) [51]. ANUDEM is a method of morphological interpolation from scattered surface elevation points, elevation contours, streamlines, and others. This model is considered as hydrologically correct to produce DEMs [52,53]. Thus, the slope raster is calculated from the maximum rate in elevation over distance, considering the DEM cell and its eight neighbors. The slope value is calculated by applying the average maximum technique [54]. According to the Brazilian Institute of Geography and Statistics (IBGE), Pato Branco has an estimated population of 82,881 inhabitants in 2019 [50], with 94.09% residing in the urban area and 5.91% in the rural area. Therefore, the municipality has an approximate demographic density of 134.25 inhabitants/km 2 . This region has a temperate climate of the Cfb type, with well-distributed rains throughout the year. According to the database of the Instituto de Águas do Paraná (Paraná Water Institute), Pato Branco has an average annual rainfall of 2089 mm (considering 30 years of historical data).

Identification of Potential Areas
The methodology proposed here is based on the multicriteria analysis in a Geographic Information System (GIS) environment. For this, the ArcGIS 10 software was used. In the GIS environment, all the necessary information was organized in layers, sharing a single database. To ensure geospatial compatibility, all information layers were spatially referenced according to the Brazilian Geodetic System, which adopts the Reference System for the Americas-2000(SIRGAS 2000. The slope map (Figure 3b) of the study area was generated from the elevation contour every 5 m. In this step, the Digital Elevation Model (DEM) (Figure 3a) was generated from the interpolation by the topo to raster method (ANUDEM) [51]. ANUDEM is a method of morphological interpolation from scattered surface elevation points, elevation contours, streamlines, and others. This model is considered as hydrologically correct to produce DEMs [52,53]. Thus, the slope raster is calculated from the maximum rate in elevation over distance, considering the DEM cell and its eight neighbors. The slope value is calculated by applying the average maximum technique [54].
To obtain the updated cartographic base of the study area, we used multispectral images from the WorldView satellite with a spatial resolution of 50 cm. In this step, the images were vectorized to extract the features of interest according to each LID model. Sidewalks, parking, and squares have been identified as potential areas for permeable pavements. Buildings with slabs were identified for the implementation of green roofs. Finally, parks and areas with vegetation are of interest for the installation of rain gardens. Eng 2020, 1, FOR PEER REVIEW 5 To obtain the updated cartographic base of the study area, we used multispectral images from the WorldView satellite with a spatial resolution of 50 cm. In this step, the images were vectorized to extract the features of interest according to each LID model. Sidewalks, parking, and squares have been identified as potential areas for permeable pavements. Buildings with slabs were identified for the implementation of green roofs. Finally, parks and areas with vegetation are of interest for the installation of rain gardens.

Identification of Conditions for the Implementation of LIDs
The constraints related to the implantation of LIDs in urban areas were defined according to the data observed in the literature [33,35,38,[55][56][57] (Table 1). In the table, the "required area size" indicates the criteria for implementing the LID in terms of the required area extension. Slope restrictions are established by design criteria, where low slope areas will ensure lower investment costs for surface modification, higher natural storage capacities, and longer infiltration times due to low flow velocities compared to the high slope areas [58]. Land use refers to the type of cover to which each LID is associated.

Multicriteria Analysis
After the identification of the areas that meet the requirements for LIDs implementation, we proceeded with the exclusion of areas that were not following the adopted premises (Table 1). For this, it was applied to the map algebra between the map of the possible place of LIDs application and the restrictions maps to verify the areas where the LID can be applied.

Identification of Conditions for the Implementation of LIDs
The constraints related to the implantation of LIDs in urban areas were defined according to the data observed in the literature [33,35,38,[55][56][57] (Table 1). In the table, the "required area size" indicates the criteria for implementing the LID in terms of the required area extension. Slope restrictions are established by design criteria, where low slope areas will ensure lower investment costs for surface modification, higher natural storage capacities, and longer infiltration times due to low flow velocities compared to the high slope areas [58]. Land use refers to the type of cover to which each LID is associated.

Multicriteria Analysis
After the identification of the areas that meet the requirements for LIDs implementation, we proceeded with the exclusion of areas that were not following the adopted premises (Table 1). For this, it was applied to the map algebra between the map of the possible place of LIDs application and the restrictions maps to verify the areas where the LID can be applied.
Thus, the areas suitable for the permeable pavements consist of sidewalks, parking lots, and squares with a slope of less than 5%. For rain garden adoption, we selected parks and areas with vegetation and slope of less than 12%. Regarding green roofs, the only limiting factor for their use is the type of roof, which must be of the flat slab type.
In this work, we have integrated all information levels based on the GIS environment. Thus, the geospatial data and technical criteria have been normalized using the fuzzy technique, as described in [59]. Fuzzy logic has been extensively applied in various studies, primarily in the treatment of spatial data and real-world representation [58,[60][61][62]. Areas unsuitable for the installation of the LID were excluded by applying Boolean logic. The Boolean model involves a logical combination of binary maps using conditional operators of type A AND B (returns the intersection A ∩ B), A NOT B (returns the elements contained only in A), A OR B (returns the elements of A and B except those at the intersection A ∩ B). The result is expressed in binary form, 0 for not satisfied, and 1 for satisfied.

Flood Susceptibility Map
To contextualize the potential areas for the implementation of LIDs concerning the occurrence of floods or critical areas for runoff, we also generated the map of susceptibility to flooding events. The overall workflow for this analysis is represented in Figure 4.
squares with a slope of less than 5%. For rain garden adoption, we selected parks and areas with vegetation and slope of less than 12%. Regarding green roofs, the only limiting factor for their use is the type of roof, which must be of the flat slab type.
In this work, we have integrated all information levels based on the GIS environment. Thus, the geospatial data and technical criteria have been normalized using the fuzzy technique, as described in [59]. Fuzzy logic has been extensively applied in various studies, primarily in the treatment of spatial data and real-world representation [58,[60][61][62]. Areas unsuitable for the installation of the LID were excluded by applying Boolean logic. The Boolean model involves a logical combination of binary maps using conditional operators of type A AND B (returns the intersection A ∩ B), A NOT B (returns the elements contained only in A), A OR B (returns the elements of A and B except those at the intersection A ∩ B). The result is expressed in binary form, 0 for not satisfied, and 1 for satisfied.

Flood Susceptibility Map
To contextualize the potential areas for the implementation of LIDs concerning the occurrence of floods or critical areas for runoff, we also generated the map of susceptibility to flooding events. The overall workflow for this analysis is represented in Figure 4. In this case, for the multicriteria analysis, the AHP process was used [63]. The AHP technique is a compensatory mathematical model [64], which was used to define the weights of the criteria. Once weights were defined, the flood susceptibility map was generated from map algebra using a weighted linear combination (WLC). Thus, the map algebra to generate the flood susceptibility map (FSM) was applied according to Equation (1)

Results and Discussion
Initially, to identify the areas with potential for the use of LID in the central area of Pato Branco, land use was observed and related to the selection criteria for each model (Figures 5-7). Thus, note that the area for permeable pavements is larger in relation to the others, followed by the rain garden and, finally, the green roof.  In this case, for the multicriteria analysis, the AHP process was used [63]. The AHP technique is a compensatory mathematical model [64], which was used to define the weights of the criteria. Once weights were defined, the flood susceptibility map was generated from map algebra using a weighted linear combination (WLC). Thus, the map algebra to generate the flood susceptibility map (FSM) was applied according to Equation (1):

Results and Discussion
Initially, to identify the areas with potential for the use of LID in the central area of Pato Branco, land use was observed and related to the selection criteria for each model (Figures 5-7). Thus, note that the area for permeable pavements is larger in relation to the others, followed by the rain garden and, finally, the green roof.
The predominance of areas for permeable sidewalks and rain gardens is explained by the prevalence of parking lots, sidewalks, parks, and vacant lots in the area. Thus, it is interesting that these systems are adopted if there is any further modification in the areas that can be built on.
Melville-Shreeve et al. [65], who was interested in understanding the use of LID in the United Kingdom, also identified the wide use of permeable pavement concerning other models. The authors identified that one of the reasons for the results is the fact that in this system, water can be stored or, still, be used as a way of improving water quality. Besides, Jato-Espino et al. [66] also comment in their study that permeable pavements are the easiest systems to implement in a consolidated urban area, due to their multifunctionality, since they act in the infiltration and storage of rainwater and allow the control of runoff.     The predominance of areas for permeable sidewalks and rain gardens is explained by the prevalence of parking lots, sidewalks, parks, and vacant lots in the area. Thus, it is interesting that these systems are adopted if there is any further modification in the areas that can be built on.
Melville-Shreeve et al. [65], who was interested in understanding the use of LID in the United Kingdom, also identified the wide use of permeable pavement concerning other models. The authors identified that one of the reasons for the results is the fact that in this system, water can be stored or, still, be used as a way of improving water quality. Besides, Jato-Espino et al. [66] also comment in their study that permeable pavements are the easiest systems to implement in a consolidated urban area, due to their multifunctionality, since they act in the infiltration and storage of rainwater and allow the control of runoff.
Regarding green roofs, the criterion of the existence of buildings with flat slab coverage was considered. Due to flat slabs occurrences being higher in buildings with more than four floors, green roofs provide an alternative to mitigate floods in areas with a higher rate of vertical urbanization, as well as to contribute to the reduction of heat islands [67]. Thus, there is a predominance of singlestory buildings or with less than four floors in the area, which would require major adaptations to receive a green roof.
For geomorphology (Figure 8), the study area has a rugged relief. The slope is characterized by the predominance of slopes between 0% and 20%, with some regions reaching values close to 70%. The slope level is classified as nearly slope (<3%), gentle slope (3%-5%), sloping slope (5%-12%), moderate slope (12%-20%), moderate steep slope (20%-45%), and steep (>45%). Regarding green roofs, the criterion of the existence of buildings with flat slab coverage was considered. Due to flat slabs occurrences being higher in buildings with more than four floors, green roofs provide an alternative to mitigate floods in areas with a higher rate of vertical urbanization, as well as to contribute to the reduction of heat islands [67]. Thus, there is a predominance of single-story buildings or with less than four floors in the area, which would require major adaptations to receive a green roof.
We observed that the central area has a predominant slope of 0% to 5% (Figure 8). With consolidated urbanization, this area has the potential for retrofit, mainly with the transformation of paved areas into permeable pavements. Rain gardens also require land with slopes up to 12%, however, the areas identified for its implementation are located mainly on the sides of the study area, which have a higher incidence of moderate to accentuated slopes.
Considering the slope of the site (Figure 8), it is visible that the favorable area for permeable pavements and rain gardens decreases from the area originally expected (Figures 9 and 10). This is explained by the imposed slope limitation (5% for permeable floors and 12% for rain gardens). In greater inclinations, the implementation of LID requires adaptations and studies of its own, which would demand a higher cost and may make its execution unfeasible.
Concerning the area destined for each LID, with the use of the soil as the main condition, it is noted that of the total study area (1,857,630.05 m 2 ), 25.87% is available or fits the selection criteria ( Table 2). After applying the conditions, the possibility of using permeable pavements and rain gardens decreases (59.59% and 51.98%, respectively), resulting in a total of 13.30% of the area. Eng 2020, 1, FOR PEER REVIEW 9 Figure 8. Slope map.
We observed that the central area has a predominant slope of 0% to 5% (Figure 8). With consolidated urbanization, this area has the potential for retrofit, mainly with the transformation of paved areas into permeable pavements. Rain gardens also require land with slopes up to 12%, however, the areas identified for its implementation are located mainly on the sides of the study area, which have a higher incidence of moderate to accentuated slopes.
Considering the slope of the site (Figure 8), it is visible that the favorable area for permeable pavements and rain gardens decreases from the area originally expected (Figures 9 and 10). This is explained by the imposed slope limitation (5% for permeable floors and 12% for rain gardens). In greater inclinations, the implementation of LID requires adaptations and studies of its own, which would demand a higher cost and may make its execution unfeasible.   Concerning the area destined for each LID, with the use of the soil as the main condition, it is noted that of the total study area (1,857,630.05 m 2 ), 25.87% is available or fits the selection criteria  Thus, the results (Table 2) show that permeable pavements are recommended for 6.30% of the entire study area, green roofs can be implanted in 3.97% of the area and 3.03% of this can receive gardens of rain. In total, 13.30% of the central region of Pato Branco can receive LID.
As for rainwater harvesting, there are no restrictions on their use. In addition, a variation of tank sizes is possible, according to the demand or built area. For the analysis of the application of rainwater harvesting, a survey of the construction areas in the study area was carried out, where it was found that there are 1536 properties (Table 3). To evaluate its application, three tank sizes were defined according to the built area of each building. Thus, tanks of 3, 5, and 10 m 3 will be associated with homes with less than 100 m 2 , from 100 to 300 m 2 and above 300 m 2 , respectively. Thus, it was observed that the use of rainwater harvesting could store up to 11,145,000 L of rainwater. Steffen et al. [68], in their study, they observed that rainwater storage tanks have a high potential for saving drinking water, especially in regions with annual rainfall above 762 mm. However, the main obstacles that may exist in the adoption of rainwater harvesting are the high capital costs, mainly in single-family houses, the maintenance costs in multi-family buildings, and the lack of knowledge of the system by the population [69].
Through the data obtained, it is noticeable that rainwater harvesting is a LID method that can be widely used due to the availability of locations. However, infiltration-based typologies, such as permeable pavements, have a better performance in reducing runoff [70]. Figure 11 shows a summary of all areas suitable for the application of LID. Observe a higher concentration of areas suitable for the implementation of green roofs and permeable pavements in regions with greater urbanization and soil drainage. The rain gardens are more notable in areas far from the central town. Besides, the central core area is characterized by slight inclination, which allows the implementation of systems that have the slope as a limiting factor.
Eng 2020, 1, FOR PEER REVIEW 12 Figure 11. Areas indicated for the use of LID models. Figure 12 shows the flood susceptibility map. As already mentioned, this map considers factors such as land use, rainfall, slope, and geology. Analyzing this map, note that the central region of the study area tends to receive the largest flow of runoff and, consequently, presents the highest risk of flooding. In fact, from population reports, floods occur annually in the high susceptibility area shown on the map.  Figure 12 shows the flood susceptibility map. As already mentioned, this map considers factors such as land use, rainfall, slope, and geology. Analyzing this map, note that the central region of the study area tends to receive the largest flow of runoff and, consequently, presents the highest risk of flooding. In fact, from population reports, floods occur annually in the high susceptibility area shown on the map. Analyzing Figures 11 and 12 it is noticeable that the permeable pavement is one of the types of LIDs indicated for that region. Currently, this area does not have any type of LID that allows the permeability of the surface.
The GIS-based multicriteria analysis proved to be efficient at choosing locations with potential for LIDs installation. In this investigation, we have presented a method to identify areas that meet all the established requirements, which helps companies and government officials in the decisionmaking process.
Moreover, observe that for the use of LIDs in environments with consolidated urbanization should be considered the land use and the slope as factors for choosing the most appropriate models to the reality of the local. In this sense, the use of LIDs, combined with the conventional urban drainage systems that already exist, contributes to the natural water cycle in the urban environment.
In this study, no distinctions were made between public and private properties. One of the biggest problems related to the implantation of sustainable urban drainage systems is that, for its effectiveness, there is a dependence on the adhesion of private owners, which makes the development of legal instruments essential to enable the adoption of these practices. Another problem, mentioned by Baptista and Nascimento [71], is related to the search, by individuals, for the most economical solution to the detriment of the most ecological.
Although rainwater management is expensive, the costs of damage resulting from rainwater runoff can be even higher. The costs of rainwater damage, including damage from property flooding, degradation of water quality and loss of natural resources, reveal that the potential benefits outweigh the cost of rainwater management.
Urban drainage management and planning must start at the level of individual property, moving forward with linking these activities to the neighborhood level and, finally, scaling them up to the level of the hydrographic basin. Moreover, public education on the need for stormwater mitigation and the benefits of effective stormwater management can motivate the public and authorities to increase efforts and investments in stormwater programs [72,73]. The effective Analyzing Figures 11 and 12 it is noticeable that the permeable pavement is one of the types of LIDs indicated for that region. Currently, this area does not have any type of LID that allows the permeability of the surface.
The GIS-based multicriteria analysis proved to be efficient at choosing locations with potential for LIDs installation. In this investigation, we have presented a method to identify areas that meet all the established requirements, which helps companies and government officials in the decision-making process.
Moreover, observe that for the use of LIDs in environments with consolidated urbanization should be considered the land use and the slope as factors for choosing the most appropriate models to the reality of the local. In this sense, the use of LIDs, combined with the conventional urban drainage systems that already exist, contributes to the natural water cycle in the urban environment.
In this study, no distinctions were made between public and private properties. One of the biggest problems related to the implantation of sustainable urban drainage systems is that, for its effectiveness, there is a dependence on the adhesion of private owners, which makes the development of legal instruments essential to enable the adoption of these practices. Another problem, mentioned by Baptista and Nascimento [71], is related to the search, by individuals, for the most economical solution to the detriment of the most ecological.
Although rainwater management is expensive, the costs of damage resulting from rainwater runoff can be even higher. The costs of rainwater damage, including damage from property flooding, degradation of water quality and loss of natural resources, reveal that the potential benefits outweigh the cost of rainwater management.
Urban drainage management and planning must start at the level of individual property, moving forward with linking these activities to the neighborhood level and, finally, scaling them up to the level of the hydrographic basin. Moreover, public education on the need for stormwater mitigation and the benefits of effective stormwater management can motivate the public and authorities to increase efforts and investments in stormwater programs [72,73]. The effective management of urban waters must also be attentive to the projected changes to the climate, as they should lead to an increase in temperatures, frequency, and intensity of extreme rain events [74].
In this context, the effective implementation of sustainable urban drainage systems must be a mutual effort between the population and urban planners. Thus, the knowledge of the retrofit techniques appropriate for each region, as well as the scale and to what extent they can be applied, helps managers and decision-makers to plan urban drainage systems sustainably and effectively. As for property owners, the adoption of sustainable techniques can cause, in addition to the aesthetic benefit, an increase in the economic value of the property.

Conclusions
Stimulated by the challenge of adapting cities to climate change, rainwater management is one of the priorities for sustainable development. In this sense, LIDs aim to restore the hydrological flow regime of urban areas, promoting infiltration, retention and evapotranspiration in the place. However, the effectiveness of LIDs is influenced not only by their design, but also by their location.
In a general context, the methodology used can identify areas with potential applications of LIDs to improve the management of rainwater. In this sense, the methods analyzed are alternatives for the management of urban drainage and, consequently, can minimize the impacts caused by the urbanization process and contributes to the natural water cycle in the urban environment.
Finally, we can conclude that the use of geoprocessing and multicriteria analysis tools are effective in the management of the urban territory. The proposed methodology assists individuals who may be interested in using one of the LIDs techniques in their homes, as well as companies and government agencies in aid decision-making and urban planning with a focus on sustainable urban drainage. Thus, the adoption of LID can bring numerous benefits, acting in an integrated manner with existing drainage structures. Even without the elimination of flood events, the adoption of sustainable drainage techniques contributes to the mildness of the climate, increased biodiversity, and the quality of the water that reaches the rivers.
The proposed model can be applied to other regions that are also susceptible to flooding, with adaptation to local characteristics and scenery. Besides, the model can be improved by adding criteria such as social and economic characteristics, which can be incorporated into the modeling proposed here. For practical use of the results, modeling and field studies must be considered to take into account all the complex watershed processes necessary to define the project specifications before implementation. In further investigations, the complexity of the model will be increased when an adaptive decision tree will be implemented, based on social, economic criteria, public facilities, and other dimensions that can impact the concept of the smart city.