Influence of Grid Resolution and Assimilation Window Size on Simulating Storm Surge Levels

: Grid resolution and assimilation window size play significant roles in storm surge models. In the Bohai Sea, Yellow Sea, and East China Sea, the influence of grid resolution and assimilation window size on simulating storm surge levels was investigated during Typhoon 7203. In order to employ a more realistic wind stress drag coefficient that varies with time and space, we corrected the storm surge model using the spatial distribution of the wind stress drag coefficient, which was inverted using the data assimilation method based on the linear expression C d = (a + b × U 10 ) × 10 − 3 . Initially, two grid resolutions of 5 ′ × 5 ′ and 10 ′ × 10 ′ were applied to the numerical storm surge model and adjoint assimilation model. It was found that the influence of different grid resolutions on the numerical model is almost negligible. But in the adjoint assimilation model, the root mean square (RMS) errors between the simulated and observed storm surge levels under 5 ′ × 5 ′ and 10 ′ × 10 ′ grid resolutions were 11.6 cm and 15.6 cm, and the average PCC and WSS values for 10 tidal stations changed from 89% and 92% in E3 to 93% and 96% in E4, respectively. The results indicate that the finer grid resolution can yield a closer consistency between the simulation and observations. Subsequently, the effects of assimilation window sizes of 6 h, 3 h, 2 h, and 1 h on simulated storm surge levels were evaluated in an adjoint assimilation model with a 5 ′ × 5 ′ grid resolution. The results show that the average RMS errors were 11.6 cm, 10.6 cm, 9.6 cm, and 9.3 cm under four assimilation window sizes. In particular, the RMS errors for the assimilation window sizes of 1 h and 6 h at RuShan station were 3.9 cm and 10.2 cm, a reduction of 61.76%. The PCC and WSS values from RuShan station in E4 and E7 separately showed significant increases, from 85% to 98% and from 92% to 99%. These results demonstrate that when the assimilation window size is smaller, the simulated storm surge level is closer to the observation. Further, the results show that the simulated storm surge levels are closer to the observation when using the wind stress drag coefficient with a finer grid resolution and smaller temporal resolution.


