Implementation of a 3D ocean model to understand upland lake wind-driven circulation

A community numerical ocean model is used to extend the understanding of wind-driven circulation in small upland lakes. A 3D model of a case study lake (Llyn Conwy, Wales, UK) is calibrated against measured velocity profiles via adjustment of the bottom roughness coefficient. Validation against a separate set of measured velocity profiles confirms the ability of the model to resolve key features of the flow field. Sensitivity analysis shows that the velocity field responds rapidly to changes in the wind forcing. Analysis of the gross circulation using Empirical Orthogonal Functions reveals a persistent two-gyre circulation pattern in the upper half layer of the water column driven by the interaction of wind and bathymetry. At the bottom, the flow is characterised by locally strong currents and analysis of vertical circulation over short time scales shows strong currents in the deepest parts of the lake basin and the responsiveness of the water column to changes in wind speed and direction. Even in small lakes, the assumption of uniform wind stress across the water surface is not always justified and topographic sheltering or other catchment roughness effects give rise to heterogeneity in the wind field. An idealized experiment for the case study lake shows that differences in circulation emerge if the wind stress is allowed to vary across the lake. Energetic wind forcing in upland areas can drive an energetic lake circulation that has important implications for mixing and sediment dynamics. 3D numerical modelling of wind-driven circulation should be more widely used to provide insights into physical limnology to support a wide range of ecological, biogeochemical and palaeoenvironmental studies.