Hydraulic Planning in Insular Urban Territories: The Case of Madeira Island—São João Stream, Funchal, Portugal

: This study’s primary goal was to conduct an analysis regarding the ﬂ ood susceptibility of the main watercourse of the São João (Funchal) drainage basin. In addition, if proven necessary, we also aimed to suggest mitigation measures, such as sizing a detention basin and promoting adjustments of the riverbed’s roughness coe ﬃ cient. This study also resorted to geomorphological data—obtained during the watershed characterization process—that were then utilized in the SIG ArcGIS software, in order to estimate the expected peak ﬂ ow rate, considering a return period of 100 years using the Gumbel distribution. Finally, the Manning–Strickler equation was utilized to determine the river discharge point’s drainage capacity; the reason for that was to verify whether its drainage capacity was su ﬃ cient to drain the entire volume of rainwater associated with an extreme ﬂ ood event. In summary, the results obtained by this study indicate that the drainage capacity of the river discharge point of the São João watershed (Funchal) is insu ﬃ cient when considering an extreme ﬂ ood event, for a return period of 100 years. Hence, it became necessary to explore the two aforementioned mitigation measures: ﬁ rst, regarding the detention basin, its sizing was calculated through both the Dutch method and the simpli ﬁ ed triangular hydrograph method; second, aiming to increase the drainage capacity of the river discharge point, it is suggested that the roughness coe ﬃ cient should also be modi ﬁ ed.


