1. Introduction
A storm surge involves local sea level oscillation or a nonperiodic abnormal rise or fall in sea level caused by strong winds and sudden changes in the air pressure associated with the passage of storms such as tropical cyclones, temperate weather systems, and oceanic squall lines. In China, the frequency of storm surge disasters triggered by tropical cyclones, especially low-lying coastal area inundation caused by extreme storm surge events, has threatened societies and economies in coastal areas [
1]. In the Pearl River Estuarine area of Guangdong Province, China, the super typhoons Hato and Mangkhut successively made landfall along the coast in 2017 and 2018, leading to widespread storm surge inundations in the coastal cities [
2,
3], causing tremendous economic losses and threatening the sustainable development of the coastal areas.
In the context of climate change, global sea levels have been rising significantly since the late 19th century; global sea levels rose by approximately 14 cm in the 20th century [
4,
5,
6]. According to long-term observation data from coastal tide gauge stations, sea levels along the coast of China have also exhibited an upward trend over the past 40 years [
7]. Based on simulations using the Coupled Model Intercomparison Project Phase 6 (CMIP6) model, simulations have demonstrated that since the 1980s, climate change has led to a noticeable seasonal shift in the global occurrence of super typhoons [
8]. Some studies have shown that under high-emission scenarios, low-lying coastal areas will experience more frequent storm surge disasters due to the rapid rise in sea levels [
9]. According to the statistical results from the coastal tide data in China, estuarine areas, such as the Pearl River Estuary, have a high probability of experiencing strong storm surges and are among the regions with the highest risk for storm surge disasters in coastal China [
10]. Accordingly, conducting simulation studies and predicting storm surge inundations in estuarine areas have great significance for disaster prevention.
Numerical models for storm surges have matured in terms of grid technology, dynamic frameworks, and parameterization schemes [
3,
11]. In the past ten years, with the rapid economic and population development in estuarine areas, simulations of storm surge inundation in estuarine areas have attracted the attention of many scholars, who have established fine grids for small regions (50–300 km) along the coast and developed storm surge inundation simulation and forecast models to support storm surge disaster prevention in coastal areas [
12,
13,
14,
15]. With enhanced computational capabilities, the land resolution in storm surge inundation models has increased from several hundred meters to 10 m, enabling block-level resolution and providing more detailed storm surge inundation simulation results [
16,
17]. However, considerable computational resources are required.
At the same time, the study of multiple influencing factors of storm surges in estuarine areas has also attracted the attention of many scholars. A numerical model (Delft3D), previously calibrated and validated, was used to simulate the lagoon hydrodynamics under different scenarios combining the mean sea level (MSL) rise and frequent storm surge events; moreover, the numerical results demonstrated that an MSL rise could change vertical zonation and threaten local habitats [
18]. When considering spring tides, storm surges, and river inflows under sea-level increases of 0.4 m and 0.9 m, the floodplains in mature estuarine systems may become submerged and experience a considerable increase in inundation depths once a certain threshold in elevation has been exceeded; conversely, immature estuarine systems can experience increases in the relative inundation extent and substantial changes in hydrodynamics such as tidal range and current velocity [
19]. A hydrodynamic model (Delft3D-FM) was used to study the flooding from both the coastal and fluvial influences in the Napa River and its interactions with the San Francisco Bay, and the results showed that large tidal amplitudes diminished the storm surge propagation upstream and that phasing between peak fluvial discharge and high tide was important for predicting the locations and occurrences of the highest water levels; additionally, the interactions between tides, river discharges, and storm surges are complex [
20]. A high-resolution, unstructured-grid, surge-wave fully coupled model and a hybrid typhoon wind model was used to assesses the highest storm tide hazard level along the coast of Taiwan, and this assessment is important to establishing an early warning of riverine inundation if the occurrence of the highest storm surge coincides with the highest astronomical tide [
21]. The Advanced Circulation (ADCIRC) model coupled with the Precipitation-Runoff Modeling System (PRMS) was used to simulate the impact from the three different versions of hurricanes on the Fox Point Hurricane Barrier (FPHB) in the Rhode Island and Narragansett Bay area; moreover, the storm surge resulting from the strongest version during the first landfall exceeded 7 m; this value exceeded the height of the FPHB, and the results all showed that northern Narragansett Bay could be particularly vulnerable to both storm surge- and rainfall-driven flooding, especially if the FPHB suffered a power outage [
22]. A range of meteorological, storm track direction, and storm bearing parameters that produced the highest sustained wind speeds were introduced into meteorological conditions, coupled with storm surges, tides, and waves, for the risk assessment of special infrastructure in Delaware Bay, such as nuclear power plants; the simulations resulted in a maximum stillwater elevation and wave height of 7.5 m and 2.5 m, respectively, which demonstrated that the estimated elevation probability had an annual exceedance of less than 10
−4 [
23].
The above research indicates that multiple factors need to be considered in the study of storm surge inundation in estuarine areas. Due to constraints in the basic data, weather forecasts, and computational capabilities, storm surge inundation forecasting in estuarine areas will become an important focus. To carry out storm surge inundation forecasts in estuarine areas, statistical analyses of the historical storm surge data, disaster data, and typhoon information in the study area need to be carried out and the results used as oceanographic and meteorological background parameters. In addition, the factors influencing storm surge inundations include typhoon landfall intensity, landfall location, landfall direction, typhoon speed, astronomical tide, nearshore waves, and nearshore topography. The existence of many influencing factors can lead to interference in storm surge inundation forecasts. Accordingly, the main focus of this paper includes the establishment of a reasonable method for storm surge inundation experiments and the analysis of the simulation results to highlight the significant factors. The purpose of this study is to reduce the number of factors influencing storm surge inundation through a numerical simulation assessment, establish a lightweight storm surge inundation database in follow-up work, and carry out rapid storm surge inundation forecasting and disaster risk estimation; the results from this study can provide technical support for sustainable development in estuarine areas.
3. Results
To carry out research on storm surge inundations in small regions, it is necessary to conduct statistical analyses of historical storm surge data, disaster data, and typhoon information of the study area, and their results can be used as oceanographic and meteorological background parameters. In addition, the factors affecting storm surge inundations include typhoon landfall intensity, landfall location, landfall direction, typhoon speed, astronomical tides, nearshore waves, and nearshore topography. Therefore, reasonable experiments need to be set up, and the simulation results need to be analyzed.
3.1. Statistics of Historical Typhoons and Storm Surges
In terms of typhoon landfall intensity, from 1949 to 2022, there were 209 typhoons passing within 300 km of the study area, an average of approximately 2.9 typhoons per year [
35,
36]. The typhoons that caused the top ten historical tide levels at the Sanzao tide gauge station were typhoons 8908, 9302, 9316, 9615, 0307, 0814, 0915, 1415, 1713, and 1822. Among them, there were four at the typhoon level, five at the strong typhoon level, and one at the super typhoon level. The typhoons that caused the top two historical tide levels in the study area were typhoons 1822 and 0814; both of these had strong landfall intensities. Evidently, when the maximum wind speed in the typhoon center is greater than 35 m/s, storm surge disasters are likely to occur in the study area.
In terms of typhoon landfall, between 1949 and 2022, there were 183 typhoons that passed and made landfall within 300 km of the study area [
35,
36], accounting for 87.6% of the total typhoons. In terms of the typhoon landfall direction, the primary direction was in the northwest quadrant (W, WNW, NW, NNW, and N), accounting for 84.7% of all landfall typhoons, and the typhoons that caused the top ten historical tide levels in the study area were all in this direction (
Figure 4). In terms of typhoon speed, approximately 52.8% of the typhoons had a speed ranging from 15 to 25 km/h, with an average speed of 19.8 km/h. In terms of typhoon impact time, typhoons influence the study area from April to December each year, with summer (July to September) being the predominant period of influence, accounting for 67.9% of all typhoons.
3.2. Settings for the Numerical Simulations of the Multifactor Storm Surge Inundation
According to the historical storm surge data, forecast experience, and topographic features, the typhoon landfall intensity, landfall direction, landfall location, and astronomical tide are the main factors influencing storm surge inundations, and the measured tide levels at coastal tide gauge stations are an important basis for recording storm surge events. Regardless of the combination of and changes in influencing factors, the measured tide levels at the tide gauge stations can likely be used to determine the storm surge inundation state. This deduction needs to be verified in the subsequent research to this study. Based on the above historical typhoon statistics, we set the typhoon speed to 20 km/h and superimposed the astronomical tide curve from the average value of astronomical high tides for July to September at the Sanzao tide gauge station to establish three experimental schemes of the factors influencing storm surge inundation as follows:
Scheme A: Storm surge inundation simulation with different landfall locations.
The typhoon moves in the west-northwest (WNW) direction with a landfall intensity of 930 hPa. Three landfall location scenarios are set up such that the total tide levels at the Sanzao tide gauge station are similar to each other to assess the differences in storm surge inundation (see
Figure 5a for the typhoon paths).
Scenario A1: The typhoon passes at a distance of 40 km from the south of the study area, and the maximum storm surge is superimposed on the astronomical low tide;
Scenario A2: The typhoon passes at a distance of 70 km from the south of the study area, and the maximum storm surge is superimposed on the astronomical slack tide;
Scenario A3: The typhoon passes at a distance of 110 km from the south of the study area, and the maximum storm surge is superimposed on the astronomical high tide.
Scheme B: Storm surge inundation simulation with different landfall intensities.
The typhoon moves in the west-northwest (WNW) direction and passes 40 km south of the study area. Three typhoon intensity scenarios are set up such that the total tide levels at the Sanzao tide gauge station are similar to each other to assess the differences in storm surge inundation (see
Figure 5b for the typhoon paths).
Scenario B1: The landfall intensity is 930 hPa, and the maximum storm surge is superimposed on the astronomical low tide;
Scenario B2: The landfall intensity is 945 hPa, and the maximum storm surge is superimposed on the astronomical slack tide;
Scenario B3: The landfall intensity is 960 hPa, and the maximum storm surge is superimposed on the astronomical high tide.
Scheme C: Storm surge inundation simulation with different landfall directions.
The typhoon intensity is 945 hPa. Three typhoon landfall direction (west, northwest, and north) scenarios are set up such that the total tide levels at the Sanzao tide gauge station are similar to each other to assess the differences in storm surge inundation (see typhoon tracks in
Figure 5c).
Scenario C1: The landfall direction is westward;
Scenario C2: The landfall direction is northwest;
Scenario C3: The landfall direction is northward.
3.3. Multifactor Numerical Simulation Analysis of Storm Surge Inundation
Based on the above experimental settings, the storm surge inundation model is driven by the wind field of the Holland typhoon model [
40].
Figure 6 and
Table 4 show the diagrams and statistics, respectively, of the storm surge inundation area and depth under Scenarios A1, A2, and A3. When the typhoon landfall intensity and direction remain unchanged, the landfall location is variable, and the highest tide levels at the Sanzao tide gauge station are similar, based on the distribution of the maximum storm surge inundation field (
Figure 6). With the exception of a small difference in the inundation range in the eastern part of the Jitimen estuary, the inundation ranges and depths are similar among the three scenarios. From the difference statistics (
Table 4) under the three scenarios, there is an approximately 5% difference in the inundation area and a less than 5% difference in the average inundation depth; these results indicate minimal variations.
Figure 7 and
Table 5 show the images and statistics of the storm surge inundation areas and depths under scenarios B1, B2, and B3, respectively. When the typhoon landfall location and direction remain unchanged, the landfall intensity varies, and the highest tide levels at the Sanzao tide gauge station are similar, according to the distribution of the maximum storm surge inundation field (
Figure 7); with the exception of a very minor difference in the inundation range in the eastern part of the Jitimen estuary, the inundation areas and depths are very similar under the three scenarios. As for the difference statistics (
Table 5), under the three scenarios, the difference in the inundation area is less than 5%, and the difference in the average inundation depth is less than 5%, indicating minimal variations.
Figure 8 and
Table 6, respectively, show the images and statistics of the storm surge inundation areas and depths under scenarios C1, C2, and C3. When the typhoon landfall location and intensity remain unchanged, the landfall direction is variable, and the highest tide levels at the Sanzao tide gauge station are similar, based on the distribution of the maximum storm surge inundation field (
Figure 8). However, notable differences are observed in the inundation areas and depths under the three scenarios. The inundation extent of the land area to the east of the Jitimen estuary widely varies, with the largest inundation area and average inundation depth occurring under scenario C3, followed by those under scenario C2, and the smallest under scenario C1. According to the difference statistics (
Table 6), the inundation areas under scenarios C2 and C3 are 16.8% and 29.8% larger than that under scenario C1, respectively, and the average inundation depths under scenarios C2 and C3 are 5.8% and 14.7% larger than that under scenario C1, respectively.
The above numerical experiments indicate that the study area has been flooded when the tide level exceeds 3 m at the Sanzao station. Based on
Table 4,
Table 5 and
Table 6, the differences in storm surge inundation among the three scenarios are analyzed in terms of the duration of high tide levels (exceeding 3 m) at the Sanzao station. In scenarios A1, A2, and A3, the durations of high tide levels exceeding 3 m are 3.13, 3.25, and 3.34 h, respectively, with differences of less than 7%. In scenarios B1, B2, and B3, the durations of the high tide levels are 2.79, 2.88, and 3.00 h, respectively, with differences of less than 8%. Evidently, under conditions where the highest tide levels at the Sanzao station and the duration of high tide levels are similar, the differences in the storm surge inundation area and average inundation depth among the three scenarios for typhoon landfall locations (A1, A2, and A3) are very small. Similarly, the differences in storm surge inundation area and average inundation depth among the three scenarios for typhoon landfall intensities (B1, B2, and B3) are also very small. However, in scenarios C1, C2, and C3, the durations of high tide levels exceeding 3 m are 2.83, 3.24, and 3.71 h, respectively. C2 and C3 are 13.1% and 31.1% larger than C1, respectively, showing significant differences. This difference is attributed to the topographic features and the direction of the typhoon landfall.
The Jitimen waterway is a north–south trumpet-shaped estuary, with the mouth facing south. When the typhoon landfall location and intensity remain unchanged and the landfall direction is westward (scenario C1), the study area is mainly affected by easterly winds. The Jitimen estuary is not conducive to the accumulation of storm surges, with a high tide duration of 2.83 h; this duration is the shortest among the three scenarios. When the typhoon landfall direction is northwestward (scenario C2), the study area is mainly affected by southeasterly winds. Compared to that in scenario C1, the Jitimen estuary is more conducive to the accumulation of storm surges, with a high tide duration of 3.24 h; this duration is 13.1% longer than that in scenario C1. When the typhoon landfall direction is northward (scenario C3), the study area is mainly affected by southerly winds; compared to those in scenarios C1 and C2, C3 is most conducive to the accumulation of storm surges, with a high tide duration of 3.71 h; this duration is 33.1% longer than that in scenario C1. In the typhoon landfall direction scheme, a longer high tide duration results in a greater total amount of seawater entering the Jitimen estuary area during the storm surge, leading to a larger inundation area and a correspondingly greater average inundation depth.
4. Discussion and Conclusions
With global climate change, the risk of extreme storms and the resulting storm surges in estuarine areas is increasing. Our numerous research findings indicate that the factors influencing storm surges in estuarine areas are diverse and complex, involving multiple aspects. Storm surges can interact with fragile estuarine ecosystems, river runoff, and seawalls [
18,
20,
22] and can also affect the hydrodynamics and geomorphic features of estuaries [
19]. Research in this field is extensive and beneficial. The factors influencing storm surges in estuarine areas can be divided into two categories: (1) different hydrodynamic influencing factors in estuaries, such as tides, runoff, waves, nearshore structures, and ecosystems, and (2) different meteorological influencing factors, mainly involving typhoon landfall intensity, location, and direction. Scholars have conducted extensive research on the first category of influencing factors and have achieved abundant research results; however, research on storm surge inundation under different meteorological influencing factors is very limited. Research in the latter area has the potential to improve the forecast level of storm surge inundation. The circulation range of typhoons is generally hundreds of kilometers or even larger, while the range of small-scale estuaries is generally less than 200 km. The circulation range of typhoons is much larger than that of small-scale estuaries, and the influencing factors of typhoons can alter the characteristics of storm surges in small-scale estuarine areas. From the perspective of defending against storm surge disasters in estuarine areas, accurate storm surge forecasts are needed before disasters occur, reasonable storm surge decision-making information is needed when disasters occur, and detailed storm surge investigations and assessments are needed after disasters occur. Clearly, forecasting and researching storm surge inundation in estuarine areas are extremely important tasks.
With the advancement of computer technology and the development of storm surge models, numerical forecasting models for storm surges provide favorable support for predicting storm surge inundation in estuarine areas. Storm surge models require high-resolution computational grids (in the order of meters or tens of meters) to characterize the complex terrain of estuarine areas. The computation of a storm surge inundation event takes a long time (from tens of minutes to several hours); thus, storm surge decision-making is negatively affected due to the time-consuming computations with numerical forecasting models. Additionally, typhoon forecast information is updated every few hours, posing a challenge for storm surge inundation model calculations. In contrast to numerical forecasting systems, a multi-scenario storm surge inundation database constructed based on historical typhoon information can be used for storm surge inundation forecasting; this is a method closely related to the present study. Before the typhoon’s impact, inputting the typhoon forecast information and forecasted tide levels for estuarine areas into the database can yield inundation forecasting results within minutes, facilitating the dissemination of storm surge warnings and decision-making information; this method is quicker than numerical forecasting models. However, this method is only applicable to small-scale estuarine areas, and the storm surge inundation databases for each estuarine area may differ, requiring continuous research support.
In this paper, using the Jitimen area as a study area, we established and validated a numerical model for storm surge inundation with refined grids in small regions, conducted a multifactor, multi-scenario numerical simulation of storm surge inundation based on the statistical characteristics of historical data affecting the meteorology and occurrence of storm surges in the study area, analyzed the differences in the numerical results, and summarized the significant factors influencing storm surge inundation. The main conclusions of this paper are as follows:
- (1)
High-resolution basic geographic information data for the study area were obtained, and a numerical model for storm surge inundation based on unstructured triangular grids was established, with a refined grid resolution of 30–40 m for the land region. The 2017 super typhoon Hato and the 2018 super typhoon Mangkhut were selected as validation case examples. The errors between simulated and measured maximum storm surges were 3.2% and 5.0%, respectively, and the errors between simulated and measured highest tide levels were 1.2% and 2.4%, respectively, at the Sanzao tide gauge station, indicating excellent model performance.
- (2)
Based on the typhoon data from 1949 to 2022 and historical data from the Sanzao tide gauge station in the study area, the main characteristics of the typhoons affecting the study area were calculated. It was found that the typhoon intensity must reach at least the typhoon level for storm surge disasters to occur, the main typhoon landfall direction was found to occur in the northwest quadrant (W, WNW, NW, NNW, and N), the average typhoon speed was approximately 20 km/h, and the main impact occurred during the summer months (July to September) each year.
- (3)
Based on the statistical data of historical typhoons, using a typhoon speed of 20 km/h, and superimposing the average value of astronomical high tides in July to September, three experimental schemes with different typhoon landfall locations, intensities, and directions that influence storm surge inundation were set up to assess the differences in the numerical results. The results showed that when the typhoon landfall direction remained unchanged and the highest tide levels at the Sanzao tide gauge station were similar, the differences between the numerical results for the typhoon landfall location and typhoon intensity schemes were less than 5%, and the inundation characteristics were similar. However, when the typhoon location and intensity were unchanged and the highest tide levels at the Sanzao tide gauge station were similar, the numerical results for the typhoon landfall direction scheme differed significantly. The inundation areas under the north and northwest direction scenarios were 16.8% and 29.8% larger than those under the west direction scenario, respectively, and the average inundation depths under the north and northwest direction scenarios were 5.8% and 14.7% greater than those under the west direction scenario, respectively, indicating large differences in inundation characteristics.
The reasons for the differences in the storm surge inundation among the three scenarios were analyzed using the duration of the high tide levels (exceeding 3 m) at the Sanzao tide gauge station. Under the conditions with relatively small differences in the maximum tide levels and the high tide duration at the Sanzao station, the differences in the storm surge inundation area and average inundation depth among the three scenarios of typhoon landfall location (A1, A2, and A3) and the three scenarios of typhoon landfall intensity (B1, B2, and B3) are minimal. However, a significant difference in the duration of the high tide levels occurs among the three scenarios of the typhoon landfall direction (C1, C2, and C3), and differences in the storm surge inundation characteristics among them are also evident. The numerical results revealed that when the highest tide levels are similar, compared with the intensity and location of the typhoon, the typhoon landfall direction and the topographic characteristics of the estuary have a greater impact on storm surge inundation.