Introduction
A storm surge is an unusual change in sea level induced by tropical cyclones (typhoons, hurricanes) or temperate cyclones (cold tides).For a long time, the damage caused by storm surges has been a major threat to human life, financial status, and infrastructure in coastal areas [1][2][3][4][5][6][7][8].Southeastern coastal areas in China are often affected by tropical cyclones from the Northwest Pacific Ocean.These disastrous typhoons could bring about serious economic losses and threaten people's lives in low-lying coastal areas [9][10][11][12][13][14][15].The disasters and losses brought by storm surges are difficult to estimate.In 1953, a huge storm surge caused tremendous disaster in Britain, Germany, the Netherlands, and other European countries, and resulted in thousands of deaths in the Netherlands [1].In 2005, Hurricane Katrina produced a storm surge level of 4-7 m; at least 1833 individuals were injured or killed and the direct economic loss amounted to USD 108 billion [5].Therefore, some factors of the storm surge model need to be thoroughly investigated in order to effectively prevent and assess storm surge disasters.
Grid resolution and assimilation window size are two important factors that can affect a storm surge model.In recent years, many scholars have conducted a lot of research on these factors.Jing et al. [16] utilized the Weather Research and Forecasting (WRF) model to explore the large eddy simulation (LES) of the high-wind area near the maximum wind radius of "Tiantu" (2021) using various grid resolutions.They found that altering the grid resolution could significantly affect the local turbulence structure of LES within a limited area, resulting in notable variations in eddy structure and intensity under different grid resolution conditions.Kerr et al. [17] used the Simulating Waves Nearshore (SWAN) model and an Advanced Circulation model (ADCIRC) to simulate the tidal harmonic components and hurricane wave results for Hurricane "Ike" (2008).Their study found that lower resolutions compromised simulation accuracy in coastal areas due to incorrect transmission or lateral attenuation.Moon et al. [18] highlighted that higher-resolution models for complex terrains and coastlines tend to produce higher average surges with a better simulation performance.Dukhovskoy et al. established a storm surge model with high resolution and used it in the Appalachian Gulf in the northeastern Gulf of Mexico during Hurricane Dennis (2005) [19].Through accurate-resolution simulations of coastal areas and waterways with intricate geometries, they unveiled the unexpected high storm surge processes in this region.Mentaschi et al. [20] developed a wave and storm surge prediction model with high resolution, which exhibited a good predictive ability for both sea level and effective wave height compared with satellite altimeters, tide gauges, and buoys as well as notable improvements compared to previous studies with lower resolutions regarding the reproduction of nearshore dynamics.Garzon et al. [21] found that the North American Mesoscale Forecast System (NAM) and European Center for Medium-Range Weather Forecasts (ECMWF) systems have the highest vertical and horizontal resolution, and indeed displayed the best root mean square deviation (RMSD) and correlation coefficient (CC), while the simulations of water levels based on the weather forecast systems with a higher horizontal resolution obtained better results.Makris et al. [22] used phase-resolving models with a fine resolution for the prediction of average ocean currents at sea level and depth in coastal areas, which are affected by atmospheric forcing and astronomical tides.The results address the significant needs of port authorities, ship pilots, and navigators in battling the problems of vessels impacting the harbor bed during mooring, towage, and berth operations using high-resolution and short-term sea-state forecasting.Fernández-Montblanc et al. [1] designed an unstructured hydrodynamic storm surge and tidal model for Europe.The tidal surge model accounts for the atmospheric pressure, wind, and astronomical tide.It was found that increasing the resolution of atmospheric forcing also improves the predictive ability, most extremely in shallow areas where wind is the main driver of surge production.Mohanty et al. [23] made an attempt to improve the storm surge prediction with a longer lead time using high-resolution mesoscale model outputs.Their findings suggest that the early warnings of tropical cyclones (TCs) obtained by the India Meteorological Department (IMD) should include the surge predictions from these highly reliable mesoscale model products with a 96-72 h lead time in order to mitigate the catastrophic loss associated with storm surges.Based on the research conducted by relevant scholars on grid resolution, it is evident that the effect of grid resolution on simulations of storm surge cannot be disregarded and is an important factor in storm surge simulation.
When analyzing the diurnal variation and hot wave phenomenon of Mars, Zhao et al. introduced different assimilation window lengths into a 4D local ensemble transform Kalman filter (4D-LETKF) to eliminate the artificially caused resonance phenomenon [24].The study found that the short assimilation window length can effectively eliminate false resonance.Wang et al. [25] used the ensemble Kalman filter (EnKF) method to compare the influence of assimilating T-TREC-retrieved winds (VTREC) versus radial velocity (Vr) on analyzing and forecasting Typhoon "Qaxi" (2008).The study found that the different assimilation windows had different effects in terms of data assimilation on the inversion wind field and radial wind speed.Zheng et al. [26] used a data assimilation scheme based on the 4DVar method to improve the prediction ability of an existing storm surge model in the North Sea of Germany.By diminishing the assimilation window's size, they found that the prediction accuracy was enhanced.Based on the ensemble Kalman filter, Kim et al. used a carbon tracker inverse simulation system to study the influence of assimilation window size on the estimation of surface carbon dioxide flux in Asia [27].The study found that when the assimilation window is shorter, the uncertainty of the optimized surface carbon dioxide flux is greater.DiNapoli et al. [28] applied a preoperative 4-day storm surge ensemble prediction system, called Model for Storm Surge Simulations (M3S), to the southwestern Atlantic continental shelf (SWACS) region.The system assimilated tidal level and elevation data using the four-dimensional ensemble square root filter (4DEnSRF) method.Their study showed that the first 2 days of the 4-day prediction depended on initial conditions, while the last 2 days were influenced by external forcing.Optimal initial conditions were obtained with a 12 h assimilation window size.Khan et al. [29] used numerically efficient hydrodynamics-waves coupled modeling and presented its practical real-time computational set-up; the results show that along this landfalling coastal section, the standard error in the maximum surge level amounts to 2.06, 1.73, and 0.66 m for the T-60 h, T-36 h, and T-12 h forecasts, respectively.Madsen et al. [30] demonstrated the positive impact of coastal altimetry observations when used in a statistical blending method together with tide gauge observations.A positive impact was demonstrated when the blended product was assimilated into a hydrodynamic model of the North Sea and Baltic Sea, showing that the simplified, computational cost-effective assimilation method improves the modeled sea level field.The aforementioned studies suggest that further investigation is warranted to explore the impact of assimilation window size on simulating storm surge levels.
In order to assess the influences of grid resolution and assimilation window size on simulating the storm surge levels, seven experiments were designed in this paper.In order to produce a more realistic wind stress drag coefficient, we have corrected the adjoint assimilation model using a finer grid resolution and smaller assimilation window size.Firstly, we designed four experiments (E1-E4) to examine the influence of different grid resolutions on both the numerical storm surge model and the adjoint assimilation model.Subsequently, we also discussed the impact of assimilation window size on the adjoint assimilation model of storm surges in E4-E7.
This paper consists of four sections.The second section introduces the numerical storm surge model and adjoint assimilation model, as well as the experimental design.A detailed investigation into how the grid resolution and assimilation window size affect the simulated levels of storm surge is provided in the third section.Conclusions are presented in the last part.

