Topographical Analysis of the 2013 Typhoon Haiyan Storm Surge Flooding by Combining the JMA Storm Surge Model and the FLO-2D Flood Inundation Model

The floods associated with the effects of an incoming tropical cyclone have an immense effect in the Philippines, especially with respect to agriculture, industry, livelihood, and public safety. Knowledge of how such storm surge flooding can affect the community is therefore of great importance. In this study, the mechanisms behind Typhoon Haiyan’s anomalous storm surge flooding in 2013, which resulted in more than 6300 casualties and 2.86 billion USD worth of damage in the Philippines, were investigated. The Japan Meteorological Agency (JMA) storm surge model and the FLO-2D flood model were used to simulate Typhoon Haiyan’s storm surge height and the extent of inundation, respectively. The storm surge input data were obtained from JMA typhoon data, and the digital terrain models used were gathered from the airborne interferometric synthetic aperture radar data. The model’s accuracy was also validated using field validation data of the extent of the observed storm surge in affected coastal areas. Topographical analysis of the inundated regions showed the effects of coastal shape, elevation, and position relative to the typhoon’s approach angle on storm surge flow depth and velocity. Storm surge maximum velocity appears to increase as the fluid flows to an increasingly elevated area. Observing fluid velocity in a coastal area with uniform storm surge discharge from all directions also showed that flow velocity tends to increase at the center. Greater flood depths were experienced in areas with lower coastal elevation and not directly located at the coast, compared to higher elevation coastal areas. Greater extents of storm surge flooding are expected in coastal areas that have a concave shape, as fluid is more likely to be dispersed when hitting a convex coast. Extents are likewise observed to be greater in coastal regions that are located perpendicular to the direction of the typhoon. The research also validated the option of using a combination of typhoon and flood models to simulate the inundation flooding caused by extreme weather events.


Introduction
The natural hazards generated by a tropical cyclone impose significant threats, especially to coastal communities.These hazards, such as rainfall-induced flooding and landslides, and the storm surges generated by the strong winds, have been known to result in substantial loss of life and damage to properties in the Philippines.One particular cyclone, Typhoon Haiyan in 2013, affected 16 million people and resulted in about 6300 deaths, making it the deadliest typhoon to hit the country in recent history [1].The number of casualties was primarily associated with the anomalous storm surge flooding generated by the typhoon [2].
However, modeling storm surges is not as common and well known as flood and debris flow modeling.This is mainly due to the rarity of storm-surge-related hazards compared with the hazards associated with flooding events.In recent years, however, several storm surge models have been developed by atmospheric and coastal researchers.These storm surge models include the Sea Lake and Overland Surges from Hurricanes model by the National Oceanic and Atmospheric Administration (NOAA) and National Hurricane Center [3,4], the Mike21 and Mike Flood models [5] developed by the DHI Water and Environment Corporation, the Delft3D coastal surge model [6] by Deltares, and the Japan Meteorological Agency (JMA) storm surge model [7] from the Japan Meteorological Agency.
The 2013 Typhoon Haiyan storm surge has been evaluated by several researchers using different models, these topics include: ensemble forecasting [8,9], boulder and sediment transport [10,11], storm surge simulations [2,12,13,14], local amplifications [15], and spatial variation of damages [16].Most of these research efforts focused on simulation of the storm surge of the event, but very few of those presented the inland flooding generated by the storm surge.[2] presented a simulated storm surge inundation and field validation results using the same models as in the present study.However, the paper showed only one affected area (Tacloban City) and focused more on the social aspects of the event rather than the scientific factors behind.
To determine the factors leading to the extreme flooding caused by the Haiyan storm surge, it is important to analyze the meteorological and geographical conditions present during the event that caused the aggravation of storm surge flooding.In other words, the topographical characteristics of the areas that made them more susceptible to inundation compared with the other areas that experienced the same typhoon intensity must be considered.
The objective of this paper is to analyze the geographical characteristics of the most affected areas during the typhoon storm surge event.The motivation in doing so is the lack of related literature discussing the vulnerability of coastal areas to storm surges as a result of their topographical characteristics.These analyses include the effects of coastal shape, elevation, the position relative to the typhoon's angle of approach to the storm surge depth, extent, and flow velocity.However, the meteorological conditions that instigated such intensity, which in turn caused the anomalous storm surge, are not covered in this research, but will be discussed in greater detail in our future research.
This research utilized the JMA storm surge model for storm surge simulation.The JMA model uses a two-dimensional shallow water equation and is used to model typhoon-induced storm surges that form in the northwestern Pacific basin.While there are a number of storm surge models that can simulate the same effect with much finer resolution, the JMA model was chosen for this research due to its higher advantage in terms of computational speed.The JMA model's capability of simulating reliable storm surge timeseries with less computational time is ideal for simulating larger domains, while not losing the objective of the research.The inland storm surge flooding, however, which the model lacks, was supplemented using the FLO2D flood model, a flood routing model that simulates overland flow over complex topography.In this study, field validation data are also provided to show the accuracy of the models used.The structure of the paper is arranged as follows: Section 2 describes the data and the data sources used, as well as the governing equations and procedures behind the storm surge and flood models.Section 3 shows the results of the inundation model, comparison with field observations, and the analysis of the storm surge depth and velocity in response to the difference in topography.Finally, Section 4 provides a summary and final conclusion for the research.

