A notable increase in failure of managed European honeybee Apis mellifera L. colonies has been reported in various regions in recent years. Although the underlying causes remain unclear, it is likely that a combination of stressors act together, particularly varroa mites and other pathogens, forage availability and potentially pesticides. It is experimentally challenging to address causality at the colony scale when multiple factors interact. In silico experiments offer a fast and cost-effective way to begin to address these challenges and inform experiments. However, none of the published bee models combine colony dynamics with foraging patterns and varroa dynamics.We have developed a honeybee model, BEEHAVE, which integrates colony dynamics, population dynamics of the varroa mite, epidemiology of varroa-transmitted viruses and allows foragers in an agent-based foraging model to collect food from a representation of a spatially explicit landscape.We describe the model, which is freely available online (www.beehave-model.net). Extensive sensitivity analyses and tests illustrate the model's robustness and realism. Simulation experiments with various combinations of stressors demonstrate, in simplified landscape settings, the model's potential: predicting colony dynamics and potential losses with and without varroa mites under different foraging conditions and under pesticide application. We also show how mitigation measures can be tested.Synthesis and applications. BEEHAVE offers a valuable tool for researchers to design and focus field experiments, for regulators to explore the relative importance of stressors to devise management and policy advice and for beekeepers to understand and predict varroa dynamics and effects of management interventions. We expect that scientists and stakeholders will find a variety of applications for BEEHAVE, stimulating further model development and the possible inclusion of other stressors of potential importance to honeybee colony dynamics.

The health of managed and wild honeybee colonies appears to have declined substantially in Europe and the United States over the last decade. Sustainability of honeybee colonies is important not only for honey production, but also for pollination of crops and wild plants alongside other insect pollinators. A combination of causal factors, including parasites, pathogens, land use changes and pesticide usage, are cited as responsible for the increased colony mortality.However, despite detailed knowledge of the behaviour of honeybees and their colonies, there are no suitable tools to explore the resilience mechanisms of this complex system under stress. Empirically testing all combinations of stressors in a systematic fashion is not feasible. We therefore suggest a cross-level systems approach, based on mechanistic modelling, to investigate the impacts of (and interactions between) colony and land management.We review existing honeybee models that are relevant to examining the effects of different stressors on colony growth and survival. Most of these models describe honeybee colony dynamics, foraging behaviour or honeybee – varroa mite – virus interactions.We found that many, but not all, processes within honeybee colonies, epidemiology and foraging are well understood and described in the models, but there is no model that couples in-hive dynamics and pathology with foraging dynamics in realistic landscapes.Synthesis and applications. We describe how a new integrated model could be built to simulate multifactorial impacts on the honeybee colony system, using building blocks from the reviewed models. The development of such a tool would not only highlight empirical research priorities but also provide an important forecasting tool for policy makers and beekeepers, and we list examples of relevant applications to bee disease and landscape management decisions.

World‐wide declines in pollinators, including bumblebees, are attributed to a multitude of stressors such as habitat loss, resource availability, emerging viruses and parasites, exposure to pesticides, and climate change, operating at various spatial and temporal scales. Disentangling individual and interacting effects of these stressors, and understanding their impact at the individual, colony and population level are a challenge for systems ecology. Empirical testing of all combinations and contexts is not feasible. A mechanistic multilevel systems model (individual‐colony‐population‐community) is required to explore resilience mechanisms of populations and communities under stress.We present a model which can simulate the growth, behaviour and survival of six UK bumblebee species living in any mapped landscape. Bumble‐BEEHAVE simulates, in an agent‐based approach, the colony development of bumblebees in a realistic landscape to study how multiple stressors affect bee numbers and population dynamics. We provide extensive documentation, including sensitivity analysis and validation, based on data from literature. The model is freely available, has flexible settings and includes a user manual to ensure it can be used by researchers, farmers, policy‐makers, NGOs or other interested parties.Model outcomes compare well with empirical data for individual foraging behaviour, colony growth and reproduction, and estimated nest densities.Simulating the impact of reproductive depression caused by pesticide exposure shows that the complex feedback mechanisms captured in this model predict higher colony resilience to stress than suggested by a previous, simpler model. Synthesis and applications. The Bumble‐BEEHAVE model represents a significant step towards predicting bumblebee population dynamics in a spatially explicit way. It enables researchers to understand the individual and interacting effects of the multiple stressors affecting bumblebee survival and the feedback mechanisms that may buffer a colony against environmental stress, or indeed lead to spiralling colony collapse. The model can be used to aid the design of field experiments, for risk assessments, to inform conservation and farming decisions and for assigning bespoke management recommendations at a landscape scale.