Typhoon and Stations
Typhoon 7203 formed in the Micronesia region on 5 July 1972 and evolved into a tropical storm on 7 July, peaking early on 11 July.Subsequently, Typhoon 7203 moved northwestward and passed into both the East China Sea and Yellow Sea before making landfall in Tanggu, Tianjin, on 27 July.After traversing Yanshan Mountain, it dissipated at the periphery of the Inner Mongolia Plateau.As a result of its impact, Beijing, Tianjin, Hebei Province, Shandong Province, and the Liaodong Peninsula experienced heavy to torrential rainfall accompanied by winds ranging in force from 9 to 12 as well as strong storm surges in the Bohai Sea.Notably, Qinhuangdao witnessed its highest recorded tidal level at approximately 2.48 m during this event.From 20:00 on 27 July 1972, heavy rainfall occurred in the jujube forest located in the suburban areas of Beijing.Yanshan's peculiar topographic features enhance wind convergence, which occurs when horizontal winds cause a net inflow of air into a region (more air entering a vertical column than leaving it).The maximum precipitation recorded within a span of 24 h reached an unprecedented level of 479.2 mm, which has held the record for precipitation in Beijing for a long time.Consequently, mudslides ensued in the mountainous regions of Huairou and Yanqing, resulting in extensive structural damage and claiming approximately 50 lives.In this study, Figure 1 shows the typhoon's path and the locations of 10 tidal stations, utilized as inputs into the storm surge model.These data are from the typhoon database website ("http://www.wztf121.com/analyse/wzs",accessed on 2 July 2024).torrential rainfall accompanied by winds ranging in force from 9 to 12 as well as strong storm surges in the Bohai Sea.Notably, Qinhuangdao witnessed its highest recorded tidal level at approximately 2.48 m during this event.From 20:00 on 27 July 1972, heavy rainfall occurred in the jujube forest located in the suburban areas of Beijing.Yanshan's peculiar topographic features enhance wind convergence, which occurs when horizontal winds cause a net inflow of air into a region (more air entering a vertical column than leaving it).
The maximum precipitation recorded within a span of 24 h reached an unprecedented level of 479.2 mm, which has held the record for precipitation in Beijing for a long time.
Consequently, mudslides ensued in the mountainous regions of Huairou and Yanqing, resulting in extensive structural damage and claiming approximately 50 lives.In this study, Figure 1 shows the typhoon's path and the locations of 10 tidal stations, utilized as inputs into the storm surge model.These data are from the typhoon database website ("http://www.wztf121.com/analyse/wzs",accessed on 2 July 2024).

Numerical Storm Surge Model and Adjoint Assimilation Model
In the present study, a depth-averaged shallow-water equation (SWE) model was utilized to calculate the storm surge level.In the Cartesian coordinate system, the primary equations of the storm surge model are composed of the vertically integrated continuum equation and the momentum equation, as follows [31][32][33][34][35][36]:

Numerical Storm Surge Model and Adjoint Assimilation Model
In the present study, a depth-averaged shallow-water equation (SWE) model was utilized to calculate the storm surge level.In the Cartesian coordinate system, the primary equations of the storm surge model are composed of the vertically integrated continuum equation and the momentum equation, as follows [31][32][33][34][35][36]: where Generally speaking, there are some errors in the simulated results caused by approximate assumptions, parameters, etc.In order to decrease the error between the simulated storm surge levels and the observed sea level, we introduced the Lagrange multiplier method and constructed the following form of objective function: where ζ is the simulation, ζ obs is the observation, and K ζ is the constant.The corresponding Lagrange function is constructed in the following form: By means of the derivation approach provided by He et al. [37], the corresponding adjoint equations are achieved: where u α , ζ α and v α are the corresponding adjoint variables of ζ, u and v.
The pressure and wind fields used in this model are derived from field observations of meteorological parameters interpolated over the storm surge model grid, based on the work of Jelesnianski et al. [38].The following is the pressure field: where P 0 is the typhoon's central pressure, P α is the sea's surface pressure at r, and P ∞ is ambient pressure.The wind field is given by the following: where R = radius of maximum wind speed W R ; V ox and V oy = migration velocities of the typhoon's center; → i and → j = unit vector in x and y axes;

Experimental Design
In order to investigate the effects of grid resolution and assimilation window size on the storm surge level, seven experiments are designed in this section, as presented in Table 1.The wind stress drag coefficient from Wu [39] is used in experiments E1 and E2, with grid resolutions of 10 ′ × 10 ′ and 5 ′ × 5 ′ , respectively.In experiments E3 (10 ′ × 10 ′ ) and E4 (5 ′ × 5 ′ ), the adjoint assimilation method using a linear formula of wind stress drag coefficient is employed to assess the impacts of different grid resolutions on the simulated storm surge levels.The assimilation window size for these two experiments is set as 6 h.In experiments E5, E6, and E7, with a grid resolution of 5 ′ × 5 ′ , the assimilation window sizes are adjusted to durations of 3 h, 2 h, and 1 h, respectively.Subsequently, in conjunction with experiment E4, we analyze the influence of assimilation window size on simulating storm surge levels.

Experiments
Grid Resolution Assimilation Window

Model Setup
In this study, the areas applied in the numerical model and adjoint assimilation model of the storm surge are 110-130 • E and 24-44 • N, covering the Bohai Sea, Yellow Sea, and East China Sea.The bottom friction coefficient in those areas is set at a constant of 0.0016.The initial flow velocity and sea surface level are set to 0 in the model.The open boundary conditions are the first island chain and the Taiwan Strait.There is no water flowing in or out of the coastal area on the closed boundary.The bathymetric data are derived from ETOPO1.The barometric pressure and wind stress are taken as the driving forces, and the time interval is 60 s.The model adopts an Arakawa C grid.The observation data used for the experiment are taken from the storm surge levels at 10 tidal stations.

