Study of the Effect of an Environmentally Friendly Flood Risk Reduction Approach on the Oman Coastlines during the Gonu Tropical Cyclone (Case Study: The Coastline of Sur)

: Tropical cyclones may be destructive in the coastal region, such as the Gonu tropical cyclone, which affected the Arabian Peninsula and parts of southern Iran in 2007. In this study, a coupled MIKE 21/3 HD/SW (hydrodynamic/spectral wave) model was used to simulate the inland ﬂooding inside the Sur port during the Gonu tropical cyclone. The MIKE 21 Cyclone Wind Generation (CWG) tool was utilized to generate the cyclone’s wind and pressure ﬁeld. The required input data were obtained from the International Best Track Archive for Climate Stewardship (IBTrACS) and imported into the CWG tool. In this study, the wind and pressure ﬁelds were compared between the analytical vortex model and European Centre for Medium-Range Weather Forecasts (ECMWF) data during the Gonu cyclone passage. Moreover, by developing a new model, artiﬁcial Mangroves’ effect on inland ﬂooding was investigated. The results show that, contrary to the ECMWF data, the analytical vortex models well captured the storm event’s wind and pressure ﬁeld. Furthermore, the ﬂood hazard is calculated based on the inundation depth, ﬂow velocity, and area’s vulnerability. The ﬂood hazard map shows that 5% of the coast is at high-risk, 49% is at medium-risk, and 46% is at low-risk class in the Sur port. By applying Mangroves as ﬂood risk reduction, the high-risk area is almost completely removed. However, medium and low-risk zones increase by 50% and 50%, respectively. This information could be helpful in disaster risk reduction and coastal management in the future.


Introduction
The storm hazards impact the coastal zones. These extreme conditions may cause damage to the coastal facilities and human life [1,2]. Many communities are located in the coastal area and comprise a large amount of population. Hence, these communities are commercial hubs and industrial activities centers and are valuable areas [3,4].
Due to the increasing global temperature, the severity and the frequency of tropical cyclones in the oceans have increased since 1970 [5]. Based on the statistics, there are over 80 tropical cyclones worldwide annually [6]. The tropical cyclones in the Atlantic and Eastern Pacific Oceans are known as Hurricane, the Western Pacific Ocean as Typhoon and the Indian Ocean and the Arabian Sea as Cyclone [7]. After forming the cyclones in the Arabian Sea, they mostly move toward the south and southeast Oman, west India, and south Pakistan and rarely enter the Oman Sea [8]. However, these cyclones can move different paths, affecting the coasts of countries such as Iran and even the United Arab Most of the tropical cyclones in the Arabian Sea are formed in pre-monsoon (May) and post-monsoon seasons (October and November) [20]. However, some cyclones are formed between June to early September (Kharif season). The Kharif season monthly means of daily average wind data for the Sur port station were obtained from European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis of the ERA-interim data from 1981 to 2019 ( Figure 2b). As shown in the Kharif season, the wind speed increased and reached up to 4 m/s, and the northeast is the dominant wind direction in the Sur port.

The Governing Equation of MIKE21 SW
The basis of the MIKE 21 SW model for wave prediction is to solve the energy transfer Khaniki et al., simulated the characteristic wave height of six tropical storms in the Indian Ocean [11]. Dibajnia et al., evaluated the maximum wave height simulation from the 1889's Gonu cyclone in the Arabian Sea [8]. Fritz et al. assessed the depth of the flood by Gonu cyclone in some affected ports of Oman including the Sur port, as a field study [10]. Mashhadi et al., used a third-generation SWAN model for modelling the Gonu cyclone and compared the numerical results through a buoy in the Chabahar port [12]. Bakhtiari et al., simulated the recent storms that occurred in the Arabian Sea, such as the Gonu cyclone, through the two-dimensional MIKE 21 model and compared the results of the hydrodynamic and wave models with the measured data in the Oman Gulf [13]. Other research studies have been conducted to evaluate the oceanographic responses to Gonu tropical cyclone and Wave pattern generated [14][15][16][17].
The previous studies and modelling concentrated more on evaluation and predicting the effects of offshore storms. However, increasing the wave height may affect coasts and land areas as well as offshore locations. The Gonu cyclone caused some flooding by the storm surge in some parts of the Oman coast, including Muscat, Qirat, Sur, and Ras al-Hadd. The flood damage in Sur port was significant in which approximately 4294 residential homes were damaged [18]. The flooding hazard assessment for Sur port and understanding how the Gonu cyclone proved to be such a destructive event are the purposes of this study. This information could be helpful in disaster risk reduction. In this study, the flood caused by the Gonu cyclone in the coastal part of Sur port was modelled by MIKE 21 modelling software. Meanwhile, artificial Mangroves' effect as an environmentally friendly flood-risk-reduction approach on the coastlines of Sur was investigated during the Gonu tropical cyclone. For this purpose, the impact of Mangroves on coastal inland flooding was considered. Further, the results of a bigger model related to the monitoring and modelling of the Makran Coastline project [19] were used to define the existing smaller model's boundary conditions. This study's main objectives are to investigate the risk of flood due to the Gonu cyclone and identify vulnerable areas in the Sur port. This information can help the managers and governments for coastal planning and management in the future. Moreover, the effect of using an environmentally friendly method to mitigate the destructive influences of cyclones is another goal of this study.

