Abstract
Extreme surface winds and wave heights of tropical cyclones (TCs)—pose serious threats to coastal community, infrastructure and environments. In recent decades, progress in numerical wave modeling has significantly enhanced the ability to reconstruct and predict wave behavior. This review offers an in-depth overview of TC-related wave modeling utilizing different computational schemes, with a special attention to WAVEWATCH III (WW3) and Simulating Waves Nearshore (SWAN). Due to the complex air–sea interactions during TCs, it is challenging to obtain accurate wind input data and optimize the parameterizations. Substantial spatial and temporal variations in water levels and current patterns occurs when coastal circulation is modulated by varying underwater topography. To explore their influence on waves, this study employs a coupled SWAN and Finite-Volume Community Ocean Model (FVCOM) modeling approach. Additionally, the interplay between wave and sea surface temperature (SST) is investigated by incorporating four key wave-induced forcing through breaking and non-breaking waves, radiation stress, and Stokes drift from WW3 into the Stony Brook Parallel Ocean Model (sbPOM). 20 TC events were analyzed to evaluate the performance of the selected parameterizations of external forcings in WW3 and SWAN. Among different nonlinear wave interaction schemes, Generalized Multiple Discrete Interaction Approximation (GMD) Discrete Interaction Approximation (DIA) and the computationally expensive Wave-Ray Tracing (WRT) A refined drag coefficient (Cd) equation, applied within an upgraded ST6 configuration, reduce significant wave height (SWH) prediction errors and the root mean square error (RMSE) for both SWAN and WW3 wave models. Surface currents and sea level variations notably altered the wave energy and wave height distributions, especially in the area with strong TC-induced oceanic current. Finally, coupling four wave-induced forcings into sbPOM enhanced SST simulation by refining heat flux estimates and promoting vertical mixing. Validation against Argo data showed that the updated sbPOM model achieved an RMSE as low as 1.39 m, with correlation coefficients nearing 0.9881.