A numerical investigation of the impact of turbulence on the feeding rates of Oithona davisae

Individual based numerical simulations of the copepod, Oithona davisae, feeding on motile prey, Oxyrrhis marina, under variable turbulent conditions are performed. These simulations correspond to laboratory observations conducted by Saiz et al. [Saiz, E., Calbet, A., and Broglio, E., 2003. Effects of small-scale turbulence on copepods: the case of Oithona Davisae. Limnol. Oceanogr., 48:1304–1311.]. The flow field in the simulation is reconstructed by a kinematic simulation whose characteristic scales are derived from the grid mesh and the dissipation rates of the laboratory experiments. The kinematic simulation provides a simplified model, which while not fully realistic, captures the basic relevant feature of turbulence. A hop and sink swimming behaviour is prescribed for O. davisae, while O. marina moves along helical paths with random changes of directions. Three possible effects are tested: the existence of a time threshold in the duration of the contacts between predator and prey, a progressive reduction of the perceptive distance with increasing turbulence level and an abrupt reduction in feeding of O. davisae when the flow speed, in relation to the copepod position, is higher than a prescribed threshold. This last approach introduces an intermittency in the feeding which depends on the variations of velocity both in space and time within the numerical box. The introduction of the time threshold causes a dome-shaped relationship between the simulated enhancement factor and the dissipation rate, while with the other two effects, a monotonic decrease in the enhancement factor is observed, with values reasonably close to the ones observed in the laboratory experiment. In all the cases, the use of realistic values of biological parameters (e.g. swimming behaviour) reproduces response curves in the range of the observations.