Simulation of Typhoon-Induced Storm Tides and Wind Waves for the Northeastern Coast of Taiwan Using a Tide–Surge–Wave Coupled Model

: The storm tide is a combination of the astronomical tide and storm surge, which is the actual sea water level leading to ﬂooding in low-lying coastal areas. A full coupled modeling system (Semi-implicit Eulerian-Lagrangian Finite-Element model coupled with Wind Wave Model II, SELFE-WWM-II) for simulating the interaction of tide, surge and waves based on an unstructured grid is applied to simulate the storm tide and wind waves for the northeastern coast of Taiwan. The coupled model was driven by the astronomical tide and consisted of main eight tidal constituents and the meteorological forcings (air pressure and wind stress) of typhoons. SELFE computes the depth-averaged current and water surface elevation passed to WWM-II, while WWM-II passes the radiation stress to SELFE by solving the wave action equation. Hindcasts of wind waves and storm tides for ﬁve typhoon events were developed to validate the coupled model. The detailed comparisons generally show good agreement between the simulations and measurements. The contributions of surge induced by wave and meteorological forcings to the storm tide were investigated for Typhoon Soudelor (2015) at three tide gauge stations. The results reveal that the surge contributed by wave radiation stress was 0.55 m at Suao Port due to the giant offshore wind wave (exceeding 16.0 m) caused by Typhoon Soudelor (2015) and the steep sea-bottom slope. The air pressure resulted in a 0.6 m surge at Hualien Port because of an inverted barometer effect. The wind stress effect was only slightly signiﬁcant at Keelung Port, contributing 0.22 m to the storm tide. We conclude that wind waves should not be neglected when modeling typhoon-induced storm tides, especially in regions with steep sea-bottom slopes. In addition, accurate tidal and meteorological forces are also required for storm tide modeling.


Introduction
Typhoon-induced storm surge and wind waves are major forces that pose a potential hazard in the form of coastal inundation and to shipping routes. Storm surge is classified as long gravity waves, whereas wind wave are short-period waves. Only long gravity waves can ascend to land, causing inundation that leads to severe loss of life and substantial damage to infrastructure in low-lying areas near the coast. Although short-period wind waves could be dangerous over the ocean, they would break at shallow waters and cannot climb land. Storm tide is the actual level of the water in the sea and adjacent tidal waterways during a typhoon event. The storm tide level results from the additive combination of the normal astronomical tide and the increase in water level due to the storm surge. The total surge depends on the interaction of wave forcing, wind stress, inverted barometer

Description of Study Area
The computational domain in the present study covers the region within a longitude of 114 • to 130 • E and latitude of 19 • to 29 • N. This area is composed of the western Pacific Ocean, the Taiwan Strait, the South China Sea and the East China Sea, which are located in the east, the west, the south and the north of Taiwan, respectively (Figure 1a). The bathymetric data used in the present study were obtained from coastal digital elevation models (DEMs) and ETOPO1 [23], a global relief model. ETOPO1 integrates land topography and the ocean bathymetry of Earth's surface with a 1-arc-minute resolution. It is quite clear in Figure 1a  wave-induced radiation stress was usually neglected in early storm surge modeling [8][9][10][11][12][13][14][15][16]. The contributions of wind waves to the momentum equation were replaced by adding an additional drag coefficient to the wind stress equation. Although this simple method is convenient in storm surge modeling, the momentum transfer rate from wave breaking is highly dependent on the slope and depth of the sea-bottom [17]. Fully coupled wave-current models based on structured or unstructured grids have been recently developed for storm tide modeling despite the high computational cost [6,7,[18][19][20][21][22].
To better understand the wave effects on storm tide modeling, the tide-surge-wave coupled modeling system based on unstructured grids was implemented for the Taiwan Waters in this study, and fully coupled simulations were conducted for storm tides and wind waves using an unstructured-grid finite-element model. This paper is organized as follows. The study site and unstructured grid system for the numerical model are described in Section 2. A brief outline of models used in the present study is given in Section 3. Section 4 describes the model validation for typhoon-induced storm tides and waves. In Section 5, numerical experiments were conducted to investigate the contributions of wave-, wind stress-and air pressure-induced surge to storm tides. The conclusions are summarized in the last section.

