Classification and dynamic simulation of the vertical density structure of lakes1

Abstract

Field data from two lakes of widely differing geometry and size are analyzed in terms of four nondimensional numbers which allow the principal mixing processes in each lake to be identified. The numbers are based on basin geometry, density stratification, wind stress, and rates of inflow and outflow. The procedure highlights the differences in the dynamics of the two lakes and allows assessment of the validity of the assumption of one‐dimensionality. The result is that both lakes were dominated by one‐dimensional, but different processes. The dynamics of the epilimnion of the smaller lake were dominated by stirring from surface wind and cooling, whereas shear at the pycnocline was also significant in the larger lake. In neither case did the effects of the earth’s rotation, inflow, or outflow generate significant horizontal gradients. A one‐dimensional numerical model (DYRESM) was used to simulate the vertical temperature and salinity structures of both lakes over lengthy periods, with good results. The model is based on the parameterization of the important physical processes in a framework of horizontal layers of variable thickness and was applied in both lakes without alteration. The interpretative power of the model is demonstrated by examination of the formation and erosion of a thermal inversion in the larger lake.