Introduction
It has become a consensus that climate change is directly connected with the increase in the periodicity of natural disasters. The effects related to climate change are projected to intensify over the course of the 21st century and even beyond, with all of the potential associated negative impacts [1]. The frequency of extreme weather events-e.g., strong rainstorms that cause floods across the world-has been increasing, mostly due to climate change. Because of the rising temperatures that result from global warming, climate change has become a significant menace to the hydrological cycle of river catchments. Indeed, global warming is related to disturbances of precipitation's intensity and frequency, considering certain climatic conditions [2,3]. Moreover, climate change also sustainable flood management finds one of its most important bases in reliable design flood estimation [42].
Thus, considering this framework, meteorological and climatic events and forces have to be seen as pertinent issues. Indeed, the significant levels of instability that exist in the hydrogeological field pose a serious problem-namely, in terms of the effects that they might end up having on the population, socioeconomic activities, infrastructure, and productivity [13]. However, the fact is that floods-like most natural phenomena-cannot be prevented. On the contrary, certain activities carried out by people can, in fact, enhance the impacts of these natural disasters; if one adds climate change to those activities, then not would only the impacts of floods increase, but they would also become more recurrent events. When floods occur, rivers tend to occupy areas that are inherent to their natural course but have been encroached upon by human settlements over the decades. Hence, it is usual for floods to be misclassified as natural disasters, whereas the most accurate classification would be a "calamitous event", given that it affects the population and their activities. The impact that this type of phenomenon has in Madeira, Portugal, is significantly high, which then constitutes a matter of social significance due to the number of human victims and the substantial levels of damage that floods cause to buildings, infrastructure, and industries. In order to ensure that proper management of the land is accomplished, it becomes necessary to conduct a hydrogeological risk assessment, since this allows for the estimation of losses generated by an extreme event, given a certain return period. It is even worth mentioning that this assessment can be expressed either in terms of negative impacts or through the probability of such an event occurring [13,43]. When assessing risk, three aspects need to be considered: hazard (H), vulnerability (V), and exposure (E). With hazard, we point out how likely it is for flooding to occur within a specified timeframe and location. With vulnerability, we refer to the extent of the impact that will affect one or more elements at risk, considering a phenomenon of a certain level of intensity. Finally, with exposure, we refer to the level of harm that could affect multiple elements at risk-e.g., individuals, goods, and real estate. The final result from the combination of these three components can be considered a risk assessment [44].
According to the European Parliament's Directive 2007/60/EC issued on 23 October 2007 [45], regarding the evaluation and control of flood risk, hydraulic studies are necessary for risk assessment. Therefore, there is a requirement to resort to accurate predictive instruments and analytical methodologies in order to identify the real conditions of risk related to the incidence of a given flood event. Regarding urban flood studies, these must provide tools for establishing suitable estimation and prevention strategies. Additionally, flood investigations in urban areas have to consider various aspects, such as the process of gathering topographical data, the characterization of phenomenological processes related to flood currents and their effects on buildings and other structures, and the algorithm selection to solve model equations [46]. Ultimately, these studies must result in graphical outputs-that is, hazard and risk maps. Hence, opting for the correct perspective might be crucial, since graphical results often depend on the author's position when observing reality. When selecting a "point of view", it is important to take multiple factors into account -for instance, the synthesis of the multidimensionality of the features to be represented, the various methodologies that are able to provide a characterization of them, and the representation's purpose. It is worth mentioning that the information from this perspective is particularly important when dealing with environmental emergencies. Moreover, effective risk communication is key to engaging all stakeholders and guaranteeing appropriate management strategiesespecially in emergency scenarios [47]. Therefore, presenting the results through virtual scenarios is relevant, as this allows a suitable option to illustrate the conditions that a given region might end up experiencing when a flood takes place. In this case, methodologies must be based on models that accurately depict the hydraulic phenomena, as the images that are generated have to precisely express the computed water surfaces [48]. Multiple mitigation measures could be mentioned when analyzing floods. Those measures can be split into two categories: structural and non-structural. On the one hand, the latter category encompasses mitigation strategies related to the prediction, prevention, and adaptation processes; on the other hand, protective measures tend to be associated with structural measures, as they often involve engineering structures. However, currently, structural measures are being questioned due to the fact that these structures might end up amplifying the flooding risk and its impact if a failure or breach was to occur; additionally, these measures promote higher levels of vulnerability to these phenomena in flood-prone areas. Hence, non-structural measures are often considered as a viable alternative to lessen flooding's impacts. Examples of non-structural measures can be found in the meteorological and hydrological monitoring processes, land-use planning and regulation, flood hazard and risk mapping, and flood risk and emergency management [45]. A decrease in natural hazards entails various actions, ranging from prevention to emergency management [45,49].
Hence, it is imperative to develop regional studies focused on the flood vulnerability at the basin scale, so as to establish strategies to monitor and manage these events [50].
This study's focus was on conducting a hydrological analysis of the São João watershed (Funchal), with the objective of determining the anticipated peak flow rate, considering a recurrence time of 100 years, and comparing that value with the stream discharge point's drainage capacity. After demonstrating that the discharge point's hydraulic features are not sufficient to ensure that the drainage process is successfully concluded, considering the estimated peak flow rate, it became necessary to estimate the dimensions of a detention basin, aiming to guarantee that the entire flow would be drained and that the stream discharge point would operate normally. Additionally, one of this study's goals was to suggest a structural intervention in this watershed's discharge point-namely, modifying the roughness coefficient associated with the channel's riverbed and wallswithout generating significant urban impacts. As a result, the stream's minimum features were assessed to enhance its capacity of drainage without requiring any alterations to its dimensions.

Area of Study
The focus of this analysis is the São João watershed (Funchal), which is situated in the Madeira Archipelago-an island group in the North Atlantic, more specifically between latitudes 30°01′ N and 33°08′ N and longitudes 15°51′ W and 17°16′ W [51,52]. More precisely, as presented in Figure 1, the watershed under study belongs to the Funchal municipality and acts as an area of precipitation catchment that ultimately supplies one of the most important streams of this same municipality.
Moreover, this watershed has suffered from significant losses due to floods in recent years-namely, in 2010 and 2013. Considering that it is located in an extensively urbanized region, the level of soil sealing is considerably high. Thus, the presence of buildings and pavements in the area of study ends up enhancing the flooding problems that this region experiences [43,53,54]. Furthermore, Figure 2 shows that this watershed's stream discharge point contains sedimentation and vegetation, which contributes to a lower drainage capacity.  Regarding the level of conservation of this watershed's mainstream, it is mostly uniform in the urbanized areas, which could be simply confirmed in situ. However, the significant amounts of vegetation and sedimentation that are accumulated in the stream might be justified by the stream's reduced slope, which obstructs a rapid drainage process and hampers its ability to drag larger sediments.

