Covariation between mean vole density and variability drives the numerical response of storks to vole prey

Hušek et al. (Popul Ecol 55:363–375, 2013 ) showed that the numerical response of storks to vole prey was stronger in regions where variability in vole density was higher. This finding is, at first sight, in contradiction with the predictions of life-history theory in stochastic environments. Since the stork productivity-vole density relationship is concave, theory predicts a negative association between the temporal variability in vole density and stork productivity. Here, we illustrate this negative effect of vole variability on stork productivity with a simple mathematical model relating expected stork productivity to vole dynamics. When comparing model simulations to the observed mean density and variability of thirteen Czech and Polish vole populations, we find that the observed positive effect of vole variability on stork numerical response is most likely due to an unusual positive correlation between mean and variability of vole density.