Data and Methods
To simulate the storm surge flooding generated by Typhoon Haiyan 2013, the research utilized the tropical cyclone best track data [17] obtained from the JMA typhoon archive.
The study domain elevation and topographical characteristics were represented using the airborne interferometric synthetic aperture radar-derived digital terrain model (DTM) with a 5 m spatial resolution.The bathymetric data used for the domain were the 2 min Global Gridded Elevation data (ETOPO2) gathered from the NOAA.
To simulate the Haiyan-induced storm surge height and the extent of inundation for the most affected coastal areas, the JMA storm surge model and the FLO-2D flood routing model were used together.

Storm Surge Model (JMA)
To simulate Typhoon Haiyan's storm surge height, the JMA storm surge model, a twodimensional numerical model that uses a vertically integrated momentum equation and a continuity equation [7], was employed.This model computes storm surge height caused by wind and inverted barometric effect, and its inputs are the typhoon track and bathymetric data.
The storm surge model's governing equations are as follows: ) where U and V are mass fluxes in the x and y directions, defined as where u and v are the x and y direction velocities, z is the vertical length, D is the water depth below mean sea level,  is the surface elevation, 0  is the water column height corresponding to the inverse barometer effect, f is the Coriolis parameter, g is the acceleration due to gravity,  is water density, ( ) where is the wind speed, and is the wind velocity.The drag coefficient is set using the results of [18] and [19]: We solved these equations by numerical integration, using the explicit finite difference method.For more details on the JMA storm surge model, refer to [20].
Typhoon track data were composed of meteorological descriptions of Typhoon Haiyan, such as its formation and dissipation dates, eye center location (latitude and longitude), maximum sustained wind speed in kts, and central pressure in hPa.Using these input data, which are essential for running the JMA storm surge model, we ran a storm surge simulation using three domains, as shown in Figure 1 At the sea-land interface, the normal component of velocity is zero.At open boundaries, the sea surface elevation is given by the inverted barometer effect; a gravity wave radiation condition was used in this study.Tides can be simulated but are excluded, since the major object of this simulation is for surge heights Our simulations produced storm surge height time series for each domain.We modified the storm surge height to incorporate astronomical tides computed using WxTide, a Windows-based tide and current prediction program.We determined the coastal provinces targeted for flood simulation from the results of our JMA storm surge simulation.The storm tide (storm surge depth + astronomical tide) time series from the JMA model was used as the input file for our storm surge inundation simulation.
The JMA storm surge model simulation identified the provinces of Iloilo, Leyte, Palawan, Samar, and Eastern Samar as having experienced the highest storm surge during Typhoon Haiyan.The location of the provinces can be viewed in the coastal inundation model results.The greatest (3) estimated storm surge height reached about 3.9 m in Batad, Iloilo (Figure 2).We then used the catchments of these provinces for our flood inundation simulation.