Schematic of the Methodology
As Figure 3 demonstrates, this study's methodology can be divided into 6 different steps. In order to conduct an accurate morphometric and hydrological characterization of the watershed under study, it was important to undertake a thorough review of the existing literature in order to ensure that all of the pertinent information would be gathered. After this first step, methodologies suggested by multiple authors were taken into account, aiming to guarantee a satisfactory level of accuracy and reliability for this study. In terms of the remaining steps mentioned in Figure 3, these are described in the next section.

Morphometric Characterization of the Watershed
Regarding the morphometric characterization of the São João watershed (Funchal), it became necessary to resort to the parameters presented below [55][56][57][58]:  Gravelius index-KC: The relationship between the perimeters of the basin under study and a perfectly circular one-both with identical areas-was utilized to estimate the level of similarity of the watershed's geometric shape to a perfect circle [57]. This parameter can be obtained through Equation (1). As this is a dimensionless parameter, a watershed with a value closer to "1" will have a shape similar to a perfect circle, regardless of its dimensions, which would translate into higher levels of flood propensity [57].
where P = the watershed's perimeter (km); A = the area of the watershed (km 2 ).  Elongation factor-KL: The relationship between the watershed being analyzed and a rectangle-both with identical areas-was used to estimate the elongation of the watershed, regardless of its dimensions. This can be determined using Equation (2). If the value obtained for this parameter is higher than "2", then the watershed can be characterized as an elongated one [57].
where LE = equivalent length (km); lE = equivalent width (km); KC = Gravelius index-a dimensionless parameter; A = the area of the watershed (km 2 ).  Shape factor-KF: This relates the watershed's average length and width. This parameter can be obtained through Equation (3). Lower values are associated with more elongated watersheds and, therefore, with watersheds where the risk of flooding is lower-regardless of their size. Moreover, if the value obtained is "1", then the basin will have a square format.
where A = the area of the watershed (km 2 ); LB = the watershed's length (km). The length of a given watershed might be estimated using the gap that exists between the stream's discharge point and its furthest point. Nonetheless, it is worth mentioning that a watershed's length does not necessarily have to be equal to its main watercourse's length. Some variations between these two values might arise as a result of the larger size that the main watercourse tends to have, mostly because of its sinuosity. By resorting to the MDE file, which was developed by LREC-RAM (the Regional Civil Engineering Laboratory of the Autonomous Region of Madeira), it became feasible to morphologically characterize both the São João watershed (Funchal) and its main channel. In order to avert restrictions associated with using a single method, the data that were gathered during this study were utilized in the equations provided by various authors.
First, to conduct a morphometric analysis, it was necessary to establish a hierarchy based on the order and magnitude of the watercourses; for that reason, both the Strahler and Shreve classifications were utilized [54]. Indeed, these two classifications can be estimated by conducting an analysis based on the DEM file-a process that involves obtaining the "flow direction" and "flow accumulation" rasters through a tool named "flow order" [54]. Additionally, studies have pointed out that the Strahler classification is highly connected with a watershed's ratio of branching/bifurcation. Equation (4) allows for the estimation of each bifurcation degree [43,54,55,57].
where Ni = the number of watercourses labeled as "i"-a dimensionless parameter; Ni+1 = the number of watercourses labeled as "i + 1"-a dimensionless parameter.
To obtain this, it is necessary to divide the number of channels of a certain order by the number of channels encompassed in the order immediately above, regardless of their dimensions. Moreover, the average level of bifurcation can be obtained based on Equation (5).
where Ni = the number of watercourses labeled as "i"-a dimensionless parameter; Ni+1 = the number of watercourses labeled as "i + 1"-a dimensionless parameter; N1 = the number of watercourses of order "1". As this parameter only denotes the arithmetic mean of the bifurcation ratios, it is also dimensionless. Additionally, a key aspect for an accurate morphometric characterization of any watershed is its concentration time. This parameter indicates the time that the watershed's total area needs to contribute for the drainage process that will culminate in the stream's discharge point [43,54,55,57].
Considering that the methods utilized to calculate the concentration time can be classified as empirical, multiple methodologies might end up producing varying results for the same parameter. Therefore, in order to avoid extreme results, it is advisable to calculate the arithmetic mean. In this case, the arithmetic mean was calculated using the results derived from the methodologies of Kirpich, Témez, and Giandotti (Equations (6)-(8), respectively) [43].
where tC = the concentration time (hours); L = the main watercourse's length (km); i = the main watercourse's slope (m/m).
where tC = the concentration time (hours); A = the area of the watershed (km 2 ); L = the main watercourse's length (km); HM = the watershed's average height (m).

