Numerical Modeling of Wave-Enhanced Turbulence in the Oceanic Upper Layer

A coupled model of air-wave-sea interaction is modified based on a new roughness formulation and the latest data. The model parameters for aerodynamic roughness from below (ARB) and wave-dependent roughness from above (ARA, z 0a ) are assumed equal. The combined roughness is assumed to be a function of friction velocity, gravity, air and seawater densities, and wave age (c w ). The model is used in a study of wave-enhanced turbulence under breaking waves to predict turbulent dissipation ( e), ARA, and drag coefficient (C d ). Both waves and shear production are considered as sources of ocean turbulent energy. The atmospheric part of the model is used only to specify a correct condition at the interface. Numerical experiments are performed to study the e-distribution, z 0a and C d , and to compare with data. The major achievement is model verification using all available data. The first full application of this model is in conjunction with an ocean circulation model in a coupled circulationwave system. Simulations show that the e-distribution is strongly dependent on local wind-forced wave heights. For each wind and wave state there is a particular wavedependent depth that is verified by data. The comparison shows that the model predicted e agrees well with the observed e of the z ‐4 law distribution of Gargett (1989). Simulations also show that waves have an important role in causing e to differ from the classical wall-layer theory and z 0a , with a value of 0.30 for the empirical constant a a . The model-predicted e, z 0a , C d and C gd agree well with data.