Flood Risk Assessment Index for Urban Mobility with the Aid of Quasi-2d Flood Model Applied to an Industrial Park in S ã o Paulo, Brazil

: Risk can be deﬁned as the relationship between the likelihood of a hazard causing a potential disaster and its consequences. This study aims to assess the likelihood that a new industrial region, located in the state of S ã o Paulo (Brazil), will be ﬂooded, causing the disruption of the mobility system and local economic activities. To fulﬁll this aim, a new approach is proposed by combining the vector information of the highway network that serves the region with the result of a quasi 2-D raster ﬂood model, generating a set of interpreting rules for classifying the safety of routes. The model called MODCEL is a quasi-2D hydrodynamic model that represents the watershed using compartments called cells, and it was adapted to work using a raster ﬁle format in which each pixel is represented as a ﬂow cell connected to its surroundings by the Saint-Venant equations without the inertia terms. Therefore, this study proposes an assessment framework that can be replicated for similar problems of ﬂood risks to mobility. The possible effects of ﬂood events on the accessibility to areas of interest are determined, indicating a possible disruption to economic activities and transportation and allowing for planning alternatives in advance.


Introduction
There is a strong correlation between urban development, flooding, and flood protection measures. Protection works have been triggered by development, but also development has been boosted by the existence of protection works [1]. This means that the sense of safety introduced by flood defense works tends to favor urban development in originally hazardous areas.
This may be a controversial trend, as it can increase the residual risks as a consequence. Therefore, flood risks may increase, in fact, even in "protected areas", since the exposure of goods and people can be intensified by the false perception of complete safeness [2]. That is, the territory put in safe conditions by hydraulic works is only safe under a certain threshold, since every project is defined for a certain reference event related to a determined return period.
In fact, this trend can be seen in many countries, where a large number of cities have been built or developed in disaster-prone regions [3]. Hence, recognizing the range of possible threatening impacts of floods is important for successfully establishing an effective flood risk management process [4].
It is noteworthy that flooding may cause serious problems in the cities' infrastructure systems and, as a consequence, affect essential products and service provisions. In the state of São Paulo, where the largest Brazilian industrial park is located, many companies suffer serious financial and competitiveness losses due to excessive rainfall and floods. According to the Research Department of the Industrial Federation of São Paulo State (FIESP) [5], 59% of the state industries are negatively affected by flooding, which result in delays in the supply chain, due to failure in the transportation service, increasing the information demand for planning and studies to mitigate losses.
Therefore, the contribution of this study lays on a proposal of a flood risk assessment framework, based on hydrodynamic modeling, that aims to evaluate the potential damage in the mobility system of a new industrial park located in the state of São Paulo (Brazil). To fulfill this aim, a new approach is proposed by combining the vector information of the highway network that serves the region with the result of a quasi 2-D raster flood model, generating a set of interpreting rules for classifying the safety of routes.
This work will help decision-makers assess the practicability of creating alternative routes, or to increase its flood defense, to find an optimal strategy that significantly reduces flood-induced delays.
The proposed framework, presented in more detail in the next section, was applied to a case study, but it can be replicated for similar problems of flood risks to mobility in any other region.

Risk Assessment Definitions
Risk, hazard, and related definitions are found in several documents, such as [6,7]. For the United Nations [8], risk is constituted by the probability of having negative consequences or expected losses (to the individual property, goods, interruption of economic activities, and environmental damage), resulting from interactions between natural or man-induced hazards and the conditions of the vulnerability/capacity of the system.
Hazard, instead, refers to "a dangerous phenomenon, substance, human activity or condition that may cause loss of life, injury or other health impacts, property damage, loss of livelihoods and services, social and economic disruption, or environmental damage" [8]. Natural hazards usually cannot be controlled or reduced. The risk, in turn, is liable to be managed by minimizing its impacts on the interest system [9].
Vulnerability can be understood, in a general way, as the degree to which a system is likely to experience harm, due to exposure to a hazard [10]. At last, resilience refers to the way in which the population and organizations manage to respond positively to the adverse conditions of an event [11].
In urban flood studies, hazard is triggered by the occurrence of intense rainfalls that interact with the watershed to produce discharges and, eventually, flooding. Then, it will be a function of the probability of flood occurrence, while its consequences will depend on the flood properties, such as flood depths, duration, and flow velocity. On the other side, local vulnerability will define how critical these flood properties are, depending on physical, social, economic, and environmental factors, which can intensify the susceptibility of a given community to the impacts of a given hazard [12].