Precipitation Analysis
The precipitation analysis that was conducted in this research was based on a probabilistic analysis regarding intense short-term extreme events. To enable this analysis, data were gathered from public sources-including information about precipitation automatically recorded by the National Water Resources Information System (SNIRH). The Gumbel distribution was selected here because it was the probabilistic methodology that would best fit the already-acquired data and the anticipated forecasts for the watersheds located in Madeira [43,54]. Hence, Equation (9) can be utilized to estimate the annual maximum daily precipitation.
where PEST = the estimated maximum daily precipitation, in annual terms (mm); PM = the average precipitation, in annual terms (mm); S′ = the standard deviation of the sample (mm); KT = the frequency factor-a dimensionless parameter. where: where Xi = the sample value (mm); XM = the sample mean (mm); n′ = the number of samples.
where TR = the return period (years). A posteriori, given a certain duration, the intensity of the precipitation can be determined utilizing Equation (12): where I = the intensity of the precipitation (mm/h); PEST = the estimated maximum daily precipitation, in annual terms (mm); tC = the concentration time (hours); k = the coefficient of time distribution-a dimensionless parameter. where: where tC = the concentration time (hours). Given that the annual maximum daily precipitation is applicable only to events that last an entire day, the coefficient of time distribution assumes a key importance. Thus, since a watershed's concentration time is equal to the duration of the precipitation event, if one was to utilize the total level of daily precipitation, it would ultimately result in oversized hydraulic structures [55,57].

