Conceptual and Analytical Framework as Flood Risk Mapping Subsidy

: There are still gaps in deﬁning values and category classiﬁcations of exposed items in quantitative damage analysis. This paper proposes a framework that reﬁnes the development of ﬂood risk analysis at a local scale. This study ﬁrst performs a quantitative risk analysis, based mainly on secondary data; it then attempts to communicate the results graphically, aiming to reduce the ﬁnancial and human resources required. We propose an easily standardized database in a GIS environment, analyzing the inﬂuence of a reservoir for ﬂood control and the construction of replicable local-scale risk curves. Hydrological (HEC-HMS) and 2D hydrodynamic (HEC-RAS) models were used to simulate hydrographs considering different return periods. For damage estimation, the processing included vectorization of lots, building use deﬁnition with Google Street View, classiﬁcation of standard designs, and a ﬁeld survey to validate those classes. In monetary value, this study calculated the effect of the construction of a reservoir for damage reduction, showing the potential to determine the effectiveness of measures adopted to mitigate ﬂood impacts. In addition, for each simulated return period, exposure, hazard, and damage maps can be established, making it possible to perform a complete risk analysis. of the still Synthesis, a relevant element in scientiﬁc work, is associated with locating, sorting, evaluating, and combining information. It can usually be achieved through transparent, reproducible, and organized procedures for summarizing such information [21]. This study is comprised of hydrological and hydrodynamic simulations, damage estimation, and risk map generation. Figure 2 shows the steps required to achieve these goals.


Introduction
In the current global context of climate emergencies, there is a need for awareness and ownership of the concept of risk so that a paradigm shift can be carried out [1]. The goal is to leave the fixed idea of trying to control extreme natural events and move towards an idea of adapting our reality into a more sustainable one. Move from resistance to resilience, from control to governance. The focus is no longer just on managing the disaster and the after-effects, but also more emphatically on risk management, that which precedes, anticipates, and prevents. These ideas have been proposed in the Hyogo Framework for Action and the Sendai Framework for Disaster Risk Reduction [2]. The latter are documents adopted by members of the United Nations at conferences seeking to promote integrated and inclusive measures. The Sendai Framework highlights the need to understand the specific local characteristics and to collaborate at various scales and between different sectors. The four priorities defined are understanding risk, strengthening governance, investing in resilience, and increasing preparedness for more effective responses ("Build Back Better"; recovery, rehabilitation, and reconstruction).
In this context, some of the challenges in management are bureaucracy, discontinuity of public policies, and lack of cooperation between stakeholders [3], as well as having to deal with anthropogenic changes, such as environmental degradation and accelerated urbanization [4]. In this scenario, risk analyses have been increasingly recognized as necessary to develop more efficient management [5]. In Brazil, the work of civil defenses against disaster scenarios related to natural phenomena at the municipal level stands out. However, disaster management is more applicable to the country than risk management. Following the terrain profile, other physical characteristics vary in the longitudina direction. The coastal area is hot and humid, with high rates of rainfall, especially in June and Atlantic tropical evergreen forests, which have been extensively modified because o the cultivation of sugar cane, mangroves, and palm trees. In the westernmost part of th basin, the vegetation cover is similar to that of a semiarid region, with xerophilic and deciduous species, especially cacti and bromeliads [13].
The precipitation regime tends to have high indices in the mid and lower course reducing in the river′s upper course. With an area of approximately 2900 km 2 , the latte has March as the wettest month, with an average annual precipitation between 600 and 800 mm, intermittent flow, and prominent flood peaks with medium recurrence. With approximately 2100 km 2 , the mid-course registers the most significant events between May and August, with an annual average rainfall of 1500 mm. With approximately 170 km 2 , the lower course registers an average yearly rainfall of 2000 mm. It is still possibl for rain, originating from atmospheric systems, to reach the entire basin. The large slope in the local topography also contribute to the high energy flows and the subsequen damage caused to riverside cities [14].
The municipality of Palmares (shown in Figure 1), our main damage analysis area, i within the Una River basin perimeter, a region with a comprehensive history of floods. I has an area of 339.3 km 2 and an estimated 63,500 inhabitants in 2020 [15]. Palmares stand out for its strategic location, close to capital cities and important inland cities. With a solid commercial and service sector, one of its main economic activities is the sugar industry. The municipality of Palmares (shown in Figure 1), our main damage analysis area, is within the Una River basin perimeter, a region with a comprehensive history of floods. It has an area of 339.3 km 2 and an estimated 63,500 inhabitants in 2020 [15]. Palmares stands out for its strategic location, close to capital cities and important inland cities. With a solid commercial and service sector, one of its main economic activities is the sugar industry.
In 2010, the Una River basin experienced the largest hydrological events ever recorded. On 17 June, a weather phenomenon called the "Easterly Wave" arrived in Pernambuco, intensified by the Atlantic Ocean temperature and strong trade winds. This combination resulted in a large concentration of rain in a short period and, consequently, in flood waves in the Una River basin. The precipitation of one day was equivalent to 70% of what was expected for the entire month. The flood event affected 67 cities in Pernambuco state and 19 in Alagoas state, with approximately 280,000 people affected and 79 deaths. In total, 102,420 people were displaced. In addition, 14,136 houses, 403 schools, 4478 km of roads, 142 bridges, two state hospitals, 85 health posts, and four municipal hospitals were damaged or destroyed. Furthermore, 27 municipalities declared a state of emergency and 17 declared calamity [16][17][18][19].
The 2010 flood event motivated the development of several preventive measures, such as the Serro Azul Dam (see Figure 1). Located in the municipality of Palmares, its construction dates back to 2017, and its drainage area covers 3299 km 2 . With a damming capacity of 303 million cubic meter, the dam can still regulate a flow of 850 L/s for other uses, including human supply, irrigation, leisure, and fishing [20].