Coastal Inundation Model
One of the JMA storm surge model's limitations is that it does not simulate storm surge inundation.The extent of storm surge flooding predictions is crucial for pre-emptive evacuation and long-term urban planning.Storm surge inundation maps were employed in this study using the FLO2D model, a flood routing model that simulates channel, overland, and street flow over complex topography.We chose the FLO2D model to conduct our flood simulation because of its compatibility with the JMA storm surge model.The FLO2D model's limitations, such as the absence of wave coefficients and receding parameters, were not deemed significant enough to drastically affect this study's objective and results.FLO2D's governing equations are as follows: where t represents time; h is the flow depth; hV is the depth-averaged velocity in one of eight flow directions; x , e r is the excess rainfall intensity; f S is the friction slope component based on Manning's equation; 0 S is the bed slope pressure gradient; and g is the acceleration due to gravity.This equation represents a one-dimensional, depth-averaged channel flow.For the floodplain, although FLO-2D is a multi-direction flow model, FLO-2D's equations of motion are applied by computing average flow velocity across a grid element boundary one direction at a time [21].
The catchments to simulate for flooding were chosen using the previous storm surge simulation; we selected the top five provinces with the highest simulated storm surge.We applied boundary grids to the province shapefiles to create a catchment.Data with a grid resolution of 20 m × 20 m were used through the provided DTM data.The flood simulation runs smoothly if the storm surge's estimated peak discharge, Q , divided by the grid element surface area, A , does not exceed 3 m 3 s −1 .After considering the simulation time, we set the grid size to 20 m. Figure 3 shows a sample of inundation results, such as the maximum water surface depth and maximum flow depth, as shown in the FLO-2D model.We used the calculated hydro-curve produced using the JMA storm surge model as input data for the FLO2D model, and then conducted inundation modeling.

Results of Storm Surge Coastal Inundation Model kilometer
We overlaid storm surge inundation results from the FLO-2D model on a satellite image for validation.Validation maps show an inundation color scale with ranges of 0.01-0.5 m (light yellow), 0.51-1.00m (yellow), 1.01-1.50m (red), 1.51-2.00m (violet), and 2.01 m and above (blue).Validation points shown as blue circles observed flooding, and points represented as red circles experienced no inundation.
Figure 4 shows a comparison of observed inundation and simulated flooding for Tolosa, Leyte.Observed inundation and simulation results generally agreed.Flooding was observed (blue circles) in areas where the model had simulated storm surge inundation, and no flooding was observed (red circles) in areas with no simulated inundation.The same was true for the storm surge validation maps of El Nido, Palawan (Figure 5), Carles, Iloilo (Figure 6), and Marabut, Samar (Figure 7).However, the same cannot be said about the storm surge validation map of Guiuan, Eastern Samar (Figure 8).Flooding was observed at several validation points but showed no simulated storm surge, with the flooded area located 2000 km away from the coast (Figure 8).This is because the flooding that interviewees experienced may have been caused by rainfall, as Eastern Samar experienced about 120 mm of precipitation (Figure 8, lower left image) during Typhoon Haiyan.

Model Results Analysis
After establishing the reliability of the models that were used, it was necessary to determine why some coastal areas experiences worse storm surge compared to other areas.Based on the research by [22], Iloilo, Samar, Leyte, Palawan, and Eastern Samar provinces were among the top 25 Philippine coastal provinces most vulnerable to storm surge.A hypothetical track of a typhoon with Haiyan's wind strength would have a storm surge height as high as 7.45 m.This demonstrates that these provinces might possess geographical characteristics that make them more susceptible to damaging storm surges.

Effect of Elevation on Storm Surge Flow Depth
Varying storm surge flooding behavior can be seen in different elevation topographies.Figures 9 and 10 show cross-sectional and x height storm surge flooding in some of the simulated catchments.Figure 9 shows the effect of elevation on storm surge height, with lower elevation areas having greater storm surge height, as seen in the A to A′ and B to B′ transects.However, Figure 10 shows that coastal areas are not automatically vulnerable to storm surges if the area's elevation exceeds the maximum storm surge height.