Drainage Capacity of the River Discharge Point and Peak Flow Rate
The Manning-Strickler equation, presented in Equation (14), was utilized to calculate the stream discharge point's drainage capacity; then, we established a comparison between the obtained value and the projected flow considering an extreme event, for a return period of 100 years. Moreover, to estimate the projected flow, multiple methodologies with a significant level of support among researchers were utilized-namely, Forti, Rational, Giandotti, and Mockus (Equations (16)- (19), respectively).
where QM = the stream discharge point's drainage capacity (m 3 /s). AM = the area of the river discharge point's cross-section (m 2 ); R = the hydraulic radius (m); i = the river discharge point's slope, in average terms (m/m); n = the coefficient of roughness of the walls and riverbed (m −1/3 s) (Table A1). where: where B = the stream discharge point runoff section's width (m); h = the stream discharge point runoff section's height (m); AM = the area of the stream discharge point's cross-section (m 2 ). It is worth mentioning that previous studies that have focused on this area were used as the main sources to gather information regarding aspects such as the stream's height and width in the discharge point area [57]. In fact, the confirmation of this first parameter was possible due to the utilization of the georeferencing process.
where QForti = the peak flow rate according to Forti (m 3 /s); A = the area of the watershed (km 2 ); b = for this parameter, the value was "2.35" if the maximum daily precipitation level stayed below 200 mm, and "3.25" for levels above 200 mm; c = for this parameter, the value was "0.5" if the maximum daily precipitation level stayed below 200 mm, and "1" for levels above 200 mm.
where QRational = the peak flow rate according to the Rational methodology (m 3 /s); C = the coefficient of surface runoff (Table A2); I = the intensity of the precipitation (mm/h); A = the area of the watershed (km 2 ).
where QGiandotti = the peak flow rate according to Giandotti (m 3 /s); λ = the reduction coefficient (Table A3); A = the area of the watershed (km 2 ); PMAX = the height of precipitation, considering that the event has a duration identical to the concentration time (mm); tC = the concentration time (hours).
where QMockus = the peak flow rate according to Mockus (m 3 /s); A = the area of the watershed (km 2 ); PEST = the estimated level of precipitation (cm); C = the coefficient of surface runoff (Table A2); tC = the concentration time (hours).
In order to guarantee that the population is secure, the dimensions of hydraulic structures have to consider a fill rate below 85% [55,57]. Hence, the implementation of mechanisms that enable the regulation of the runoff-for instance, spillways-assumes a significant level of importance.
As mentioned above, Equation (20) was used to calculate the fill rate. If the discharge point has insufficient drainage capacity to deal with the level of rainflow that exists in the watershed and cannot ensure that the safety margin is accomplished, it becomes necessary to estimate the dimensions of accurate structures of mitigation, such as detention basins. The fill rate is related to a given section's drainage capacity assuming a certain flow. Therefore, in a scenario where this parameter is greater than 100%, the section is not capable of dealing with such a high level of water, which ultimately leads to overflow [57].

Detention Basin Sizing
In scenarios where the discharge point is not capable of handling the volume of rainwater, it is necessary to design a spillway in order to guarantee the flow's normalizationan important aspect, considering that this flow will ultimately reach the stream's discharge point. In this case, a spillway of the Cipolletti type was picked due to its capacity for facilitating the runoff and reducing turbulence in the areas where there is contact with water [56,59]. Its dimensions were calculated using Equation (21).
Once the flow to be drained to the discharge point has been established and regulated, it becomes possible to calculate the level of water that the detention basin will retain. To determine this level of water, two approaches were considered: the Dutch method, and the simplified triangular hydrograph (STH) (Equations (22) and (23), respectively).
where QS = the flow drained by the spillway (m 3 /s); LSD = the sill's width (m 3 /s); HD = the height of the waterline above the sill (m).
where VA = the volume of storage (m 3 ); QP = each methodology's peak flow rate (m 3 /s); QS = the flow drained by the spillway (m 3 /s); tC = the concentration time (hours). It should be noted that the basis for Equation (23) can be found in the geometric examination that was conducted based on the STH approach ( Figure A1). Indeed, this equation was established by taking into account an event that lasts at least twice as long as a watershed's concentration time. Considering that the last particle of rainwater to reach the stream's discharge point would come from the farthest point and would be generated in the last moment of the precipitation event, it becomes clear that it would be necessary to consider the value of the concentration time for the volume drained by the river's discharge point [57].
Given that the Dutch method fails to account for the damping and delay of the precipitation hydrograph, the hydraulic structures whose dimensions are estimated considering this methodology might end up being oversized [60], as demonstrated in Figure 4, where qs is the spillway's runoff capacity, tc is the concentration time, tMAX is the maximum duration of precipitation (base), td is the delay associated with the beginning of the process of accumulation of water in the detention basin, Ha,MAX is the maximum capacity of storage, and i(tMAX) is the intensity of precipitation associated with the maximum duration. Thus, it was demonstrated that, when resorting to the Dutch method, the storage process and precipitation initiate at the same time, which is an unrealistic scenario, given that storage will not begin until the moment when the downstream-drained flow exceeds the runoff capacity of the spillway.