Urban Transport Disruption
Critical infrastructures provide essential products and services to cities, including transportation, drinking water, water treatment, energy, telecommunications, information technology, food, banking and finance, public health and healthcare, commercial facilities, and emergency services [13][14][15].
It is noteworthy that these infrastructure systems, in modern cities, are connected directly or indirectly through one or more intervening components. Therefore, disruption or failure of one component can propagate through the network and subsequently affect the other components [14,15].
In this sense, disruptions to such systems can cause drastic socioeconomic consequences and impede the development of sustainable cities [16]. Therefore, reliable and secure operations of infrastructure systems are indispensable to national well-being and development [15].
In this sense, extensive research has been conducted to evaluate the impacts of flooding on communities and economies [17]. However, fewer studies have been focused on the vulnerability of the infrastructure system to flooding events and the consequences on the corresponding services [15,18,19].
Regarding the weather-related disruption of the transport sector, flooding is considered its prevalent cause [20,21]. When flooding closes a certain road, not only are the cars directly caught by water affected, but a greater part of the transportation network may be affected by the consequent traffic jam. As the transportation network enables the movement of people and goods and is the basis of economic activities, its disruption may have drastic socio-economics consequences [16]. According to the World Bank [22], approximately 73% of global expected annual damage (EAD) is caused by surface and river floods. Additionally, the average global EAD for transport infrastructure assets is $14.6 billion and about $8 billion in low and middle-income countries.
The functionality of a transportation network is based on many factors, but, in theory, if vehicles traveling on roads can drive near the roads' speed limit, there would be no notable congestion, except at peak traffic hours. However, in the case of flood events, some roads may have diminished their level of service, resulting in a decreased transportation-network effectiveness [23].
In this sense, some researchers have developed methodologies for flood risk assessment, considering the transport system disruption. It is known that flooding on a road does not necessarily preclude people from driving along it [20]. Therefore, in order to transition from a binary view of a flooded road being considered 'open' or 'closed', a more complete approach considering the relationship between the flood depths and the disruption of the road transport system can be established. Usually, this relationship considers the flood depth and the impact on the vehicle's speed and on the road network connectivity [20,24,25].
In the [20] study, the function that describes the relationship between the depth of standing water and vehicle speed was constructed by fitting a curve to video analysis supplemented by a range of quantitative data that has been extracted from existing studies and other safety literature. The main advantage of this methodology is that it is simple to incorporate this proposed function into existing transport models to produce better estimates of flood induced delays.
Another interesting study was carried out by [15]. The authors proposed a flood risk assessment framework for interdependent infrastructure systems based on the strongest path method (SPM) and 2D hydrodynamic modelling. The model investigates the impacts of direct and indirect connections among infrastructure systems, explores their relative importance in the network, and prioritizes the associated risks for more efficient resource allocation.
In the study carried out by [26] the vulnerability is determined by a susceptibleimpacted-susceptible (SIS) diffusion model to capture the impact of flooding on the road network connectivity. With model results, the authors were able to estimate a relationship between vehicle speed and flood depth on road networks. Some other authors have developed this same relationship vehicle speed and flood depth but using hydrodynamic models. In the [25] study, the authors calculated the flooding level in urban areas, caused by rainfall, by using the spatial runoff assessment tool (S-RAT) and the flood inundation model (FLO-2D model). Using modelling results, they combined a flood-vehicle Speed curve and a rainfall-flood depth curve to prepare a traffic disruption map.
Another similar approach was undertaken by [24], where the flood depth was modeled using a static and hydrodynamic model and then assigned to each road segment. After that, flood depth was turned into one of the factors that causes an increase in travel time.
Also using hydrodynamic modelling, ref. [27] proposed an integrated framework linking meteorological information, land use functions, and hydrodynamic models with a safety speed function to relate flood depth to a reduction in speed in order to determine road network vulnerability.
Considering the consulted references, this study intends to present a flood risk assessment framework that also correlates the flood depths with the disruption of the road transport system. The proposed methodology is based on the combination of vector information of the highway network that serves the region with the result of a quasi 2-D raster flood model, generating a set of interpreting rules for classifying the safety of routes.