Description of Study Area
The computational domain in the present study covers the region within a longitude of 114° to 130° E and latitude of 19° to 29° N. This area is composed of the western Pacific Ocean, the Taiwan Strait, the South China Sea and the East China Sea, which are located in the east, the west, the south and the north of Taiwan, respectively (Figure 1a). The bathymetric data used in the present study were obtained from coastal digital elevation models (DEMs) and ETOPO1 [23], a global relief model. ETOPO1 integrates land topography and the ocean bathymetry of Earth's surface with a 1-arc-minute resolution. It is quite clear in Figure 1a that the sea-bottom elevations along the eastern coast of Taiwan are very steep. The bottom elevations vary rapidly from tens of meters at the shoreline to several thousand meters at 10 km off the coast. With respect to the north coast, the slopes of the sea-bottom are relatively flatter. A total of 278,630 non-overlapping unstructured triangular cells and 142,041 nodes were used in the horizontal to fit the complex shoreline of Taiwan and its adjacent small islands. The resolutions of the mesh range from 30 km to 300 m. The fine meshes are along the coastline of Taiwan and small offshore islands, while the coarse meshes are along the open ocean boundaries (Figure 1b).

Methodology
The effects of the wave forcing and meteorological forcing induced by typhoons on storm tide modeling are taken into account using the coupled modeling system for wind waves, tide and storm surge. The numerical tide-surge-wave modeling system used in the present study consists of a hydrodynamic model (Semi-implicit Eulerian-Lagrangian Finite-Element model, SELFE [24]) and a third-generation spectral wind wave model (Wind Wave Model II, WWM-II [25]). The same computational domain and unstructured-grid system are shared in the coupled model SELFE-WWM-II [26], which avoids the need for data interpolation between the two models and ensures parallel efficiency. SELFE computes the depth-averaged current and water surface elevation passed to WWM-II, while WWM-II passes the radiation stress to SELFE by solving the wave action equation. Figure 2 shows a schematic diagram of the coupled modeling system. More details about the processes of model coupling and model validation with analytical solutions can be found in [26].

The Hydrodynamic Model
The hydrodynamic model, SELFE, coupled with WWM-II, is used to compute tides and surges. SELFE is a semi-implicit Eulerian-Lagrangian finite-element model that has been widely used for

Methodology
The effects of the wave forcing and meteorological forcing induced by typhoons on storm tide modeling are taken into account using the coupled modeling system for wind waves, tide and storm surge. The numerical tide-surge-wave modeling system used in the present study consists of a hydrodynamic model (Semi-implicit Eulerian-Lagrangian Finite-Element model, SELFE [24]) and a third-generation spectral wind wave model (Wind Wave Model II, WWM-II [25]). The same computational domain and unstructured-grid system are shared in the coupled model SELFE-WWM-II [26], which avoids the need for data interpolation between the two models and ensures parallel efficiency. SELFE computes the depth-averaged current and water surface elevation passed to WWM-II, while WWM-II passes the radiation stress to SELFE by solving the wave action equation. Figure 2 shows a schematic diagram of the coupled modeling system. More details about the processes of model coupling and model validation with analytical solutions can be found in [26].

Methodology
The effects of the wave forcing and meteorological forcing induced by typhoons on storm tide modeling are taken into account using the coupled modeling system for wind waves, tide and storm surge. The numerical tide-surge-wave modeling system used in the present study consists of a hydrodynamic model (Semi-implicit Eulerian-Lagrangian Finite-Element model, SELFE [24]) and a third-generation spectral wind wave model (Wind Wave Model II, WWM-II [25]). The same computational domain and unstructured-grid system are shared in the coupled model SELFE-WWM-II [26], which avoids the need for data interpolation between the two models and ensures parallel efficiency. SELFE computes the depth-averaged current and water surface elevation passed to WWM-II, while WWM-II passes the radiation stress to SELFE by solving the wave action equation. Figure 2 shows a schematic diagram of the coupled modeling system. More details about the processes of model coupling and model validation with analytical solutions can be found in [26].

