Development of Prediction Model for Storm Surge Hazards in the Developing Countries

: Bangladesh is one of the most vulnerable countries in the world with around 718,000 deaths in the past fifty years. This country is especially in danger for cyclones because of its location in the triangular-shaped Bay of Bengal. The scientific scenario suggests that enlarged sea surface temperature will intensify cyclone movement. Tropical cyclones generate storm surges. Storm surges severely change the coastal environment, damage coastal structures, destroy forests and crops, inun-date the coastline with saltwater and cause loss of lives. Due to overcrowding in the mainland in Bangladesh, poor and landless people live on the coast and they face frequent cyclones and associated surges. They affect food and drinking water and there is danger of transmission risks of infectious diseases, such as diarrhea, malaria, eye infections, skin diseases, etc. Some problems following a cyclone are usually created for their low literacy rate and poor knowledge of the environment. The tangible monitoring and warning of cyclones and associated surges should be given more priority for the region. The main objectives of this paper are to highlight the existing activities of the model in storm surges and related areas in the Bay of Bengal. We will explain the progress of a location-specific real-time standpoint prediction system for providing effective and timely surge forecasts. We will also introduce a model through numerical experiments using a severe cyclone in April 1991 to predict the storm surges that will be used to reduce economic losses and the number of deaths during a strong storm surge in the coastal area of Bangladesh.


Introduction
Storm surges are associated with severe weather such as tropical cyclones that constitute the world's most catastrophic natural disasters [1].Among the most threatening calamities, storm surges stand out as the most damaging and undoubtedly as a cause of death and ruin as massive as that of earthquakes and tsunamis [2].In history, low-lying deltaic regions of Bangladesh and Myanmar are heavily subjected to storm surge hazards.Fortunately, in most recent storm surge cases, the death toll is not as high as previously since the storm surge warning systems in fact worked well.However, we keep in mind that storm surges could be a cause of high death rates in the future while the number of residents is rising in coastal regions [3,4].The key components contributing to calamitous surges in Bangladesh are as follows [5]:  Convergence of the Bay  High astronomical tides  Positive cyclone track  Shallow coastal water  Densely populated low-lying island.An updated and accurate storm surge prediction system is required to mitigate and prevent coastal disasters.It is essential to study local allocation characteristics exclusively along with previous storm surge heights in order to update the numerical forecast systems.Various ranges of motion must be investigated through numerical discrete structure of the governing equations to forecast surges in coastal areas exactly.It is also indispensable to examine long-term observation data perfectly in predicting the areas that storm surges might affect with precision [6][7][8].As the landward boundary is approached, substantial grid refinement is typically required to resolve important schemes and prevent energy from aliasing.In providing an adequate resolution in the near shore region without increasing the size of the discrete problem, a numerical method must be used that permits a very high degree of grid flexibility.

Model Equation
The following depth-averaged vertically integrated form of the mass conservation equation, and the x and y components of the momentum equation, respectively, are used in investigating the dynamical process in the sea [6]: In Equations ( 1)-( , where Ω is the angular speed of the earth rotation about its own axis, and φ represents the latitude of a position of interest) is the Coriolis parameter; g (=9.81 m/s 2 ) is the acceleration due to local gravity; ρ is the sea water den- sity, assumed to be uniform; f C (=0.0026) is the dimensionless bottom friction coefficient; x T and y T are the and y components, respectively, of the wind stress acting on the sea surface.The wind stress components mentioned above are derived following Ali [5] as To derive the components of the wind stress, wind field is required.As in [7], the wind field is generated through the following formula: In Equations ( 4) and ( 5), a u and a v stand for the x and y components of the sur- face wind, respectively, a ρ is the air density, D C (= 0.0028) is the surface drag coeffi- cient, 0 V is the maximum sustained wind at the maximum radial distance R and a r is the distance between the cyclone center and the point at which the wind field is desired. x

Boundary Conditions
The study domain has three open boundaries and one closed boundary.At the closed boundary (coastal or island boundary), the normal component of the depth averaged velocity is set to zero.Following Paul et al. [6], the western, eastern and southern open boundary conditions are, respectively, taken into account.These lead to ( ) where a, T , and ϕ stand for, respectively, the amplitude, period and phase of the tidal constituent under consideration.Meghna River mouth is considered as open boundary and the freshwater discharge through the river is considered following Roy [7].This leads to where Q denotes the volume of freshwater river discharge in a unit time and B represents its width.