Materials and Methods
Flood risks to mobility assessment must consider the available routes, land use, occupation patterns, and flood characteristics, establishing reference indicators to assess lock down conditions [28]. In this context, a specific flood risk assessment framework for urban mobility was developed.
The proposed framework (as seen in Figure 1) was elaborated in four stages: (1) the identification of previous studies, (2) the vulnerability to flood assessment, (3) the flood hazard assessment, and (4) the flood risk assessment. More details regarding each phase are presented in the following.

Study Area
The state of São Paulo, southeast region of Brazil, has a total area of 248,193.1 km 2 (~2.9% of the entire country) and approximately 43.7 million inhabitants (~22% of Brazilian population). Its urbanization level is about 96.37%, which shows its high urban expressiveness [29].
This state also has the highest Brazilian GDP, corresponding to 32.5% of the national GDP. The two most significant economic activities undertaken inside the state are services (77.02% of GDP) and industries (21.43%) [29].
Considering the importance of industries in state economics, this research presents a case study that consists of the evaluation of alternative entry and exit routes of a new industrial park located in the state of São Paulo, as shown in Figure 2.

Phase 1: Previous Studies
The first phase of the proposed methodology consisted of raising information of the studies carried out at three levels of government: federal, state, and municipal.
Previous studies were gathered to better understand the region where the interest site is located and to identify the most vulnerable areas to flood events. These studies will help to assess the site's vulnerability to flooding.

Phase 2: The Vulnerability to Floods Assessment
Streets and highways connect regions and contribute to their economic development. Flood events not only injury the traffic but can also damage the structures of bridges and crossings and affect economic activities. This situation particularly affects industrial areas where interruptions in the pathways can lead these sites to have severe supply chain crises, isolation, and even affect their operations.
In this sense, the vulnerability to the flood indicator should point to how vulnerable to floods the mobility systems and the interest site itself are. In this study, the flood vulnerability assessment follows a general qualitative classification, as shown in Table 1. The flood vulnerability assessment was made according to the following steps: 1.
Indicating the possible routes between the site and the nearest highway.

2.
Highlighting the existing bridges and culverts on each route.

3.
Defining the watershed area of each river crossing (bridges and culverts).

4.
Calculating the hydrological vulnerability of each river crossing based on the normalization of its watershed area with values between 0 and 1, with a higher value being more vulnerable (Equation (1)). where: wa-watershed areas of each river crossing (m 2 ); hv-hydrological vulnerability of a river crossing; f norm -area normalization function. The area equivalent to the density value for the third quartile of the entire sample is taken as a reference value. Any watershed with an area greater than this reference value is at the maximum normalized value, that is equal to 1. For watersheds with areas smaller than this reference, the normalized value is based on a linear distribution. This procedure was made to reduce possible "flattening" of the evaluation scale, due to the presence of extreme values.

5.
Computing the route vulnerability. The route vulnerability is then given by the highest value among two options (OP 1 and OP 2). This configuration aims to eliminate the potential bias of routes with several river crossings, but only a few of them with high vulnerability, in a situation that could lead to a false low value of vulnerability. The adopted options were: • OP 1: the route vulnerability is given by the simple average of the hydrological vulnerability of all existing river crossings on the route (Equation (2)); where: RV OP 1 -route vulnerability option 1; hv-hydrological vulnerability of existing river crossings on the route; n-number of all of the existing river crossings on the route.
• OP 2: the route vulnerability is given by the normalization of the average watershed area of all river crossings on its way. Firstly, the average area of the watershed of all river crossings on each route is taken (Equation (3)). Then this average area value is normalized with the same normalization function previously presented (Equation (4)). where: AWA-average of watershed areas of all river crossings on the route (m 2 ); wa-watershed areas of each river crossing on the route (m 2 ); n-number of all of the existing river crossings on the route. RV OP 2 -route vulnerability option 2;