The Hydrodynamic Model
The hydrodynamic model, SELFE, coupled with WWM-II, is used to compute tides and surges. SELFE is a semi-implicit Eulerian-Lagrangian finite-element model that has been widely used for

The Hydrodynamic Model
The hydrodynamic model, SELFE, coupled with WWM-II, is used to compute tides and surges. SELFE is a semi-implicit Eulerian-Lagrangian finite-element model that has been widely used for the simulation of tsunami propagation [27], water quality and ecosystem dynamics [28,29], oil spill dynamics [30], generating inundation maps [31,32] and modeling inundation induced by extreme river flows, typhoons and hurricanes [33,34]. The SELFE model used in the present study is a two-dimensional vertically integrated version with a barotropic mode (SELFE-2DH). The governing equations in two-dimensional form are given as: where η(x, y, t) is the free-surface elevation; H = η + h, h is the bathymetric depth; u(x, y, t) and v(x, y, t) are the horizontal velocity in the x, y direction, respectively; f is the Coriolis factor; g is the acceleration due to gravity;ψ is the earth's tidal potential; α is the effective earth elasticity factor; ρ 0 is the reference density of water; and P A (x, y, t) is the atmospheric pressure at the free surface. τ sx and τ sy are the wind stress in the x, y direction, which can be expressed as: where C s is the wind drag coefficient; ρ a is the air density; and W x , W y are the wind speed in x, y components. Although C s is often modeled as an increasing function of wind velocity, Powell et al. [35] suggested that C s should be limited at high wind speed. The formula for calculating the C s in SELFE-2HD is as: (0.61 + 0.063 × 6.0), W < 6.0 (0.61 + 0.063W), 6.0 ≤ W ≤ 50.0 (0.61 + 0.063 × 50.0), W > 50.0 (6) where W is the wind speed at 10 m above the sea surface. τ rx and τ ry are the wave-induced stresses in the x, y direction, which can be computed following [36]: where S xx , S xy and S yy are the wave radiation stresses and are represented according to [37]: Water 2017, 9, 549 5 of 24 where N is the wave action, σ is the wave relative angular frequency, θ is the wave direction, and C g and C p are the wave group velocity and wave phase velocity, respectively. τ bx and τ by are the bottom shear stress in the x, y direction and are computed by following formula: where C d is the bottom drag coefficient, and can be parameter where n is the Manning coefficient and will be determined through the model validation. A time step of 120 s was used for the hydrodynamic model in simulations with no signs of numerical instability for the present model grid system. The time step of WWM-II is 600 s, i.e., SELFE and WWM-II exchange information every five hydrodynamic time steps. Chen et al. [38] used the coupled ADCRIC+SWAN model to predict the storm surges and wind waves in Mobile Bay, Alabama. They exchanged the hydrodynamic model information with the wind wave model every 30 min and suggested that model results are insensitive to increases in the exchange frequency.

The Spectral Wind Wave Model
The governing equation for the third-generation spectral wind wave model, WWM-II, is the wave action equation, given as:

∂N ∂t
where N is the wave action, C gx and C gy are the wave group velocity in the x, y direction, u and v are the horizontal velocity in the x, y direction, σ is the wave relative angular frequency, θ is the wave direction, C σ and C θ are the propagation velocity in σ, θ space, and S tot is the sum of the source terms for wave variance. The maximum wave direction in WWM-II is 360 • , and this measure is discretized into 36 bins. The low-and high-frequency limits of the discrete wave period are 0.03 and 1.0 Hz, respectively, and are divided into 36 bins. The bottom friction is set equal to 0.067 based on the formulation of JONSWAP [39]. The wave breaking in shallow water area is computed in WWM-II using the method presented by [40] with a constant wave breaking coefficient of 0.78.