Input
In our study, we need several types of input.The meteorological input, namely maximum sustained wind radius, maximum sustained wind speed and storm track information are obtained from the Bangladesh Meteorological Department (BMD) through personal communication.The time varying positions of the storm and its nature are presented in Table 1 for a better understanding.The study also needs bathymetry data which is collected from the British Admiralty Chart.Shepard interpolation is used to supply water depth at the grid points of the three schemes, as we will see later, representing water.Further, the study needs tidal constants to generate a tidal response on the area of interest, which are taken from the study due to Paul et al. [6].Following Paul et al. [9], the freshwater discharge through the river per unit time is taken as kg/s 5100 = Q and the remaining parameters have been assumed to have their standard value.

Set-Up of Nested Scheme
In order to incorporate the coastal complexity with minimum cost, a nested grid technique is used in this study.A high-resolution fine grid model, referred to as FMS, is nested into a scheme with relatively low resolution, which is referred to as CMS.It is of interest to note here that the FMS is designed to incorporate all the major offshore islands along the coast of Bangladesh.Now to incorporate properly the land-sea interface and bottom topographic detail of the Meghna estuarine area, which is referred to as the world's most vulnerable zone, a scheme with very high resolution is nested into the FMS.The innermost scheme is referred to as VFMS.To have a clear idea about the schemes, Figure 1 and Table 2 are inserted for a better perspective.It is to be noted here that Paul et al. [6] first used the of 10 MOL to solve shallow water equations in predicting water levels due to the tide-surge interaction along the coast of Bangladesh, where they used a nested grid technique without high resolution of grids.Thus, this study is an improvement on that of Paul et al. [6].The aim of the study is to solve the governing equations with a numerical method of lines (MOL).It is of interest to note here that the method is efficient over the standard finite difference method due to having some benefits, especially, in computational coast, stability criterion, and simplicity in solving partial differential equations.Thus, in this regard, the equations and the boundary conditions are discretized only for spatial derivatives/variables by means of a semi-implicit three point central finite difference technique.We consider the discrete points in the xy-plane by defining If any dependent variable ) , , ( t y x ξ at a grid point ) , ( , then with the aid of the notations For every grid point ) , ( , where Again, for every point ) , ( can be written as ( )  , Equation ( 3) can be written as ( ) The boundary conditions specified by Equation ( 6), the elevations at  in the following manner:

Working Procedure and Model Run
Discretized form of Equations ( 8)-(10) are available in [6].However, briefly they are inserted here for a better understanding.First Equation ( 8) is solved with RK(4,4) method for the elevation ζ at the internal (even, odd) grid points representing water of the CMS.Discretized BCs given by Equations ( 11)-( 13) are then used to estimate ζ at the bound- ary (even, odd) grid points.An interpolation is then adopted to obtain ζ at the leftover grid points representing water and the land-sea interface.The wind field is then generated using Equations ( 4) and (5).Then Equation ( 9) is used for estimating u at the interior (odd, odd) grid points representing water by inserting the parametric values involved and finally v is evaluated in a similar fashion at the interior (even, even) grid points solving Equation (10).After computing ζ , u , and v in the CMS, the scheme is coupled with the FMS following a process found in [6] to have the boundary values from the CMS to run the FMS.In a similar manner, the VFMS is run with the boundary information from the FMS.Along the northern open boundary segment, b U is taken into account for the VFMS using Equation ( 14).This process is repeated over time providing the updated values of ζ , u and v as initial conditions for computing WLs ( ζ ) owing to tide, surge and nonlinear interaction of tide and surge.For computing tide, surge and total water levels (water levels due to tide-surge interaction) all the models, namely tide, surge, and tide-surge interaction were run considering the time step as A stable tidal oscillation was obtained after four tidal cycles of integration that provided the profile of the sea surface if the cyclone is not taken into account.In getting the pure surge level, the surge model (in absence of astronomical tide) was also run from the cold start.In achieving total WLs (due to the dynamic interaction of tide and surge), the tide model was run first.After having a stable tidal regime, the surge model was then made to run over it.

