Assessing the performance of formulations for nonlinear feedback of surface gravity waves on ocean currents over coastal waters

This study presents applications of a two-way coupled wave-circulation modelling system over coastal waters, with a special emphasis of performance assessments of two different methods for nonlinear feedback of ocean surface gravity waves on three-dimensional (3D) ocean currents. These two methods are the vortex force (VF) formulation suggested by Bennis et al. (2011) and the latest version of radiation stress (RS) formulation suggested by Mellor (2015). The coupled modelling system is first applied to two idealized test cases of surf-zone scales to validate implementations of these two methods in the coupled wave-circulation system. Model results show that the latest version of RS has difficulties in producing the undertow over the surf zone. The coupled system is then applied to Lunenburg Bay (LB) of Nova Scotia during Hurricane Juan in 2003. The coupled system using both the VF and RS formulations generates much stronger and more realistic 3D circulation in the Bay during Hurricane Juan than the circulation-only model, demonstrating the importance of surface wave forces to the 3D ocean circulation over coastal waters. However, the RS formulation generates some weak unphysical currents outside the wave breaking zone due to a less reasonable representation for the vertical distribution of the RS gradients over a slopping bottom. These weak unphysical currents are significantly magnified in a two-way coupled system when interacting with large surface waves, degrading the model performance in simulating currents at one observation site. Our results demonstrate that the VF formulation with an appropriate parameterization of wave breaking effects is able to produce reasonable results for applications over coastal waters during extreme weather events. The RS formulation requires a complex wave theory rather than the linear wave theory for the approximation of a vertical RS term to improve its performance under both breaking and non-breaking wave conditions.