Global Model for the Prediction of Ocean Tides
For many practical applications of modeling in coastal environments, accurate predictions of tidal currents and elevations are always indispensable. Due to the necessity of simulating the interaction between astronomical tide and storm surge, the driving forces at the open boundaries of the tide-surge model are tidal elevations. In the present study, a global ocean tidal prediction model, TPXO, developed by Oregon State University, was adopted to extract the tidal harmonic constants and then specified at the open boundaries of SELFE-2HD for simulating tidal propagation. TPXO uses inverse theory and measured data assimilated from tide gauges and the TOPEX/Poseidon satellite. This model finds an optimum balance between observations and hydrodynamics. The extracted tidal harmonic constants provided complex amplitudes of the earth-relative sea-surface elevation for eight primary (M 2 , S 2 , N 2 , K 2 , K 1 , O 1 , P 1 and Q 1 ) tidal constituents. The detailed methodologies used to compute the tides in TPXO can be found in [41,42].

Parametric Cyclone Model
The meteorological boundary conditions for storm surge calculation consist of the wind speed and atmospheric pressure of typhoons. A common way to generate wind and air pressure fields for typhoons is to reconstruct them using the analytical parametric cyclone model. Over the years, many parametric cyclone models have been developed to provide meteorological information for storm surge modeling [43][44][45][46][47][48]. Jakobsen and Madsen [49] investigated and compared parametric cyclone models based on cyclone position, central pressure, maximum wind speed and radius to maximum wind speed. They found that the analytical models provide very similar air pressure and tangential wind speed distributions for a cyclone. Therefore, a parametric cyclone model presented by Holland [45] was employed: where P A is the air pressure; P n is the ambient pressure; P c the central air pressure of the typhoon; R max is the radius to the maximum wind speed; r is the radial distance from the typhoon center; W H is the wind speed; ρ a is the air density; f is the Coriolis factor; and B is a parameter that characterizes the scale of the typhoon. The wind field was modified from gradient height to 10 m above sea surface with a wind-reducing coefficient of 0.85 following [50]. The formula of B = 1.5 + (980 − P c )/120 presented by Hubbert et al. [51] was adopted. With regard to the estimation of R max , the formula used to simulate the storm surge in the northwest Pacific Ocean [52] was adopted in this study: where ∆P = P n − P c ; c 0 , c 1 and ε were set to be 5.0377, −0.0232 and 0.4502, respectively [52].
Holland model is an axisymmetric model, the asymmetric structure of a typhoon cannot be represented by the model. Highly asymmetric structures in a landfalling typhoon often lead to large errors in storm surge forecasting. To overcome this limitation, Georgiou [53] introduced a hurricane translation speed (W t ) to Holland model and was used in the present study: W is the total wind speed; W t is the typhoon translation speed; W x , W y are the wind speed in x, y components; θ 1 is the angle from the direction of the hurricane movement; θ 2 is the azimuthal angle with respect to the typhoon's eye; and θ 3 is the inflow angle at the surface. Although the local wind and surge estimates are less sensitive to inflow angle [50], the inflow angle θ 3 was considered in the present study, and can be described as a function of r [54]: Table 1 lists the parameters used in SELFE-WWM-II and parametric cyclone model.

Indicators of Model Performance
Three criteria are adopted to assess the model performance for water level simulation, including the mean absolute error (MAE), the root mean square error (RMSE), and percent bias (PBIAS) [55]. The optimal value of PBIAS is 0, with low-magnitude values indicating accurate model simulations. Positive values of PBIAS indicate that the model was overestimating while negative values indicate that the model was underestimating [56]. The equations for these three criteria are shown as follows: η s i is the simulations and η m i is the measurements.