Results and Discussion
Water levels due to tide, surge and their interaction are calculated at the stations shown in Figure 1, but for the sake of brevity, the peak water levels due to surge and interaction of tide and surge are presented in Table 3 and the computed time series of water levels due to the tide, surge and the interaction of tide and surge are presented in Figure 2 with observed data.It is seen from Table 3 that our computed peak surge levels (water levels due to meteorological forcing only) along the region of interest vary between 2.97-6.51m with 5.29 m at Chittagong.Paul et al. [6] predicted 3.52-6.70m surge along the coast of Bangladesh.Further, according to BMD, the peak surge at Chittagong at the time of landfall was 5.50 m.Thus, our computed surge level agrees well with the reported data by BMD at Chittagong and the simulated surge levels by Paul et al. [6].Again, our computed peak total water levels were found to vary between 3.86 m (Kuakata) and 7.28 m (Companigonj), which agree fairly well with the results obtained in Paul et al. [6] (see Table 3).Our computed time series of water levels with respect to the MSL due to tide, surge and the interaction of tide and surge associated with the storm in April 1991 at Hiron Point, Char Chenga (Hatiya) and Chittagong are displayed in Figure 2. The corresponding observed data collected from Bangladesh Inland Water Transport Authority (BIWTA) are also presented in Figure 2. It is seen from Figure 2 that when the surge is away from the coast, the tide dominates the surge level; whether the opposite characteristics can be found when the surge is nearer to the coast as is expected.However, it is perceived from Figure 2 that our computed water levels due to the interaction of tide and surge are in good agreement with the data obtained from BIWTA.

Figure. 2
Computed water levels with respect to the mean sea level due to tide, surge and their interaction with observed data at Hiron Point, Char Chenga and Chittagong.In each subplot, a black dotted curve represents the configuration for tide, black solid curve for surge, red dotted curve for tide-surge interaction, and a circle represents observed data.For validation of our simulated results, the root mean square analysis is carried out between the results attained in the study and observed data from BIWTA.The that came out in this regard are presented in Table 4.The root mean square error (RMSE) values obtained in this regard due to Paul et al. [9] are also presented in the same table for comparison.It is seen from Table 4 that the results attained by the study are considerable and comparable with those presented in [9].

Conclusions
In this study, the MOL in coordination with the RK (4,4) method is used to solve vertically integrated shallow water equations in Cartesian coordinates for simulating water levels along the coast of Bangladesh.The water levels due to the non-linear interaction of tide and surge associated with the cyclone in April 1991 are found to be in reasonable agreement with observed data from BIWTA on the basis of RMSE values.The outcome of this study thus can be utilized in practical forecasting.

Citation:
Naher, H.; Paul, G.C. Development of Prediction Model for Storm Surge Hazards in the Developing Countries.Environ.Sci.Proc.2021, 5, 7. https://doi.org/10.3390/IECG2020-08630 Academic Editor: Jesus Martinez Frias Published: 30 November 2020 Publisher's Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.Copyright: © 2020 by the authors.Licensee MDPI, Basel, Switzerland.This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Figure 1 .
Figure 1.Domains of the three schemes CMS, FMS and VFMS and actual coastal and island boundaries along with ten representative locations at which computed results are presented (after Paul et al. [9]).
discretize the equations of interest as follows.

2 M
60 = Δt s to ensure Courant-Friedrich-Lewy (CFL) stability criterion.It is noteworthy that our model can run individually and simultaneously.The first tide model was run from the cold start 0levels with respect to the mean sea levels (MSL) considering the effect of tidal constituent along the southern open boundary of the CMS in the absence of meteorological forcing.

Table 1 .
Time series for the positions and the nature of the April 1991 cyclone (Source: BMD).

Table 3 .
[6]puted peak water levels simulated by the present study with respect to the mean sea level (MSL) and those obtained in[6].

Table 4 .
Estimated RMSEs in meters.The errors have been estimated between computed and observed water levels from 02.00 UTC of April 28 to 02.00 UTC of April 30 for the storm in April 1991.The observed data were obtained from BIWTA.