Materials and Methods
Synthesis, a relevant element in scientific work, is associated with locating, sorting, evaluating, and combining information. It can usually be achieved through transparent, reproducible, and organized procedures for summarizing such information [21]. This study is comprised of hydrological and hydrodynamic simulations, damage estimation, and risk map generation. Figure 2 shows the steps required to achieve these goals. municipal hospitals were damaged or destroyed. Furthermore, 27 municipalities d a state of emergency and 17 declared calamity [16][17][18][19].
The 2010 flood event motivated the development of several preventive m such as the Serro Azul Dam (see Figure 1). Located in the municipality of Palm construction dates back to 2017, and its drainage area covers 3299 km 2 . With a d capacity of 303 million cubic meter, the dam can still regulate a flow of 850 L/s f uses, including human supply, irrigation, leisure, and fishing [20].

Materials and Methods
Synthesis, a relevant element in scientific work, is associated with locating, evaluating, and combining information. It can usually be achieved through tran reproducible, and organized procedures for summarizing such information [2 study is comprised of hydrological and hydrodynamic simulations, damage est and risk map generation. Figure 2 shows the steps required to achieve these goal

Hydrological and Hydrodynamic Modeling
Technological advances have expanded the possibilities regarding modelin are several models currently available. Li et al. (2019) [22] proposed a coupled called TPMF from the TOPMODEL and MIKEFLOOD models, in a study ap China. TOPMODEL is a semi-distributed top hydrological model and MIKEFLO hydrodynamic model that simulates 1D flow for channels and 2D flow for floodp this specific work, a 1D-2D coupled model was used to simulate flood events. software is the European LISFLOOD, a physical and distributed hydrologica created for flood prediction and for the analysis of climate change in transnationa Gai et al. (2019) [23] applied this model in a study developed in China. The highlighted LISFLOOD as an appropriate hydrological model for water re management planning at the basin scale.

Hydrological and Hydrodynamic Modeling
Technological advances have expanded the possibilities regarding modeling. There are several models currently available. Li et al. (2019) [22] proposed a coupled model called TPMF from the TOPMODEL and MIKEFLOOD models, in a study applied in China. TOPMODEL is a semi-distributed top hydrological model and MIKEFLOOD is a hydrodynamic model that simulates 1D flow for channels and 2D flow for floodplains. In this specific work, a 1D-2D coupled model was used to simulate flood events. Another software is the European LISFLOOD, a physical and distributed hydrological model created for flood prediction and for the analysis of climate change in transnational basins. Gai et al. (2019) [23] applied this model in a study developed in China. The authors highlighted LISFLOOD as an appropriate hydrological model for water resources management planning at the basin scale.
In this study, the floodplains were simulated using the Hydrologic Engineering Center (HEC) models from the US Army Corps of Engineers. HEC models are free, have wide applicability, constant updates and technical support. The Hydrologic Modeling System (HEC-HMS) has the versatility to set up a simulation scheme according to the region s characteristics and data availability. Other applications of these software can be found in Devi et al. (2019) [24], with projections of future scenarios, in Nharo et al. (2019) [25], with flood mapping even in the face of scarcity of data, and in Ribeiro Neto et al. (2015) [14], with analysis of the spatial and temporal variation of floods. As in Hadimlioglu et al. (2020) [26], a 2D model was sought to better represent the flow dynamics, and the progression of the flood, spatially and temporally. Another relevant aspect was the easy connectivity with Geographic Information Systems, both for pre-and post-processing. The model was set up for this application, as shown in Figure 3.
In this study, the floodplains were simulated using the Hydrologic Engineering Center (HEC) models from the US Army Corps of Engineers. HEC models are free, have wide applicability, constant updates and technical support. The Hydrologic Modeling System (HEC-HMS) has the versatility to set up a simulation scheme according to the region′s characteristics and data availability. Other applications of these software can be found in Devi et al. (2019) [24], with projections of future scenarios, in Nharo et al. (2019) [25], with flood mapping even in the face of scarcity of data, and in Ribeiro Neto et al. (2015) [14], with analysis of the spatial and temporal variation of floods. As in Hadimlioglu et al. (2020) [26], a 2D model was sought to better represent the flow dynamics, and the progression of the flood, spatially and temporally. Another relevant aspect was the easy connectivity with Geographic Information Systems, both for pre-and post-processing. The model was set up for this application, as shown in Figure 3. The river analysis system (HEC-RAS) model simulates the routing flow of rivers and canals. The scheme defined for this study considered the streamflow simulated with the HEC-HMS as input in the HEC-RAS model (upstream boundary condition), as shown in Figure 3, and the normal depth for the downstream boundary condition. The calculation of the normal depth requires the energy slope, but can be replaced by the bottom channel slope, as carried out in this study. The 2D version used in the simulations has an algorithm of solution based on finite volume that can handle different flow regimes with more accuracy, facilitating the modeling of steep rivers. This study used version 5.0.7 of the HEC-RAS, which calculates unsteady flow using the complete Saint Venant equations [27]. The 2D grid was created by defining a polygon that covered the entire region of interest (the urban area of the city of Palmares). Additional information helped in the delineation process, including the digital elevation model (DEM) and the image of the RAS Mapper (interface for geospatial data processing and visualization). The DEM had a spatial resolution of 1 m and accuracy of 25 cm.
The Serro Azul reservoir was represented in the HEC-HMS with the aim of simulating the effect of flood control. The geometry and dimensions of the hydraulic devices were informed in the HEC-HMS model, resulting in the output discharge calculated upstream of the Palmares with, and without, the reservoir. Details about parameter estimation and the calibration process of the models in the Una River basin have been described in previous studies [14,28]. The models satisfactorily represented the The river analysis system (HEC-RAS) model simulates the routing flow of rivers and canals. The scheme defined for this study considered the streamflow simulated with the HEC-HMS as input in the HEC-RAS model (upstream boundary condition), as shown in Figure 3, and the normal depth for the downstream boundary condition. The calculation of the normal depth requires the energy slope, but can be replaced by the bottom channel slope, as carried out in this study. The 2D version used in the simulations has an algorithm of solution based on finite volume that can handle different flow regimes with more accuracy, facilitating the modeling of steep rivers. This study used version 5.0.7 of the HEC-RAS, which calculates unsteady flow using the complete Saint Venant equations [27]. The 2D grid was created by defining a polygon that covered the entire region of interest (the urban area of the city of Palmares). Additional information helped in the delineation process, including the digital elevation model (DEM) and the image of the RAS Mapper (interface for geospatial data processing and visualization). The DEM had a spatial resolution of 1 m and accuracy of 25 cm.
The Serro Azul reservoir was represented in the HEC-HMS with the aim of simulating the effect of flood control. The geometry and dimensions of the hydraulic devices were informed in the HEC-HMS model, resulting in the output discharge calculated upstream of the Palmares with, and without, the reservoir. Details about parameter estimation and the calibration process of the models in the Una River basin have been described in previous studies [14,28]. The models satisfactorily represented the rainfall-runoff processes and the streamflow routing, providing the required information for the risk analysis.
The alternating block method [29] calculated the design storm hyetograph needed to simulate flood events for different return periods. The application of this method requires an intensity-duration-frequency (idf) curve, the specification of the return period, and the total duration of precipitation. Two idf curves, corresponding to two stations located in the municipalities of Palmares (Equation (1)) and São Benedito do Sul (Equation (2)), were used to represent extreme precipitation in the Una River basin [30]. The station at Palmares had a time series of 83 years  and São Benedito do Sul 25 years (1993Sul 25 years ( -2017. The maximum precipitation data was fitted to the Gumbel probability density function and its adherence was verified using the Kolmogorov-Smirnov test with significance level of 5%. Figure 1 shows the region of influence of each idf equation used to calculate the design storm in the sub-basins for six return periods (5, 10, 25, 50, 100, and 200 years). Finally, the duration of the hyetographs was defined as twice the concentration time of the sub-basin. The Kirpich method, based on the length of the longest river and the elevation difference between the farthest point of the longest river and its mouth, was used to estimate the concentration time.
where i is the intensity of the precipitation (mm·h −1 ), T is the return period in years, and t is the duration of the precipitation in minutes. The HEC-HMS simulation was compared to the maximum streamflow obtained by Dantas [17] using a daily time series of 39 years at the Palmares stream gauge station (shown in Figure 1). The annual maximum streamflow series fit the generalized extreme value distribution, which allowed the discharge calculation for the six return periods (5, 10, 25, 50, 100, and 200 years).

Damage Analysis
Vulnerability, intrinsically linked to the exposed community and its respective demographic conditions, is often classified as damage resulting from flooding, especially in the case of quantitative analysis [31]. If model validation is already difficult for a hazard analysis because of the lack of data, it is much harder for a damage analysis. The absence of validation data is a common problem, resulting mainly from high financial and human demands. However, a standardized, reliable, comparable, and consistent database is essential to study the behavior, variability, and, consequently, the relationship between damage and the characteristics of extreme events [10]. Thus, the initial idea was to develop criteria for estimating damage at a local scale, based mostly on secondary data.

Components of Risk
Hazard indicators characterize extreme events and influence the damage caused by floods. Depth, velocity, and energy are examples of indicators summarized by Kreibich et al. [32] and Ribeiro Neto et al. [33]. In this study, water depth and extent were used to estimate the damage associated with the structural aspects of residential buildings. Exposure and vulnerability are the other two components of risk. The demographic conditions of the community influence vulnerability, and it is estimated as a function of flood damage. Depth-damage curves are the most used method, where it is possible to assess damage according to the water depth reached in the flood event. When the return period of the event is known, the depth-damage curves can become probability-damage curves, which is the engineering flood risk [34].
Penning-Rowsell and Chatterton [35] proposed a damage classification for structural and content losses. Through discussions with contractors and architectural personnel, structural damages were defined, resulting in a relationship between water depth and loss as a percentage of a dwelling s value. The content losses were estimated based on the depreciation value of the furniture from current prices. Many studies conducted in several countries have adopted this concept [34,36,37]. For example, in Brazil, the relationship between water depth and the percentage of damaged buildings was estimated by Salgado [38]. Equation (3) summarizes the damage to buildings [39].
CRE-cost of damages to buildings CUB-basic unit cost of construction PED-percentage of damaged building AIC-inundated built area (m 2 ) The damage analysis was based on a standard single-family dwelling, adopting the Brazilian Standard NBR 12721/2005 and the CUB to characterize the buildings. The NBR 12721/2005 defines standard designs for different types of buildings (e.g., houses, buildings, and sheds) and uses (residential, commercial, and service). The CUB calculation method, which considers the cost per square meter for each type and use, is provided monthly by the unions of civil construction industry in all Brazilian states. The values used in this study are presented in Table 1. Some elements used to categorize the dwellings are the composition of the building (rooms and their quantity), number of floors, specifications of finishing, and building area. The Supplementary Material describes the classes in Table 1 according to NBR 12721/2005. In Equation (3), the CUB is adjusted by a physical depreciation factor of 50%.
The PED is determined for intervals of water depth and calculated based on the standard designs and an index of quantification of the damage extension supported by specialists in civil construction [39]. The damages are calculated for depths greater than 0.50 m because, in general, buildings have a base above the curb. Table 2 lists the PED values as a function of the standard designs and classes of depth. The inundated built area (AIC) was calculated by intersecting the vectorized lots with the flooded area calculated by the HEC-RAS. The analysis included an estimation of inventory, that is, the contents of a building. The value used was from a report on losses and damages resulting from the 2010 event prepared by the World Bank with support from the Pernambuco government [40]. The value was corrected for January 2021, considering the inflation of the period, and corresponded to the complete furniture of the housing unit. Assuming that the estimation of the World Bank considers an intermediate residential class, the values of the items refer to the standard design R1-N. Other personal items not listed were considered to add 15% to the final deal.

Buildings Dataset
All lots comprising the urban area affected by floods in Palmares were vectorized manually, taking as base map a high-resolution image from an aerial survey with spatial resolution of 20 cm and accuracy of 25 cm. In ArcGIS/ESRI, besides the calculation of the area of the lots, the vectors were associated with information (street name, building number, type, use, number of floors, roof area, and standard design), creating a georeferenced database. In addition, Google Street View, an application that provides panoramic views on the ground level, allowed the visualization of the frontage of the buildings. In Palmares, panorama photos were taken in 2012, two years before the high-resolution images.
According to Brazilian regulations, identifying each house building pattern is essential for damage estimation. Thus, the first parameter used was the built area of each lot (m 2 ) calculated in a GIS environment. Subsequently, field visits were carried out, and a second classifying element was surveyed for validation, which was the identification and number of rooms in each visited residence. This information was then used to define the thresholds of area of the lots corresponding to each standard design. Field visits were conducted between July and August 2019. With 150 interviews in total, the sample had a confidence level of 90% and a margin of error of 6.46%. Data collection was performed with the support of the Android DataScope app, a collaborative platform that allows the creation and sharing of forms, in addition to offline data acquisition. The app also allowed the taking of photographs and exporting of the results in a CSV file.
The questionnaires aimed to validate residential pattern classes. We considered selecting houses to visit, keeping a balanced number of samples for each category. However, difficulties that were encountered, such as residents absence or refusal to answer, led us to select houses randomly. In addition, a different area of the city was covered in each visit to represent differences in urban morphology.

Flood Simulation
The first tests with the HEC-HMS had streamflow overestimation at the Palmares station, probably due to the use of the same return period for all sub-basins. The alternative was to apply the tool flow ratio available in the HEC-HMS. The flow ratio allows the user to control the flow in the sub-basins, defining a factor that regulates the hydrograph, by multiplying the factor by each ordinate.
The main characteristics of the hydrodynamic model were a computational grid with cells measuring 40 m, a simulation time step of 1 min, an output time step of 30 min, and a Manning roughness coefficient of 0.035. After some tests, the best value was found for the flow ratio (0.25), resulting in the discharges listed in Table 3 (with and without application of the flow ratio). The HEC-RAS simulation tended to produce a more significant discharge reduction because of the flow routing processes represented in the hydrodynamic model. Simulations for different conditions were performed to test the sensitivity of 2D modeling in HEC-RAS. The analyzed parameters were the size of the 2D mesh cells, the Manning coefficient, and the use of break lines. The latter is a model device that forces the mesh cells to align and control the flow direction. Break lines can be used along any higher ground, levees, roads, or even the entire length of the banks of the main channel [41]. The values initially adopted for the mesh and the roughness coefficient were obtained from previous studies for the same region [14,18,42] [43].
Incorporating the Serro Azul reservoir in simulations produced discharge reductions that varied from 45.4% to 55% in HEC-HMS and from 43.5% to 54.6% in HEC-RAS. This reduction was directly proportional to the return periods, with the 200-year simulation presenting the most significant reduction percentage. Figure 4 illustrates the evolution of the water extent, considering all return periods with, and without, the effect of the Serro Azul reservoir.

Damage Estimation
Based on the limits of the flood extension registered in previous events (e.g., 2000 2010, and 2017), the blocks and lots situated inside the inundation perimeter wer vectorized (see Figure 5). There were 2733 buildings, of which 1852 were residential, 24 commercial, 185 mixed, 173 services, 86 abandoned, and 190 non-identified. Based on th field information, the buildings were classified according to Table S1 and then compare to the roof area, establishing the relation proposed in Table 4.

Damage Estimation
Based on the limits of the flood extension registered in previous events (e.g., 2000, 2010, and 2017), the blocks and lots situated inside the inundation perimeter were vectorized (see Figure 5). There were 2733 buildings, of which 1852 were residential, 247 commercial, 185 mixed, 173 services, 86 abandoned, and 190 non-identified. Based on the field information, the buildings were classified according to Table S1 and then compared to the roof area, establishing the relation proposed in Table 4.  For the city of Palmares, the time difference between the surveyed images provided by Google Street View (2012) and the scene captured in the field survey (2019) did not show significant differences, with the largest difference being the painting of frontage ( Figure 6). The standard design classification was based on the roof area of vectorized lots and the list of rooms in each residence. According to Brazilian regulations, these are also two of the most prominent elements in the differentiation of each pattern [44]. Standards also use finishing items; however, they showed less relevance for classification and were more inaccessible.  For the city of Palmares, the time difference between the surveyed images provided by Google Street View (2012) and the scene captured in the field survey (2019) did not show significant differences, with the largest difference being the painting of frontage ( Figure 6). The standard design classification was based on the roof area of vectorized lots and the list of rooms in each residence. According to Brazilian regulations, these are also two of the most prominent elements in the differentiation of each pattern [44]. Standards also use finishing items; however, they showed less relevance for classification and were more inaccessible. The built area was vital for classification, as it was easily affordable information, gained mainly through remote sensing. Among all the parameters present in the dataset, the built-up area stood out as the most relevant in the damage calculation. Google Street View, in turn, proved to be a very viable tool for carrying out analysis of this type, The built area was vital for classification, as it was easily affordable information, gained mainly through remote sensing. Among all the parameters present in the dataset, the builtup area stood out as the most relevant in the damage calculation. Google Street View, in turn, proved to be a very viable tool for carrying out analysis of this type, especially in the context of scarcity of financial, human, and data resources. Arrighi et al. [9] reinforced this idea, highlighting the surface area of the affected buildings as the first metric in the exposure analysis, in the case of their research, which was also carried out in a GIS environment at the lot scale. The number of floors is another parameter mentioned. Google Street View was the tool chosen to gather information and fill in the attribute table, which proved extremely effective.
Approximately 45% of the vectorized dwellings were in class R1-B, more basic homes with two bedrooms, a living room, a bathroom, a kitchen, and an area for a tank. Almost 34% of the houses were in the R1-N class, 12% in the RP1Q class, only 5% in the R1-A-class, and 4% were unidentified. There was also a predominance of type R1-A and R1-N residences in flatter parts of the city, planned neighborhoods with greater infrastructure and space. Standard R1-B and RP1Q dwellings, on the other hand, were primarily in areas with greater density, where the city design is more organic, such as on the banks of the river.
The World Bank [40] study estimated the damages caused by the 2010 flood event, which totaled BRL 3.4 billion in Pernambuco state (4% of the state s GDP, 60% corresponding to direct costs, and 40% to indirect costs). Residential buildings registered a cost of BRL 10,250.00 per unit. Considering the 1852 residential buildings vectorized, the damage due to dwellings would be BRL 18,983,000 according to the World Bank [40] information. Applying Equations (1) and (2), the thresholds of Table 4, and the values of CUB, the damage estimation was BRL 15,292,572.56. These values demonstrated coherence in the developed methodology.

Risk Curve for Estimation of Benefits
The analysis of flood control measures must consider the probability of events and avoid overestimating the benefits [34]. The probability of an event with the respective damage allows building damage-probability curves (also called risk curves) to be developed, whose integral is the expected annual damage. The sequence of graphs in Figure 7 illustrates the damage estimation for each return period at Palmares. Given the streamflow, the depth was estimated ( Figure 7a); then, for the same depth, the damage was calculated using Equations (1) and (2) (Figure 7b). The graph in Figure 7c shows the calculated probability associated with the streamflow, and the combination of the probability (Figure 7c) and the damage (Figure 7b) results in the relation shown in Figure 7d. The benefits of the control measures could be estimated by subtracting the areas under the curves in Figure 7d.
In the graphs of Figure 7, the solid line represents the scenario before the construction of the Serro Azul dam. In contrast, the dashed line corresponds to the conditions after the insertion of the dam. Before the dam, the calculated water depth values ranged from 6.42 to 8.78 m with a flow varying from 535 to 1606 m 3 /s. The depths lowered after installing the dam, varying between 5.05 and 7.10 m, with flows from 302 to 730 m 3 /s. The construction of the Serro Azul reservoir led to a significant damage reduction according to the simulations of the return periods. Table 5 lists the damage caused by the return period, considering simulations with, and without, the Serro Azul reservoir. The total reduction, including structure and inventory, was 88.0%, 86.7%, 76.6%, 67.5%, 59.4%, and 51.7% for 5, 10, 25, 50, 100, and 200 years, respectively, revealing an inverse relationship between the percentage of reduction and the frequency of flood events. The annual benefits given by the Serro Azul reservoir would be BRL 3.5 million, considering only the residential sector. The global benefit certainly increased when considering other cities located downstream and other sectors, such as services and commerce. streamflow, the depth was estimated (Figure 7a); then, for the same depth, the dam was calculated using Equations (1) and (2) (Figure 7b). The graph in Figure 7c shows calculated probability associated with the streamflow, and the combination of probability (Figure 7c) and the damage (Figure 7b) results in the relation shown in Fig  7d. The benefits of the control measures could be estimated by subtracting the areas un the curves in Figure 7d. In the graphs of Figure 7, the solid line represents the scenario before the construc of the Serro Azul dam. In contrast, the dashed line corresponds to the conditions after insertion of the dam. Before the dam, the calculated water depth values ranged from to 8.78 m with a flow varying from 535 to 1606 m³/s. The depths lowered after instal the dam, varying between 5.05 and 7.10 m, with flows from 302 to 730 m³/s. The construction of the Serro Azul reservoir led to a significant damage reduc according to the simulations of the return periods. Table 5 lists the damage caused by return period, considering simulations with, and without, the Serro Azul reservoir. total reduction, including structure and inventory, was 88.0%, 86.7%, 76.6%, 67.5%, 59 and 51.7% for 5, 10, 25, 50, 100, and 200 years, respectively, revealing an inv relationship between the percentage of reduction and the frequency of flood events. annual benefits given by the Serro Azul reservoir would be BRL 3.5 million, conside only the residential sector. The global benefit certainly increased when considering o cities located downstream and other sectors, such as services and commerce.

Maps of Risk Elements
The most significant hazard was concentrated in the channel and its immediate vicinity in all simulated scenarios, as in Amarnath et al. [45]. Similar to Andrade and Szlafsztein [46], exposure was presented here as flood extents, identified by the elements within the events scopes ( Figure 8). The classes of damage in Figure 8 were used for the quartiles of the monetary values: low (0-25%), medium (25-50%), high (50-70%), and very high (75-100%). Damage analysis and the exposure of buildings for each event were evaluated through overlapping floodplains and built elements, as in Tincu et al. [47]. For the events of 5, 10, 25, 50, 100, and 200 years, 210, 343, 758, 1080, 1286, and 1431 exposed buildings were identified, respectively. In engineering, risk can be interpreted by studying different events, magnitudes, and respective damages. Thus, the values forming the constructed probability-damage curves were presented in a spatially distributed manner.
When considering the social dimension of vulnerability, the characteristics related to the society under study cannot be neglected. In the case of residences, for example, income and family education level are some of the factors that will influence their quality and, consequently, their response to extreme phenomena. Poorer families are likely to have more fragile housing than those with a better income (GOERL et al., 2012) [48]. This scenario was confirmed in the study area.
high (75-100%). Damage analysis and the exposure of buildings for each event were evaluated through overlapping floodplains and built elements, as in Tincu et al. [47]. For the events of 5, 10, 25, 50, 100, and 200 years, 210, 343, 758, 1080, 1286, and 1431 exposed buildings were identified, respectively. In engineering, risk can be interpreted by studying different events, magnitudes, and respective damages. Thus, the values forming the constructed probability-damage curves were presented in a spatially distributed manner. When considering the social dimension of vulnerability, the characteristics related to the society under study cannot be neglected. In the case of residences, for example, income and family education level are some of the factors that will influence their quality and, consequently, their response to extreme phenomena. Poorer families are likely to have more fragile housing than those with a better income (GOERL et al., 2012) [48]. This scenario was confirmed in the study area.
When analyzing indicators from the Brazilian Institute of Geography and Statistics (IBGE) Demographic Census (household income, age and total number of residents in private and collective households) there was a clear difference amongst areas, almost defining two very distinct sectors, one further to the south and the other to the east. In the field visits, it was possible to have a greater understanding of the local dynamics. Of the 150 residents interviewed, the vast majority have lived in Palmares for many years. When asked about the damage caused by floods, the answer was almost unanimous, everything was lost, repeatedly, especially at the big event in 2010. In this latter, besides all personal property, the damage caused to homes was diverse. Damaged ceilings, floors, finishes and doors. Covers and walls washed away by the water. Lots of mud, dirt and stench left behind. When analyzing indicators from the Brazilian Institute of Geography and Statistics (IBGE) Demographic Census (household income, age and total number of residents in private and collective households) there was a clear difference amongst areas, almost defining two very distinct sectors, one further to the south and the other to the east. In the field visits, it was possible to have a greater understanding of the local dynamics. Of the 150 residents interviewed, the vast majority have lived in Palmares for many years. When asked about the damage caused by floods, the answer was almost unanimous, everything was lost, repeatedly, especially at the big event in 2010. In this latter, besides all personal property, the damage caused to homes was diverse. Damaged ceilings, floors, finishes and doors. Covers and walls washed away by the water. Lots of mud, dirt and stench left behind.
With great unevenness in the terrain, Palmares still preserves the traces of its creation, high density on the banks of the Una River, but already showing expansion in areas of higher elevation. Even with the construction of the Serro Azul dam, and its proven importance for damage reduction, its combination with other measures, such as the removal of buildings located on river banks, as well as the association of structural and non-structural measures, led to the most favorable conditions of risk management, as highlighted by Amarnath et al. (2015) [45].

Discussion
The evaluation of the exposure (water extent) and hazard (water depth) for each return period was successful, showing the efficiency of model application (hydrologic and hydrodynamic) as being better than other procedures adopted under limited data conditions [47]. The relationship between the magnitudes of the events and their respective consequences was developed, enabling the estimation of the average annual damage after integration of the area below the flood risk curve or damage-probability curve. The damage reduction percentage due to the flood control installed in the Una River (between 51% and 88%) was close to that obtained by Shrestha and Kawasaki [49] (40-60%). Through the methodology of this research, the annual benefits provided by the Serro Azul dam are estimated at BRL 3.5 million per year for the residential sector, not counting isolated events likely to happen within this period. When considering other municipalities downstream of the dam and other sectors (service and commerce), there would be an increase in this number.
Fadel et al. [34] discussed how apparent non-action is the least costly measure. Any intervention will always have an initial cost, and benefits will only be provided later. The authors explained that other flood events are likely to occur in the medium and long term. As a result, the cost of accumulated damage will increase, easily overcoming a scenario in which no measures are implemented. In addition to being designed to protect the cities downstream from flood events, with the end of work on its pipeline, the Serro Azul dam will have expanded use, acting to supply municipalities in the Agreste region of Pernambuco. In Palmares, besides the effects of the Serro Azul, other measures were adopted, such as removing several buildings located on river banks, which could result in more benefits that were not considered in this study. The combination of structural and non-structural measures is the most favorable scenario for risk management in the scientific literature on this topic.
Finally, regarding the performance of risk analysis, it is important to highlight that uncertainty is inherent and uncertainties are present at significant levels throughout the process. The reduction of uncertainties depends on the existence of more works on the theme developed here. Thinking beyond the observed data enables comparisons to be made regarding the performance of risk curves. Although this is a current topic, it is not a new topic, and even so, it still has many gaps and few works actually trying to reach the end of the risk analysis. Collaborative initiatives have gained increasing prominence on this topic, with concepts such as collaborative mapping, crowdsourcing, voluntary geographic information and citizen science. In this research, for example, there were many sources of uncertainties, such as topographical data, flow characteristics, or model parameters. It is extremely important for decision-makers to take this limitation into account when developing public policies. A great effort to reduce uncertainties is found in validating hydrologic, hydrodynamic, and damage modeling. Validation techniques used here included comparing the results simulated by the model and observed data, comparing the performance of different models, and analyzing results by experts in the field. However, it is reaffirmed that validation still faces several challenges, such as scarcity, non-standardization, or non-availability of data in sufficient quantity and quality. Thus, damage models are often transferred from one region to another without adaptation, which is not recommended. The introduction of specific local data contributes to reducing these uncertainties. The hope is that the multiplication of studies on the subject can help create a database for comparison and validation of results [10].

Conclusions
This study was based on local data, applying damage estimation at the lot scale. The virtual classification of buildings allowed coverage of practically the entire constructed territory, without excessive demand for financial and human resources. The lot scale allowed for a local analysis of buildings, which was information not present in the IBGE Census data which covers social demographic information. We sought to minimize generalizations regarding urban morphology.
The registration of the exposed elements was performed satisfactorily. Building classification was completed, with the roof area being the main parameter in the methodology, where the use of affordable and secondary data was prioritized. The classification of residential standard designs, an essential step in the damage analysis, was possible due to the information obtained in the field survey, contributing to the understanding of the diversity present in loco. The results obtained were validated owing to the availability of damage estimation for the 2010 flood event.
Damage estimation methods allow for the consideration of different scenarios and contribute to risk analysis. In monetary value, this study calculated the effect of the construction of a reservoir for damage reduction, showing the potential to determine the effectiveness of measures adopted to mitigate flood impacts. In addition, the maps of exposure, hazard, and damage were elaborated for each simulated return period, making possible a complete risk analysis.
The understanding of risk is established and strengthened through quality information, accessible analysis methodologies, data spatialization, decision support, and the establishment of communication between the different participants. With its continental territory, Brazil has already been a site for several disastrous flood events. The search for combining methodologies and creating more robust risk maps for local regions aims to reduce costs and damages generated in the public and private sectors, especially to the population. It is important to highlight the possibility of applying the methods used in this research to other regions and cities, as almost all the data can be obtained secondarily.