Hydraulic Planning in Insular Urban Territories: The Case of Madeira Island—São Vicente

This study aims to examine the flood propensity of the main watercourse of São Vicente drainage basin and, if relevant, to propose two methodologies to alleviate the impacts, i.e., detention basin sizing and riverbed roughness coefficient adjustment. Geomorphological data were obtained from the watershed characterization process and used through the SIG ArcGIS software for the flood propensity assessment and then for the calculation of the expected peak flow rate for a return period of 100 years through the Gumbel Distribution. Subsequently, the drainage capacity of the river mouth was verified using the Manning-Strickler equation, in order to establish whether the river mouth of the watershed has the capacity to drain the entire volume of rainwater in a severe flood event. In summary, it was possible to conclude that São Vicente’s watershed river mouth is not able to completely drain the rain flow for the established return period. Thus, its drainage capacity was guaranteed by modifying the walls and streambed roughness coefficient and by sizing the detention basin using the Dutch and the Simplified Triangular Hydrograph methods.


Introduction
Global warming with its increasing variability leads to an increased risk of both floods and drought [1]. While increases in temperatures are in question for all seasons of the year, precipitation may decrease in one season while increasing in the other season. According to some strong findings, this variability in precipitation will increase even more in the future [2]. This, in turn, will lead to a discontinuity of precipitation throughout the year and, consequently, to an increase in sudden and torrential rains. The most natural consequence of this situation is that floods occur in many basins. Floods, depending on the size of the flow in the surrounding area, affect settlements and agriculture by damaging their areas, lower and upper structures, facilities, and living things, and they interrupt human life and socio-economic activities. Sociological effects on humans from floods, psychological disorders, and the like, are also seriously affected.
The source of water, in addition to determining the amount of water falling on the surface in terms of vegetation characteristics, influences the amount of water on the floor, underground infiltration of ground in terms of soil properties, plant and underground water leaking from the amount of the residual flow to the understanding of the causes of floods and geomorphology relationship is extremely important. There are also incorrect  Likewise, the island capital's watersheds-i.e., Funchal, history has recurrently shown that São Vicente's basin suffers with major flooding events, like the ones that took place in 2010 and 2013, resulting in both civil and assets catastrophic losses. São Vicente also suffers from anthropogenic pressure, much like any other urban municipality, particularly felt by a considerable soil sealing index as a result of an urban sprawl over a semi-rural area [10,11]. Furthermore, its main watershed river mouth is covered in abundant vegetation and a vast track of sediments, that considerably reduce the drainage capacity of this water channel, as shown in Figure 2.  The state of conservation of the stream is virtually the same throughout its length within the urban perimeter, which can be clearly confirmed in situ. Abundant vegetation and sediments can be found along the watercourse, mostly due to the low slope of the streambed that reduces the flowrate's velocity and subsequently the ability to expel sediments through the watershed's river mouth. The state of conservation of the stream is virtually the same throughout its length within the urban perimeter, which can be clearly confirmed in situ. Abundant vegetation and sediments can be found along the watercourse, mostly due to the low slope of the streambed that reduces the flowrate's velocity and subsequently the ability to expel sediments through the watershed's river mouth.

Schematic of the Methodology
The methodology adopted can be summarized in 6 phases, as illustrated in Figure 3.

Schematic of the Methodology
The methodology adopted can be summarized in 6 phases, as illustrated in Figure 3. This case study followed a methodology with six major stages with the first being a deep literature review with the purpose of getting the much needed hydrological and morphometric analysis of São Vicente's watershed. To present a valid contribution for the study of flood-prone watercourses, different methodologies were also taken into consideration suggested by prestigious authors, thus ensuring the liability of this study's outcome. The following steps are further explained above. This case study followed a methodology with six major stages with the first being a deep literature review with the purpose of getting the much needed hydrological and morphometric analysis of São Vicente's watershed. To present a valid contribution for the study of flood-prone watercourses, different methodologies were also taken into consideration suggested by prestigious authors, thus ensuring the liability of this study's outcome. The following steps are further explained above.