Influence of the Grid Resolution on Simulating Storm Surge Levels
The influence of grid resolution on the storm surge level is evaluated using the numerical storm surge model and adjoint assimilation model in this section.Tables 2 and 3 display the RMS errors and absolute mean difference (AMD) errors between simulations and observations at 10 tidal stations in E1-E4 at the end of each stage.However, the results presented in Table 4 are the RMSEs between simulations and observations at each tidal station at the end of the model's operation.In other words, one in Table 3 is temporal and the other in Table 4 is spatial.It can be seen that the RMS and AMD errors are nearly identical in E1 and E2 using the numerical storm surge model.These results indicate that the influence of grid resolution on the numerical storm surge model is almost negligible.The RMS errors in E3 and E4 from the adjoint assimilation model are 15.6 cm and 11.6 cm, and the AMD errors are 12.1 cm and 8.9 cm, respectively.The results reveal that the behavior of the adjoint assimilation model in E4 is better than that in E3, indicating that the finer grid resolution can lead to a closer consistency between simulated and observed storm surge levels.In addition, Table 4 presents the RMS errors, the Pearson correlation coefficient (PCC) [40], and the Willmott skill scores (WSS) [41] between simulation and observation at 10 tidal stations during Typhoon 7203.Taking the Yantai station as an example, the RMS errors in E3 and E4 are 24.5 cm and 18.6 cm, respectively, indicating a reduction of 24.08%.The average RMS errors have decreased by 27.49%, from 17.1 cm in E3 to 12.4 cm in E4.The average PCC values for 10 tidal stations in E1-E4 are 81%, 80%, 89%, and 93%, respectively.The average WSS values in E1-E4 are 82%, 82%, 92%, and 96%.It can be seen that the PCC and WSS values in E1 and E2 are almost identical, while these two indices in E3 and E4 significantly increase; that is, the PCC and WSS in the adjoint assimilation model exceed those in the numerical storm surge model.These results in Tables 3-5 demonstrate that regardless of the values of each RMS error, the PCC, the WSS, or their average values, the finer the grid resolution is, the closer the simulation is to the observation.The storm surge levels and the differences between simulated and observed levels at Yingkou, Huludao, and Yantai tidal stations in E1-E4 are depicted in Figures 2-4.In E3-E4, the storm surge levels are simulated via the spatial distribution of the wind stress drag coefficient, which is inverted via the data assimilation method based on the linear expression C d = (a + b × U 10 ) × 10 −3 .According to the time series graphs of simulations and observations, it is apparent that the results from E1 and E2 are largely different from the observed levels, and the simulations of storm surge levels in E3 and E4 more closely approach the observed values.This result shows that the adjoint assimilation model is superior to the numerical storm surge model.At the same time, comparing the results of E3 and E4, it can be found that the result in E4 is slightly superior to that in E3, showing that simulating storm surge levels with finer grid resolutions yields results that are closer to the observations.coefficient, which is inverted via the data assimilation method based on the linear expression Cd = (a + b × U10) × 10 −3 .According to the time series graphs of simulations and observations, it is apparent that the results from E1 and E2 are largely different from the observed levels, and the simulations of storm surge levels in E3 and E4 more closely approach the observed values.This result shows that the adjoint assimilation model is superior to the numerical storm surge model.At the same time, comparing the results of E3 and E4, it can be found that the result in E4 is slightly superior to that in E3, showing that simulating storm surge levels with finer grid resolutions yields results that are closer to the observations.coefficient, which is inverted via the data assimilation method based on the linear expression Cd = (a + b × U10) × 10 −3 .According to the time series graphs of simulations and observations, it is apparent that the results from E1 and E2 are largely different from the observed levels, and the simulations of storm surge levels in E3 and E4 more closely approach the observed values.This result shows that the adjoint assimilation model is superior to the numerical storm surge model.At the same time, comparing the results of E3 and E4, it can be found that the result in E4 is slightly superior to that in E3, showing that simulating storm surge levels with finer grid resolutions yields results that are closer to the observations.The peak values of the storm surge level between the simulations and observations in E1-E4 are exhibited in Figure 5.At YingKou station, the simulated peaks of storm surge levels in E1, E2, and E4 were, respectively, underestimated by 46.35 cm, 47.69 cm, and 0.74 cm, and the peak in E3 was overestimated by 34.3 cm compared to the observed peak surge.Similarly, at HuLuDao station, the peaks of surges in E1-E4 were underestimated by 60.17 cm, 56.17 cm, 28.4 cm, and 16.4 cm compared to the observation, respectively.The simulated peaks of surge in E3 and E4 approached more closely to the observations than those in E1 and E2, indicating that the adjoint assimilation model performs better than the numerical model.Meanwhile, the peak simulated in E4 is closer to the observation than that in E3, indicating that a finer grid resolution can result in better simulation results.The peak values of the storm surge level between the simulations and observations in E1-E4 are exhibited in Figure 5.At YingKou station, the simulated peaks of storm surge levels in E1, E2, and E4 were, respectively, underestimated by 46.35 cm, 47.69 cm, and 0.74 cm, and the peak in E3 was overestimated by 34.3 cm compared to the observed peak surge.Similarly, at HuLuDao station, the peaks of surges in E1-E4 were underestimated by 60.17 cm, 56.17 cm, 28.4 cm, and 16.4 cm compared to the observation, respectively.The simulated peaks of surge in E3 and E4 approached more closely to the observations than those in E1 and E2, indicating that the adjoint assimilation model performs better than the numerical model.Meanwhile, the peak simulated in E4 is closer to the observation than that in E3, indicating that a finer grid resolution can result in better simulation results.The peak values of the storm surge level between the simulations and observations in E1-E4 are exhibited in Figure 5.At YingKou station, the simulated peaks of storm surge levels in E1, E2, and E4 were, respectively, underestimated by 46.35 cm, 47.69 cm, and 0.74 cm, and the peak in E3 was overestimated by 34.3 cm compared to the observed peak surge.Similarly, at HuLuDao station, the peaks of surges in E1-E4 were underestimated by 60.17 cm, 56.17 cm, 28.4 cm, and 16.4 cm compared to the observation, respectively.The simulated peaks of surge in E3 and E4 approached more closely to the observations than those in E1 and E2, indicating that the adjoint assimilation model performs better than the numerical model.Meanwhile, the peak simulated in E4 is closer to the observation than that in E3, indicating that a finer grid resolution can result in better simulation results.