Model Validation
The hindcasts of wave height and storm tide induced by five typhoon events, including Typhoon Saola (2012), Typhoon Soulik (2013), Typhoon Matmo (2014), Typhoon Soudelor (2015) and Typhoon Dujuan (2015), were conducted to validate the coupled model. The paths of these five typhoons were crossed or passed Taiwan from the east to the west or from the southeast to the northwest. Besides, the 10 min average wind speed near the typhoon center were in the range of 38 m/s to 51 m/s, and were classified as the moderate to the severe typhoon (according to the categories of Central Weather Bureau in Taiwan). Therefore, the surges induced by these five typhoons would be significant along the northeastern of Taiwan.
The tracks of these five typhoons and their corresponding central air pressures are shown in Figure 3. Three wave buoys (Keelung buoy, Suao buoy and Hualien buoy) and three tide gauge stations (Keelung Port, Suao Fish Port and Hualien Port) located along the northeastern coast of Taiwan were selected to compare wave heights and sea water levels between the simulations and observations. Figure 3 also illustrates the locations of wave buoys and tide gauge stations. The historical data used to construct the wind and air pressure fields for the five typhoon events were obtained from The Regional Specialized Meteorological Center (RSMC) Tokyo-Typhoon Center Best Track database. Hourly data for wave height and sea water level measured at three wave buoys and three tide gauge stations were provided by the Harbor and Marine Technology Center and the Central Weather Bureau in Taiwan.
historical data used to construct the wind and air pressure fields for the five typhoon events were obtained from The Regional Specialized Meteorological Center (RSMC) Tokyo-Typhoon Center Best Track database. Hourly data for wave height and sea water level measured at three wave buoys and three tide gauge stations were provided by the Harbor and Marine Technology Center and the Central Weather Bureau in Taiwan.  Figure 4 presents the comparisons of the modeled and observed water levels at Keelung Port, Suao Fish Port and Hualien Port for the five typhoon events. Although the water level hydrographs show that the simulations have good agreements with the observations in both astronomical tides and storm tides, the model results sometimes underestimate or overestimate the water level. These discrepancies may be due to the inaccuracy of the wind field and air pressure field created by the parametric cyclone model. The parametric cyclone model has difficulty representing the accurate structure of typhoons at certain locations and is not expected to reproduce the wind conditions far away from the center [54]. Feng et al. [57] also suffered from a similar problem when Jelesnianski's circular hurricane model was employed to simulate storm surge at Tianjin in China. The statistical errors for the differences between the simulated and observed water levels at three tidal gauge stations can be found in   Figure 4 presents the comparisons of the modeled and observed water levels at Keelung Port, Suao Fish Port and Hualien Port for the five typhoon events. Although the water level hydrographs show that the simulations have good agreements with the observations in both astronomical tides and storm tides, the model results sometimes underestimate or overestimate the water level. These discrepancies may be due to the inaccuracy of the wind field and air pressure field created by the parametric cyclone model. The parametric cyclone model has difficulty representing the accurate structure of typhoons at certain locations and is not expected to reproduce the wind conditions far away from the center [54]. Feng et al. [57] also suffered from a similar problem when Jelesnianski's circular hurricane model was employed to simulate storm surge at Tianjin in China. The statistical errors for the differences between the simulated and observed water levels at three tidal gauge stations can be found in

Validation of Wave Height and Wave Period
Due to the scarcity of measured wave direction date at the Keelung buoy, Suao buoy and Hualien buoy and measured wave period data at the Suao buoy and Hualien Buoy, model-observation comparisons of wave parameters were only conducted for wave height at the Keelung buoy, Suao buoy and Hualien buoy and for wave period at the Keelung buoy. It should be noted that the simulated and measured wave height and wave period are significant wave height and peak period. These two wave parameters are often specified in the measurement of sea states.
The comparisons of wave height between the simulation and observations at the Keelung buoy, Suao buoy and Hualien buoy for the five typhoon events are shown in Figure 5. The wave heights increased as the typhoon approached and decreased when typhoon was far away. The largest measured wave height exceeded 16 m at the Suao buoy during Typhoon Soudelor (2015) and Typhoon Dujuan (2015) (Figure 5d,e). The wave height measured at the Keelung buoy was almost 10 m during Typhoon Soulik (2013), while wave heights of 7.5 m were recorded during Typhoon Matmo (2014) and Typhoon Soudelor (2015) at Hualien buoy. The coupled model can only catch the peak values well but fails to mimic the measured wave heights around peak waves because the wind fields are not very accurate. Figure 6 represents the results of comparing the simulated and measured wave period at the Keelung buoy. As the typhoon passed, the wave periods also gradually increased. The maximum wave period is approximately 13 s during the period of Typhoon Soudelor (2015). Table 2

