Mathematical Analysis for Honeybee Dynamics Under the Influence of Seasonality

In this paper, we studied a mathematical model for honeybee population diseases under the influence of seasonal environments on the long-term dynamics of the disease. The model describes the dynamics of two different beehives sharing a common space. We computed the basic reproduction number of the system as the spectral radius of either the next generation matrix for the autonomous system or as the spectral radius of a linear integral operator for the non-autonomous system, and we deduced that if the reproduction number is less than unity, then the disease dies out in the honeybee population. However, if the basic reproduction number is greater than unity, then the disease persists. Finally, we provide several numerical tests that confirm the theoretical findings.