Modification of the Roughness Coefficient
Additionally, from a structural perspective, as a mitigation strategy, the alteration of the roughness coefficient of the watercourse's riverbed and walls was considered. One of the most significant advantages associated with this measure is the fact that it enhances the drainage capacity by diminishing the friction level. In terms of the Manning-Strickler equation, this measure can be seen as a modification of the value associated with the "n" parameter, aiming to enhance a certain watercourse's flow, which might be accomplished by changing the material that covers the stream's riverbed and walls [57].

Results
Regarding the results presented in this section, it is relevant to note that these values were obtained by resorting to the aforementioned formulae. Furthermore, an analysis of each of the indicators presented in Table 1 was conducted in order to study the characteristics of this watershed's main watercourse, from a morphometric perspective. These indicators were later correlated with reference values from different sources. First, the analysis focused on the watershed's area, as this parameter has significant importance in studying the volume drained to the stream's discharge point. Considering its area, a watershed might be characterized as follows: very large > 20 km 2 ; large > 10 km 2 ; medium > 1 km 2 ; small < 1 km 2 [61]. In this case, the São João watershed (Funchal) can be classified as "Large", which can be translated into a higher recurrence of floods. Nonetheless, it is worth noting that the reference values presented here are arbitrary and, therefore, might end up differing depending on the type of study that is being undertaken [61].
When comparing the altitudes of this watershed's borders and central region, as illustrated in Figure 5, one could point out that the borders generally present higher altitudes. Hence, this characteristic contributes to a rapid supply of water to the main watercourse, which means that the flow in the stream and, ultimately, in the discharge point will be larger. As Figure 6 demonstrates, this watershed has various low-and medium-order watercourses supplying the main one, indicating a high drainage capacity. Indeed, the hydrological behavior of a given region can be studied through the hydric density index. This parameter expresses a region's tendency to generate new channels. Thus, watersheds that possess higher levels of hydric density are prone to having numerous ephemeral channels [57].
Considering the data that were made available by the SNIRH, it became possible to study the precipitation [62]. In fact, we considered sample data that encompassed a period of 16 years, as shown in Table A4 and Figure A2. The results listed in Table 2 were obtained by resorting to a Gumbel-distribution-based probabilistic processing.
Then, when the intensity of precipitation had already been estimated, considering a 100-year return period, the peak flow rates were calculated based on the methods presented in the previous section, as displayed in Table 3. Given that the study region is located in a peripheral area with a significant presence of commercial buildings, in the Rational methodology, the coefficient associated with surface drainage was found to be 0.500 (Table 4). This coefficient represents the proportion of water that is expected to be drained superficially-in this case, approximately half of the volume of water is expected to drain superficially.   Table 4. Surface drainage coefficient adopted (source: [63]).

Urban Areas Occupation of the Land Coefficient of Surface Drainage
Commercial Area City Center 0.700-0.950 Peripheral Areas 0.500-0.700 In terms of the reduction coefficient (λ), the value that was chosen to estimate the flow through the methodology suggested by Giandotti can be seen in Table 5. Table 5. Adopted Giandotti reduction coefficient (source: [64]).

