Modeling bumble bee population dynamics with delay differential equations