  The honey bees play a role of unquestioned relevance in nature and the
comprehension of the mechanisms affecting their population dynamic is of
fundamental importance. As experimentally documented, the proper development of
a colony is related to the nest temperature, whose value is maintained around
the optimal value if the colony population is sufficiently large. Then, the
environmental temperature, the way in which this influence the nest temperature
and the colony population size, are variables closely linked to each other and
deserve to be taken into account in a model that aims to describe the
population dynamics. In the present study, as first step, the continuous-time
autonomous system proposed by Khoury, Myerscoug and Barron (KMB) in 2011 was
approximated by means a Nonstandard finite difference (NSFD) scheme in order to
obtain a set of autonomous difference equations. Subsequently, with the aim to
introduce the seasonal effects, a nonautonomous version (NAKMB) was proposed
and formulated in discrete-time domain via a NSFD scheme, by introducing a
time-dependent formulation for the queen bee laying rate and the recruitment
rate coefficients. By means the phase-plane analysis was possible to deduce
that, with an appropriate choice of the parameters, the NAKMB model admits both
a limit cycle at nonzero population size and an equilibrium point marking the
colony collapse, depending on the initial population size.