Phase 3: The Flood Hazard Assessment
The flood hazard assessment is based on the maximum water depths resulting from a hydrodynamic simulation of intense rainfall, with 100 and 500 years of return period (H 100 e H 500 , respectively), over the watershed. The hazard is then classified in a qualitative range, according to Table 2. This classification is adapted from [20], who proposed a depth break function for road transport. To map flood depths, hydrologic and hydrodynamic mathematical modeling tools were used, as shown in the following.
The hydrological study aims to elaborate the design rainfall events, related to the chosen return periods, as inputs to the hydrodynamic studies.
The duration of rainfall is inversely proportional to rainfall intensity. Thus, the shorter the rainfall duration, the greater its intensity. Short-term rainfall tends to cause higher runoff peaks in the drainage network, while longer rainfalls result in higher water volumes in the system. To evaluate the worst scenario, a design rainfall was taken with a duration equivalent to the time of concentration of the watershed, highlighting the effect of higher flood peak values.
The necessary hydrological information refers to the historical rainfall series measured in the region, which enables the elaboration of IDF (intensity-duration-frequency) equations to local rainfall. The design rainfall events, with a critical duration for the watershed and related to 100 and 500 years of return period, were elaborated based on these equations.
The resulting hydrographs, obtained by the application of the synthetic unit hydrograph [31], were used as boundary conditions for hydrodynamic modeling.
For the hydrodynamic studies, the modeling tool called MODCEL was used as a flow-cellular model for urban watershed representation. MODCEL is originally a quasi-2D model [32] that represents the urban space using a set of homogeneous compartments, called cells. The concept of flow cells was initially developed by [33] and enshrined by [34].
The cells cover the whole space of the watershed in an interconnected way, forming a flow network, linked by one-dimensional equations. Thus, MODCEL can represent the twodimensional characteristics of the watershed in a simplified way, only using 1D equations. MODCEL is also integrated with a hydrological module that performs the rainfall-runoff transformation in each cell, using a modified rational method and considering a storage capacity in each cell.
The variation of the water volume in a cell i, in a time interval t, is given by the mass balance in this cell. Thus, in differential terms, we have the following continuity equation (Equation (5)): where: Q i,k -discharge between neighboring cells i and k; Z i -water level in the center of the cell i; A si -the surface area of the water mirror in the cell i; P i -discharge produced by the rainfall occurring on the cell i; t-independent time variable.
In the MODCEL, the temporal discretization of the differential equation (Equation (5)) is made by numerically linearizing all nonlinear terms, so that there is no need for an iterative solution procedure in order to simplify the numerical model. The discharges are expanded in the Taylor series. The cell links are defined by the modeler according to the characteristics observed in the topographic and urban fabric analysis of the region to be modelled.
MODCEL is in continuous improvement, and it is used as an alternative for modeling complex urban areas. MODCEL can be found in detail in [35] and some of its applications can be found in [36][37][38][39][40].
In this study, however, MODCEL was used as a raster model to represent the region of interest in an automatic way. Raster models are those that transform each point of a digital elevation model (DEM) or digital terrain model (DTM) in a "pixel". In this case, the structure of MODCEL was adapted so that each pixel represents a flow cell that is connected to the other cells in its surroundings by the Saint-Venant dynamic equation, without the inertia terms, as represented schematically in Figure 3. This type of model has increasingly gained credibility for floodplain modeling [41][42][43][44][45][46][47][48][49][50]. The main advantage of these models is the processing speed [51]. However, it is necessary to find the equilibrium between the best resolution of the modelled area and the computational cost [52], without disregarding the cost of surveying.