Case Study
Historically, the Sur is known for being an important destination point for sailors. Today, the sea still plays an integral part in Sur's life. The port has a lagoon, and it is protected by rubble mound breakwater ( Figure 1). In this port, a considerable portion of the population inhabits the coastal areas or is engaged in commercial and recreational activities. Therefore, the storm can cause a lot of loss of life and economic damage to the port's inhabitants. The considered area is shown in Figure 2a, and the dominant direction of the waves in this area during the cyclone is from the east and south-west, and the offshore wave height could reach up to 9 m. Most of the tropical cyclones in the Arabian Sea are formed in pre-monsoon (May) and post-monsoon seasons (October and November) [20]. However, some cyclones are formed between June to early September (Kharif season). The Kharif season monthly means of daily average wind data for the Sur port station were obtained from European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis of the ERA-interim data from 1981 to 2019 ( Figure 2b). As shown in the Kharif season, the wind speed increased and reached up to 4 m/s, and the northeast is the dominant wind direction in the Sur port.

The Governing Equation of MIKE21 SW
The basis of the MIKE 21 SW model for wave prediction is to solve the energy transfer equation along with the source and sink terms. The energy transfer equation is considered in its spectral form to account for the random sea wave nature.
The wave model's governing equation is based on a wave action density spectrum N(σ,θ), defined in Equation (1).
where σ is the independent phase parameter (the relative angular frequency), and θ is the wave propagation direction. E is the energy density [21]. Most of the tropical cyclones in the Arabian Sea are formed in pre-monsoon (May) and post-monsoon seasons (October and November) [20]. However, some cyclones are formed between June to early September (Kharif season). The Kharif season monthly means of daily average wind data for the Sur port station were obtained from European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis of the ERA-interim data from 1981 to 2019 ( Figure 2b). As shown in the Kharif season, the wind speed increased and reached up to 4 m/s, and the northeast is the dominant wind direction in the Sur port.

The Governing Equation of MIKE21 SW
The basis of the MIKE 21 SW model for wave prediction is to solve the energy transfer equation along with the source and sink terms. The energy transfer equation is considered in its spectral form to account for the random sea wave nature.
The wave model's governing equation is based on a wave action density spectrum N(σ,θ), defined in Equation (1).
where σ is the independent phase parameter (the relative angular frequency), and θ is the wave propagation direction. E is the energy density [21].
Eng 2021, 2 144 MIKE 21 SW includes directional decoupling and full spectral of formulations. Parameterisation is performed in a frequency range by inputting the zero moment and the wave action spectrum's first momentum as parameters. The wave action balance equation is formulated in Cartesian or spherical coordinates (Equation (2)).
where N(x, σ, θ, t) is the action density, t is the time, x = (x, y) are the Cartesian coordinates, and v C x , C y , C σ , C θ is the propagation velocity of a wave group in the four-dimensional phase space of x, σ, and θ. S is the source term for the energy balance equation, and ∇ is the four-dimensional differential operator in the x, σ, and θ space. In this present study, the full spectral formulation was chosen.

The Governing Equation of MIKE21 FM
The hydrodynamic model of MIKE 21 FM is based on the numerical solution of the two dimensional incompressible Reynolds averaged Navier-Stokes equations summoning the assumptions of Boussinesq and hydrostatic pressure [22].
The continuity equation in Cartesian coordinates is in the form of Equation (3).
∂h ∂t The horizontal momentum shallow water equations in Cartesian coordinates is expressed as Equations (4) and (5).