Effect of Elevation on Storm Surge Velocity
Determining flood velocity is also important for early warning efforts, because even low-level fast flowing floods are likely to cause significant damage.In addition to flow depth, the FLO-2D model can also map flow velocity, inundation duration, time to maximum depth, impact force, static pressure, specific energy, and hazard level.
Figure 11 shows a noticeable flow velocity increase that can be seen in the circled area.The circled area's elevation is shown in the boxed graph.The elevation graph shows that maximum flow velocity increased in the raised section in the middle, which can be explained by fluid behavior in that flow velocity is greater in areas that would cause fluids to compress (i.e., a corresponding velocity increase in smaller areas with a constant fluid volume).

Effect of Coastal Shape on Storm Surge Flow Depth
Coastal formation characteristics can also affect an incoming storm surge's extent and depth.Figures 12 and 13 show the coastal shape's effect on storm surge flow depth and extent.Storm surge flooding extents are farther away in coastal areas with concave coasts than those with straight-line coasts (circled areas), even at the same elevation.Fluid is more likely to be dispersed when hitting a convex coast, whereas fluid hitting a concave coast is likely to accumulate in the center, leading to a greater extent of flooding.

Effect of Coastal Shape on Storm Surge Velocity
We also found storm surge differed in response to the shape of the coast.Figure 14 shows velocity behavior differences in areas with dissimilar coastal characteristics and demonstrates that the maximum velocity increases in the center of catchments with equal distance from storm surge inflows.Vector velocity grids in the boxed areas of Figure 14 show uniform storm surge discharge from all directions, which then accumulate in the center.Thus, the flow velocity increases, because the fluid now has more momentum to flow in a single direction.

Effect of Angle of Approach on Flow Depth
Differences in the typhoon's angle of approach to the coast can also affect the extent of storm surge.Figure 15 shows that extent of inundation was farther in areas that were directly hit by the typhoon.Typhoons that move onshore perpendicular to the coast tend to produce higher storm surges than those moving parallel to or at an oblique angle to the coast [23].[24] also conducted research regarding the response in storm surge height to wind direction for cities exposed at the end of a long inlet, as compared to cities exposed in straight coastline.Changes in storm surge heights and extents were observed by simulating changes in the angle of directional wind.The same methodology was applied by [25] to this event, supporting the results that the orientation of coastal area relative to the typhoon's wind direction is a major factor in simulating expected storm surge.

Scale in meters
Figure 15.Effect of typhoon's approach angle on storm surge flow depth.The area delineated by the red box shows coasts directly perpendicular to the typhoon's wind direction, whereas the blue boxed area shows the coast that is parallel to the wind direction.

Conclusions
We simulated Typhoon Haiyan's storm surge height and inundation on the most affected coastal provinces in the Philippines, using the JMA storm surge model and the FLO-2D flood routing model, respectively.The reliability of the models used was proven through validation using actual storm surge records for Typhoon Haiyan.Note that we only conducted qualitative analysis for the effect of quantitative analysis, because of lack of validation data.
Storm surge model results identified Iloilo, Samar, Leyte, Palawan, and Eastern Samar provinces as having the greatest storm surge heights during Typhoon Haiyan, with the maximum storm surge height as high as 3.9 m.We subsequently used these province catchments for FLO-2D flood modeling.
We analyzed the aggravation of storm surge depth and extent in terms of the coast's topographical characteristics.Analysis shows significant effects of elevation, coastal shape, and inland approach angle on flood depth and storm surge velocity.The detailed findings are as follows:  Flow depth analysis in terms of elevation indicated that deeper flooding events were experienced in areas with lower ground elevations (estuaries, low elevation sand lines), whereas coastal areas with higher elevations were not inundated even when the area was directly located on the coast. Flow velocity analysis in terms of elevation indicated that the maximum velocity increased in the center of catchments with higher elevation and equal distance from storm surge inflows; floods flowed faster in areas that would drive the fluid to compress. Flow depth analysis in terms of coastal shape indicated that some areas are flooded more than others, even with the same elevation, due to differences in coastal shape.This is due to the fact that fluid is more likely to be dispersed when hitting a convex coast, whereas fluid hitting a concave coast is likely to accumulate in the center, leading to a greater extent of flooding.
 Flow velocity analysis in terms of coastal shape revealed that in some regions with same level of elevation, flood velocity appeared to be faster in the central areas of catchment that had inflow parameters in all directions.
 Flow depth analysis in terms of a typhoon's angle of approach indicated that extents are farther in the areas that are directly hit by the typhoon; coastlines that are perpendicular to the to the typhoon's directional approach displayed a greater tendency to produce storm surges to a greater extent than those that are parallel to the coast.
