Unravelling the impact of climate change on honey bees: An ensemble modelling approach to predict shifts in habitat suitability in Queensland, Australia

Abstract

Honey bees play a vital role in providing essential ecosystem services and contributing to global agriculture. However, the potential effect of climate change on honey bee distribution is still not well understood. This study aims to identify the most influential bioclimatic and environmental variables, assess their impact on honey bee distribution, and predict future distribution. An ensemble modelling approach using the biomod2 package in R was employed to develop three models: a climate-only model, an environment-only model, and a combined climate and environment model. By utilising bioclimatic data (radiation of the wettest and driest quarters and temperature seasonality) from 1990 to 2009, combined with observed honey bee presence and pseudo absence data, this model predicted suitable locations for honey bee apiaries for two future time spans: 2020-2039 and 2060-2079. The climate-only model exhibited a true skill statistic (TSS) value of 0.85, underscoring the pivotal role of radiation and temperature seasonality in shaping honey bee distribution. The environment-only model, incorporating proximity to floral resources, foliage projective cover, and elevation, demonstrated strong predictive performance, with a TSS of 0.88, emphasising the significance of environmental variables in determining habitat suitability for honey bees. The combined model had a higher TSS of 0.96, indicating that the combination of climate and environmental variables enhances the model's performance. By the 2020-2039 period, approximately 88% of highly suitable habitats for honey bees are projected to transition from their current state to become moderate (14.84%) to marginally suitable (13.46%) areas. Predictions for the 2060-2079 period reveal a concerning trend: 100% of highly suitable land transitions into moderately (0.54%), marginally (17.56%), or not suitable areas (81.9%) for honey bees. These results emphasise the critical need for targeted conservation efforts and the implementation of policies aimed at safeguarding honey bees and the vital apiary industry.

Keywords: Apis mellifera; biomod2; climate change; ensemble modelling; honey bees; species distribution modelling.