This paper investigates the interaction between stratification and turbulence by means of turbulence models. The standard and the advanced turbulent kinetic energy -dissipation (k-e) model are derived theoretically, including algebraic stress relations. It is shown that a certain empirical constant in the standard model turns out to be a complicated imphcit function in the advanced model, namely, a function of the turbulent shear number, the turbulent buoyancy number, and a wall correction. For a better understanding and physical interpretation of the k-e models, an analysis is carried out for a simplified case where diffusive fluxes are neglected. For this ideahzation it is shown that (1) the flux Richardson number R] has a certain lower bound R• due to the estabhshment of convection, (2) a steady state flux Richardson number R} t (which is defined here for this purpose) labels the borderline between the tendency of turbulence to decrease or collapse (R] > R} t) or to increase (R] < R}t), and (3) the wen-known upper limit for turbulent shear flow, R• • 0.25, fits our theory. Using the standard model, the advanced model, a modified version of the level-2 model of Mellor and Yamada and a modified version of Kochergin's model, the evolution of thermal stratification in the northern North Sea during the Fladen Ground Experiment (FLEX'76) is simulated numerically and compared with the measurements. In this specific application, the two k-e models performed best. 1. Introduction Ocean-atmosphere interactions are strongly determined by the dynamics of the upper mixed layer and the evolution and erosion of thermoclines. In the water column we find a complex interplay of several processes forced by different fluxes through the boundaries: momentum flux at the surface due to wind and at the bottom due to friction, and buoyancy flux through the surface due to heat flux, evaporation, precipitation or sea ice dynamics and through the bottom due to erosion and deposition of sediment and suspended particulate matter (for recent observations of highly complex phenomena in the upper mixed layer see Mourn and Caldwell [1994]). There are strong motivations to deal with the upper mixed layer of the ocean, including fishery, primary production and eutrophication, and the occurrence of toxic algae (red tide). Finally, our climate dynamics is closely coupled to the storage of heat and CO2 in the upper mixed layer of the world ocean. In shelf seas and estuaries, we often find a lower mixed layer at the bottom, where important geomorphological and transport processes take place. This layer can be separated from the upper mixed layer by a region of stable stratification. If such a stratified region is absent, both mixed layers condense into only one. Physical and numerical models are an important tool for understanding and predicting the complicated processes in these regions of the ocean. Since the first comprehensive review of upper mixed-layer modelling by Kraus [1977], a large number of turbulence closure models were presented. They ...