Evolution Wave Condition Using WAVEWATCH III for Island Sheltered Area in the South China Sea

: Wave conditions around islands in the South China Sea (SCS) are of signiﬁcant interest due to their importance for marine operations and coastal engineering. Understanding and accurately predicting wave characteristics in this region are crucial. In this study, the third-generation wave model WAVEWATCH III is employed to examine wave conditions around islands in the SCS. Accord-ing to the water depth and signiﬁcant wave height, the sea state around the island was classiﬁed into two categories: typhoon sea state and moderate sea state. Several popular wind input–dissipation source terms (ST2, ST4 and ST6) are used to assess the typhoon sea state and the moderate sea state separately. The results are validated by ﬁeld wave data. ST4 and ST6 show good performance in signiﬁcant wave height for moderate sea states, while ST2 is good at the mean wave period. For the typhoon sea state, ST2 gives the best results in signiﬁcant wave height with larger correlation coefﬁcients and a smaller RMSE. The above results provide valuable insights into the effects of different source terms on the accuracy of wave simulations for different sea states. The spatial distribution of the significant wave heights is also demonstrated with ST2, which may be useful for assessing the wave conditions of marine structures from the large scale of the SCS to the island scale of the Yongle Atoll.


Introduction
Wave conditions are very important for the safety of marine structures. The accurate assessment of wave conditions can optimize the design parameters and enhance the performance of marine structures [1]. Marine structures, such as very large floating structures, fish farming cages and plate forms, are very important for island development, especially in remote areas away from the mainland, providing spaces for living and production. However, the wave conditions around islands and reefs are highly dependent on the complex submarine topography. Due to the specific bathymetry around islands in the South China Sea, the physical processes involved in wave evolution are different between deep-sea areas and island areas. Therefore, the use of the WWIII model has been extended from deep-sea and continental shelf areas to shallow-water nearshore areas. The applicability of the WWIII model in near-island areas has been extensively verified [2][3][4][5]. Moreover, several studies have confirmed that the WWIII model provides generally good results at different wind forces in deep water by comparing it with other models such as SWAN, WAM, WAVAD and WWIII [6][7][8][9].
The WWIII model offers a diverse range of source terms, particularly for wind inputdissipation, and scholars have conducted numerous studies to determine their applicability and accuracy [10][11][12]. These include the WAM Cycle 3 (ST1) package [13,14], the Tolman and Chalikov (1996) package [15] (ST2), the WAM Cycle 4 (ST3) package [16,17], the 2 of 19 Ardhuin et al. (2010) package [18] (ST4), and the Zieger et al. (2015) [19] source term (ST6). Several scholars have examined the performance of different wind input-dissipation source terms in different regions of the globe under different wind forces. By comparing the performance of ST2, ST3, ST4 and ST6, Liu et al. [20] found that ST2 was the least accurate formula for simulating Hurricane Ivan, while ST4 results were in better agreement with the measurements. The same hurricane simulation results [21] show that the most appropriate formula is not the same for different areas (Gulf of Mexico). For the northern Tyrrhenian Sea and off the Mediterranean Spanish coast, Mentaschi et al. [22] showed that the source term introduced by ST4 is the best overall choice. However, Kalantzi et al. [23] found that neither ST1 nor ST2 can properly simulate waves in the presence of swells. Foli et al. [24] demonstrated that ST2STAB (a variant of ST2) and ST6 were suitable for simulating wave height and direction, respectively, in the entire West African region. In particular, the ST2 has gained wide acceptance in typhoons and rough sea conditions [25][26][27][28][29].
In general, waves passing around islands will experience refraction, diffraction, and deformation [30,31]. Moreover, the South China Sea (SCS) is known to have abundant marine resources, particularly in the vicinity of reef islands [32,33]. The measurement location in our research has a depth of 30 m. Due to the significant differences in wave conditions, we have categorized them into a moderate sea state and a typhoon sea state based on the characteristic wave length. In the following sections, we aim to present a representative wave forecasting scheme for the Yongle Atoll area under different oceanic conditions, namely, the typhoon sea state and moderate sea state. Section 2 describes the study area, data and wave model. The numerical wave conditions are validated using field wave data, and the best parameter choices and plan are selected for the moderate sea state and the typhoon sea state in Section 3. Conclusions are given in Section 4.

