Simulation of Varroa mite control in honey bee colonies without synthetic acaricides: Demonstration of Good Beekeeping Practice for Germany in the BEEHAVE model.

The BEEHAVE model simulates the population dynamics and foraging activity of a single honey bee colony (Apis mellifera) in great detail. Although it still makes numerous simplifying assumptions, it appears to capture a wide range of empirical observations. It could, therefore, in principle, also be used as a tool in beekeeper education, as it allows the implementation and comparison of different management options. Here, we focus on treatments aimed at controlling the mite Varroa destructor. However, since BEEHAVE was developed in the UK, mite treatment includes the use of a synthetic acaricide, which is not part of Good Beekeeping Practice in Germany. A practice that consists of drone brood removal from April to June, treatment with formic acid in August/September, and treatment with oxalic acid in November/December. We implemented these measures, focusing on the timing, frequency, and spacing between drone brood removals. The effect of drone brood removal and acid treatment, individually or in combination, on a mite-infested colony was examined. We quantify the efficacy of Varroa mite control as the reduction of mites in treated bee colonies compared to untreated bee colonies. We found that drone brood removal was very effective, reducing mites by 90% at the end of the first simulation year after the introduction of mites. This value was significantly higher than the 50-67% reduction expected by bee experts and confirmed by empirical studies. However, literature reports varying percent reductions in mite numbers from 10 to 85% after drone brood removal. The discrepancy between model results, empirical data, and expert estimates indicate that these three sources should be reviewed and refined, as all are based on simplifying assumptions. These results and the adaptation of BEEHAVE to the Good Beekeeping Practice are a decisive step forward for the future use of BEEHAVE in beekeeper education in Germany and anywhere where organic acids and drone brood removal are utilized.

Long-term effects of antibiotic treatments on honeybee colony fitness: A modelling approach.

Gut microbiome disequilibrium is increasingly implicated in host fitness reductions, including for the economically important and disease-challenged western honey bee Apis mellifera. In laboratory experiments, the antibiotic tetracycline, which is used to prevent American Foulbrood Disease in countries including the US, elevates honey bee mortality by disturbing the microbiome. It is unclear, however, how elevated individual mortality affects colony-level fitness.We used an agent-based model (BEEHAVE) and empirical data to assess colony-level effects of antibiotic-induced worker bee mortality, by measuring colony size. We investigated the relationship between the duration that the antibiotic-induced mortality probability is imposed for and colony size.We found that when simulating antibiotic-induced mortality of worker bees from just 60 days per year, up to a permanent effect, the colony is reduced such that tetracycline treatment would not meet the European Food Safety Authority's (EFSA) honey bee protection goals. When antibiotic mortality was imposed for the hypothetical minimal exposure time, which assumes that antibiotics only impact the bee's fitness during the recommended treatment period of 15 days in both spring and autumn, the colony fitness reduction was only marginally under the EFSA's threshold. Synthesis and Applications. Modelling colony-level impacts of antibiotic treatment shows that individual honey bee worker mortality can lead to colony mortality. To assess the full impact, the persistence of antibiotic-induced mortality in honey bees must be determined experimentally, in vivo. We caution that as the domestication of new insect species increases, maintaining healthy gut microbiomes is of paramount importance to insect health and commercial productivity. The recommendation from this work is to limit prophylactic use of antibiotics and to not exceed recommended treatment strategies for domesticated insects. This is especially important for highly social insects as excess antibiotic use will likely decrease colony growth and an increase in colony mortality.

Honeybee colonies compensate for pesticide-induced effects on royal jelly composition and brood survival with increased brood production.

Sublethal doses of pesticides affect individual honeybees, but colony-level effects are less well understood and it is unclear how the two levels integrate. We studied the effect of the neonicotinoid pesticide clothianidin at field realistic concentrations on small colonies. We found that exposure to clothianidin affected worker jelly production of individual workers and created a strong dose-dependent increase in mortality of individual larvae, but strikingly the population size of capped brood remained stable. Thus, hives exhibited short-term resilience. Using a demographic matrix model, we found that the basis of resilience in dosed colonies was a substantive increase in brood initiation rate to compensate for increased brood mortality. However, computer simulation of full size colonies revealed that the increase in brood initiation led to severe reductions in colony reproduction (swarming) and long-term survival. This experiment reveals social regulatory mechanisms on colony-level that enable honeybees to partly compensate for effects on individual level.

Honey bee colony performance affected by crop diversity and farmland structure: a modeling framework.

Forage availability has been suggested as one driver of the observed decline in honey bees. However, little is known about the effects of its spatiotemporal variation on colony success. We present a modeling framework for assessing honey bee colony viability in cropping systems. Based on two real farmland structures, we developed a landscape generator to design cropping systems varying in crop species identity, diversity, and relative abundance. The landscape scenarios generated were evaluated using the existing honey bee colony model BEEHAVE, which links foraging to in-hive dynamics. We thereby explored how different cropping systems determine spatiotemporal forage availability and, in turn, honey bee colony viability (e.g., time to extinction, TTE) and resilience (indicated by, e.g., brood mortality). To assess overall colony viability, we developed metrics, PH and PP, which quantified how much nectar and pollen provided by a cropping system per year was converted into a colony's adult worker population. Both crop species identity and diversity determined the temporal continuity in nectar and pollen supply and thus colony viability. Overall farmland structure and relative crop abundance were less important, but details mattered. For monocultures and for four-crop species systems composed of cereals, oilseed rape, maize, and sunflower, PH and PP were below the viability threshold. Such cropping systems showed frequent, badly timed, and prolonged forage gaps leading to detrimental cascading effects on life stages and in-hive work force, which critically reduced colony resilience. Four-crop systems composed of rye-grass-dandelion pasture, trefoil-grass pasture, sunflower, and phacelia ensured continuous nectar and pollen supply resulting in TTE > 5 yr, and PH (269.5 kg) and PP (108 kg) being above viability thresholds for 5 yr. Overall, trefoil-grass pasture, oilseed rape, buckwheat, and phacelia improved the temporal continuity in forage supply and colony's viability. Our results are hypothetical as they are obtained from simplified landscape settings, but they nevertheless match empirical observations, in particular the viability threshold. Our framework can be used to assess the effects of cropping systems on honey bee viability and to develop land-use strategies that help maintain pollination services by avoiding prolonged and badly timed forage gaps.