Morphometric Characterization of the Watershed
The key parameters used for the morphometric characterization of a watershed are [4,5,[12][13][14][15][16][17]: • Gravelius Index-K C : This parameter ratio between the perimeter of the watershed with a hypothetical circular watershed with the same area, telling how close to a circular shape the watershed is [13,16]. The Gravelius index can therefore be calculated using Equation (1), being a dimensionless parameter and characterizing the watershed as flood-prone the closer it is to 1 [16].
where: P = Perimeter of the watershed, km; A = Watershed area, km 2 .
• Elongation Factor-K L : The elongation factor determines the ratio between the watershed shape and a rectangle with the same area and can be calculated using Equation (2). It is dimensionless and characterizes the watershed as elongated if the result is higher than 2 [16]. where: L E = Equivalent length, km; l E = Equivalent width, km; K C = Gravelius Index, dimensionless; A = Watershed area, km 2 .
• Shape Factor-K F : The Shape Factor describes the ratio between the average width of the watershed with its length. It can be calculated using Equation (3) and is also dimensionless, indicating the basin's elongation degree. The higher is the Shape Factor, the less flood-prone the watershed is and more elongated. where: A = Watershed area, km 2 ; L B = Length of the watershed, km.
The watershed's length can be found measuring the distance between the farthest point to the watershed's river mouth. The watershed's length must not be confused with its main watercourse length, as the last is commonly bigger due to the stream's sinuosity. Using the MDE file provided by the Regional Civil Engineering Laboratory of the Autonomous Region of Madeira (LREC-RAM), it was possible to morphometrically characterize São Vicente's watershed and therefore its main watercourse. The morphometrical data gathered in this study was then processed using different author's methodologies with their very own specific restraints.
A morphometric analysis consists of establishing a hierarchy between the watershed's watercourses-i.e., Strahler or Shreve's hierarchies-according with their order or magnitude [17]. Both types of hierarchy can be conducted following a hydrological analysis of the DEM file, thus obtaining the "flow accumulation" and the "flow direction" rasters using the "flow order" tool [15]. Nevertheless, the Strahler's hierarchy is deeply associated to a given watershed's bifurcation ratio, with the many degrees of bifurcation being calculated using the Equation (4) [7,12,[14][15][16][17][18]. where: N i = Number of watercourses classified as "i"; dimensionless; N i+1 = Number of watercourses classified as "i + 1", dimensionless.
This dimensionless coefficient is obtained by the ratio of the number of watercourses of a given order by the number of watercourses of the immediately higher order. The average bifurcation value can be calculated based on Equation (5). where: N i = Number of watercourses classified as "i"; dimensionless; N i+1 = Number of watercourses classified as "i + 1", dimensionless; N 1 = Number of first-order watercourses.
The bifurcation ratio is also a dimensionless parameter as it merely represents an arithmetic mean of all bifurcation ratios. Moreover, the concentration time of a watershed reveals itself as a key-factor towards the good morphometric characterization of a watershed. It determines the required amount of time needed for all the watershed area contribute to the process of rainfall drainage until it finally crosses the river mouth [12,14,17,18].
Since the equations used to determine the concentration time of a watershed are empirical, each methodology presents different values for the same parameters, and therefore it is advised to use an arithmetic mean of them all, avoiding extremes. In this study, it was used to gather values from the Kirpich (Equation (6)), Témez (Equation (7)), and Giandotti (Equation (8)) methodologies [15].  where: t C = Time of concentration, hours; L = Length of the main watercourse, km; i = Slope of the main watercourse, m/m.
where: t C = Time of concentration, hours; A = Watershed area, km 2 ; L = Length of the main watercourse, km; H M = Average height of the watershed, m.

Precipitation Analysis
The hydrological study is performed based on a probabilistic analysis of extreme highintensity and short-duration events that took place in São Vicente's watershed throughout history. This data was obtained using National Weather Resources Information System (SNIRH), which also publishes precipitation data recorded automatically in many different stations through the island. Regarding the implemented probabilistic methodology, the Gumbel's Distribution was adopted as it seemed more appropriate to process the obtained data and meet Madeira Island's watersheds projections [19]. Hence, the maximum annual daily precipitation can be calculated using Equation (9).
where: P EST = Estimated annual maximum daily precipitation, mm; P M = Average annual precipitation, mm; S = Sample standard deviation, mm; K T = Frequency Factor, dimensionless. where: where: X i = Sample value, mm; X M = Sample mean, mm; n = Number of samples. where: T R = Return period, years.
After establishing the daily precipitation for an extreme phenomenon, the precipitation intensity with a particular duration can be obtained using Equation (12).
where: I = Precipitation intensity, mm/h; P EST = Estimated annual maximum daily precipitation, mm; t C = Time of concentration, hours; k = Time distribution coefficient, dimensionless. where: where: t C = Time of concentration, hours. The time distribution coefficient is a primary parameter since the annual maximum daily precipitation is only valid for events lasting 24 h. Thus, as the duration of precipitation is equal to the time of concentration of the watershed, using the total amount of daily precipitation in the hydrologic analysis would lead to oversized hydraulic structures [15,20].

Drainage Capacity of the River Mouth and Peak Flow Rate
A river mouth's drainage capacity can be calculated using the Manning-Strickler equation (Equation (14)) and compared to the expected flow for an extreme event with a return period of 100 years. On the other hand, the expected flow for these type of events can also be calculated using the commonly known and used Forti (Equation (16)); Rational (Equation (17)); Giandotti (Equation (18)); and Mockus' (Equation (19)). equations. where: Q M = Drainage capacity of the river mouth, m 3 /s. A M = Area of the river mouth cross-section, m 2 ; R = Hydraulic radius, m; i = Average slope of the river mouth region, m/m; n = Roughness coefficient of the riverbed and walls, m −1/3 s, Table A1. where: where: B = Width of the river mouth runoff section, m; h = Height of the river mouth runoff section, m; A M = Area of the river mouth cross-section, m 2 .
The width and height of the stream in the region of the mouth were obtained through previous studies in the region [15], and the first parameter was confirmed via the georeferencing process. where: where: Q Rational = Peak flow rate by the rational methodology, m 3 /s; C = Surface runoff coefficient, Table A2; I = Precipitation intensity, mm/h; A = Watershed area, km 2 .
One of the most important design criteria for hydraulic infrastructures is to determine a Fill Rate value lower than 85%, therefore considering a safety margin thus ensuring the safety of the population and their assets [15,21]. Additionally, to control the river mouth's flowrate, it is also required to implement runoff restraining mechanisms, namely spillways.
As mentioned before, the Fill Rate value can be calculated using Equation (20) and if the river mouth's runoff capacity proves to be insufficient to drain the rain flow of the given watershed and respect the proposed safety margin, then it must be designed a flood mitigation mechanism, such as a detention basin.
where: FR = Fill Rate, %; Q P = Peak flow rate of each methodology, m 3 /s; Q M = Drainage capacity of the river mouth, m 3 /s.
The Fill Rate parameter refers to the ability of a drainage section to drain a particular flow. Therefore, if the Fill Rate value is greater than 100%, the section is no longer able to drain the total volume of water without overflowing [15].

Detention Basin Sizing
As stated before, if the river mouth cross section is insufficient to drain all the rain flow collected along the watershed, it must be designed some type of spillway to restrain the flowrate, keeping it under the downstream expected limit. In this study, it was adopted a Cipolleti type of spillway as this is a type has features that reduce turbulence regions of contact with the water, thus making it easier to drain all the stream's runoff [5,6,20]. The design of this type of spillways can be performed using the Equation (21).
Knowing the amount of flowrate that must be drained by the stream's river mouth, it is also possible to estimate the volume of water that needs to and will be retained by the detention basin. In this sense, two different methodologies were adopted, namely the Dutch Method (Equation (22)) and the Simplified Triangular Hydrograph (STH; Equation (23)).
where: Q S = Flow drained by spillway, m 3 /s; L SD = Width of the sill, m 3 /s; H D = Height of the waterline above the sill, m. where: V A = Storage Volume, m 3 ; Q P = Peak flow rate of each methodology, m 3 /s; Q S = Flow drained by the spillway, m 3 /s; t C = Concentration time, hours.
Based on the STH geometric analysis ( Figure A1), the Equation (23) was then formulated, considering an event with a duration of at least twice the concentration time of the given watershed. This, took into consideration that the last rain particle to reach the river mouth came from the farthest region and that it also would take place at the last instant of precipitation, indicating that it would need to be equal to the time of concentration to be considered as drained by the river mouth [15].
These methodologies were selected since the Dutch Method does not consider the delay and damping of the precipitation hydrograph, ultimately resulting in the overdesign of the infrastructure [22], as illustrated in Figure 4, where Q s : states for the runoff capacity of the spillway; t c : meaning the concentration time; t MAX : being the maximum precipitation duration (standard); t d : being the time delay until the beginning of water accumulation in the detention basin; H a,MAX : being the maximum storage capacity; and i(t MAX ): meaning the precipitation intensity for the maximum duration.
where: QS = Flow drained by spillway, m 3 /s; LSD = Width of the sill, m 3 /s; HD = Height of the waterline above the sill, m.
where: VA = Storage Volume, m 3 ; QP = Peak flow rate of each methodology, m 3 /s; QS = Flow drained by the spillway, m 3 /s; tC = Concentration time, hours. Based on the STH geometric analysis (Figure A1), the Equation (23) was then formulated, considering an event with a duration of at least twice the concentration time of the given watershed. This, took into consideration that the last rain particle to reach the river mouth came from the farthest region and that it also would take place at the last instant of precipitation, indicating that it would need to be equal to the time of concentration to be considered as drained by the river mouth [15].
These methodologies were selected since the Dutch Method does not consider the delay and damping of the precipitation hydrograph, ultimately resulting in the overdesign of the infrastructure [22], as illustrated in Figure 4, where qs: states for the runoff capacity of the spillway; tc: meaning the concentration time; tMAX: being the maximum precipitation duration (standard); td: being the time delay until the beginning of water accumulation in the detention basin; Ha,MAX: being the maximum storage capacity; and i(tMAX): meaning the precipitation intensity for the maximum duration. Thus, it was confirmed that, in the Dutch Method, storage begins at the same time with precipitation, which does not correspond to reality as storage will only begin when the flow drained downstream is greater than the spillway's runoff capacity.

Modification of the Roughness Coefficient
Another structural measure that was taken into consideration was the modification of the roughness coefficient of the walls and streambed of the main watercourse, as this would avoid friction between the water and the channel, therefore increasing its drainage Thus, it was confirmed that, in the Dutch Method, storage begins at the same time with precipitation, which does not correspond to reality as storage will only begin when the flow drained downstream is greater than the spillway's runoff capacity.

Modification of the Roughness Coefficient
Another structural measure that was taken into consideration was the modification of the roughness coefficient of the walls and streambed of the main watercourse, as this would avoid friction between the water and the channel, therefore increasing its drainage capacity. This methodology is based on changing the value of the parameter "n" in the Manning-Strickler equation to improve the flowrate of a given watercourse by considering another type of material or at least its conservation status for its walls and streambed [15].

Results
The results shown here correspond to the data obtained by applying the formulas already described. Thus, to evaluate the morphometric features of the main watercourse of São Vicente, an individual analysis of each parameter listed in Table 1 was conducted, correlating them with reference values proposed in various bibliographies.
The first parameter refers to the watershed's area, a key-factor for the volume of water that needs to be drained through the river mouth. It can be classified as: Very Large > 20 km 2 ; Large > 10 km 2 ; Medium > 1 km 2 and Small < 1 km 2 [23]. This way, as illustrated in Table 1, the considered watershed can be classified as "Very Large", and more flood prone compared to smaller watersheds. Nevertheless, the standard values are arbitrary and may vary according to the analysis conducted throughout the study [23] and with the flood-prone character of the basin. As shown in Figure 5, the São Vicente watershed's borders are considerably higher than the central area, indicating a steep slope and subsequently a very fast supply of the main watercourse, thus increasing the volume of water flowing through the stream that will end in its river mouth. Regarding this watershed's singularities, Figure 6 shows a higher number of streams which also suggests that it has a higher drainage capacity-there are many low and medium order tributary streams that supply the main watercourse. In addition, the drainage system is an index that translates to the hydrographic tendency that a watershed has to create new streams. Thus, basins with higher hydric densities tend to present more tributary streams, this happening as a consequence of the ability to generate new streams [15,17].
Once again, this analysis was only possible thanks to the data available on the National Information System on Water Resources (SNIRH) [24], assessing data samples from a considered period of sixteen years as shown in Table A4 and Figure A2. Therefore, using the Gumbel Distribution's probabilistic process, it was possible to obtain the values shown in Table 2. Regarding this watershed's singularities, Figure 6 shows a higher number of streams which also suggests that it has a higher drainage capacity-there are many low and medium order tributary streams that supply the main watercourse. In addition, the drainage system is an index that translates to the hydrographic tendency that a watershed has to create new streams. Thus, basins with higher hydric densities tend to present more tributary streams, this happening as a consequence of the ability to generate new streams [15,17].   Table 4. Surface drainage coefficient adopted (Source: [25]).

Urban Areas Land Occupation Surface Drainage Coefficient
Commercial Area City Center 0.700-0.950 Peripheral Areas 0.500-0.700 The value of the reduction coefficient (λ) used in the calculation of the flow through Giandotti's methodology is presented in Table 5. Once again, this analysis was only possible thanks to the data available on the National Information System on Water Resources (SNIRH) [24], assessing data samples from a considered period of sixteen years as shown in Table A4 and Figure A2. Therefore, using the Gumbel Distribution's probabilistic process, it was possible to obtain the values shown in Table 2. After obtaining the precipitation intensity index estimated for a return period of 100 years, it was calculated peak flowrates, as presented in Table 3, using the aforementioned methodologies and equations. The surface runoff coefficient particularly used in the rational methodology was 0.500 (Table 4) since the area under study is a peripheral region with commercial buildings. This value translates essentially to the ratio of water that tends to run on the stream's surface, i.e., 50% of the total precipitation.  Table 4. Surface drainage coefficient adopted (Source: [25]).

Land Occupation Surface Drainage Coefficient
Commercial Area City Center 0.700-0.950 Peripheral Areas 0.500-0.700 The value of the reduction coefficient (λ) used in the calculation of the flow through Giandotti's methodology is presented in Table 5. Table 5. Adopted Giandotti's reduction coefficient (Source: [26]).

Area (km 2 )
λ Equivalent "C" <300 0.346 1.250 Regarding the drainage capacity of the considered watershed's river mouth, the Manning-Strickler equation ultimately confirmed the need to further implement a structural flood mitigation measure like the detention basin, where the values obtained are summarized in Table 6. It must also be noted that the walls and the streambed have different roughness coefficients. Consequently, the drainage capacity of the river mouth was calculated through the weighted arithmetic mean of the corresponding coefficients, considering that the stone and mortar walls are in good condition (n = 0.020) and the streambed is made of a rocky surface with abundant vegetation in poor condition (n = 0.040). Another key-factor that must be taken into consideration is the very low slope of the river mouth, tending to reduce the flowrate velocity and consequently the drainage capacity of the cross section. To simulate and model a critical situation, it was then considered a slope of 0.01 m/m in the reference section. As presented in Table 6, the Fill Rate is higher than the established limit of 69% for the Rational, Giandotti, and Mockus methods, yet again clearly indicating the need to implement mitigation and flowrate control measures for the river mouth section. Considering this, it was designed a detention basin with the flowrates calculated of the methodologies, affected by the spatial restrains and the anthropogenic pressures of an urban area, namely the already existing infrastructures nearby the watercourse.
Since the design of detention basins depend on the exceeding limits of the flowrate for the watershed's river mouth, a Cipolleti's trapezoid spillway was also designed to restrain and control the flowrate that will end downstream. This type of spillway's features can be found in Table 7.
Afterwards, the detention basins were also designed through the Dutch and the STH Methods, which are merely simplified methodologies that do not take into consideration many key-factors and consequently result in the overdesign of the hydraulic infrastructure. Also, the width and heigh of the detention basin cross section were both fixed with same values of the existing one to reduce this measure's environmental impacts on an urban area.
Thereupon, the only geometric parameter that may vary was its length, being limited to the main watercourse length. Using all the methodologies and both methods made possible to present the following length values in Table 8. Finally, changing the roughness coefficient of the streambed and walls was also considered as an alternative flood mitigation structural measure towards preventing its impacts whilst maintaining the same amount and features of the streambed vegetation. This way, Table 9 values were particularly chosen to clearly improve the conservation status of the streambed, thus reducing friction between the drained water and the material covering the watercourse, subsequently increasing its drainage capacity. At last, the modified roughness coefficients of the stream walls were considered to have its surface covered in concrete in a good condition status, although the streambed would remain with the same rocky and abundant vegetation features, nevertheless in good condition. The values used for these coefficients are summed in Table 10.

Discussion
As this study's main goal was to check if it was needed to put into action flood mitigation measures to further prevent major impacts in São Vicente's watershed, the use of a detention basin revealed itself as valid and useful structural measure towards controlling its river mouth's flowrate [27]. At first, the Fill Rate was 98%, 135% and 99%, respectively for Rational, Giandotti, and Mocku's methodologies, which ultimately decreased to only 61% after adopting the detention basin measure. This structural measure's outcome is clear evidence that it may enable the river mouth to work below 85% of its full capacity. Moreover, this proves the accuracy of the Regional Directorate for Territorial Management and Environment (DROTA) prediction, as presented in Table 11. Table 11. Watersheds with high flood risk. (Source: [28]).

São Vicente São Vicente
This study aimed to cause the least possible impact over the considered watercourse and its surroundings since it is believed that the presence of natural elements in cities present itself and act as a vital condition for the environmental recovery of the urban territory [29]. Additionally, a nature and urban systems symbiosis is typically found as a key-factor or goal to further achieve a territory or city's sustainability [30,31]. Nevertheless, uncontrolled urban sprawl is something that can take place especially in rural areas, thus creating urban voids [32].
As it was not made any change to the stream's cross section, namely its height and width, the only variable parameter was its length. It was based on this concept that the Dutch Method presented an abnormal oversize of the detention basin's length when compared to the watershed's main course's length. Therefore, according to this method it would be needed to change one or both cross section dimensions and so it cannot be considered valid for the aforementioned urban design settings.
The exact same conditions were imposed for the STH method, with it showing a different and this time valid approach since the detention basin's length was shorter than the watershed's main watercourse length.
As for the change of the stream bed and walls roughness coefficient, it was decided to remain with the abundant vegetation and sediments along the watercourse but improving its conservation status by performing a correct maintenance as this would outcome in a cheaper process with less wasting of time and resources. Also, there's no need to perform maintenance on the stream walls frequently, since the mechanical abrasion only happens in alluvial events that tend to result in the drainage of higher volumes of water and large sediments.
Despite being a relatively simple structural measure, the change of the roughness coefficient of this stream, resulted in a significant way, enabling its river mouth to work below the Fill Rate limit, that itself, emphasizing that both methodologies-i.e., the STH method and the changing of the stream's roughness coefficient-can be implemented together, to optimize and reduce the required detention basin's length.
As a final remark, it should be noted that the methodologies were simplified and therefore do not consider local peculiarities. Thus, this may result in oversized hydraulic infrastructures because of conservative considerations and inputs.

Conclusions
This study revealed how flood-prone São Vicente's watershed is in the event of extreme rainfall occurrence, as it was already predicted by DROTA's own Flood Risk Report. The watershed's drainage capacity is highly decreased by the presence of abundant vegetation and a huge number of sediments throughout the watercourse, resulting in a lower flowrate in an already low-slope stream and river mouth. The insufficient drainage capacity of the river mouth was verified through 3 of 4 methodologies used in this study namely: Rational, Giandotti, and Mocku's.
Regarding the two methods used during this study, the Dutch Method did not present coherent results as it indicates the need for very long detention basins in relation to the watershed's main watercourse. On the other hand, the Simplified Hydrograph Method presented not only satisfactory results but also to be easier to implement as there is no need for change either the stream's height or width.
Even though changing the watershed's streambed and walls roughness coefficient may seem a relatively simple and unworthy measure to consider, it surely proved to mitigate the flood's impact, fulfilling its main goal by preserving infrastructures and people's assets.
Afterall, this study leaves a clear open-door to others that may complement its contents and methodologies by optimizing its techniques. To improve the often-complex urban hydraulic system and demand, it is also expected that new studies take notice of the need to reduce sediment deposition as it seems to make a huge long-term impact over the watercourse's drainage capacity and ultimately to prevent a major flood impact [33]. On the other hand, mechanical abrasion of this stream's walls and the amount of time that often takes to local public authorities to perform any type of maintenance have been two strong reasons for how degraded the main course tributaries are and subsequently by the lower water quality discharged [34,35] and therefore also need to be studied and improved. Furthermore, the urban growth ratio projected for the municipality of São Vicente and how it may impact soil waterproofing and ultimately surface run-off should be a top concern and studied, complementing this work's outcome.
Generally, this study enhances the methodologies and techniques used in similar case studies as valid and appropriate towards scientific development based on flood scenarios modelling and simulations [36,37].

Data Availability Statement:
The data presented in this study are openly available. Also, it is possible to contact one of the study authors.