THEORETICAL AND NUMERICAL STUDY OF WAVE-CURRENT INTERACTION IN STRONGLY-SHEARED FLOWS

The application of wave-current interaction theory in ocean circulation models has been extensively developed over the past decade, with formulations extended to three dimensions and based either on radiation stress formulations or on the Craik-Leibovich formulation. However, few of these studies consider the interaction of waves with relatively strongly sheared current, in which current shear can affect linear wave dynamics at leading order. The problem arises from the study of the evolution of highly concentrated sediment plumes developing at the mouth of small mountainous rivers. Although the annually averaged discharge of these small mountainous rivers is trivial compared to large rivers, during the extreme flooding events triggered by typhoon or tropic cyclones, these rivers, most of which located at tectonically active mountain belts, can carry highly concentrated sediment ( up to several g/l in the river plume) into the ocean. The magnitude of river discharge velocity at the river mouth may reach several m/s, comparable to the wave phase speed in coastal water. In addition, these flooding events usually coincide with very energetic wave conditions induced by the storms. Therefore, the interaction of waves with strongly sheared current becomes a very important dynamic process at this kind of river plumes. In our study, we establish a new framework to describe the interaction of small amplitude surface gravity waves and strongly sheared currents, where shear can exist in both vertical and horizontal directions. To begin with, we limit the derivation to the case of a narrow-banded slowly varying wave train propagating shoreward in the coastal ocean outside of the surf zone. Accordingly, our problem is assumed to be finite depth without wave breaking. Later we can extend the formulation to describe a spectrum of surface waves and include wave energy dissipation. In contrast to existing formulations, where waves at most feel a weighted depth-average current which follows from a weak-current, weak-shear approximation, the present formulation allows for an arbitrary degree of vertical shear, leading to a description of the vertical structure of waves in terms of solutions to the Rayleigh stability equation. The resulting formulation leads to a conservation law for wave action, and forcing terms for the description of mean flow using the Craik-Leibovich vortex force formulation. This new framework of wave-current interaction can be applied to numerical model based on ROMS/SWAN to study dynamics in coastal waters.