Phase 4: The Risk Assessment
The flood risk assessment indicates how prone to suffering the impacts from floods are on the mobility systems and the site itself, which is estimated through the product between the vulnerability and the hazard. The general equation used to calculate the flood risk is given by Equation (6): where: R-flood risk; V-flood vulnerability; H 100 -Flood hazard for the 100-year storm. Equivalent to the maximum flood depth resulting from the 100-year storm; H 500 -Flood hazard for the 500-year storm. Equivalent to the maximum flood depth resulting from the 500-year storm.
The choice of a product allows for cancelling the final risk if one of the parcels is equal to zero. Thus, if there is no vulnerability, there is no risk, because there is no exposure, or yet, if there is no hazard, there is no risk, because there is no impact cause.
The weight system applied to the hazard indicators (H 100 e H 500 ) considered the probability of the occurrence of each event. The flooding event resulting from the 100-year storm has a probability of occurrence equal to 1% in a given year, while one from the 500-year storm has a probability equal to 0.2%. Thereby, the 100-year storm has a chance of occurrence that is five times greater than the 500-year storm. This result was used to define the weight of the indicators, rounded to two decimal places.
The equation is weighted to give more importance to the hazard that characterizes the potential impact of the flooding event. This decision was the result of a modeler decision and can be revised in each particular application.
Finally, the flood risk classification received a qualitative range, as shown in Table 3. The flood risk assessment of a mobility system considers the flood vulnerability of each potential route, accounting for the hazard indicator encountered along the route path.
To make it possible, it detached all flood hazard values from flow cells that are superimposed by route path. Each cell receives a risk value (Equation (6)) considering the flood depths in the cell and the vulnerability of the route. Then, the risk along the route is recalculated by the moving average of three consecutive cells (~180 m). This definition is proposed to guarantee that there is a systematic flooded reach, and not only a small local spot in one cell/pixel. In this way, the maximum risk will be achieved only when the route presents a continuous stretch with at least 180 m (probably more) in length with a maximum flood risk (R = 1). This configuration also indirectly measures the recession time of flooding, once long stretches of flooded areas tend to have a slower flood recession.

Phase 1: Previous Studies
The information obtained from previous studies allows an overview about flooding characteristics on the industrial park surrounding areas. Studies were carried out at three levels of government: federal, state, and municipal.
The region of the industrial park site does not show any riverine paths prone to floods, accordingly to ANA (National Water Agency) [53], but a CPRM (Geological Survey of Brazil) study [54] indicates several susceptible areas inside the municipality territory. The municipality administration indicates a large flooding zone, and a part of this flooding zone reaches the industrial park neighborhood. Figure 4 shows the existing flood risk information at the industrial park region, consolidated in a single map composed of flood vulnerability [54] and susceptibility [53], intending to offer a comparison with the obtained results in the present study. The flooded area mapped in the drainage management plan of the municipality is also plotted in the same figure. Generally, from the information found in the existing studies, it is possible to notice that the industrial park is not directly exposed to common flood events. The site is located on a plateau, several meters above the valleys of the neighboring rivers. However, flood events have a high chance of indirectly affecting the site of interest, blocking transportation or disrupting public services.

Phase 2: The Flood Vulnerability Assessment
The vulnerability assessment of the industrial park aimed to classify its propensity to suffer critical impacts from flood events, which could result in the disruption of mobility systems.
As mentioned before, the flood vulnerability assessment is given by the number of accesses to the site, considering the closer highway, as well as the hydrological vulnerability of the various river crossings existing in the path.
The methodology adopted for the flood vulnerability assessment starts by indicating the possible routes between the site and the nearest highway, and also identifying the existing bridges and culverts on each route.
The main accesses to the industrial park are highways A and B, both located on the west side of the site. The access to these highways can be achieved through four routes: Route 1, Route 2, Route 3, all of them with three river crossings, and Route 4 with six river crossings. Figure 5 shows all four routes and their river crossings, called T2 until T13. In sequence, the watershed area for each river crossing must be determined. Table 4 shows, in its second column, the watershed areas of all river crossings and the refence value up to maximum vulnerability, taken as equal to the third quartile of the watershed area distribution. The function used in the area normalization and to calculate the hydrological vulnerability of river crossings is show in Figure 5. Once calculated and once the hydrological vulnerability of each river crossing is identified, it is possible to calculate the route vulnerability. Its value is given by the highest value among two options (OP 1 and OP 2) previously described.
Finally, Table 4 and Figure 6 also show the vulnerability evaluation of potential routes of the industrial park site.    Upstream areas were modeled using a hydrologic model, with the objective of estimating input hydrographs as boundary conditions to the hydrodynamic modeling domain. The hydrodynamic modeling domain covers different types of land use and cover, which were analyzed to compose a weighted runoff coefficient, based on [55]. The final CN coefficient for the North and South watersheds is shown in Table 6. In Figure 7B, the two boundary conditions, the North and South watersheds, are represented, as well as the existing types of land use of the watershed.