Area (km 2 )
Equivalent "C" <300 0.346 1.250 The Manning-Strickler equation was utilized to analyze the stream discharge point's drainage capacity; the results derived from this equation can be found in Table 6. Moreover, this equation allows for a more accurate analysis of whether this watershed's discharge point would require the implementation of a detention basin. Given that the stream's walls and bed presented distinct roughness coefficients, it became necessary to resort to a weighted mean considering these two coefficients. On the one hand, the walls are in good condition-and, therefore, n = 0.020; on the other hand, the riverbed's surface has the presence of sediments and vegetation, which implies that, in comparison to the walls, it presents a poorer state of conservation (n = 0.040). Additionally, it is worth noting that the slope in the region of the stream's discharge point is significantly low, and this ends up contributing to a decrease in terms of the flow's velocity and the section's drainage capacity. Thus, to consider a critical scenario, a slope of 0.01 m/m was considered for the reference section. Considering that the Fill Rate surpassed the predefined level of 85% for three of the four methodologies adopted-Giandotti, Mockus, and Rational-it was found to be necessary to implement mitigation measures, namely, a detention basin. Thus, the dimensions of the detention basin were estimated taking into account the flows previously obtained, in addition to the existing spatial and urban limitations due to the presence of infrastructure in the stream's surroundings.
To calculate the detention basin's appropriate dimensions, first, it was necessary to size a Cipolletti trapezoid spillway, since this would enable the flow's regularization and control. As a result, the spillway's characteristics are listed in Table 7. Then, in a second stage, both the Dutch method and the STH method were utilized to size the detention basin. However, due to the fact that these two methods are relatively simplified, using them might lead to oversized structures. In addition, this structure's height and width were defined as identical to the existing cross-section, aiming to diminish the impacts of the implementation of a detention basin from both the environmental and urban perspectives. The sole remaining variable was its length, which should not be larger than the main channel's length. Table 8 presents the dimensions obtained through the Dutch method and the STH method regarding this structure. Lastly, altering the roughness coefficient without changing the riverbed's vegetation was seen as a substitute approach to mitigate the impact of these events. As such, the values presented in Table 9 are related to the enhancement of the riverbed's level of conservation, which is expected to lead to a decrease in terms of friction and, consequently, to an increase in the stream's drainage capacity. The wall's revised roughness coefficient refers to a concrete-finished surface that presents a satisfying level of conservation, while the riverbed's sedimentation and vegetation are maintained. Table 10 presents a summary of the coefficients utilized.

Discussion
This study's main goal was to evaluate the need for simplified mitigation measures regarding flooding events in the São João watershed (Funchal). When it was demonstrated that this watershed required such measures, the implementation of a detention basin-a structural measure-was efficient in regulating the volume of water at the stream's discharge point [57]. In fact, while the fill rate obtained through three of the four methodologies utilized-Giandotti, Mockus, and Rational-was initially above 85% (107%, 96%, and 110%, respectively), it fell to a far more satisfying value as a result of the detention basin (54%-a value that is in line with the safety margin that had been pre-established). Additionally, this study's findings are in line with the flood risk assessments carried out by the Regional Directorate for Territorial Ordering and Environment (DROTA), as demonstrated in Table 11, thereby demonstrating that this study presents a satisfying level of accuracy. This study also focused on provoking the fewest possible alterations in the surrounding areas of the stream, as preserving natural elements and values in urban regions in currently considered to be a relevant requirement for the environmental revitalization of these areas [66]. Furthermore, ensuring that urban and natural systems harmoniously coexist is key to accomplishing the goal of sustainability [67,68]. Conversely, an unorganized process of urbanization might create urban voids [69].
It was therefore concluded that both the width and height of the streams should not be altered, which ultimately resulted in length being the sole existing variable. Given this assumption, despite the Dutch method presenting efficient results in terms of the regularization of the flow, it cannot be applied because it would require an alteration of either the height or width, given that this structure's length was larger than the main waterway's.
Identical restraints were considered with regard to the STH method. In this case, when comparing the two lengths, the detention basin's value did not surpass the one presented by the main waterway, meaning that the methodology was applicable.
Regarding the modification that was suggested for the stream's roughness coefficient, it was determined to only improve its conservation level, which means that the riverbed's vegetated and stony characteristics would be maintained, as completely removing all of the sediments and vegetation would be a complex, recurrent, and costly process. In terms of the walls, frequent maintenance should not be needed, considering that abrasioninduced wear would only take place in alluvial channels that drain substantial volumes of water and sizable sediments.
Even though it was a simplified measure, modifying the stream's roughness coefficient can be seen as an effective strategy, since it allows the discharge point to operate below the predetermined filling limit. Additionally, it is worth noting that this measure and the STH method can be applied simultaneously, resulting in a detention basin that would not demand such a large length.
Still, it is important to note that the methodologies mentioned above are relatively simplified-that is, they do not consider local specificities. As a result, the outcomes have a significant safety margin incorporated, leading to oversized structures.