Spectral Balancing Equation
Third-generation wave models, such as WAVEWATCH III, are based on the wave spectrum balance equation and widely applied in both theoretical and practical applications [34]. In general, WWIII simulates the propagation of waves by solving the wave ac-

Spectral Balancing Equation
Third-generation wave models, such as WAVEWATCH III, are based on the wave spectrum balance equation and widely applied in both theoretical and practical applications [34]. In general, WWIII simulates the propagation of waves by solving the wave action balance equation: where F is the fundamental spectrum; N is the wave action density spectrum; S represents the net effect of sources and sinks for the spectrum F. k, θ, and σ represent the wavenumber, wave direction, and intrinsic frequency, respectively, and x and t are variables that represent space and time [35]. Then, the balance Equation (1) is extended to the following conservative form: .
where c g = (c g sin θ, c g cos θ) denotes the group velocity; d is the mean water depth; and U is the current velocity (depth-and time-averaged over the scales of individual waves). S is a coordinate in the direction θ, and m is a coordinate perpendicular to S. Equation (3) is used for Cartesian coordinates. For large-scale calculations, this equation is usually applied to spherical coordinates, defined by longitude λ and latitude φ [36].
where R is the radius of the earth, and U φ , and U λ are current components. θ = 0 corresponds to waves traveling from west to east. In WWIII wave models, energy is conserved in the absence of currents. However, when currents are present, the energy of the spectrum is no longer conserved because the currents shift the average momentum acting on the waves [37,38]. Therefore, the wave action density spectrum is the preferred choice in the WWIII model.
The net source term S in fact consists of four main components: an atmosphere-wave interaction term S in , a nonlinear wave-wave interaction term S nl , the dissipation term S ds , and the wave-bottom interaction term S bot . However, in this study, we need to consider the effect of shallow water on wave evolution around the Yongle Atoll. Specifically, depthinduced breaking (S db ) is crucial for numerical simulations in shallow water. To provide a more realistic initial wave growth, a linear input term S ln can also be added. In summary, the main source terms considered in the calculation process are presented in the following equation [26,39]: S = S ln + S in + S nl + S ds + S bot + S db (11)

Spatial Grids
To balance computational efficiency with accuracy, an unstructured triangular grid consisting of 42,445 nodes and 84,961 cells was selected for our calculations, as shown in Figure 2 [40]. Although our primary focus is on the Yongle Atoll area in the SCS, this grid covers the entire expansive computational domain of the SCS. To ensure computational efficiency, a coarser grid is utilized outside of the Yongle Atoll. The grid resolution becomes progressively finer closer to the Yongle Atoll, enhancing the accuracy of our calculations.

Wave Model Setup
The study utilized WAVEWATCH III version 5.16 to investigate the wind input-dissipation source terms for ocean wave simulations. The code for different source items is modular, providing multiple options for the same source term [41]. Wind input term in S and dissipation terms ds S are two different processes, but the balance of these two source terms determines the overall growth characteristics of wave energy and is considered interrelated. In this study, we focus on the effect of different input-dissipation source term packages, which are listed below.  [19].
The wave energy balance equation [39] is solved using the same source terms as described in Table 1. The spectrum is discretized using 30 frequencies ranging from 0.04118 to 0.6532 Hz with an increment factor of 1.1 and 36 directions with a 10-degree directional increment. To improve computational efficiency, WWIII utilizes four types of time steps: a maximum global time step and a maximum CFL time step for x-y of 20 s; a minimum source term time step of 20 s; and a maximum CFL time step for x-y of 10 s. Other parameter schemes remain the same as the default values of the WWIII model, and more details can be found in the user manual [35].