Design Storm
The design rainfall was estimated through the intensity-duration-frequency relationship. The time of concentration was estimated through the equation of George Ribeiro [56], defined as follows (Equation (7)): where: L-Length of main watershed valley (km); S-Average slope of the basin (M/m); p-Coefficient of vegetation cover of the basin.
The time of concentration of the watershed taking the industrial park site as reference is estimated at 1114.14 min. To use a critical rainfall for the modeled area, a design rainfall of 1440 min, equivalent to the time of concentration of the watershed, was calculated, allowing for the evaluation of a critical flood peak.
From this rainfall duration, the intensity and precipitation height were calculated for the simulation scenario, using two rainfall equations (IDF) for two rain gauges. The first IDF was the equation from the DAEE (Department of Water and Electric Energy of São Paulo) [57] for the rain gauge 1, defined by: where: i-rainfall intensity, corresponding to the duration t and recurrence time T, in mm/min; T-recurrence time in years; t-duration of rainfall in minutes.
The Thiessen polygon method was used to define the influence areas of each considered gauge, as shown in Figure 8. Table 7 shows the total height for both rainfall design events in both rain gauges.  The watersheds of the boundary conditions are entirely located within the area of influence of the rain gauge 1.

Hydrodynamic Studies
The modeled domain in the studied region is 216 km 2 and was discretized in a grid of 59,985 cells with 3600 m 2 each. The map of Figure 9 shows the runoff used in the hydrodynamic studies for the domain area where the industrial park lies. The manning and the runoff coefficient used for each land cover are presented in Table 8. The simulation results of the flood depths for both events of RP100 years and RP500 years are presented in Figure 10.  Although the location of the industrial park has no direct contact with flooding waters, the region has a very critical situation. The main river valley localized close to the A and B highways suffers from severe flood depths.
Additionally, the flood discharges of rivers at bridges and culverts were evaluated to verify potential impacts on drainage structures. These points are shown in Figure 11 Table 9 shows the result of the hydraulic evaluation of river crossings in the industrial park region, highlighting the estimated flood discharges and hydraulic capacity of the drainage structures. The scale of input data gives this bias to simulated flood discharges, but it does not affect the flood mapping. The evaluation of the flood peak must be carried out with more detailed data.

The Flood Hazard Assessment
The results of the flood hazard evaluation considering the storm with return periods of 100 years and 500 years are presented in Figure 12.
The areas with very high hazard classifications are concentrated along the valleys of rivers, showing two critical zones: between the A and B highways and in the basin outfall region. This last region presents a large and occupied area impacted by flooding that has high and very high hazard indicators.
As it was expected, the area closest to the industrial park site is not prone to a flood hazard, but the general situation of the watershed, within the modeling domain, shows large areas with high and very high hazard classification, that can affect the normal operation of transportation and public services.

Phase 4: The Risk Assessment
In the case of this study, the flood risk expresses the probability of an area to become flooded and cause difficulties in the access to the industrial park.
For the evaluation of the flood risk of the mobility system, the risk of each route was calculated considering the relation between the hazard indicators encountered along the route path and the route vulnerability (see the explanation in item "4.2. Phase 2: The Flood Vulnerability Assessment" and "4.3. Phase 3: The Flood Hazard Assessment").
The industrial park was classified as having a very high flood risk with a final average classification of 0.86 (in a scale from 0.00 to 1.00), as shown in Table 10.