Influence of Assimilation Window Size on Simulated Storm Surge Levels
The influence of the assimilation window size on the simulated levels of storm surge is addressed by setting 6 h, 3 h, 2 h, and 1 h sizes under a grid resolution of 5 ′ × 5 ′ .The RMS errors and AMD errors between the simulations and observations in E4-E7 are listed in Table 5.The RMS errors in E4-E7 are, on average, 11.6 cm, 10.6 cm, 9.6 cm, and 9.3 cm, respectively, and the mean RMS errors in E5, E6, and E7 are 8.62%, 17.24%, and 19.83% lower when compared to that in E4.Additionally, the RMS errors, the PCC, the WSS between the simulations and observations at 10 tidal stations, and their averages in E4-E7 at 10 tidal stations are listed in Table 6.Taking RuShan station as an example, the RMS errors in E4-E7 are 10.2 cm, 7.7 cm, 6.1 cm, and 3.9 cm, respectively.The RMS errors in E5, E6, and E7 are 24.51%,40.20%, and 61.76% lower when compared with that in E4.The average PCC values for 10 tidal stations in E4-E7 are 93%, 95%, 97%, and 97%, respectively.The average WSS values in E4-E7 are 96%, 97%, 98%, and 98%.These two evaluation indicators increase as the assimilation window size gradually decreases, and the accuracy of the storm surge models in E6 and E7 shows little difference.The above results indicate that the simulations of storm surge levels are closer to the observed levels when the assimilation window size is smaller.Via comprehensive comparison, we can see that the assimilation window size of 1 h or 2 h is applicable in the adjoint assimilation model of storm surge when simulating storm surge levels.The storm surge levels and the differences between simulations and observations at YingKou, LongKou, and RuShan tidal stations in E4-E7 are depicted in Figures 6-8.
In particular, at RuShan station, the PCC and WSS values in E4 and E7 show significant increases from 85% to 98% and from 92% to 99%, respectively.Comparatively, the simulated values in E4 and E5 are far from the observed values, and the simulated values in E6 and E7 are closer to the observations.These findings all suggest that the storm surge levels simulated with smaller assimilation window sizes are closer to the observed levels.