Data and Modeling
The data are used to model the coastal flood caused by the Gono cyclone in the Sur port including wind, bathymetry, and water level. The data were obtained using the available international databases, which are given in Table 2. One of the principal constituents for the modelling in storm condition is the precise wind and pressure field because the wind and pressure cause shear stresses, which made surge current and wave [23,24]. If the available wind fields cannot capture the cyclone wind fields or wind data are not available at all, then the charts of atmospheric surface pressure can be used [25].
In this study, the MIKE 21 Cyclone Wind Generation (CWG) tool (Danish Hydraulic Institute (India), dhigroup.com, Headquarters location: Hørsholm, Denmark, Founded: 1964, Parent organization: DHI Water & Environment) was utilized to generate the cyclone's wind and pressure field. The Young and Sobey (1981), Holland (1980), and Rankine (1872) vortex models were used to generate the wind and pressure fields during the passage of the Gonu cyclone [26][27][28]. Cyclone position, central pressure, maximum sustained wind speed, and radius to the maximum wind (R Vmax ) are the model's requirement input data [25]. International Best Track Archive for Climate Stewardship (IBTrACS) global tropical cyclone database has been used to obtain these data. The results were compared with observed wind and pressure data obtained from 3 hourly Jask meteorology synoptic station and the ECMWF reanalysis ERA5 data. ERA5 has a special resolution of 0.25 • and a temporal resolution of an hour. The wind and pressure field modelled using the CWG tool was validated using available data from the Jask meteorological synoptic Station. Due to the Gonu cyclone track's proximity to the Jask station (about 16 km), the ECMWF ERA5 could not capture the cyclone wind and pressure fields. However, ERA5 is accurate enough in non-cyclonic weather conditions. The analytical vortex models well captured the wind and pressure field during cyclone passage from the study area ( Figure 3).
For the statistical analysis of the dimensionless goodness-of-fit indicator, the Nash and Sutcliffe coefficient of efficiency is calculated as follows [29]: where, O i and P i represent the sample (of size N) containing the observations and the simulation, respectively. O is the mean of the observed values SD represents the standard deviation of the observation data, and the root mean square error (RMSE) is given by Equation (7): Eng 2021, 2 146 Eng 2021, 2, FOR PEER REVIEW 6 (a) (b) For the statistical analysis of the dimensionless goodness-of-fit indicator, the Nash and Sutcliffe coefficient of efficiency is calculated as follows [29]: where, Oi and Pi represent the sample (of size N) containing the observations and the simulation, respectively. Ō is the mean of the observed values SD represents the standard deviation of the observation data, and the root mean square error (RMSE) is given by Equation (7): Ritter and Munoz-Carpena proposed four model performance classes based on NSE ranges (Table 3), denoted as Unsatisfactory, Acceptable, Good, and Very good [30]. During the Gonu cyclone passage from the Jask synoptic station (7 June), the NSE values are 0.68, 0.66, and 0.72 for the Young and Sobey, Holland, and Rankine vortex model results, respectively. In contrast, the NSE of the ECMWF ERA5 is 0.57. According to the criteria for the goodness-of-fit evaluation (Table 3), the Young and Sobey, Holland and Rankine vortex model results' performance is "Acceptable." However, the ECMWF ERA5 data performance rating is "Unsatisfactory," while, in the normal condition, the ECMWF ERA5 data performance is "Good" (NES = 0.81).
Wind fields were extracted from the MIKE 21 CWG a grid extending from the latitudes of 22.15-24.15° N and longitudes of 59-60.45° E, from the 5-6 June 2007 ( Figure 4). As shown in Figure 4, during the typhoon's passage, the MIKE 21 CWG data represented the wind field accurately within the Sur port's vicinity. Ritter and Munoz-Carpena proposed four model performance classes based on NSE ranges (Table 3), denoted as Unsatisfactory, Acceptable, Good, and Very good [30]. During the Gonu cyclone passage from the Jask synoptic station (7 June), the NSE values are 0.68, 0.66, and 0.72 for the Young and Sobey, Holland, and Rankine vortex model results, respectively. In contrast, the NSE of the ECMWF ERA5 is 0.57. According to the criteria for the goodness-of-fit evaluation (Table 3), the Young and Sobey, Holland and Rankine vortex model results' performance is "Acceptable." However, the ECMWF ERA5 data performance rating is "Unsatisfactory," while, in the normal condition, the ECMWF ERA5 data performance is "Good" (NES = 0.81).
Wind fields were extracted from the MIKE 21 CWG a grid extending from the latitudes of 22.15-24.15 • N and longitudes of 59-60.45 • E, from the 5-6 June 2007 ( Figure 4). As shown in Figure 4, during the typhoon's passage, the MIKE 21 CWG data represented the wind field accurately within the Sur port's vicinity.
The topographic data (bathymetry and hypsometry) have been generated by merging two different datasets. The Digital Elevation Model (DEM) with 30 m spatial resolution, which could be interpolated onto the model, was created for the shore zone as hypsometry data. The bathymetry dataset was complemented by the General Bathymetric Chart of the Ocean (GEBCO). However, due to the importance of sea depth data near the shore, the study area's shore profile was modified [31]. Figure 5 shows the modification beach profile. In addition, a dfs2 file was prepared for importing into the MIKE 21 through ArcGIS in the land area. For this purpose, the raster map of the area was created through the IDW tool in ArcGIS at the elevation points. Further, the limited DEM file was available, which can be imported into MIKE 21 as region-specific data through converting it to a dfs2 file.  shore, the study area's shore profile was modified [31]. Figure 5 shows the modification beach profile. In addition, a dfs2 file was prepared for importing into the MIKE 21 through ArcGIS in the land area. For this purpose, the raster map of the area was created through the IDW tool in ArcGIS at the elevation points. Further, the limited DEM file was available, which can be imported into MIKE 21 as region-specific data through converting it to a dfs2 file.  Figure 7 shows a comparison between the HYCOM GOFS 3.1 data and the simulated water level. It can be seen that the numerical result reasonably matches the HYCOM GOFS 3.1 data. The NSE value in this comparison is 0.91, and the performance rate of fit is "Very Good" ( Table 3 Eng 2021, 2, FOR PEER REVIEW 9 coefficient must be considered to select suitable values. The most crucial factors considered in the previous studies are surface roughness, vegetation, seasonality, level of water, surface irregularities, and obstructions [32]. In this study, the Manning's coefficient based on land cover data is used. The model was calibrated based on HYCOM GOFS 3.1 water level data to achieve a suitable local Manning coefficient for both open water and land area in Sur port.     Figure 7 shows a comparison between the HYCOM GOFS 3.1 data and the simulated water level. It can be seen that the numerical result reasonably matches the HYCOM GOFS 3.1 data. The NSE value in this comparison is 0.91, and the performance rate of fit is "Very Good" ( Table 3). Note that the HYCOM GOFS 3.1 water level data does not exist between the 4 June 2007 12:00 a.m. and 6 June 2007 12:00 p.m.   The Manning coefficient's determination is a challenge because of its empirical nature, and judgment and experience are needed. First of all, the parameter affecting the coefficient must be considered to select suitable values. The most crucial factors considered in the previous studies are surface roughness, vegetation, seasonality, level of water, surface irregularities, and obstructions [32]. In this study, the Manning's coefficient based on land cover data is used. The model was calibrated based on HYCOM GOFS 3.1 water level data to achieve a suitable local Manning coefficient for both open water and land area in Sur port. Table 4 shows the final Manning's coefficient values for the study areas. Meanwhile, Manning's coefficient value (n) in Table 4

Results and Discussion
The flood model results caused by the Gono cyclone in the Sur port are presented in two sections. First of all, the inundation and current speed model prediction, and the other, Mangrove for flood mitigation scenario model. The mentioned sections are presented in the following.

Inundation Model Results
Gonu cyclone generated destructive storm surge at the Sur port. Figure 8 shows the Storm Surge in the Sur port during the Gonu cyclone event. As shown in Figure 8, when the cyclone passed by the considering area, the storm surge near the Sur shore achieved its peak value. Meanwhile, it seems that there are no significant differences among the results of storm surge from different analytical vortex models.

Inundation Model Results
Gonu cyclone generated destructive storm surge at the Sur port. Figure 8 shows the Storm Surge in the Sur port during the Gonu cyclone event. As shown in Figure 8, when the cyclone passed by the considering area, the storm surge near the Sur shore achieved its peak value. Meanwhile, it seems that there are no significant differences among the results of storm surge from different analytical vortex models.   Figure 9b. The maximum flow velocity varies all over the study area, with its maximum magnitude being more than 0.9 m/s.   Figure 9b. The maximum flow velocity varies all over the study area, with its maximum magnitude being more than 0.9 m/s. In this study, the flood hazard is calculated through a matrix (Table 5) based on the inundation depth, flow velocity, and area's vulnerability. Vulnerability or, on the other hand, nature of the area is based on the type of buildings, the construction methods, and land-use at three levels: Low (Multi-story apartments), Medium (Typical residential area (2-storey), commercial and industrial properties), and High (Bungalows, mobile homes, busy roads, parks, single-story schools, campsites, etc.) vulnerabilities [34]. Table 5. Flood hazard matrix thresholds as a function of inundation depth and flow velocity and nature of the area.

Nature of Area Inundation Depth × Flow
Velocity (m 2 /s)

Low Vulnerability
Medium Vulnerability High Vulnerability In this study, the flood hazard is calculated through a matrix (Table 5) based on the inundation depth, flow velocity, and area's vulnerability. Vulnerability or, on the other hand, nature of the area is based on the type of buildings, the construction methods, and land-use at three levels: Low (Multi-story apartments), Medium (Typical residential area (2-storey), commercial and industrial properties), and High (Bungalows, mobile homes, busy roads, parks, single-story schools, campsites, etc.) vulnerabilities [34]. Table 5. Flood hazard matrix thresholds as a function of inundation depth and flow velocity and nature of the area. The flood hazard map based on the literature discussed above and according to Table 5 is reported in Figure 10. As shown in Figure 10, approximately 5% of the coast is assigned to the high-risk class, 49% of Medium, and 46% of low-risk. The area with a low-risk to life is located mostly in the central and eastern parts of the model in which the breakwater is located. While the area with high-risk to life is shown in the western part of the Sur coast. The medium-risk to life zone exists in both the eastern and western parts of the coast.

Mangroves for Mitigation Scenario
Mangroves are the first line of defense against coastal flooding in many tropical and subtropical coastal areas because of their strong structures. Mangroves reduce waves and storm surges, and the roots of Mangroves protect the coastline against erosion and stabilizing the soil [35]. Furthermore, Mangroves survive in extreme environmental conditions and salty water. Because of all these benefits and the suitable Omani coastlines' favorable situation to grow Mangroves forests, it seems Mangroves forests are an environmentally friendly approach to protect the coastal zones against inland flooding, especially during tropical cyclones. The Sur port coastal flood characteristics have evaluated in two scenarios, including with and without Mangroves.
In this research, a Mangroves forest along the Sur port Shoreline was considered by determining a local Manning coefficient and the Mangroves forest's width. Figure 11 shows the assuming location of the Mangroves along the Sur port shoreline in the simulations. The Manning's coefficient (M) for open Mangroves in the coupled model was considered to be 20 m 1/3 /s 1 [36], and the width of Mangroves was assigned as 90 m.

Mangroves for Mitigation Scenario
Mangroves are the first line of defense against coastal flooding in many tropical and subtropical coastal areas because of their strong structures. Mangroves reduce waves and storm surges, and the roots of Mangroves protect the coastline against erosion and stabilizing the soil [35]. Furthermore, Mangroves survive in extreme environmental conditions and salty water. Because of all these benefits and the suitable Omani coastlines' favorable situation to grow Mangroves forests, it seems Mangroves forests are an environmentally friendly approach to protect the coastal zones against inland flooding, especially during tropical cyclones. The Sur port coastal flood characteristics have evaluated in two scenarios, including with and without Mangroves.
In this research, a Mangroves forest along the Sur port Shoreline was considered by determining a local Manning coefficient and the Mangroves forest's width. Figure 11 shows the assuming location of the Mangroves along the Sur port shoreline in the simulations. The Manning's coefficient (M) for open Mangroves in the coupled model was considered to be 20 m 1/3 /s 1 [36], and the width of Mangroves was assigned as 90 m. Figure 12a shows the spatial map of flood characteristic for inundation depth for the mitigation flood scenario with Mangroves. The maximum inundation all over the computational domain is almost less than 1 m in this scenario. The flow velocity pattern for the storm event is presented in Figure 12b. The maximum flow velocity for the mitigation flood scenario with Mangroves varies all over the study area, with its maximum magnitude being less than 0.5 m/s. In this research, a Mangroves forest along the Sur port Shoreline was considered by determining a local Manning coefficient and the Mangroves forest's width. Figure 11 shows the assuming location of the Mangroves along the Sur port shoreline in the simulations. The Manning's coefficient (M) for open Mangroves in the coupled model was considered to be 20 m 1/3 /s 1 [36], and the width of Mangroves was assigned as 90 m.  The flood hazard map, according to the Table 5, are plotted in Figure 13 for the mitigation flood scenario model with Mangroves. As shown in Figure 13, for this scenario, approximately 50% of the coast is assigned to the medium-risk class and 50% to the lowrisk. By assuming the effect of Mangroves in the study area, the high-risk area almost removes completely. However, both medium and low-risk zones increase in the study area. The area with a low-risk to life is located mostly in the central and eastern parts of the model where the rubble mound breakwater is situated. The medium-risk to life areas exist mostly in the western parts of the coast, while some parts of the east coast are in medium-risk to life class.  The flood hazard map, according to the Table 5, are plotted in Figure 13 for the mitigation flood scenario model with Mangroves. As shown in Figure 13, for this scenario, approximately 50% of the coast is assigned to the medium-risk class and 50% to the low-risk. By assuming the effect of Mangroves in the study area, the high-risk area almost removes completely. However, both medium and low-risk zones increase in the study area. The area with a low-risk to life is located mostly in the central and eastern parts of the model where the rubble mound breakwater is situated. The medium-risk to life areas exist mostly in the western parts of the coast, while some parts of the east coast are in medium-risk to life class. risk. By assuming the effect of Mangroves in the study area, the high-risk area almost removes completely. However, both medium and low-risk zones increase in the study area. The area with a low-risk to life is located mostly in the central and eastern parts of the model where the rubble mound breakwater is situated. The medium-risk to life areas exist mostly in the western parts of the coast, while some parts of the east coast are in medium-risk to life class.

Conclusions
Tropical cyclones are a serious threat to human lives and properties in coastal areas, which cause high destructive effects in these zones. Evaluating new approaches to protect and reduce cyclone's destructive impacts on habitats and properties is an important issue in coastal engineering and management. This study investigates the risk of floods due to

Conclusions
Tropical cyclones are a serious threat to human lives and properties in coastal areas, which cause high destructive effects in these zones. Evaluating new approaches to protect and reduce cyclone's destructive impacts on habitats and properties is an important issue in coastal engineering and management. This study investigates the risk of floods due to the Gonu cyclone in the Sur port. Moreover, the impact of an assumption Mangroves forest as a flood risk reduction approach on the flood characteristics (inundation and flow velocity) is evaluated considering two scenarios, including with and without Mangroves forest along the shoreline. For this purpose, the coupled MIKE 21/3 HD/SW (hydrodynamic/spectral wave) model was simulated. Furthermore, the Young and Sobey, Holland and Rankine analytical vortex models were used to generate the wind and pressure fields during the Gonu cyclone passage. The following conclusions can be drawn:

•
The analytical vortex models well captured the wind and pressure field during the storm event. However, the ECMWF ERA5 data did not accurately capture the wind and pressure field due to the proximity of the Gonu cyclone track to the study area. However, the ECMWF ERA5 data are accurate enough in non-cyclonic weather conditions. • The storm surge near the Sur shore achieved its peak value of 0.44 m (by the Rankine vortex model). While the maximum surge in the normal condition was almost less than 0.1 m. This indicates that the Gonu cyclone led to a significant rise in seawater level that caused inland flooding. Meanwhile, it seems that there are no significant differences among the results of storm surge from different analytical vortex models. The maximum inundation all over the computational domain is almost less than 1 m for the flood scenario with Mangroves.

•
The maximum flow velocity varies all over the study area with its maximum magnitude of more than 0.9 m/s for the flood scenario without Mangroves. While for the flood scenario with Mangroves, the maximum flow velocity is less than 0.5 m/s.

•
The flood hazard map shows that 5% of the coast is at high-risk, 49% is at medium-risk, and 46% is at low-risk class for the flood scenario without Mangroves. • By applying Mangroves as flood risk reduction, the flood hazard map shows that the high-risk area almost removes completely. However, medium and low-risk zones increase by 50% and 50%, respectively.

•
The eastern part of the coast of Sur port seems to be well protected by the breakwater. While most parts of the west coast are unprotected and exposed to damage from storms and cyclones.
• Impact of large flood in urban and critical facilities and infrastructure can be reduced by timely and correct action before a disastrous storm, so managers and governments should be considered appropriate methods to mitigate cyclones' destructive effects. The use of artificial mangrove forests as an environmentally friendly method of protecting the coast can be an effective way to protect the Sur port against climate change and extreme environmental conditions.