Based on the equations of motion and the assumption that ocean turbulence is of isotropy or quasi-isotropy, we derived the closure equations of the second-order moments and the variation equations for characteristic quantities, which describe the mechanisms of advection transport and shear instability by the sum of wave-like and eddy-like motions and circulation. Given that ocean turbulence generated by wave breaking is dominant at the ocean surface, we presented the boundary conditions of the turbulence kinetic energy and its dissipation rate, which are determined by energy loss from wave breaking and entrainment depth respectively. According to the equilibrium solution of the variation equations and available data of the dissipation rate, we obtained an analytical estimation of the characteristic quantities of surface-wave-generated turbulence in the upper ocean and its related mixing coefficient. The derived kinetic dissipation rate was validated by field measurements qualitatively and quantitatively, and the mixing coefficient had fairly good consistency with previous results based on the Prandtl mixing length theory.