Through these analyses, we can conclude that not all coastal areas are vulnerable to storm surges.Some regions may be more inundated than others depending on their physical inland characteristics, which can be used to improve early warning systems and long-term coastal urban planning.


are the x and y components of wind stress on the sea surface, and bx  by  are the x and y components of bottom friction stress.Wind stresses are expressed as( , ) : (1) 114° E-135° E, 5° N-25° N; (2) 118° E-128° E, 8° N-24° N; and (3) 115° E-160° E, 5° N-25° N. Domain 1 covers the storm surge rise and fall over the Philippine Area of Responsibility (region in the northwestern pacific where the Philippines' national meteorological agency monitors significant weather disturbances), domain 2 identifies the regions that recorded the highest sea level rise during Typhoon Haiyan, and domain 3 shows Typhoon Haiyan's life span from cyclogenesis to dissipation.

Figure 1 .
Figure 1.Domains used for simulating the storm surge height using JMA (Japan Meteorological Agency) storm surge model.

Figure 2 .
Figure 2. Most affected coastal areas (upper panel) and their respective simulated storm surge time series results (lower panel).

Figure 3 .
Figure 3. Maximum water surface depth and maximum flow depth in m as shown in the FLO2D model.

Figure 4 .
Figure 4. Storm surge inundation and validation points of Tanauan and Tolosa, Leyte.Validation points: blue circle (observed flooding) and red circle (no flooding).

Figure 5 .
Figure 5. Storm surge inundation and validation points of El Nido, Palawan.Validation points: blue circle (observed flooding) and red circle (no flooding).

Figure 6 .
Figure 6.Storm surge inundation and validation points of Carles, Iloilo.Validation points: blue circle (observed flooding) and red circle (no flooding).

Figure 7 .
Figure 7. Storm surge inundation and validation points of Marabut, Samar.Validation points: blue circle (observed flooding) and red circle (no flooding).

Figure 8 .
Figure 8. Storm surge inundation and validation points of Guiuan, Eastern Samar.The upper left figure shows the transect elevation for a cross-sectional transect (blue arrow) in the storm surge inundation map (right image).The lower left image shows accumulated rainfall during the event in mm/day, the area circled in purple is the area shown on the right panel.

Figure 9 .Figure 10 .
Figure 9. Storm surge inundation (left panel) and cross-sectional x height elevation (red line in right panel); storm surge flooding (black line) in terms of elevation for transects A to A′ and B to B′ (blue solid lines in the left panel,) showing effect of elevation on storm surge height.

Figure 11 .
Figure 11.Map of maximum storm surge velocity (ms −1 ).Storm surge inland flow direction is shown with blue arrows.Inset panel shows surface elevation versus storm surge height, with velocity of transect (red arrow) increasing as ∆H decreases.

Figure 12 .
Figure 12.Effect of coastal shape on flow depth and extent in areas with the same elevation (Marabut, Samar).The right panel shows inland inundation, whereas the left panel shows the elevation contour (pink line is flooding extent), proving that the inundated regions (circled areas) are low-lying coasts.

Figure 13 .
Figure 13.Effect of coastal shape on flow depth and extent in areas with the same elevation (Ajuy, Iloilo).The right panel shows inland inundation, whereas the left panel shows elevation contour (pink line is flooding extent), proving that inundated regions are low-lying coasts.

Figure 14 .
Figure 14.Maximum velocities (ms −1 ) at constant elevation, blue arrows showing inland flow direction.Lower panels display zoomed-in regions showing directional velocity of flow, green and yellow areas having higher velocity flow compared with the blue ones.