A simple model of plankton population dynamics coupled with a LES of the surface mixed layer

The concentration of phytoplankton in the sea is affected by biological processes, such as growth/mortality rates, predatory zooplankton concentrations and nutrient levels. Phytoplankton concentrations are also influenced by physical processes, in particular the mixing properties of the local fluid environment. On planktonic scales (∼10-1000μm) one can assume the local turbulent flow is isotropic, with no distinction between horizontal and vertical mixing. However, agglomerations of phytoplankton into patches are observed on larger scales of up to hundreds of metres, whose formation will be influenced by the anisotropic advection/mixing properties and large-eddy structures prevalent in the surface mixed layer. This paper presents the results of the coupling of a large-eddy simulation (LES) model of the mixed layer with an advection-diffusion system of coupled equations for nitrate-phytoplankton-zooplankton (NPZ) concentration, incorporating sub-grid parameterizations of the biological processes. Typically these include phytoplankton growth due to light levels and ambient nitrate concentration, offset by grazing losses due to the presence of zooplankton. The primary goal of this work is to investigate how the characteristics of the mixed layer turbulence influence the observed distribution of phytoplankton. One novel feature is the incorporation of a 'vortex-force' term in the LES code in order to generate Langmuir circulations. It has been speculated that the enhanced mixing rates associated with 'Langmuir turbulence' play a significant role in regulating planktonic activity. Results derived from the coupled LES-NPZ model, run with and without the presence of Langmuir circulations, are presented in order to investigate these ideas.