Wave Model Setup
The study utilized WAVEWATCH III version 5.16 to investigate the wind inputdissipation source terms for ocean wave simulations. The code for different source items is modular, providing multiple options for the same source term [41]. Wind input term S in and dissipation terms S ds are two different processes, but the balance of these two source terms determines the overall growth characteristics of wave energy and is considered interrelated. In this study, we focus on the effect of different input-dissipation source term packages, which are listed below.  [19].
The wave energy balance equation [39] is solved using the same source terms as described in Table 1. The spectrum is discretized using 30 frequencies ranging from 0.04118 to 0.6532 Hz with an increment factor of 1.1 and 36 directions with a 10-degree directional increment. To improve computational efficiency, WWIII utilizes four types of time steps: a maximum global time step and a maximum CFL time step for x-y of 20 s; a minimum source term time step of 20 s; and a maximum CFL time step for x-y of 10 s. Other parameter schemes remain the same as the default values of the WWIII model, and more details can be found in the user manual [35].

Wind Field Data
This study employs the Weather Research and Forecasting (WRF) model to perform reanalysis of wind fields. The wind field was forced with 1 h time step WRF from reanalysis dataset, with a spatial resolution of 0.1 • in both longitude and latitude, and the number of grids is 201 and 151, respectively. Various studies [39,[51][52][53] have shown that the accuracy of the wind field has a significant impact on the wave field. In order to verify the accuracy of the wind field, the correlation coefficient (Corrcoef), root mean square error (RMSE) and bias (Bias) are introduced in this paper to describe the relationship between the measured and simulated values. These statistics are defined by the following equations: where x i and y i represent the measured and simulated values, respectively, and N is the number of samples. The wind velocity simulation and observation comparison are presented in Figure 3. The results show a high correlation coefficient of 0.905, indicating good agreement. However, the simulated wind velocity is slightly higher than the measured wind velocity, with a bias of 1.233. Based on these results, the simulated wind field is suitable for use as input data for the calculation model.

Sea Condition
Moderate sea state is defined as SWH less than one meter in this paper, while severe wave condition, referred to as typhoon sea state, corresponds to SWH larger than one meter at Yongle Atoll. It should be noted that for the most common wave steepness of 1:30, one meter SWH corresponds to a wave length of 30 m, which is the wave depth of the measurement points as shown in Figure 4. Notably, the water depth of both points (Point 1 and Point 2) is 30 m.
The typhoon sea conditions in this article calculate the wave evolution from 3 September to 15 September 2020, during which the wind level in the target sea remains within level 5 (i.e., wind speed below 10.7 m/s), and the sea state is categorized as level 4.

Sea Condition
Moderate sea state is defined as SWH less than one meter in this paper, while severe wave condition, referred to as typhoon sea state, corresponds to SWH larger than one meter at Yongle Atoll. It should be noted that for the most common wave steepness of 1:30, one meter SWH corresponds to a wave length of 30 m, which is the wave depth of the measurement points as shown in Figure 4. Notably, the water depth of both points (Point 1 and Point 2) is 30 m.

Sea Condition
Moderate sea state is defined as SWH less than one meter in this paper, while severe wave condition, referred to as typhoon sea state, corresponds to SWH larger than one meter at Yongle Atoll. It should be noted that for the most common wave steepness of 1:30, one meter SWH corresponds to a wave length of 30 m, which is the wave depth of the measurement points as shown in Figure 4. Notably, the water depth of both points (Point 1 and Point 2) is 30 m.  The typhoon sea conditions in this article focus specifically on Typhoon No. 22 of 2020, known as "Vamco". "Vamco" originated in the sea east of the Philippines on 9 November 2020 and reached the eastern part of the SCS on 12 November. The path of the typhoon propagation from southeast to northwest is depicted in Figure 5.
To gain a clearer understanding of the wind velocity during the movement of the typhoon, Figure 6 displays the comparison between the measured and simulated wind velocity for "Vamco". It is evident that the simulated wind velocity is considerably higher than the measured wind velocity near the peak of the typhoon, with the simulated values being about 1.8 times greater than the measured values. The velocity vectors of the typhoon's passage are shown in Figure 7. The typhoon sea conditions in this article calculate the wave evolution from 3 September to 15 September 2020, during which the wind level in the target sea remains within level 5 (i.e., wind speed below 10.7 m/s), and the sea state is categorized as level 4.
The typhoon sea conditions in this article focus specifically on Typhoon No. 22 of 2020, known as "Vamco". "Vamco" originated in the sea east of the Philippines on 9 No-vember 2020 and reached the eastern part of the SCS on 12 November. The path of the typhoon propagation from southeast to northwest is depicted in Figure 5.
typhoon propagation from southeast to northwest is depicted in Figure 5.
To gain a clearer understanding of the wind velocity during the movement of the typhoon, Figure 6 displays the comparison between the measured and simulated wind velocity for "Vamco". It is evident that the simulated wind velocity is considerably higher than the measured wind velocity near the peak of the typhoon, with the simulated values being about 1.8 times greater than the measured values. The velocity vectors of the typhoon's passage are shown in Figure 7.  To gain a clearer understanding of the wind velocity during the movement of the ty-phoon, Figure 6 displays the comparison between the measured and simulated wind velocity for "Vamco". It is evident that the simulated wind velocity is considerably higher than the measured wind velocity near the peak of the typhoon, with the simulated values being about 1.8 times greater than the measured values. The velocity vectors of the typhoon's passage are shown in Figure 7.      During the moderate sea state, the SWHs at both point 1 and 2 did not exceed 1 m and showed fluctuations in response to changes in wind velocity. A comparison of the time history curves between the simulated SWHs and the measured values revealed the performance of the different input-dissipation source term schemes in Figure 8. All three input-dissipation source term schemes were able to adequately capture the trend of the SWHs in the moderate sea state. However, the SWHs simulated by the ST2 scheme were consistently larger than the measured values, especially after 8 September 2020. This suggests that the ST2 scheme dissipates energy more slowly than the ST4 and ST6 schemes, as evidenced by the slow decrease in SWH simulated by ST2 as wind speed decreased oscillating after 8 September 2020 [19]. Based on the scatter plot comparison in Figure 9 and the data presented in Table 2, it is evident that all of the source schemes exhibit a strong linear relationship with the measured significant wave heights (SWH). ST4 and ST6 show particularly promising results in terms of correlation coefficients, RMSE, and bias. Specifically, ST4 demonstrates the highest correlation coefficients of 0.860 and 0.861 at the respective points, indicating a strong relationship between the simulated and measured SWH. It is closely followed by ST6 with correlation coefficients of 0.845 and 0.835. Similarly, ST6 performs better in terms of RMSE and bias, which are smaller than those of the other schemes. ST4 is the second-best performing scheme, following closely behind ST6. Based on the analysis of the above results, ST4 and ST6 are the most suitable schemes for simulating the SWHs under the moderate sea state as they have the highest correlation coefficients and smaller RMSE and bias values.  At point 1, the performance of each source term indicates that ST2 has a poor correlation between the simulated and the measured data, as reflected by its high average relative error of 46.24%. In contrast, both the ST4 and ST6 schemes provide a better representation of SWH, with small and close average relative errors for both. The results at point 2 also support this conclusion. Overall, the analysis suggests that the ST4 and ST6 schemes are more suitable for simulating SWHs under moderate sea state conditions. Based on the scatter plot comparison in Figure 9 and the data presented in Table 2, it is evident that all of the source schemes exhibit a strong linear relationship with the measured significant wave heights (SWH). ST4 and ST6 show particularly promising results in terms of correlation coefficients, RMSE, and bias. Specifically, ST4 demonstrates the highest correlation coefficients of 0.860 and 0.861 at the respective points, indicating a strong relationship between the simulated and measured SWH. It is closely followed by ST6 with correlation coefficients of 0.845 and 0.835. Similarly, ST6 performs better in terms of RMSE and bias, which are smaller than those of the other schemes. ST4 is the second-best performing scheme, following closely behind ST6. Based on the analysis of the above results, ST4 and ST6 are the most suitable schemes for simulating the SWHs under the moderate sea state as they have the highest correlation coefficients and smaller RMSE and bias values.

Mean Wave Period Validation
The mean wave period was validated in a similar manner to the SWH, and the simulated and measured periods were compared ( Figure 10). The mean wave period simulated generally reflects the variation during this period, and the measured periods fluctuate up and down around 3 s. Each of the three schemes accurately describes the changes in the time period up to 9 September 2020, while the simulated mean wave periods for all three schemes after this date are larger than the measured values, with the ST4 scheme being significantly overestimated.

Mean Wave Period Validation
The mean wave period was validated in a similar manner to the SWH, and the simulated and measured periods were compared ( Figure 10). The mean wave period simulated generally reflects the variation during this period, and the measured periods fluctuate up and down around 3 s. Each of the three schemes accurately describes the changes in the time period up to 9 September 2020, while the simulated mean wave periods for all three schemes after this date are larger than the measured values, with the ST4 scheme being significantly overestimated.
In terms of the mean wave period simulation, ST2 exhibits the best performance for points 1 and 2. The average relative errors recorded for ST2 are the smallest at 14.10% and 14.29% for points 1 and 2, respectively. The ST6 is a close second in terms of performance. Therefore, the ST2 and ST6 schemes are more suitable for mean wave period simulation in moderate sea states. In terms of the mean wave period simulation, ST2 exhibits the best performance for points 1 and 2. The average relative errors recorded for ST2 are the smallest at 14.10% and 14.29% for points 1 and 2, respectively. The ST6 is a close second in terms of performance. Therefore, the ST2 and ST6 schemes are more suitable for mean wave period simulation in moderate sea states.
The input-dissipation source terms exhibit variable performance when simulated with different parameters. ST2 performs best when estimating mean wave periods, while ST4 and ST6 outperform other source terms in estimating SWHs in moderate sea states.

Significant Wave Height (SWH)
Three different schemes (referred to as ST2, ST4 and ST6) were employed to simulate the SWH trend. All three schemes provide a good description of the trend and show a peak at the same time as the measured data in Figure 11. However, the extreme values of all three simulations are overestimated. ST2 overestimates the extreme SWH by about 1.8 times the measured value, while ST4 and ST6 overestimate it by 2.3 and 2.6 times, respectively.  Figure 11. Comparison of SWH between simulated and measured data at point 1 (a) and 2 (b) during Vamco. The red dots and the three different colored lines (for ST2, ST4 or ST6 source terms) are the measured and simulated wave heights, respectively.
From the linear fit results, as shown in Table 3 and Figure 12, it can be observed that all three schemes exhibit a strong linear relationship with the measured values. However, when comparing the results in Table 3, it can be seen that the ST2 scheme has smaller RMSE and bias values, indicating a more accurate simulation of the typhoon wave evolution in the island area.
Considering the above results, it can be concluded that ST2 has a higher correlation coefficient with the measured values while maintaining a low average relative error. Therefore, compared to other source terms, ST2 is more suitable for simulating SWH in the island area under typhoon sea states. Table 3. Statistical data of SWH using different source term schemes during Vamco.

Source Scheme
Measurement Point Corrcoef RMSE Bias Figure 11. Comparison of SWH between simulated and measured data at point 1 (a) and 2 (b) during Vamco. The red dots and the three different colored lines (for ST2, ST4 or ST6 source terms) are the measured and simulated wave heights, respectively.
Comparing the SWH results, it is evident that the calculated values using ST2 are closer to the true values, as seen from the results at points 1 and 2. At point 1, the average relative errors for ST2, ST4 and ST6 are 16.95%, 17.7% and 77.19%, respectively. At point 2, the errors are 22.08%, 59.94% and 91.66%, respectively. The relative errors at both points follow the same order, indicating that ST4 and ST6 overestimate the SWH more before and after the typhoon passage, which suggests that the schemes reduce the dissipation term.
Based on the analysis of the wind velocity data and the simulation results of SWH, it can be concluded that the ST2 provides the closest estimation to the measured values of SWH during Vamco. Despite the overestimation of extreme values, the ST2 scheme still provides a good description of the peak of SWH if there is no error in the wind field. Therefore, the ST2 scheme is more accurate in describing the wind-generated wave process and fits better in the study of typhoon waves near the islands and reefs.
From the linear fit results, as shown in Table 3 and Figure 12, it can be observed that all three schemes exhibit a strong linear relationship with the measured values. However, when comparing the results in Table 3, it can be seen that the ST2 scheme has smaller RMSE and bias values, indicating a more accurate simulation of the typhoon wave evolution in the island area.   Considering the above results, it can be concluded that ST2 has a higher correlation coefficient with the measured values while maintaining a low average relative error. Therefore, compared to other source terms, ST2 is more suitable for simulating SWH in the island area under typhoon sea states.

Mean Wave Period Validation
The simulation results of the three schemes ( Figure 13) show that the mean wave period variation is relatively stable before the typhoon impact. However, after the typhoon passed, the mean wave period increased significantly and deviated from the trend of the measured values, which is inconsistent with the actual situation. None of the three source terms provide satisfactory results for the simulation of the mean wave period.

Mean Wave Period Validation
The simulation results of the three schemes ( Figure 13) show that the mean wave period variation is relatively stable before the typhoon impact. However, after the typhoon passed, the mean wave period increased significantly and deviated from the trend of the measured values, which is inconsistent with the actual situation. None of the three source terms provide satisfactory results for the simulation of the mean wave period.  After comparing the wave parameters mentioned above, it is evident that the simulation of the mean wave period is not as accurate as that of SWH, which is consistent with the previous analysis [54,55]. The simulation of typhoon-generated waves is a challenging task due to the rapid fluctuations in wind speed and direction [24]. Additionally, this marine environment is likely to be influenced by strong local currents [56,57]. The effect of current action will also be the next step to be explored.
Although the selection of source terms may depend on different regions and sea conditions, ST2 has been widely recognized for its good performance in simulating severe sea conditions, including typhoon sea states. For instance, Sheng et al. [26] found that the ST2 formula was the best option for simulating typhoon Fung-wong around the Zhoushan islands. Similarly, during Hurricane Katrina in the Gulf of Mexico, Montoya et al. [29] showed that the ST2 formula provided better results. Additionally, ST2 has also demonstrated good simulation performance under rough sea conditions [25,28,29].

Wave Spectrum
The wave spectrum is a statistical representation of wave surface elevation and is considered the most important form of wave description [58]. Numerous studies have focused on wave spectra [59][60][61]. In this section, the optimal source term determined in the previous section, ST2, is used to simulate the wave spectrum and compare the simulated results with the measured values.
During the typhoon transit, the wave spectrum (shown in Figure 14) exhibits a double peak spectrum mainly due to the presence of wind-sea states induced by local winds and waves from the open ocean around the island. By comparing the simulated wave spectrum with the field measurements, it is found that the shape and variance density of the wave spectrum simulated by ST2 agree well with the measured values, further demonstrating the accuracy of the ST2 simulations in the typhoon sea state.

Spatial Distribution of Significant Wave Heights
Using ST2, the optimal source term determined above, the spatial distribution of SWH under the typhoon sea state is studied, as shown in Figure 15. Prior to the approach of the typhoon, the SWH in the Yongle Atoll area was 1 m or less. As the typhoon approached, the SWH gradually increased throughout the region, particularly in the deep waters outside the island. The SWH peaked at 1:00 UTC on 14 November 2020, with peaks reaching 12 m in deep water. However, in the lagoon area of the Yongle Atoll, the peak occurred a few hours after the outer deep water area and did not exceed 4 m.

Spatial Distribution of Significant Wave Heights
Using ST2, the optimal source term determined above, the spatial distribution of SWH under the typhoon sea state is studied, as shown in Figure 15. Prior to the approach of the typhoon, the SWH in the Yongle Atoll area was 1 m or less. As the typhoon approached, the SWH gradually increased throughout the region, particularly in the deep waters outside the island. The SWH peaked at 1:00 UTC on 14 November 2020, with peaks reaching 12 m in deep water. However, in the lagoon area of the Yongle Atoll, the peak occurred a few hours after the outer deep water area and did not exceed 4 m.

Spatial Distribution of Significant Wave Heights
Using ST2, the optimal source term determined above, the spatial distribution of SWH under the typhoon sea state is studied, as shown in Figure 15. Prior to the approach of the typhoon, the SWH in the Yongle Atoll area was 1 m or less. As the typhoon approached, the SWH gradually increased throughout the region, particularly in the deep waters outside the island. The SWH peaked at 1:00 UTC on 14 November 2020, with peaks reaching 12 m in deep water. However, in the lagoon area of the Yongle Atoll, the peak occurred a few hours after the outer deep water area and did not exceed 4 m.

Summary and Conclusions
This study evaluates the performance of three popular wind input-dissipation source terms of the WAVEWATCH III numerical wave model in hindcasting ocean waves in the Yongle Atoll area of the SCS. As the physical processes of the moderate and typhoon sea states differ, the three source terms are compared graphically and statistically with the field wave data to identify the most suitable source term for each sea state.
In a moderate sea state with SWHs less than one meter, all three schemes show good

Summary and Conclusions
This study evaluates the performance of three popular wind input-dissipation source terms of the WAVEWATCH III numerical wave model in hindcasting ocean waves in the Yongle Atoll area of the SCS. As the physical processes of the moderate and typhoon sea states differ, the three source terms are compared graphically and statistically with the field wave data to identify the most suitable source term for each sea state.
In a moderate sea state with SWHs less than one meter, all three schemes show good results in the growth range. However, ST4 and ST6 perform better than ST2 in capturing the decaying process, indicating that dissipation is more sensitive. However, ST4 and ST6 do not perform well in simulating the mean wave period during and after the decaying process, whereas ST2 shows better results.
In the typhoon sea state, all three source terms tend to overestimate SWH, but ST2 provides the best results throughout the entire Typhoon Vamco process. The accuracy of SWH is highly influenced by wind velocity, and the shared increases between wind velocity and SWH peaks is consistent with the measured values. Therefore, obtaining more accurate wind speed and source terms is crucial for improving the accuracy of SWH predictions. ST2 is effective in simulating the shape and variance density of the wave spectrum, providing further evidence of its accuracy in typhoon sea states. However, all three source terms tend to overestimate the mean wave period, especially during the decaying process of the typhoon sea state. In sum, different input-dissipation packages are validated using measured wave data around the island, and the best packages vary for the typhoon sea state and the moderate sea state.
This study provides valuable insights into wave hindcasting in the Yongle Atoll area of the South China Sea. Valuable insights are provided for understanding wave dynamics and the impact of different source terms on wave simulation accuracy. These insights can inform and guide similar studies and wave modeling efforts to improve coastal hazard assessments and climate change impact studies.