From winter to summer and back: Lessons from the parameterization of a seasonal food web model for the Białowieża forest.

Dynamic food web models describe how species abundances change over time as a function of trophic and life-history parameters. They are essential to predicting the response of ecosystems to perturbations. However, they are notoriously difficult to parameterize, so that most models rely heavily either on allometric scaling of parameters or inverse estimation of biomass flows. The allometric approach makes species of comparable body mass have near-identical parameters which can generate extinctions within a trophic level. The biomass flow approach is more precise, but is restricted to steady-states, which is not appropriate for time-varying environments. Adequately parameterizing large food webs of temperate and arctic environments requires dealing both with many species of similar sizes and a strongly seasonal environment. Inspired by the rich empirical knowledge on the vertebrate food web of the Białowieża forest, we parameterize a bipartite food web model comprising 21 predators and 124 prey species. Our model is a non-autonomous coupled ordinary differential equations system that allows for seasonality in life-history and predation parameters. Birth and death rates, seasonal descriptions of diet for each species, food requirements and biomass information are combined into a seasonal parameterization of a dynamic food web model. Food web seasonality is implemented with time-varying intrinsic growth rate and interaction parameters, while predation is modelled with both type I and type II functional responses. All our model variants allow for >80% persistence in spite of massive apparent competition, and a quantitative match to observed (seasonal) biomasses. We also identify trade-offs between maximizing persistence, reproducing observed biomasses, and ensuring model robustness to sampling errors. Although multi-annual cycles are expected with type II functional responses, they are here prevented by a strong predator self-regulation. We discuss these results and possible improvements on the model. We provide a general workflow to parameterize dynamic food web models in seasonal environments, based on a real case study. This may help to better predict how biodiverse food webs respond to changing environments.