Validation of Wave Height and Wave Period
Due to the scarcity of measured wave direction date at the Keelung buoy, Suao buoy and Hualien buoy and measured wave period data at the Suao buoy and Hualien Buoy, model-observation comparisons of wave parameters were only conducted for wave height at the Keelung buoy, Suao buoy and Hualien buoy and for wave period at the Keelung buoy. It should be noted that the simulated and measured wave height and wave period are significant wave height and peak period. These two wave parameters are often specified in the measurement of sea states.
The comparisons of wave height between the simulation and observations at the Keelung buoy, Suao buoy and Hualien buoy for the five typhoon events are shown in Figure 5. The wave heights increased as the typhoon approached and decreased when typhoon was far away. The largest measured wave height exceeded 16 m at the Suao buoy during Typhoon Soudelor (2015) and Typhoon Dujuan (2015) (Figure 5d,e). The wave height measured at the Keelung buoy was almost 10 m during Typhoon Soulik (2013), while wave heights of 7.5 m were recorded during Typhoon Matmo (2014) and Typhoon Soudelor (2015) at Hualien buoy. The coupled model can only catch the peak values well but fails to mimic the measured wave heights around peak waves because the wind fields are not very accurate. Figure 6 represents the results of comparing the simulated and measured wave period at the Keelung buoy. As the typhoon passed, the wave periods also gradually increased. The maximum wave period is approximately 13 s during the period of Typhoon Soudelor (2015). Table 2 summarizes the statistical errors for wave height and wave period. The maximum wave period MAE, RMSE and PBIAS are 1.27 m, 1.56 m and 12.41% at the Suao buoy for wave height and 1.31 s, 1.51 s and −9.62% at the Keelung buoy, respectively. Overall, our tide-surge-wave coupled model demonstrates a great performance in simulating storm tides and wave heights. The best-simulated wave height, wave period and water level results were produced with a Manning coefficient set equal to 0.025. The Manning coefficient n was a constant in the model due to the scarcity of type of the sea-bottom material. However, the bottom drag coefficient C d is varied with H.
The spatial distributions of wave height predicted by WWM-II during Typhoon Soudelor (2015) are shown in Figure 7. The extent of wave height over 18 m increases as the typhoon moved northwestward. The maximum wave heights appear in the right forward quadrant of the typhoon center. This is due to the enhancement of wind velocity along the right side of typhoon by adding the typhoon translation speed [58].
Water 2017, 9,549 12 of 24 northwestward. The maximum wave heights appear in the right forward quadrant of the typhoon center. This is due to the enhancement of wind velocity along the right side of typhoon by adding the typhoon translation speed [58].

Results and Discussion
The storm tide is caused by several components, including the astronomical tide, the low-pressure of typhoons, the wind stress of typhoons, the wave setup and the Coriolis force [59]. In the present study, three numerical experiments (as shown in Table 3) were conducted to investigate the contribution of air pressure, wind stress and waves to the storm tide at three tide gauge stations (Figure 3) during Typhoon Soudelor (2015). The baseline is a full coupled modeling; Case 1 is a decoupled modeling in which only the wave-induced radiation stress was not included; Case 2 is a decoupled modeling in which the wave-induced radiation stress and wind stress were not included; and Case 3 is a decoupled modeling in which the wave-induced radiation stress and air pressure were not included. Table 3. Designed numerical experiments for investigating the factors contributing to surge.

Conditions Investigated Factors
Baseline Simulation with air pressure, wind stress and wave None Case 1 Simulation only with air pressure, wind stress Wave Case 2 Simulation only with air pressure Wind stress Case 3 Simulation only with wind stress Air pressure