Discussions and Conclusions
Floods, especially the ones resulting from intense precipitation, are the predominant cause of weather-related disruption of the transport sector. As the transportation network enables the movement of people and goods and is the basis of economic activity, its disruption may have drastic socio-economics consequences.
Since extreme weather events can cause direct or indirect damage to the transport sector, many researchers have developed methodologies for flood risk assessment, considering the transport system disruption. In most cases, when considering the impact of flooding on transport disruption, it is assumed that a road is either 'open' or 'closed'. However, flooding on a road does not necessarily preclude people from driving along it.
In this sense, the aim of this paper was to present a framework, based in hydrodynamic modeling, to establish a relationship between flood depths and the disruption of the road transport system. This new approach combines the vector information of the highway network that serves the region with the result of a quasi 2-D raster flood model, generating a set of interpreting rules for classifying the safety of routes.
The proposed methodology was used to assess the likelihood that a new industrial region, located in the state of São Paulo (Brazil), could be flooded, causing the disruption of the mobility system and local economic activities. It can be particularly interesting for decision-makers to plan alternative paths and strategies in advance to reduce flood-induced delays.
The proposed framework was elaborated in four phases. The first phase was the identification of previous studies that helped to better understand the region where the site is located, identify the most flood prone areas, and create the bases to assess flood vulnerability.
The second phase is the flood vulnerability assessment. Considering that the area is well above the quota of the valleys of the rivers nearby, the vulnerability and flood risk analysis was based on the possibility of the disruption of the mobility system. In this study, the vulnerability assessment considered the number of available alternative routes and the river crossings on their way. The results of the flood vulnerability assessment show that all four routes that connect the industrial park to the most important highways in the region are classified as high or very high, regarding flood vulnerability. This situation highlights the transport infrastructure vulnerability.
The third phase is the flood hazard assessment. In this study, the hazard analysis is based on the result of hydrological studies and hydrodynamic simulation of intense rainfall over the watershed. The floods resulting from hydrological events with 100 and 500 years of recurrence time were evaluated and the extent and depth of flooding in each event were determined. To understand these hydraulic phenomena in a more systemic way and to make more accurate estimates of their consequences, mathematical hydrodynamic models must be used. For the hydrodynamic studies, the modeling tool called MODCEL, a quasi-2D model [32], was used as a raster model to represent the region of interest.
The fourth and final phase is the risk analysis. The flood risk assessment is a combined evaluation of the vulnerability and hazard analyses that aims to determine the possible effects of flood events on the mobility system of the industrial park, indicating the dis-ruption of its operation. The result of risk analysis showed that the mobility system is a critical item.
The proposed methodology is adherent with some recent research carried out by many authors, as it considered that there is a relationship between flood depth and a decrease in traffic velocity [15,20,[24][25][26][27]. The obtained results indicate that the proposed methodology can achieve reliable and sound results.
However, it is noteworthy that this framework is highly dependent on physical and flood information during the elaboration and calibration of the model. These data, and especially the ones related to the river crossings main characteristics, are not easily found information. Despite this, the proposed framework proved to be satisfactory in evaluating the possible effects of flood events on the accessibility of the industrial park, indicating the possible disruption to economic activities and transportation, and allowing for planning alternatives in advance.
The results also showed that, to reduce the vulnerability of the industrial park, it is necessary to create alternative routes or to increase its flood defense. However, these mitigating measures may be infeasible, as they go beyond the private sphere, and are very costly. It is important to emphasize that characterizing the possible safe places to direct routes to is a step to be thought about calmly and confirmed with the existing and/or created flood spots.
Another relevant point is the computational performance and the modeling method adopted. It is observed that the size and detail of the area to be modeled can make modeling unfeasible. However, this modeling approach was capable of providing an equilibrium between the best resolution of the modelled area and the computational cost, allowing the creation of a relationship between the flood depths and the disruption of the road transport system.
Finally, it is important to notice that this framework was applied in a case study but can be replicated for similar problems of flood risks to mobility.