Bombus terrestris in a mass-flowering pollinator-dependent crop: A mutualistic relationship?

Bumblebees (Bombus spp.) rely on an abundant and diverse selection of floral resources to meet their nutritional requirements. In farmed landscapes, mass-flowering crops can provide an important forage resource for bumblebees, with increased visitation from bumblebees into mass-flowering crops having an additional benefit to growers who require pollination services. This study explores the mutualistic relationship between Bombus terrestris L. (buff-tailed bumblebee), a common species in European farmland, and the mass-flowering crop courgette (Cucurbita pepo L.) to see how effective B. terrestris is at pollinating courgette and in return how courgette may affect B. terrestris colony dynamics. By combining empirical data on nectar and pollen availability with model simulations using the novel bumblebee model Bumble-BEEHAVE, we were able to quantify and simulate for the first time, the importance of courgette as a mass-flowering forage resource for bumblebees. Courgette provides vast quantities of nectar to ensure a high visitation rate, which combined with abundant pollen grains, enables B. terrestris to have a high pollination potential. While B. terrestris showed a strong fidelity to courgette flowers for nectar, courgette pollen was not found in any pollen loads from returning foragers. Nonetheless, model simulations showed that early season courgette (nectar) increased the number of hibernating queens, colonies, and adult workers in the modeled landscapes. Synthesis and applications. Courgette has the potential to improve bumblebee population dynamics; however, the lack of evidence of the bees collecting courgette pollen in this study suggests that bees can only benefit from this transient nectar source if alternative floral resources, particularly pollen, are also available to fulfill bees' nutritional requirements in space and time. Therefore, providing additional forage resources could simultaneously improve pollination services and bumblebee populations.

Bumble-BEEHAVE: A systems model for exploring multifactorial causes of bumblebee decline at individual, colony, population and community level.

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.

Modeling Effects of Honeybee Behaviors on the Distribution of Pesticide in Nectar within a Hive and Resultant in-Hive Exposure.

Recently, the causes of honeybee colony losses have been intensely studied, showing that there are multiple stressors implicated in colony declines, one stressor being the exposure to pesticides. Measuring exposure of individual bees within a hive to pesticide is at least as difficult as assessing the potential exposure of foraging bees to pesticide. We present a model to explore how heterogeneity of pesticide distribution on a comb in the hive can be driven by worker behaviors. The model contains simplified behaviors to capture the extremes of possible heterogeneity of pesticide location/deposition within the hive to compare with exposure levels estimated by averaging values across the comb. When adults feed on nectar containing the average concentration of all pesticide brought into the hive on that particular day, it is likely representative of the worst-case exposure scenario. However, for larvae, clustering of pesticide in the comb can lead to higher exposure levels than taking an average concentration in some circumstances. The potential for extrapolating the model to risk assessment is discussed.

Predicting Honeybee Colony Failure: Using the BEEHAVE Model to Simulate Colony Responses to Pesticides.

To simulate effects of pesticides on different honeybee (Apis mellifera L.) life stages, we used the BEEHAVE model to explore how increased mortalities of larvae, in-hive workers, and foragers, as well as reduced egg-laying rate, could impact colony dynamics over multiple years. Stresses were applied for 30 days, both as multiples of the modeled control mortality and as set percentage daily mortalities to assess the sensitivity of the modeled colony both to small fluctuations in mortality and periods of low to very high daily mortality. These stresses simulate stylized exposure of the different life stages to nectar and pollen contaminated with pesticide for 30 days. Increasing adult bee mortality had a much greater impact on colony survival than mortality of bee larvae or reduction in egg laying rate. Importantly, the seasonal timing of the imposed mortality affected the magnitude of the impact at colony level. In line with the LD50, we propose a new index of "lethal imposed stress": the LIS50 which indicates the level of stress on individuals that results in 50% colony mortality. This (or any LISx) is a comparative index for exploring the effects of different stressors at colony level in model simulations. While colony failure is not an acceptable protection goal, this index could be used to inform the setting of future regulatory protection goals.

BEEHAVE: a systems model of honeybee colony dynamics and foraging to explore multifactorial causes of colony failure.

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.

Towards a systems approach for understanding honeybee decline: a stocktaking and synthesis of existing models.

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.

Pupal developmental temperature and behavioral specialization of honeybee workers (Apis mellifera L.).

Honeybees (Apis mellifera) are able to regulate the brood nest temperatures within a narrow range between 32 and 36 degrees C. Yet this small variation in brood temperature is sufficient to cause significant differences in the behavior of adult bees. To study the consequences of variation in pupal developmental temperature we raised honeybee brood under controlled temperature conditions (32, 34.5, 36 degrees C) and individually marked more than 4,400 bees, after emergence. We analyzed dancing, undertaking behavior, the age of first foraging flight, and forager task specialization of these workers. Animals raised under higher temperatures showed an increased probability to dance, foraged earlier in life, and were more often engaged in undertaking. Since the temperature profile in the brood nest may be an emergent property of the whole colony, we discuss how pupal developmental temperature can affect the overall organization of division of labor among the individuals in a self-organized process.