Effects of Waves on Storm Tides
The wave-induced surge (wave setup, η wave ) was computed by subtracting the results of Case 1 from that of the baseline. Figure 8 present the comparison of the simulated water level between the baseline and Case 1 at three tide gauge stations, and the time series of η wave is also depicted. The contributions of η wave to the storm tide are quite different at the three tide gauge stations. The η wave appears to be non-significant at Keelung Port and Hualien Port (Figure 8a,c), but it is nonetheless significant, reaching 0.55 m at Suao Port (Figure 8b). The spatial distributions of wave height and η wave are shown in Figure 9. The η wave can be as high as 1.0 m along the coastline in relatively shallow areas, (Figure 9b) because of the offshore giant wave (Figure 9a). Huang et al. [18] found that the wave-induced forces result in an additional 0.3-0.5 m of surge relative to surge-only simulation for an Ivan-like hurricane impacting Tampa Bay, Florida. Bertin et al. [7] concluded that the maximum wave setup exceeds 0.3 and 0.4 m in the Bay of Biscay, France during the extra-tropical storms Xynthia and Joachim. The result obtained in this study is similar to their findings and reveals that the wave-induced surge is significant and essential for storm tide modeling when large wave heights occur in a region with a steep sea-bottom slope [17,[60][61][62].

Results and Discussion
The storm tide is caused by several components, including the astronomical tide, the low-pressure of typhoons, the wind stress of typhoons, the wave setup and the Coriolis force [59]. In the present study, three numerical experiments (as shown in Table 3) were conducted to investigate the contribution of air pressure, wind stress and waves to the storm tide at three tide gauge stations ( Figure 3) during Typhoon Soudelor (2015). The baseline is a full coupled modeling; Case 1 is a decoupled modeling in which only the wave-induced radiation stress was not included; Case 2 is a decoupled modeling in which the wave-induced radiation stress and wind stress were not included; and Case 3 is a decoupled modeling in which the wave-induced radiation stress and air pressure were not included. Simulation only with wind stress Air pressure

Effects of Waves on Storm Tides
The wave-induced surge (wave setup, wave  ) was computed by subtracting the results of Case 1 from that of the baseline. Figure

Effects of Wind Stress on Storm Tides
The simulations of Case 2 were compared against those of Case 1 to evaluate the contribution of wind-induced surge ( wind  ) to the storm tide ( Figure 10

Effects of Wind Stress on Storm Tides
The simulations of Case 2 were compared against those of Case 1 to evaluate the contribution of wind-induced surge ( wind  ) to the storm tide ( Figure 10).

Effects of Wind Stress on Storm Tides
The simulations of Case 2 were compared against those of Case 1 to evaluate the contribution of wind-induced surge (η wind ) to the storm tide ( Figure 10). A larger η wind of almost 0.22 m can be found at Keelung Port (Figure 10a), while weak contributions of η wind are at Suao Port and Hualien Port (Figure 10b,c). A simple one-dimensional equation (as shown in Equation (24)) can be used to represent the η wind [57]: As the total water depth (H) becomes shallow, the wind-induced surge (η wind ) is larger for the same wind stress. The sea-bottom elevations in the region of Keelung Port are shallower compared to those around Suao Port and Hualien Port (Figure 1a). The wind direction was onshore at Keelung Port during the approach of Typhoon Soudelor (2015) (Figure 11a,b). The effect of wind stress on storm tide would be significant in the shallow water area with onshore wind. According to the results derived in Sections 5.1 and 5.2, it is concluded that the typhoon-induced wind wave would be larger along coast with the steep sea-bottom slope, for instance in the eastern coast of Taiwan, while the wind-induced surge would be smaller there.
Water 2017, 9,549 18 of 24 ( Figure 10b,c). A simple one-dimensional equation (as shown in Equation (24)) can be used to represent the wind  [57]: As the total water depth ( H ) becomes shallow, the wind-induced surge ( wind  ) is larger for the same wind stress. The sea-bottom elevations in the region of Keelung Port are shallower compared to those around Suao Port and Hualien Port (Figure 1a). The wind direction was onshore at Keelung Port during the approach of Typhoon Soudelor (2015) (Figure 11a,b). The effect of wind stress on storm tide would be significant in the shallow water area with onshore wind. According to the results derived in Sections 5.1 and 5.2, it is concluded that the typhoon-induced wind wave would be larger along coast with the steep sea-bottom slope, for instance in the eastern coast of Taiwan, while the wind-induced surge would be smaller there.     is 0.6 m at Hualien Port but is less than 0.1 m at Keelung Port and Suao Port. Hualien Port is the station closest to the center of the typhoon when Typhoon Soudelor (2015) made landfall and was thus exposed to the lowest air pressure. The air_pre.


can also be estimated using a simple formula based on the inverted barometer effect [60]:  Figure 12 illustrates the comparisons of the Case 1 and the Case 3 simulations and the surge induced by air pressure (η air_pre. ). The largest η air_pre. is 0.6 m at Hualien Port but is less than 0.1 m at Keelung Port and Suao Port. Hualien Port is the station closest to the center of the typhoon when Typhoon Soudelor (2015) made landfall and was thus exposed to the lowest air pressure. The η air_pre. can also be estimated using a simple formula based on the inverted barometer effect [60]:

Effects of Air Pressure on Storm Tides
The air pressure was approximately 955 hpa around Hualien Port (Figure 11c), and therefore the η air_pre. = 1.04 × (1013 − 955) = 60.32 cm = 0.6032 m. This answer is very close to the result from the numerical model (as shown in Figure 12c). . This answer is very close to the result from the numerical model (as shown in Figure 12c).

Future Work
This is a preliminary study for the simulation of storm tide along the coast of Taiwan using a coupled tide-surge-wave modeling system. A further application modeling the coastal inundation associated with typhoons will be performed using high-spatial resolution grids to properly represent the natural barriers and seawalls. Using a three-dimensional coupled tide-surge-wave modeling system would likely be close to reality [63] and allow a proper representation of the vertical circulation that occurs in breaking zones [7]. Therefore, a three-dimensional coupled model should be applied to forecast the astronomical tide, storm tide and wave heights along the coast of Taiwan. Besides, the parametric cyclone model is insufficient for an accurate representation of the structure of a typhoon. Therefore the numerical weather model, such as WRF (Weather Research and Forecasting) model, is essential for storm tide modeling or forecasting in the near future.

Conclusions
A two-dimensional hydrodynamic (SELFE-2DH) and wind wave (WWM-II) coupled modeling system (SELFE-WWM-II) was applied to hindcast the storm tide and waves induced by typhoons for the northeastern coast of Taiwan. The detailed comparisons generally show good agreement between the model simulations and the available data measured at tide gauge stations and wave buoys. The numerical experiments were designed to analyze the contributions of waves, wind stress

Future Work
This is a preliminary study for the simulation of storm tide along the coast of Taiwan using a coupled tide-surge-wave modeling system. A further application modeling the coastal inundation associated with typhoons will be performed using high-spatial resolution grids to properly represent the natural barriers and seawalls. Using a three-dimensional coupled tide-surge-wave modeling system would likely be close to reality [63] and allow a proper representation of the vertical circulation that occurs in breaking zones [7]. Therefore, a three-dimensional coupled model should be applied to forecast the astronomical tide, storm tide and wave heights along the coast of Taiwan. Besides, the parametric cyclone model is insufficient for an accurate representation of the structure of a typhoon. Therefore the numerical weather model, such as WRF (Weather Research and Forecasting) model, is essential for storm tide modeling or forecasting in the near future.

Conclusions
A two-dimensional hydrodynamic (SELFE-2DH) and wind wave (WWM-II) coupled modeling system (SELFE-WWM-II) was applied to hindcast the storm tide and waves induced by typhoons for the northeastern coast of Taiwan. The detailed comparisons generally show good agreement between the model simulations and the available data measured at tide gauge stations and wave buoys. The numerical experiments were designed to analyze the contributions of waves, wind stress and air pressure to the storm tide associated with Typhoon Soudelor (2015) using a well-calibrated SELFE-WWM-II model. The results reveal that the wave-induced surge could reach to 0.55 m at Suao Port due to the offshore giant wind wave (exceeding 16.0 m) and its steep sea-bottom slope. The effects of wind stress were only slightly significant at Keelung Port, contributing 0.22 m to the storm tide as a result of shallow water and onshore wind. The air pressure results in a 0.6 m surge at Hualien Port because of the inverted barometer effect. It must be noted that the influence of waves, wind stress and air pressure on storm tides is also highly dependent on the track of the typhoon. However, it is concluded that waves should not be omitted in modeling typhoon-induced storm tides. In addition, accurate tidal and meteorological forces are indispensable in suppressing inaccuracies in storm tide predictions. Moreover, more detailed wave data from offshore buoys are needed to elucidate and verify the simulated wave dynamics in the coastal waters of Taiwan.