Conclusions
In relation to Typhoon 7203, we have explored the effects of the grid resolution and assimilation window size on simulations of storm surge levels in the Bohai Sea, Yellow Sea, and East China Sea in this paper.In the adjoint assimilation model, we used the inverted spatial distribution of the wind stress drag coefficient to calculate the storm surge level via the data assimilation method based on the linear expression Cd = (a + b × U10) × 10 −3 .
In order to investigate the influences of different grid resolutions on simulations of storm surge levels, we conducted four experiments.In E1 and E2, two grid resolutions of 10′ × 10′ and 5′ × 5′ were set up in the numerical storm surge model.Similarly, two experiments were carried out on the adjoint assimilation method using grid resolutions of 10′ × 10′ and 5′ × 5′ in E3 and E4.The results of the evaluation indicators, comprising the values of RMS error, AMD error, PCC, WSS, and their average, reveal that in E1 and E2, the impact of the grid resolution was minimal and practically negligible in the numerical model of storm surges.However, in E3 and E4, when using the adjoint assimilation model, varying grid resolutions were found to significantly affect the simulation accuracy.It was found that finer grids can yield more precise simulation levels.
The influence of the assimilation window's size on simulations of storm surge levels with a 5′ × 5′ grid resolution was investigated in the present study.Using the adjoint as-

Conclusions
In relation to Typhoon 7203, we have explored the effects of the grid resolution and assimilation window size on simulations of storm surge levels in the Bohai Sea, Yellow Sea, and East China Sea in this paper.In the adjoint assimilation model, we used the inverted spatial distribution of the wind stress drag coefficient to calculate the storm surge level via the data assimilation method based on the linear expression C d = (a + b × U 10 ) × 10 −3 .
In order to investigate the influences of different grid resolutions on simulations of storm surge levels, we conducted four experiments.In E1 and E2, two grid resolutions of 10 ′ × 10 ′ and 5 ′ × 5 ′ were set up in the numerical storm surge model.Similarly, two experiments were carried out on the adjoint assimilation method using grid resolutions of 10 ′ × 10 ′ and 5 ′ × 5 ′ in E3 and E4.The results of the evaluation indicators, comprising the values of RMS error, AMD error, PCC, WSS, and their average, reveal that in E1 and E2, the impact of the grid resolution was minimal and practically negligible in the numerical model of storm surges.However, in E3 and E4, when using the adjoint assimilation model, varying grid resolutions were found to significantly affect the simulation accuracy.It was found that finer grids can yield more precise simulation levels.
The influence of the assimilation window's size on simulations of storm surge levels with a 5 ′ × 5 ′ grid resolution was investigated in the present study.Using the adjoint assimilation method in E4-E7, assimilation window sizes with intervals of 6 h, 3 h, 2 h, and 1 h were assessed, respectively.The results demonstrate that the performance in E6 and E7 was superior to that in E4 and E5, as suggested by the comparison of RMS error, AMD error, PCC, and WSS.In particular, the average PCC and WSS values reached 97% and 98% when the assimilation window sizes were 2 h and 1 h in E6-E7.Therefore, the effect of assimilation window size on storm surge levels is also very important.Smaller assimilation window size showed enhanced accuracy and higher PCC and WSS in simulating storm surge levels.
In addition to the grid resolution and assimilation window size, there are many other factors that can also affect the storm surge model, such as the bottom friction coefficient, the time interval of external forcing, the type of external forcing, and the wind stress drag coefficient.In future research, we will consider the effects of these factors on model performance.

Figure 1 .
Figure 1.Typhoon path and locations of tidal stations.Star symbols represent the tidal stations' locations.The solid line indicates the path of Typhoon 7203.Circles represent the time.

Figure 1 .
Figure 1.Typhoon path and locations of tidal stations.Star symbols represent the tidal stations' locations.The solid line indicates the path of Typhoon 7203.Circles represent the time.
wind stress drag coefficient; ρ α = air density; A = the horizontal eddy viscosity coefficient; g = gravitational acceleration; t = time; u and v = velocity in x and y directions, respectively; k = bottom friction factor; h = unperturbed water depth; W x and W y = the surface wind field; ρ ω = sea water density; P α = sea surface pressure;x and y = longitude and latitude; ζ = surface level; f = Coriolis parameter.

Figure 2 .
Figure 2. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at YingKou station during Typhoon 7203.

Figure 3 .
Figure 3. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at HuLuDao station during Typhoon 7203.

Figure 2 .
Figure 2. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at YingKou station during Typhoon 7203.

Figure 2 .
Figure 2. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at YingKou station during Typhoon 7203.

Figure 3 .
Figure 3. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at HuLuDao station during Typhoon 7203.Figure 3. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at HuLuDao station during Typhoon 7203.

