Abstract
This paper describes a theoretical model of the ocean surface mixed layer. A nonstationary one-dimensional system of equations for an Ekman ocean layer and for heat conduction is closed with a system of dynamic turbulence equations. The latter consists of the turbulent energy equation and an equation for the turbulent energy decay function. The eddy viscosity coefficient is then determined from these equations. The values of the surface and lower boundary temperature and of wind velocity as a function of time are taken from Halpern (1974) data. The computational results show that an increase of wind velocity produces deepening of the mixed layer with some time lag. In the region of the jump in density at the bottom of the layer, the increase in the surface wind generates large velocity gradients across an inner boundary layer. Comparison of the solution with experimental data shows that the model realistically simulates the dynamics of mixed layer deepening both qualitatively and quantitatively.