Conclusions
This study's findings indicate that the São João watershed (Funchal) is vulnerable to flooding events, especially when extreme levels of precipitation occur; indeed, such vulnerability is also noted in DROTA's Flood Risk Report. This scenario occurs due to the stream bed's characteristics-namely, the fact that its surface has sediments and vegetation, which end up diminishing the flow's velocity and the stream's drainage capacity. This is particularly true for areas that have a very reduced slope-namely, the stream discharge point. In fact, the deficient drainage capacity that this watershed's discharge point possesses is supported by three of the four methodologies that were used-Giandotti, Mockus, and Rational.
In terms of the mitigation strategies, the Dutch method had no applicability, as it demanded an excessive length for the detention basin-even larger than the main watercourse's length. In contrast, the STH method brings more adequate results, as it allows this structure to be implemented without modifying the stream's height or width.
Lastly, the alteration proposed regarding the roughness coefficient also proved to be efficient. Among the advantages that this strategy possesses is the fact that it is a relatively simple measure to implement.
Given the impossibility of addressing all of the aspects that constitute a more extensive and effective study in this case, further research could be conducted to complement or improve the findings presented here-for instance, studying the drainage capacity associated with the existing urban hydraulic system, aiming to diminish the volume stored in detention basins; analyzing the deposition of sediments considering the main watercourse's entrainment velocity [70]; monitoring the artificial channels' walls' deterioration levels, as a result of abrasion, and analyzing the maximum amount of time to go ahead with maintenance (for instance, desilting and silting processes); studying how the artificial channels' processes of degradation can be correlated with the tributaries' water quality [71,72]; analyzing the impact of urbanization on the increase of the flow; planning and estimating the costs associated with the adoption of the mitigation strategies suggested in the present study; studying the impact of tide levels on the drainage capacity of artificial water channels and the probability associated with flooding events in the area of the stream's discharge point; and studying the impact of this type of channel from the perspective of territorial planning and management-i.e., how they can be adapted to watersheds located in rural areas.
The findings of this study are consistent with the results and outcomes suggested in analogous studies that also considered case studies and simulations as catalysts for scientific advancement [73,74]  Funding: Through the funding provided by MÁXIMA DINÂMICA REPARAÇÕES E CON-STRUÇÕES LDA and DANIEL SILVA SANTOS UNIP LDA, we were able to obtain the necessary resources to conduct the research, analyze relevant data, and draft this article. Furthermore, their trust and investment in our work demonstrate the recognition of the value and potential of this research. We hope that this enduring partnership continues to promote research and development in the field of civil construction. Once again, we are deeply grateful to MÁXIMA DINÂMICA REPARAÇÕES E CONSTRUÇÕES LDA and DANIEL SILVA SANTOS UNIP LDA for their invaluable support.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are openly available. Moreover, it is possible to contact one of the study's authors.
Acknowledgments: Special thanks are due to MÁXIMA DINÂMICA REPARAÇÕES E CON-STRUÇÕES LDA and DANIEL SILVA SANTOS UNIP LDA, both Civil Construction Companies, for the generous financial support granted, which made this study and article publication possible. Their significant contribution was fundamental to the success of this research project. We would like to express our sincere gratitude to MÁXIMA DINÂMICA REPARAÇÕES E CONSTRUÇÕES LDA and DANIEL SILVA SANTOS UNIP LDA for their continuous support in advancing scientific knowledge within the field of civil construction. Without their sponsorship, this publication would not have been possible.