Figure 3 .
Figure 3. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at HuLuDao station during Typhoon 7203.Figure 3. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at HuLuDao station during Typhoon 7203.

Figure 4 .
Figure 4. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at YanTai station during Typhoon 7203.

Figure 5 .
Figure 5. Peaks of simulated and observed storm surge levels in E1-E4 at HuLuDao and YingKou stations.

Figure 4 .
Figure 4. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at YanTai station during Typhoon 7203.

Figure 4 .
Figure 4. Simulated storm surge levels in E1-E4 and the observed level (top), and the differences between the simulation and observation (bottom) at YanTai station during Typhoon 7203.

Figure 5 .
Figure 5. Peaks of simulated and observed storm surge levels in E1-E4 at HuLuDao and YingKou stations.

Figure 5 .
Figure 5. Peaks of simulated and observed storm surge levels in E1-E4 at HuLuDao and YingKou stations.

Figure 6 .
Figure 6.Simulated storm surge levels in E4-E7 and observed level (top), and the differences between the simulation and observation (bottom) at YingKou station during Typhoon 7203.

Figure 6 .
Figure 6.Simulated storm surge levels in E4-E7 and observed level (top), and the differences between the simulation and observation (bottom) at YingKou station during Typhoon 7203.

Figure 6 .
Figure 6.Simulated storm surge levels in E4-E7 and observed level (top), and the differences between the simulation and observation (bottom) at YingKou station during Typhoon 7203.

Figure 7 .
Figure 7. Simulated storm surge levels in E4-E7 and observed level (top), and the differences between the simulation and observation (bottom) at LongKou station during Typhoon 7203.Figure 7. Simulated storm surge levels in E4-E7 and observed level (top), and the differences between the simulation and observation (bottom) at LongKou station during Typhoon 7203.

Figure 7 . 17 Figure 8 .
Figure 7. Simulated storm surge levels in E4-E7 and observed level (top), and the differences between the simulation and observation (bottom) at LongKou station during Typhoon 7203.Figure 7. Simulated storm surge levels in E4-E7 and observed level (top), and the differences between the simulation and observation (bottom) at LongKou station during Typhoon 7203.J. Mar.Sci.Eng.2024, 12, x FOR PEER REVIEW 15 of 17

Figure 8 .
Figure 8. Simulated storm surge levels in E4-E7 and observed level (top), and the differences between the simulation and observation (bottom) at RuShan station during Typhoon 7203.

Table 2 .
Root mean square errors between simulated and observed levels in E1-E4 (unit: cm).
* The typhoon process selected in this study is 72 h.There are 12 stages during which each stage is set as 6 h in E1-E2.There are also 12 stages during which the assimilation window size is 6 h in E3-E4.These 12 points of comparison in Table2represent the average RMSEs between simulations and observations for 10 tidal stations in E1-E4 at the end of each stage.

Table 3 .
Absolute mean difference errors between simulated and observed levels in E1-E4 (unit: cm).The typhoon process selected in this study is 72 h.There are 12 stages during which each stage is set as 6 h in E1-E2.There are also 12 stages during which the assimilation window size is 6 h in E3-E4.These 12 points of comparison in Table2represent the average AMDEs between simulations and observations for 10 tidal stations in E1-E4 at the end of each stage. *

Table 4 .
Root mean square errors (unit: cm), Pearson correlation coefficient, and Willmott skill scores between the simulated and observed storm surge levels in E1-E4 at 10 tidal stations.

Table 5 .
Root mean square errors and absolute mean difference errors between simulations and observations in E4-E7 (unit: cm).

Table 5 .
Cont.The typhoon process selected in this study lasts 72 h.In stages 12, 24, 36, and 72, the assimilation window sizes are 6, 3, 2, and 1 h in E4-E7.The 12, 24, 36, and 72 points of comparison in Table5represent the RMSEs and AMDEs between simulations and observations in E4-E7 at the end of each stage. *

Table 6 .
Root mean square errors (unit: cm), Pearson correlation coefficient, and Willmott skill scores between the simulated storm surge levels and observations in E4-E7 at 10 tidal stations.