Modeling the spring blooms of ciliates in a deep lake

In contrast to the macro/mesozooplankton, microzooplankton has received much less attention in ecosystem models. In many modeling studies, microzooplankton has been either entirely neglected, or else, data were often not available for validation, or agreement between the observed and the simulated abundances was rather poor. In this study, we compare the simulation results from several alternative models considering different formulations of ciliate growth in a hydrodynamically driven 1D nutrient-phytoplankton–multiple zooplankton model, with long-term datasets from the deep, monomictic Lake Constance. We show that the parameterization of the limitation of ciliate growth with a constant specific mortality rate and/or predation by copepods leads to uncontrolled ciliate blooms. In contrast, implementation of a density-dependent mortality rate enables reproduction of algae–ciliate dynamics over a variety of environmental settings encompassed by the 14-year dataset spanning 21 years in a lake undergoing oligotrophication. Considering the numerous processes that can be responsible for the dampening of ciliate blooms, our findings suggest that employing a simple density-dependent mortality term offers a pragmatic solution for the challenge of including the microzooplankton, characterized by an overwhelming complexity of trophic interactions, in ecosystem models.