Next Article in Journal
Spatial Origin and Diversification of the Lycocerus fainanus Species Group (Coleoptera, Cantharidae), with Descriptions of Four New Species from China and Vietnam
Next Article in Special Issue
Global Warming, Advancing Bloom and Evidence for Pollinator Plasticity from Long-Term Bee Emergence Monitoring
Previous Article in Journal
Combination of Modified Atmosphere and Irradiation for the Phytosanitary Disinfestation of Trogoderma granarium Everts (Coleoptera: Dermestidae)
Previous Article in Special Issue
Joint Impacts of Drought and Habitat Fragmentation on Native Bee Assemblages in a California Biodiversity Hotspot
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 12
Issue 5
10.3390/insects12050443
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ensemble Models Predict Invasive Bee Habitat Suitability Will Expand under Future Climate Scenarios in Hawai’i

by
Jesse A. Tabor
1,2 and
Jonathan B. Koch
3,4,*
1
Department of Geography & Environmental Studies, University of Hawai’i, 200 W. Kāwili Street, Hilo, HI 96720, USA
2
Department of Biology, Utah State University, 5305 Old Main Hill, Logan, UT 84322, USA
3
Tropical Conservation Biology & Environmental Science Graduate Program, University of Hawai’i, Hilo, 200 W. Kāwili Street, Hilo, HI 96720, USA
4
Pollinating Insect—Biology, Management, and Systematics Research Unit, U.S. Department of Agriculture—Agricultural Research Service, 1410 N. 800 E., Logan, UT 84341, USA
*
Author to whom correspondence should be addressed.
Insects 2021, 12(5), 443; https://doi.org/10.3390/insects12050443
Submission received: 9 April 2021 / Revised: 28 April 2021 / Accepted: 11 May 2021 / Published: 13 May 2021
(This article belongs to the Special Issue Non-Apis Pollinators and Global Change)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Climate change exacerbates the threat of biological invasions by increasing climatically suitable regions for species to invade outside of their native range. Island ecosystems may be particularly sensitive to the synergistic effects of climate change and biological invasions. In Hawai’i there are 21 non-native bees that have the capacity to spread pathogens and compete for resources with native bees. We performed an ensemble of species distribution models (SDM) for eight non-native bee species (Hymenoptera: Anthophila) in Hawai’i to predict climatically suitable niches across current and future climate scenarios. We found a significant difference in habitat suitability between SDMs that were constructed with specimen records from their native and non-native (Hawai’i) range. Although SDMs predict expansion of suitable habitat into higher elevations under 2070 climate scenarios, species-rich areas are predicted to stay below 500 m elevation. Our models can inform decisions on the management of non-native bees in Hawai’i by assessing risk of invasion into new areas around the archipelago.

Abstract

Climate change is predicted to increase the risk of biological invasions by increasing the availability of climatically suitable regions for invasive species. Endemic species on oceanic islands are particularly sensitive to the impact of invasive species due to increased competition for shared resources and disease spread. In our study, we used an ensemble of species distribution models (SDM) to predict habitat suitability for invasive bees under current and future climate scenarios in Hawai’i. SDMs projected on the invasive range were better predicted by georeferenced records from the invasive range in comparison to invasive SDMs predicted by records from the native range. SDMs estimated that climatically suitable regions for the eight invasive bees explored in this study will expand by ~934.8% (±3.4% SE). Hotspots for the invasive bees are predicted to expand toward higher elevation regions, although suitable habitat is expected to only progress up to 500 m in elevation in 2070. Given our results, it is unlikely that invasive bees will interact directly with endemic bees found at >500 m in elevation in the future. Management and conservation plans for endemic bees may be improved by understanding how climate change may exacerbate negative interactions between invasive and endemic bee species.

1. Introduction

Biological invasions are one of the most severe threats to biodiversity and natural resources. Climate change is predicted to increase the risk of biological invasions by increasing climatically suitable regions for invasive species [1]. The intersection of climate change and invasive species is considered to be one of the main factors likely to impact bee diversity, together with land use change, exposure to pesticides, and pathogen spillover [2]. Climate constraints on invasive bee species, specifically temperature and precipitation patterns, may be reduced at their range limits, potentially allowing invasive species to expand beyond their current geographic range and into novel environments [3,4,5,6]. In addition to changes in geographic distribution of species, climate change has been shown to impact plant-pollinator phenology mismatch, bee genetic composition and body size, and species interactions [7,8,9,10].
Hawai’i is one of the most remote archipelagos in the world and has been invaded by many different invasive plants and animals, including at least 21 species of introduced bees (Table 1 and Figure 1) [11,12,13,14,15]. The invasive bees present in Hawai’i have been accidentally introduced, excluding Apis mellifera Linnaeus, 1758, which was brought to Hawai’i in 1857 primarily to deliver pollination services to non-native Prosopis pallida (Fabaceae) to support the cattle industry, and subsequently honey production in 1850 [16]. In addition to A. mellifera, two other bee species, Hylaeus leptocephalus (Morawitz, 1871) and Lithurgus scabrosus (Smith, 1859), are of European origin and are now common and widespread [13]. Eight other species, Ceratina arizonensis (Cockerell, 1898), Xylocopa sonorina (Smith, 1874), Lasioglossum imbrex Gibbs, 2010, Lasioglossum impavidum (Sandhouse, 1924), Lasioglossum microlepoides (Ellis, 1914), Lasioglossum puteulanum (Gibbs, 2009), Megachile gentilis (Cresson, 1872), and Megachile policaris (Say, 1831), are adventive from the western United States [13,14]. One species, Hylaeus albonitens (Cockerell, 1905), is from Australia [13]. Finally, nine bee species, Ceratina smaragdula (Fabricius, 1787), Ceratina dentipes (Friese, 1914), Hylaeus strenuus (Cameron, 1897), Megachile chlorura (Cockerell, 1918), Megachile diligens (Smith, 1879), Megachile fullawayi (Cockerell, 1914), Megachile lanata (Fabricius, 1775), Megachile timberlakei (Cockerell, 1920), and Megachile umbripennis (Smith, 1853), are from Southeast Asia and the South Pacific [11,13,17,18]. There is evidence that two species, Megachile umbripennis and Megachile diligens, were introduced by Polynesians in pre-contact times [13]; however, data to support this hypothesis are lacking. In contrast, the remaining 18 invasive bee species (excluding A. mellifera) have been accidentally introduced following the arrival of non-kānaka maoli to the archipelago.
Non-native bee species have the potential to become invasive pests when they cause environmental, economic, or human harm [19]. Competition for floral resources and nesting sites, alteration of pollination networks, and introductions of pathogens are all possible consequences of alien bee introductions in Hawai’i [14]. The Megachilidae family has the largest invasive bee presence on the islands, with nine species. All Megachilidae species found in Hawai’i nest in a wide array of preformed cavities, but are especially prone to nest in wood or hollow plant stems [13,20]. Most Megachile are leaf-cutters; they line the nest cavity and construct individual cells from leaves or flower petals they have cut [13,20,21]. M. timberlakei, M. umbripennis, and L. scabrosus have been suggested to have negative effects on native populations by competing for nesting sites and floral resources and altering pollination networks [20,21,22]. Five species of bees in the Apidae family, including A. mellifera, C. arizonensis, C. smaragdula, and C. dentipes, have been suggested to compete for floral resources with Hylaeus bees in Hawai’i due to the similar proximity of habitat [12,23]. Furthermore, X. sonorina is a pollinator of invasive weeds [15,24,25].
Four species from the Halictidae family have been recorded in Hawai’i. In 2013, two species, L. microlepoides and L. imbrex, were frequently documented throughout the disturbed coastal and lowland areas around O‘ahu [23]. Given their dominance in disturbed habitats, both species are suggested to be significant pollinators of invasive weeds [23]. L. impavidum recorded in Hawai’i, has been suggested to compete for floral resources with native species because of its abundant occurrence in company with Hylaeus [12,13]. In California, its presumed origin, L. impavidum is primarily a lowland species, residing in elevations below 600 m [13]. However, its Hawai’i range has been documented at coastal localities to as high as 2700 m [13]. Finally, three invasive Hylaeus bees in the Colletidae family are present in Hawai’i. H. albonitens and H. strenuus are both suggested to compete for floral resources with native Hylaeus bees [12]. H. leptocephalus, which is present in Honolulu, O‘ahu, is thought to persist in low numbers only in urban sites, without spreading into native habitat [11]. As of the writing of this manuscript, no bees from the bee families Andrenidae, Stenotridae, or Melittidae have been documented in Hawai’i.
Research on how climate change could affect bees is limited for tropical island areas [1,26,27,28,29,30,31]. However, identifying the intersection between climate change and invasive species range expansion is critical for informing the management of imperiled and/or endemic biodiversity. For example, in Hawai’i there are more than 60 species of Hylaeus endemic to one or multiple islands [12,32]. Seven species of Hylaeus have been placed under the protection of the U.S. Endangered Species Act [33]. However, nearly half of the Hylaeus species endemic to Hawai’i are threatened due to habitat loss or alteration [12]. In fact, 10 species of Hylaeus may have gone extinct, as they have not been documented for a significant amount of time [12,23]. Endemic Hylaeus are especially vulnerable to the impacts of climate-induced biological invasions because they have fewer opportunities to adapt by altering their distribution [31,33]. The smaller land area of islands generally translates into very small populations and ranges for endemic insects such as Hylaeus [31,34,35,36]. Additionally, because of the limited area, Hylaeus endemic to Hawai’i are more vulnerable to invasive species through competition, predation, and disease spread [18,36,37]. These bees may be particularly susceptible to invasion because they rely on only a few plant species from which they gather pollen [38,39,40]. Low genetic variation, small colonizing populations, and reduced species richness may limit insular bee species in their ability to adapt [41,42,43,44]. Island species have evolved with few others and have developed survival strategies based on mutualism rather than defense mechanisms against predators and competitors [14,44]. As a result, harmful effects from climate-exacerbated biological invasions can encompass the entire habitat of endemic Hylaeus more readily than a continental species habitat.
The purpose of this study was to perform a comprehensive assessment of the specific regional niches of invasive bee species in Hawai’i and assess their potential expansion in the islands based on current and future climate scenarios. Little is known about potential habitat expansion of invasive bees due to climate change in Hawai’i [17]. An ability to accurately predict the impacts of climate change on species distribution is necessary to make informed decisions for biodiversity conservation. We chose to use an ensemble of species distribution models to accomplish our research aims. Our first objective in this study was to identify the bioclimatic niche for eight invasive bees in Hawai’i based on locality records from their native range and their invasive Hawai’i range. We predicted that the invaded Hawai’i niche will predict habitat suitability differently when compared to a model that projects habitat suitability based on the respective bee’s native niche. Our hypothesis was based on previous research that suggests invasive species have a high capacity to adapt to novel environments [45]. Our second objective in this study was to determine how future climate scenarios may influence the elevational distribution of invasive bees in Hawai’i. We predicted that invasive bees would expand into higher elevation habitat in future climate models. Our hypothesis was based on research that suggests species disperse to higher elevations as the temperature warms [46].

2. Methods

2.1. Species Occurrence Data

A total of 52,511 unique locality records of invasive bees in Hawai’i were accessed from museum databases (Global Biodiversity Information Facility (GBIF; http://www.gbif.org/, Accessed: 3 June 2019) (Table S1) and Discover Life (Asher and Pickering 2011) (Table S2). The locality records, also referred to as occurrence data, were collected at different times by a diversity of collectors. All occurrence records were aggregated into 1 km2 cells corresponding to the resolution of environmental variables. Identical presence locations were removed and only 1 unique presence location was retained. Records were verified by published works and occurrence records that could not be validated by the literature or additional resources were excluded from the study. We selected 8 out of the total 21 species to construct species distribution models in Hawai’i: A. mellifera, C. dentipes, C. smaragdula, L. impavidum, L. microlepoides, L. puteulanum, M. umbripennis, and X. sonorina (Figure 2). Sufficiently digitized and georeferenced location information for the remaining 13 species was not publicly available for analysis via GBIF or other online resources. Data on the bees invasive to Hawai’i, including the timing of their invasion, are provided in Table 1.

2.2. Environmental Data

In our study, we applied 19 bioclimatic variables (derived from temperature and precipitation measures), averaged between 1970 and 2000, with a spatial resolution of 30 arcsec (~1 × 1 km), from the WorldClim 2.1 database (http://www.worldclim.org/, Accessed: 3 June 2019) [47]. The current bioclimatic variables were computed from monthly values of minimum, average, and maximum temperature and monthly precipitation. The variables used in our analysis included the following: annual mean temperature (BIO 1), mean diurnal range (mean of monthly (max temp–min temp)) (BIO 2), isothermality (BIO2/BIO7) (×100) (BIO 3), temperature seasonality (standard deviation ×100) (BIO 4), max temperature of warmest month (BIO 5), min temperature of coldest month (BIO 6), temperature annual range (BIO 5–BIO 6) (BIO 7), mean temperature of wettest quarter (BIO 8), mean temperature of driest quarter (BIO 9), mean temperature of warmest quarter (BIO 10), mean temperature of coldest quarter (BIO 11), annual precipitation (BIO 12), precipitation of wettest month (BIO 13), precipitation of driest month (BIO 14), precipitation seasonality (coefficient of variation) (BIO 15), precipitation of wettest quarter (BIO 16), precipitation of driest quarter (BIO 17), precipitation of warmest quarter (BIO 18), and precipitation of coldest quarter (BIO 19). To reduce multicollinearity among the environmental variables, a principal component analysis was conducted to highlight the relationship between the target species occurrences and the specific environmental combinations within the archipelago. Variables were chosen based on orthogonal direction and overall environmental variation following strategies implemented in the BIOMOD2 package [48].
Distributions of the 8 invasive bee species were also modeled for future climatic conditions. We used projected bioclimatic variables for the period 2070 from representative greenhouse gas concentration pathway (RCP) scenario, RCP 8.5. RCP 8.5, derived from the Coupled Model Intercomparison Project Phase 5 (CMIP5), represents a radiative forcing of +8.5 W/m for the period 2000–2100, predicted to raise the average temperature 4.3 °C by 2100 [49]. Many research groups around the world have produced different global climate models (GCMs), which have been submitted to the CMIP5. GCMs can be used to forecast climate change because they capture the processes that respond to climate forcing. The bioclimatic variables used for modeling current distribution were used to predict future distribution of the invasive bees. We downloaded future bioclimate data of 12 GCMs: BCC-CSM1-1, CCSM4, GFDL-ESM2G, GISS-E2-R, HadGEM2-AO, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, MIROCESM, MIROC5, MRI-CGCM3, and NorESM1-M from the WorldClim 1.2 database. We created an ensemble of the 12 GCMs by taking average values and used the ensemble values as predictors. Our multi-model ensemble average accounts for inherent variability among the different future climate GCMs.

2.3. Species Distribution Modeling

Species distribution models (SDM) have been used to compare species’ regional ecological niches and forecast the range shifts of species under future climate change scenarios [1,50,51,52,53]. Species distribution modeling is an approach that predicts the distribution of a species across geographic space and time using the correlation between the geographic occurrence or abundance of a species and corresponding environmental conditions to predict the most suitable habitat [52,54]. Various methods have been used in SDMs, including regression, machine learning, classification, and maximum entropy [48,55]. Discrepancies between different techniques can be large and performance can vary significantly across different algorithms. Considering the variability between algorithms, we chose to use an ensemble modeling approach. Ensemble modeling of species distributions involves simulations across more than 1 set of initial conditions, model classes, model parameters, and boundary conditions by combining individual SDMs built through different modeling algorithms [55]. By using a wide range of approaches to test the models, ensemble modeling accounts for inter-model variability and uncertainties in predictions [56]. However, prediction uncertainty may also be dependent on GCM and RCP variation [55].
The analysis was conducted in R environment v 3.6.1 (R Core Team, 2019) using the BIOMOD2 package (Grenoble, France) [48]. The algorithms used to produce an ensemble model were as follows: 3 regression methods (GAM: general additive model; GLM: general linear model; and MARS: multivariate adaptive regression splines), 3 machine learning methods (ANN: artificial neural network; GBM: generalized boosting model; and RF: random forest), 2 classification methods (CTA: classification tree analysis; FDA: flexible discriminant analysis), and 1 maximum entropy approach (MAXENT) [57]. To identify the differences between bioclimatic niche of native range and Hawai’i range, 2 independent approaches were used: 1 approach using native range occurrences to predict favorable areas in Hawai’i and 1 approach using Hawai’i range occurrences to predict favorable areas in the Hawai’i. As these models required background data and the actual absence data were unavailable, we used 10,000 pseudo-absences randomly generated in the native range environmental space and 7000 pseudo-absences randomly generated in the Hawai’i environmental space. The models were calibrated by using 80% of the occurrence points (presence and pseudo-absence) as training data and evaluated by using the remaining 20% as testing data [46]. We repeated the process of pseudo-absence generation 3 times and repeated evaluation runs 4 times per species, resulting in a total of 108 models per species (9 models, 4 evaluation runs, and 3 pseudo-absence selection procedures) under each climate scenario (i.e., current climate (1970–2000) and RCP 8.5 (2070)).
The area under the curve (AUC) of receiver operating characteristics and true skills statistics (TSS) were used to measure model validation and predictive performance. The AUC value represents the predictive power of a model [58]. According to the AUC value, the model was graded as “poor” (if AUC = 0.6–0.7), “fair” (AUC = 0.7–0.8), “good” (AUC = 0.8–0.9), or “excellent” (AUC = 0.9–1.0) [58]. TSS measure ranges from −1 to +1 where +1 indicates a perfect agreement, and a TSS value below 0.4 indicates poor model discrimination [58]. From the 108 models per species, we built ensemble models using a weighted-mean approach in which weights are awarded for each model proportionally to their evaluation metrics scores. Only the models with greater than fair predictive accuracy (TSS > 0.5) to greater than good predictive accuracy (TSS > 0.8) were used to build an ensemble from the projection outputs [4,48].

2.4. Bioclimatic Niche Analysis

Binary maps (suitable and unsuitable) were produced using the optimal threshold that maximizes the TSS score as a cut-off value using the Biomod_RangeSize function, which then converted the projected occurrence probabilities during the cross-validation procedure. These binary maps were used to measure the loss, stability, and gain of predicted suitable areas following the predicted climate scenario demonstrated with an RCP 8.5 in 2070. From these binary maps, we measured the range size of the studied invasive species as represented by the number of climatically suitable pixels across Hawai’i for the designated time period (Table 2).
We calculated niche similarity between the native and the invasive Hawai’i range using the bioclimatic variables selected from the principal component analysis [48]. Native models were calibrated using the native range and projected onto the invasive Hawai’i environment. Alternatively, invasive Hawai’i models were calibrated using the occurrence records from the invasive Hawai’i range and projected onto the Hawai’i environment. We extracted the model values and tested for differences in habitat suitability between the native and invasive Hawai’i range with a Wilcoxon test across all 8 species.

2.5. Species Richness Analysis

To determine how projected climate change will impact invasive species richness, we conducted a species richness analysis to identify the regions potentially suitable for the maximum number of invasive bees under current and future climate. We combined binary maps of climatically suitable niches for all 8 species to generate species richness (cells with a higher value indicating high species diversity) and extent maps (cells occupied by at least a single species). We calculated changes in areas of both richness and extent of potentially suitable regions under current and future climate. Using a digital elevation model accessed from the University of Hawai’i at Manoa School of Ocean and Earth Science Technology (http://www.soest.hawaii.edu/coasts/data/hawaii/dem.html, Accessed: 3 June 2019) combined with the species richness maps, we determined how future climate scenarios may influence the elevational distribution of invasive bees in Hawai’i. Finally, we tested for the effect of species richness (i.e., 1–8 species) and climate timeframe (i.e., contemporary vs. future) on the elevation distribution of species richness with a non-parametric analogue to the two-way ANOVA, the Scheirer-Ray-Hare test.

3. Results

3.1. Species Occurrence Data

As publicly available and georeferenced data was lacking for the majority of the invasive bees to Hawai’i, we were able to pursue our research objectives with 8 of 21 invasive bee species: A. mellifera, C. dentipes, C. smaragdula, L. impavidum, L. microlepoides, L. puteulanum, M. umbripennis, and X. sonorina. These species represent bees that have georeferenced data in both their native and invasive Hawai’i geographies. We summarize the data in Table 1 and provide the citation for the data in Tables S1 and S2.

3.2. Species Distribution Modeling

Based on AUC and TSS metrics, the invasive Hawai’i niche, calibrated from the Hawai’i range occurrences, predicted higher habitat suitability than records sampled from their respective native geographies (Figure 3). Due to the native occurrence model poorly predicting actual species occurrences in Hawai’i, we decided to use the invasive Hawai’i SDMs to project habitat suitability in the future. The model performance was evaluated by the scores of two (AUC and TSS) performance metrics (Table 3). The average AUC values of the eight studied species ranged from 0.937 (A. mellifera) to 0.998 (L. puteulanum and M. umbripennis), indicating that the models have excellent predictive accuracy. Likewise, the average TSS value ranged from 0.731 (A. mellifera) to 0.991 (L. puteulanum), indicating good predictive accuracy. Moreover, we only used the models with highest predictive accuracy to build an ensemble from the projection outputs (Table 4).

3.3. Bioclimatic Niche Analysis

Following our PCA approach to selecting variables to construct SDMs, we identified a combination of four from a pool of seven bioclimatic variables to construct the final SDMs for each of the eight species: BIO 1, BIO 2, BIO 7, BIO 9, BIO 12, BIO 15, and BIO 19. Based on pairwise analyses of SDMs, our results indicated little overlap between the native niche and the invasive Hawai’i niche across seven species (A. mellifera: W = 8452.5, p = 2.65 × 10−5, C. dentipes: W = 166, p = 7.80 × 10−9, L. impavidum: W = 0, p = 0.00793, L. microlepoides: W = 28, p = 5.34 × 10−7, L. puteulanum: W = 0, p = 3.18 × 10−7, M. umbripennis: W = 0, p = 1.29 × 10−8, X. sonorina: W = 849, p = 0.0193) (Figure 3). Our analyses found that the only niche overlap between the native and Hawai’i ranges occurred across the C. smaragdula distribution (W = 974.5, p = 0.09402).

3.4. Species Richness Analysis

Based on the invasive Hawai’i SDMs, areas of potentially suitable niches for the studied species mostly occupy low elevation areas (Figure 4). Out of the eight species, five had potentially suitable areas that covered less than 5% of land area under current climate in Hawai’i (Table 2). Specifically, C. dentipes, L. imbrex, L. microlepoides, L. puteulanum, and M. umbripennis are predicted to have a restricted distribution under current climate conditions, whereas A. mellifera, C. smaragdula, and X. sonorina have a wider distribution under the current climatic conditions (Table 2 and Figure 5). However, under the current climate conditions, X. sonorina is predicted to have the largest suitable habitat, covering 15% of the land area in Hawai’i. The predicted suitable bioclimatic niche for L. puteulanum under current climate conditions covered only 1% of the land area in Hawai’i. Furthermore, under current climate conditions, the SDMs predicted most of the suitable habitat under 500 m with most areas occupied by a single species (Figure 4 and Figure 6).
Our analysis identified that “hotspot areas” of high invasive species richness are distributed primarily in low elevation sites (<500 m) under current and future climate conditions. Specifically, Scheirer-Ray-Hare tests found a significant difference in elevational distribution of species richness across current and future climate scenarios (Climate (current vs. 2070): H = 803.8, df = 1, p = 0; Species Richness (SR): H = 3828.2, df = 7, p = 0; Climate: SR: H = 57.5, df = 6, p = 1.45 × 10−10) (Figure 4), with highest median species richness of invasive species to be present in habitats at low elevation sites (<500 m). While our analyses demonstrated both gains and losses of suitable habitat when comparing contemporary and future climate models (Table 2), we predict that there will be a major geographic expansion of habitat suitability for invasive bees in Hawai’i under future climate scenarios (Figure 5 and Figure 6). Our SDMs predict that climatically suitable regions would increase on average by ~934.8% (±3.4% SE) for all species except L. impavidum (Table 2). For example, a maximum increase of 2425.7% habitat suitability for L. microlepoides is estimated. This habitat suitability increase is predicted to cover 25.9% of the terrestrial habitat across the archipelago (Table 2). However, the species estimated to experience the greatest range gain, in terms of geographic coverage, is A. mellifera. This economically significant species to Hawai’i is estimated to experience high habitat suitability across 40.3% of the terrestrial habitat across the archipelago (Table 2).

4. Discussion

The models developed in our study predict that the changing climate will create additional climatically suitable areas for invasive bees in Hawai’i over the next 50 years. This research provides baseline information to aid in effective management of invasive bees by showing the areas which have suitable niches for the invasive bees. Our study of eight of the 21 invasive bee species in Hawai’i found that SDMs based on georeferenced records from their invasive range were better at predicting habitat suitability in the invasive range than records found in their native range. Furthermore, while the eight invasive bees are predicted to expand in their geographic range over the next 50 years, the expansion will likely be limited to low elevation habitats <500 m, with little evidence for an expansion up to high elevation habitats (>500 m). Of all the species studied, L. microlepoides is predicted to encounter the greatest habitat suitability expansion, again primarily at low elevation habitats across the archipelago, whereas A. mellifera will likely be present in ~40% of the terrestrial land area of the major Hawaiian islands.
The results and approach of our study may be helpful for the prevention and early detection of invasive insects, namely bees, in determining suitable niches outside of their native range. The development of SDMs is a useful approach in determining how bees and other insects will expand outside of their native niche [1,45,59]. However, it is also clear that there are limitations to SDMs in predicting the invasive spread of bees, as evidenced by data collected during ground surveys of invasive A. manicatum in northeastern North America [60]. Specifically, Graham et al. [60] found no evidence for A. manicatum in the vast majority of areas predicted by Strange et al. [45] SDMs of A. manicatum—despite standardized surveys of their study area in northeast North America. However, the absence of an invasive species in a suitable area outside their native range does not necessarily suggest poor model performance, but may be an artifact of the species not yet dispersing into the new habitats, especially if they are recent invaders [61,62]. Furthermore, other limiting factors such as nesting biology, diet breadth, and phenology may also impact detection and colonization rates of invasive species, especially on oceanic islands.
The SDMs generated in this study predict that additional suitable areas for invasive bees are expected to emerge in the higher elevation zones of Hawai’i. However, our models predict that the expansion of invasive bees across elevation in Hawai’i is limited. The creation of climatically suitable regions for invasive bees in the high-elevation regions, which are already vulnerable to climate change, may have negative consequences for endemic biodiversity in the future. Therefore, biological invasions will add pressure and increase risks to the most vulnerable ecosystems in Hawai’i. Along with climate change, anthropogenic development is considered a major driver that promotes biological invasions in island ecosystems. For example, roads play an important role in the spread of alien species by facilitating dispersal corridors into the adjacent ecosystem [63]. This may promote the dispersal of invasive bees from the coast when climate change opens up suitable regions by reducing climatic barriers for them to invade high elevation regions [3]. Therefore, monitoring and management of invasive bees in Hawai’i should account for the vulnerability posed by climate change combined with anthropogenic activities. However, the suitable regions identified may not be occupied by invasive bees due to natural dispersal barriers created by lava flows, which are predominant physical features in Hawai’i.
Geographic shifts in species range involve multiple ecological processes such as dispersal, physiology, species interactions, population interactions, and evolution operating at multiple scales [64]. SDMs do not explicitly consider these inherent biological properties, which interact with ecological processes and ultimately cascade to species persistence [52]. Furthermore, there are other potential issues such as environmental variables used in the analysis, modeling algorithm, GCMs, collinearity, model complexity, model evaluation method, and threshold values to produce binary maps that can influence model outcomes. In addition, future land use (e.g., road building) change scenarios can also alter future species distributions [1,17]. Improvements of models based on natural history and ecological information and increased availability of specimen data are crucial issues for enhancing the predictive accuracy of the models [60].
Climate change has the potential to create more suitable regions for adventive bee species in Hawai’i. Climate change alters the distribution, composition, and phenology of native species while facilitating the dispersal of invasive species by removing current climatic barriers [3,7]. In our study, it is clear that habitat suitability for invasive bees will expand considerably in comparison to their current bioclimatic niche. However, this expansion of species richness is predicted to occur primarily below 500 m. Our results suggest that in Hawai’i, cooler temperatures associated with higher elevation sites may limit the invasion of invasive bees [4]. However, climate may not be the only limiting factor for invasive bees to spread to high elevation habitats. Thus, a better understanding of species traits, dispersal pathways, and the mechanism of the natural filters that prevent colonization of invasive species are necessary. The results of our research show a diverse set of hypothetical responses by invasive bees to climate change; therefore, species-specific prioritization exercises may be helpful to better manage and monitor specific invasive bee species.
Given the endangered status of endemic Hylaeus bees in Hawai’i, it is important to monitor and predict suitable habitats of invasive bees. By highlighting the similarities and differences between the native and invasive Hawai’i bioclimatic niche, results can inform stakeholders on the invasive potential of invasive bees [17,45]. Furthermore, the SDMs in our study provide evidence of current and future risks associated with invasive bee species in Hawai’i. For example, endemic Hylaeus bees, including species protected under the U.S. Endangered Species Act, that are distributed at low elevation habitats will likely continue to interact with invasive bees over the next 50 years. However, endemic Hylaeus bees at high elevation habitats will likely experience limited interactions with invasive bees. Given these predicted ecological differences, endemic Hylaeus bees will experience different disease and resource pressures.
While the current study used an ensemble model approach to estimate future climate models and species distribution models, the accuracy of models nevertheless relies on the accuracy of the data being used. The accuracy of GCMs is essential to predictability in SDMs, and the lack of bioclimatic agreement between different databases is an important factor [65]. Furthermore, including other bioclimatic variables in the initial model building process may produce more accurate models [66]. For example, a study that combined ENVIREM variables with Bioclim variables improved model performance in 13 out of 20 species [66]. For the current study of non-native bees in Hawai’i, model performance may be improved in future modeling exercises by including topographical variables such as slope, aspect, heat load index, and terrain ruggedness.
Many other approaches to SDMs exist that may improve or alter predictability. For example, a single source modeling approach often produces biased spatial predictions [65]. However, one study found that fine-tuned individual models may sometimes perform better than ensemble models [67]. Furthermore, Zhu et al. 2021 demonstrated that using a weighted mean approach to produce ensemble models may produce models overly influenced by the extreme values of individual predictions [67]. Sample size is also a large factor of influence on the accuracy of SDMs [68,69]. Algorithms can be fine-tuned to produce models with the highest predictability [70]. However, using an ensemble modeling approach to SDMs may overlook these fine-tuning capabilities. Many different scoring metrics exist for determining accuracy and predictability of SDMs such as the Brier score, Boyce index, or the Jaccard and Sorensen indices [4,71,72,73,74]. In this study, we chose to use AUC and TSS. While this method is widely used in ecology, some research suggests that TSS can be a misleading measure of model performance because of its dependence on prevalence [73].
In our study, we presented a geographic and bioclimatic assessment of eight invasive bees in Hawai’i. However, more research and surveys are needed to document the distribution and spread of the additional 13 bee species not included in our study. Our research has important implications for the management and monitoring of biological invasions of bees in Hawai’i. Our analysis also highlights the value of using SDMs to estimate species richness under future climate scenarios. As the climate changes, new habitats will emerge that may be suitable for adventive bee species. Climate change facilitates dispersal, introduction, and naturalization of adventive species as well as reduces the resilience of local ecosystems [75]. Thus, identifying emerging pests that may pose a threat to ecosystems in Hawai’i and the agricultural economy through SDMs is a useful tool in management and conservation.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/insects12050443/s1, Table S1: Bee specimen records evaluated from the Global Biodiversity Information Facility (http://gbif.org, Accessed: 3 June 2019) and their respective DOIs across species-specific data downloads. Table S2, Bee specimen records evaluated from the American Museum of Natural History to determine their capacity for informing SDMs across their (a) native and (b) non-native (invasive) Hawai’i distribution. Not all records could be used as they did not provide publicly available georeferenced data on the website.

Author Contributions

J.A.T. and J.B.K. designed the research. J.A.T. conducted analyses and data visualization. J.A.T. and J.B.K. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by a David H. Smith Postdoctoral Research Fellowship (2018) to J.B.K. and an NSF-REU site grant (#146130) to J.A.T. The APC was funded by a voucher from MDPI Insects.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data downloaded from GBIF are described in Table S1. Scripts and data used for all analyses are available at https://github.com/jesseat12/hawaii_bee (Accessed: 8 April 2021).

Acknowledgments

Mahalo to the University of Hawai’i High Performance Computing Cluster for assistance with remote computing and Renee Bellinger for onboarding the authors. We thank Jonathan Price for earlier guidance and suggestions on this project. J.A.T. was supported by an NSF-Research experience for undergraduates site grant awarded to Becky Ostertag (PI) and Noelani Puniwai (Co-PI) (#146130). The U.S. Department of Agriculture, Agricultural Research Service (USDA-ARS) is an equal opportunity/affirmative action employer and all agency services are available without discrimination. Mention of commercial products and organizations in this manuscript is solely to provide specific information. It does not constitute endorsement by USDA-ARS over other products and organizations not mentioned.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dew, R.M.; Silva, D.P.; Rehan, S.M. Range expansion of an already widespread bee under climate change. Glob. Ecol. Conserv. 2019, 17, e00584. [Google Scholar] [CrossRef]
  2. Potts, S.G.; Imperatriz-Fonseca, V.; Ngo, H.T.; Biesmeijer, J.C.; Breeze, T.D.; Dicks, L.V.; Garibaldi, L.A.; Hill, R.; Settele, J.; Vanbergen, A.J. The Assessment Report on Pollinators, Pollination and Food Production: Summary for Policymakers; Secretariat of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services: Bonn, Germany, 2016; p. 36. ISBN 9789280735680. [Google Scholar]
  3. Hellmann, J.J.; Byers, J.E.; Bierwagen, B.G.; Dukes, J.S. Five potential consequences of climate change for invasive species. Conserv. Biol. 2008, 22, 534–543. [Google Scholar] [CrossRef]
  4. Bellard, C.; Thuiller, W.; Leroy, B.; Genovesi, P.; Bakkenes, M.; Courchamp, F. Will climate change promote future invasions? Glob. Chang. Biol. 2013, 19, 3740–3748. [Google Scholar] [CrossRef]
  5. Thomas, C.D.; Cameron, A.; Green, R.E.; Bakkenes, M.; Beaumont, L.J.; Collingham, Y.C.; Erasmus, B.F.N.; De Siqueira, M.F.; Grainger, A.; Hannah, L.; et al. Extinction risk from climate change. Nature 2004, 427, 145–148. [Google Scholar] [CrossRef]
  6. Pecl, G.T.; Araújo, M.B.; Bell, J.D.; Blanchard, J.; Bonebrake, T.C.; Chen, I.-C.; Clark, T.D.; Colwell, R.K.; Danielsen, F.; Evengård, B.; et al. Biodiversity redistribution under climate change: Impacts on ecosystems and human well-being. Science 2017, 355. [Google Scholar] [CrossRef]
  7. Walther, G.-R.; Post, E.; Convey, P.; Menzel, A.; Parmesan, C.; Beebee, T.J.C.; Fromentin, J.-M.; Hoegh-Guldberg, O.; Bairlein, F. Ecological responses to recent climate change. Nature 2002, 416, 389–395. [Google Scholar] [CrossRef]
  8. Root, T.L.; Price, J.T.; Hall, K.R.; Schneider, S.H.; Rosenzweig, C.; Pounds, J.A. Fingerprints of global warming on wild animals and plants. Nature 2003, 421, 57–60. [Google Scholar] [CrossRef]
  9. Bartomeus, I.; Ascher, J.S.; Wagner, D.; Danforth, B.N.; Colla, S.; Kornbluth, S.; Winfree, R. Climate-associated phenological advances in bee pollinators and bee-pollinated plants. Proc. Natl. Acad. Sci. USA 2011, 108, 20645–20649. [Google Scholar] [CrossRef] [Green Version]
  10. Rubidge, E.M.; Patton, J.L.; Lim, M.; Cole Burton, A.; Brashares, J.S.; Moritz, C. Climate-induced range contraction drives genetic erosion in an alpine mammal. Nat. Clim. Chang. 2012, 2, 285–288. [Google Scholar] [CrossRef]
  11. Magnacca, K.N.; Gibbs, J.; Droege, S. Notes on alien and native bees (Hymenoptera: Apoidea) from the Hawaiian Islands. Bish. Mus. Occas. Pap. 2013, 114, 61–65. [Google Scholar]
  12. Magnacca, K.N. Conservation status of the endemic bees of Hawai’i, Hylaeus (Nesoprosopis) (Hymenoptera: Colletidae). Pac. Sci. 2007, 61, 173–190. [Google Scholar] [CrossRef] [Green Version]
  13. Snelling, R.R. Bees of the Hawaiian Islands, exclusive of Hylaeus (Nesoprosopis) (Hymenoptera: Apoidea). J. Kans. Entomol. Soc. 2003, 76, 342–356. [Google Scholar]
  14. Reaser, J.K.; Meyerson, L.A.; Cronk, Q.; de Poorter, M.; Eldrege, L.G.; Green, E.; Kairo, M.; Latasi, P.; Mack, R.N.; Mauremootoo, J.; et al. Ecological and socioeconomic impacts of invasive alien species in island ecosystems. Environ. Conserv. 2007, 34, 98–111. [Google Scholar] [CrossRef] [Green Version]
  15. Barrows, E.M. Robbing of Exotic Plants by Introduced Carpenter and Honey Bees in Hawaii, with Comparative Notes. Biotropica 1980, 12, 23–29. [Google Scholar] [CrossRef]
  16. Roddy, K.M.; Arita-Tsutsumi, L. A history of honey bees in the Hawaiian Islands. J. Hawaii Pac. Agric. 1997, 8, 59–70. [Google Scholar]
  17. Shell, W.A.; Rehan, S.M. Invasive range expansion of the small carpenter bee, Ceratina dentipes (Hymenoptera: Apidae) into Hawai’i with implications for native endangered species displacement. Biol. Invasions 2019, 21, 1155–1166. [Google Scholar] [CrossRef]
  18. Shell, W.A.; Rehan, S.M. Range Expansion of the Small Carpenter Bee Ceratina smaragdula across the Hawaiian Archipelago with Potential Ecological Implications for Native Pollinator Systems. Pac. Sci. 2017, 71, 1–15. [Google Scholar] [CrossRef]
  19. Stout, J.C.; Morales, C.L. Ecological impacts of invasive alien species on bees. Apidologie 2009, 40, 388–409. [Google Scholar] [CrossRef] [Green Version]
  20. Poulsen, N.R.; Rasmussen, C. Island bees: Do wood nesting bees have better island dispersal abilities? Apidologie 2020, 51, 1006–1017. [Google Scholar] [CrossRef]
  21. Groom, S.V.C.; Hayes, S.E.; Ngo, H.T.; Stevens, M.I.; Schwarz, M.P. Recipe for disruption: Multiple recent arrivals of megachilid bees in Pacific archipelagos. J. Insect Conserv. 2014, 18, 613–622. [Google Scholar] [CrossRef]
  22. Rasmussen, C.; Carríon, A.L.; Castro-Urgal, R.; Chamorro, S.; Gonzalez, V.H.; Griswold, T.L.; Herrera, H.W.; McMullen, C.K.; Olesen, J.M.; Traveset, A. Megachile timberlakei Cockerell (Hymenoptera: Megachilidae): Yet another adventive bee species to the Galápagos Archipelago. Pan-Pac. Entomol. 2012, 88, 98–102. [Google Scholar] [CrossRef]
  23. Magnacca, K.N.; King, C. Assessing the presence and distribution of 23 Hawaiian yellow-faced bee species on lands adjacent to military installations on O’ahu and Hawai’i Island. Tech. Rep. 2013, 185, 26. [Google Scholar]
  24. Groom, S.V.C.; Stevens, M.I.; Ramage, T.; Schwarz, M.P. Origins and implications of apid bees (Hymentopera: Apidae) in French Polynesia. Entomol. Sci. 2017, 20, 65–75. [Google Scholar] [CrossRef]
  25. Manson, D.C.M. The occurrence of the carpenter bee, Xylocopa sonorina in New Zealand. Weta 1988, 11, 4–5. [Google Scholar]
  26. Imbach, P.; Fung, E.; Hannah, L.; Navarro-Racines, C.E.; Roubik, D.W.; Ricketts, T.H.; Harvey, C.A.; Donatti, C.I.; Läderach, P.; Locatelli, B.; et al. Coupling of pollination services and coffee suitability under climate change. Proc. Natl. Acad. Sci. USA 2017, 114, 10438–10442. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Le Conte, Y.; Navajas, M. Climate change: Impact on honey bee populations and diseases. Rev. Sci. Tech. 2008, 27, 485–497, 499–510. [Google Scholar] [CrossRef]
  28. Giannini, T.C.; Maia-Silva, C.; Acosta, A.L.; Jaffé, R.; Carvalho, A.T.; Martins, C.F.; Zanella, F.C.V.; Carvalho, C.A.L.; Hrncir, M.; Saraiva, A.M.; et al. Protecting a managed bee pollinator against climate change: Strategies for an area with extreme climatic conditions and socioeconomic vulnerability. Apidologie 2017, 48, 784–794. [Google Scholar] [CrossRef] [Green Version]
  29. Giannini, T.C.; Acosta, A.L.; Garófalo, C.A.; Saraiva, A.M.; Alves-dos-Santos, I.; Imperatriz-Fonseca, V.L. Pollination services at risk: Bee habitats will decrease owing to climate change in Brazil. Ecol. Model. 2012, 244, 127–131. [Google Scholar] [CrossRef]
  30. Giannini, T.C.; Costa, W.F.; Borges, R.C.; Miranda, L.; da Costa, C.P.W.; Saraiva, A.M.; Imperatriz Fonseca, V.L. Climate change in the Eastern Amazon: Crop-pollinator and occurrence-restricted bees are potentially more affected. Reg. Environ. Chang. 2020, 20, 9. [Google Scholar] [CrossRef] [Green Version]
  31. Taylor, S.; Kumar, L. Global Climate Change Impacts on Pacific Islands Terrestrial Biodiversity: A Review. Trop. Conserv. Sci. 2016, 9, 203–223. [Google Scholar] [CrossRef] [Green Version]
  32. Magnacca, K.N. Two new species of Hylaeus (Nesoprosopis) (Hymenoptera: Colletidae) from O’ahu, Hawai’i. Zootaxa 2011, 3065, 60–65. [Google Scholar] [CrossRef] [Green Version]
  33. U.S. Fish and Wildlife Service. Endangered and threatened wildlife and plants: Endangered status for 49 species from the Hawaiian islands, Proposed Rule. 50 CFR Part 17. Fed. Regist. 2015, 81, 67786–67860. [Google Scholar]
  34. Altaba, C.R. Extinction Resilience of Island Species: An Amphibian Case and a Predictive Model. Diversity 2014, 6, 43–71. [Google Scholar] [CrossRef]
  35. Harter, D.E.V.; Irl, S.D.H.; Seo, B.; Steinbauer, M.J.; Gillespie, R.; Triantis, K.A.; Fernández-Palacios, J.-M.; Beierkuhnlein, C. Impacts of global climate change on the floras of oceanic islands-Projections, implications and current knowledge. Perspect. Plant Ecol. Evol. Syst. 2015, 17, 160–183. [Google Scholar] [CrossRef]
  36. Ing, K.; Mogren, C.L. Evidence of Competition between Honey Bees and Hylaeus anthracinus (Hymenoptera: Colletidae), an Endangered Hawaiian Yellow-Faced Bee. Pac. Sci. 2020, 74, 75–85. [Google Scholar] [CrossRef]
  37. Plentovich, S.; Graham, J.R.; Haines, W.P.; King, C.B.A. Invasive ants reduce nesting success of an endangered Hawaiian yellow-faced bee, Hylaeus anthracinus. NeoBiota 2021, 64, 137–154. [Google Scholar] [CrossRef]
  38. Koch, J.B.; Sahli, H.F. Patterns of flower visitation across elevation and successional gradients in Hawai’i. Pac. Sci. 2013, 67, 253–266. [Google Scholar] [CrossRef]
  39. Sahli, H.F.; Krushelnycky, P.D.; Drake, D.R.; Taylor, A.D. Patterns of Floral Visitation to Native Hawaiian Plants in Presence and Absence of Invasive Argentine Ants1. Pac. Sci. 2016, 70, 309–322. [Google Scholar] [CrossRef]
  40. Shay, K.; Drake, D.R.; Taylor, A.D.; Sahli, H.F. Alien insects dominate the plant-pollinator network of a Hawaiian coastal ecosystem. Pacifica 2016, 70, 409–429. [Google Scholar] [CrossRef]
  41. Corlett, R.T. Invasive aliens on tropical East Asian islands. Biodivers. Conserv. 2010, 19, 411–423. [Google Scholar] [CrossRef]
  42. Jupiter, S.; Mangubhai, S.; Kingsford, R.T. Conservation of Biodiversity in the Pacific Islands of Oceania: Challenges and Opportunities. Pac. Conserv. Biol. 2014, 20, 206–220. [Google Scholar] [CrossRef]
  43. Parmesan, C.; Yohe, G. A globally coherent fingerprint of climate change impacts across natural systems. Nature 2003, 421, 37–42. [Google Scholar] [CrossRef]
  44. Cronk, Q.C.B. Islands: Stability, diversity, conservation. Biodivers. Conserv. 1997, 6, 477–493. [Google Scholar] [CrossRef]
  45. Strange, J.P.; Koch, J.B.; Gonzalez, V.H.; Nemelka, L. Global invasion by Anthidium manicatum (Linnaeus) (Hymenoptera: Megachilidae): Assessing potential distribution in North America and beyond. Biologicals 2011, 13, 2115–2133. [Google Scholar] [CrossRef]
  46. Freeman, B.G.; Scholer, M.N.; Ruiz-Gutierrez, V.; Fitzpatrick, J.W. Climate change causes upslope shifts and mountaintop extirpations in a tropical bird community. Proc. Natl. Acad. Sci. USA 2018, 115, 11982–11987. [Google Scholar] [CrossRef] [PubMed]
  47. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas: NEW CLIMATE SURFACES FOR GLOBAL LAND AREAS. Int. J. Clim. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  48. Thuiller, W.; Lafourcade, B.; Engler, R.; Araújo, M.B. BIOMOD—A platform for ensemble forecasting of species distributions. Ecography 2009, 32, 369–373. [Google Scholar] [CrossRef]
  49. Collins, M.; Knutti, R.; Arblaster, J.; Dufresne, J.-L.; Fichefet, T.; Friedlingstein, P.; Gao, X.; Gutowski, W.J.; Johns, T.; Krinner, G.; et al. Long-term climate change: Projections, commitments and irreversibility. In Climate Change 2013—The Physical Science Basis: Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2013; pp. 1029–1136. [Google Scholar]
  50. Guisan, A.; Thuiller, W. Predicting species distribution: Offering more than simple habitat models. Ecol. Lett. 2005, 8, 993–1009. [Google Scholar] [CrossRef]
  51. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 2011, 17, 43–57. [Google Scholar] [CrossRef]
  52. Elith, J.; Kearney, M.; Phillips, S. The art of modelling range-shifting species: The art of modelling range-shifting species. Methods Ecol. Evol. 2010, 1, 330–342. [Google Scholar] [CrossRef]
  53. Elith, J.; Leathwick, J.R. Species Distribution Models: Ecological Explanation and Prediction across Space and Time. Annu. Rev. Ecol. Evol. Syst. 2009, 40, 677–697. [Google Scholar] [CrossRef]
  54. Eckert, S.; Hamad, A.; Kilawe, C.J.; Linders, T.E.W.; Ng, W.; Mbaabu, P.R.; Shiferaw, H.; Witt, A.; Schaffner, U. Niche change analysis as a tool to inform management of two invasive species in Eastern Africa. Ecosphere 2020, 11. [Google Scholar] [CrossRef] [Green Version]
  55. Shrestha, U.B.; Shrestha, B.B. Climate change amplifies plant invasion hotspots in Nepal. Divers. Distrib. 2019, 10, e0134665. [Google Scholar] [CrossRef] [Green Version]
  56. Araújo, M.B.; New, M. Ensemble forecasting of species distributions. Trends Ecol. Evol. 2007, 22, 42–47. [Google Scholar] [CrossRef] [PubMed]
  57. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 2006, 190, 231–259. [Google Scholar] [CrossRef] [Green Version]
  58. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS): Assessing the accuracy of distribution models. J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  59. Koch, J.; Looney, C.; Hopkins, B.; Lichtenberg, E.M. Projected climate change will reduce habitat suitability for bumble bees in the Pacific Northwest. BioRxiv 2019. [Google Scholar] [CrossRef]
  60. Graham, K.K.; Eaton, K.; Obrien, I.; Starks, P.T. Anthidium manicatum, an invasive bee, excludes a native bumble bee, Bombus impatiens, from floral resources. Biol. Invasions 2018. [Google Scholar] [CrossRef]
  61. Zhu, G.; Gutierrez Illan, J.; Looney, C.; Crowder, D.W. Assessing the ecological niche and invasion potential of the Asian giant hornet. Proc. Natl. Acad. Sci. USA 2020, 117, 24646–24648. [Google Scholar] [CrossRef]
  62. Moo-Llanes, D.A. Inferring Distributional Shifts of Asian Giant Hornet Vespa mandarinia Smith in Climate Change Scenarios. Neotrop. Entomol. 2021. [Google Scholar] [CrossRef]
  63. McDougall, K.L.; Lembrechts, J.; Rew, L.J.; Haider, S.; Cavieres, L.A.; Kueffer, C.; Milbau, A.; Naylor, B.J.; Nuñez, M.A.; Pauchard, A.; et al. Running off the road: Roadside non-native plants invading mountain vegetation. Biol. Invasions 2018, 20, 3461–3473. [Google Scholar] [CrossRef] [Green Version]
  64. Urban, M.C.; Bocedi, G.; Hendry, A.P.; Mihoub, J.-B.; Pe’er, G.; Singer, A.; Bridle, J.R.; Crozier, L.G.; De Meester, L.; Godsoe, W.; et al. Improving the forecast for biodiversity under climate change. Science 2016, 353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Morales-Barbero, J.; Vega-Álvarez, J. Input matters matter: Bioclimatic consistency to map more reliable species distribution models. Methods Ecol. Evol. 2019, 10, 212–224. [Google Scholar] [CrossRef]
  66. Title, P.O.; Bemmels, J.B. ENVIREM: An expanded set of bioclimatic and topographic variables increases flexibility and improves performance of ecological niche modeling. Ecography 2018, 41, 291–307. [Google Scholar] [CrossRef] [Green Version]
  67. Zhu, G.; Fan, J.; Peterson, A.T. Cautions in weighting individual ecological niche models in ensemble forecasting. Ecol. Model. 2021, 448, 109502. [Google Scholar] [CrossRef]
  68. Barbet-Massin, M.; Rome, Q.; Villemant, C.; Courchamp, F. Can species distribution models really predict the expansion of invasive species? PLoS ONE 2018, 13, e0193085. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Liu, C.; Newell, G.; White, M. The effect of sample size on the accuracy of species distribution models: Considering both presences and pseudo-absences or background sites. Ecography 2019, 42, 535–548. [Google Scholar] [CrossRef] [Green Version]
  70. Muscarella, R.; Galante, P.J.; Soley-Guardia, M.; Boria, R.A.; Kass, J.M.; Uriarte, M.; Anderson, R.P. ENMeval: An R package for conducting spatially independent evaluations and estimating optimal model complexity for MAXENT ecological niche models. Methods Ecol. Evol. 2014, 5, 1198–1205. [Google Scholar] [CrossRef]
  71. Bradley, A.A.; Schwartz, S.S.; Hashino, T. Sampling Uncertainty and Confidence Intervals for the Brier Score and Brier Skill Score. Weather Forecast. 2008, 23, 992–1006. [Google Scholar] [CrossRef]
  72. Hirzel, A.H.; Le Lay, G.; Helfer, V.; Randin, C.; Guisan, A. Evaluating the ability of habitat suitability models to predict species presences. Ecol. Model. 2006, 199, 142–152. [Google Scholar] [CrossRef]
  73. Leroy, B.; Delsol, R.; Hugueny, B.; Meynard, C.N.; Barhoumi, C.; Barbet-Massin, M.; Bellard, C. Without quality presence-absence data, discrimination metrics such as TSS can be misleading measures of model performance. J. Biogeogr. 2018, 45, 1994–2002. [Google Scholar] [CrossRef]
  74. Roberts, D.R.; Bahn, V.; Ciuti, S.; Boyce, M.S.; Elith, J.; Guillera-Arroita, G.; Hauenstein, S.; Lahoz-Monfort, J.J.; Schröder, B.; Thuiller, W.; et al. Cross-validation strategies for data with temporal, spatial, hierarchical, or phylogenetic structure. Ecography 2017, 40, 913–929. [Google Scholar] [CrossRef]
  75. Walther, G.-R.; Roques, A.; Hulme, P.E.; Sykes, M.T.; Pysek, P.; Kühn, I.; Zobel, M.; Bacher, S.; Botta-Dukát, Z.; Bugmann, H.; et al. Alien species in a warmer world: Risks and opportunities. Methods Ecol. Evol. 2009, 24, 686–693. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Insects 12 00443 g001 550
Figure 1. Documentation of invasive bees to Hawai’i over time. See Table 1 for the year an invasive bee species was documented in Hawai’i.
Figure 1. Documentation of invasive bees to Hawai’i over time. See Table 1 for the year an invasive bee species was documented in Hawai’i.
Insects 12 00443 g001
Insects 12 00443 g002 550
Figure 2. Geographic distribution of 8 invasive bees to Hawai’i.
Figure 2. Geographic distribution of 8 invasive bees to Hawai’i.
Insects 12 00443 g002
Insects 12 00443 g003 550
Figure 3. Estimates of habitat suitability (HS) using native and invasive to Hawai’i georeferenced records in species distribution models. Species with asterisk represents a significant difference in HS. Species with NS represents no significant difference in HS.
Figure 3. Estimates of habitat suitability (HS) using native and invasive to Hawai’i georeferenced records in species distribution models. Species with asterisk represents a significant difference in HS. Species with NS represents no significant difference in HS.
Insects 12 00443 g003
Insects 12 00443 g004 550
Figure 4. Elevation distribution of species richness across current and future (i.e., 8.5) climate scenarios for invasive to Hawai’i models. Black line in each boxplot represents the median elevation distribution.
Figure 4. Elevation distribution of species richness across current and future (i.e., 8.5) climate scenarios for invasive to Hawai’i models. Black line in each boxplot represents the median elevation distribution.
Insects 12 00443 g004
Insects 12 00443 g005a 550Insects 12 00443 g005b 550
Figure 5. Habitat extent for all 8 invasive bee species across current (top) and future (bottom) climate scenarios for invasive to Hawai’i models.
Figure 5. Habitat extent for all 8 invasive bee species across current (top) and future (bottom) climate scenarios for invasive to Hawai’i models.
Insects 12 00443 g005aInsects 12 00443 g005b
Insects 12 00443 g006a 550Insects 12 00443 g006b 550
Figure 6. Species richness (i.e., 1–8) across current (top) and future (bottom) climate scenarios for invasive to Hawai’i models.
Figure 6. Species richness (i.e., 1–8) across current (top) and future (bottom) climate scenarios for invasive to Hawai’i models.
Insects 12 00443 g006aInsects 12 00443 g006b
Table
Table 1. Occurrence data download summary for 21 bees (Hymenoptera: Anthophila) that are invasive to Hawai’i from the Global Biodiversity Information Facility (GBIF) webpage (http://gbif.org, Accessed: 3 June 2019).
Table 1. Occurrence data download summary for 21 bees (Hymenoptera: Anthophila) that are invasive to Hawai’i from the Global Biodiversity Information Facility (GBIF) webpage (http://gbif.org, Accessed: 3 June 2019).
FamilyGenusSpeciesNo. of Records on GBIFNo. GeoRef Records on GBIFFinal No. GeoRef RecordsYear Documented in Hawai’iGeneral Native Range
ApidaeApismellifera150,293134,92650,6401857Europe (Snelling 2003)
ApidaeCeratinaarizonensis8368281011950Southwestern United States (Daly 1973)
ApidaeCeratinasmaragdula4093181071998Southeast Asia (Hirashima, 1969, Snelling 2003)
ApidaeCeratinadentipes183145701996Southeast Asia (Snelling 2003)
ApidaeXylocopasonorina1551381011874North America (Snelling 2003)
ColletidaeHylaeusleptocephalus5894902471958Europe (Snelling 2003)
ColletidaeHylaeusalbonitens924172591995Australia (Magnacca & King 2013)
ColletidaeHylaeusstrenuus1122007India (Magnacca 2011)
HalictidaeLasioglossumimbrex00492005Western North America (Gibbs 2010)
HalictidaeLasioglossumimpavidum337326281994Coastal California (Snelling 2003)
HalictidaeLasioglossummicrolepoides13,25813,2421802010Western North America and Northern Mexico (Magnacca & King 2013)
HalictidaeLasioglossumputeulanum175217491882012Eastern North America (Gibbs, 2011); USGS (https://www.usgs.gov/media/images/lasioglossum-nr-puteulanum-male-side, Accessed: 3 June 2019)
MegachilidaeLithurgusscabrosus3618141907Europe (Snelling 2003)
MegachilidaeMegachilechlorura3121988Philippines (Snelling, 2003); Southeast Asia (Rasmussen 2012)
MegachilidaeMegachilediligens50151879South Pacific according to Snelling, 2003), Southeast Asia according to Rasmussen 2012
MegachilidaeMegachilefullawayi2117121921Guam according to Cockerell, but likely brought from Asia. (Snelling 2003)
MegachilidaeMegachilegentilis7316682301899Northwestern USA, species know from southern BC (Snelling 2003); as M. palmarum Perkins
MegachilidaeMegachilelanata345302862012Southeast Asia (Gonsalez et al. 2019), India (Magnacca et al. 2013)
MegachilidaeMegachilepolicaris11148913822018Georgia and Florida, west to California and Mexico
MegachilidaeMegachiletimberlakei10881904First documented in Hawai’i, probably South Pacific region. (Snelling 2003)
MegachilidaeMegachileumbripennis177148391898Northern India and China (Timberlake, 1921)
Table
Table 2. Predicted distribution of habitat suitability across the 8 major islands of Hawai’i under contemporary and future (2070, RCP 8.5) global climate models. Estimates presented in this table approximate 1 pixel to 1 km2 (30 arc sec).
Table 2. Predicted distribution of habitat suitability across the 8 major islands of Hawai’i under contemporary and future (2070, RCP 8.5) global climate models. Estimates presented in this table approximate 1 pixel to 1 km2 (30 arc sec).
SpeciesCurrent Range (km2)Range Loss (km2) (2070, RCP 8.5)Range Gain (km2)
(2070, RCP 8.5)		% Loss
(2070, RCP 8.5)		% Gain (km2)
(2070, RCP 8.5)
A. mellifera2562 (12%) *1685538 (40.3%)6.6%216.2%
C. dentipes889 (4.4%)213987 (24.3%)2.4%448.5%
C. smaragdula748 (3.7%)03600 (21.7%)0481.3%
L. impavidum1217 (6.0%)171202 (7.0%)14.1%16.6%
L. microlepoides364 (1.8%)04827 (25.9%)02425.7%
L. puteulanum213 (1.0%)04202 (22.0%)01972.8%
M. umbripennis377 (1.8%)03308 (18.3%)01769.0%
X. sonorina2958 (14.7%)134390 (36.6%)0.4%148.4%
* Percentages in parenthesis estimate percent of suitable area for 8 invasive bees relative to total estimated area of Hawai’i in the contemporary time (1970–2000) and projected into 2070 (20,061 pixels total at 30 arc sec). Percentages without parenthesis estimate percent relative to current range size.
Table
Table 3. True skills statistic (TSS) and area under the curve (AUC) test statistics for 8 SDMs predicted the distribution of 8 invasive bees to Hawai’i.
Table 3. True skills statistic (TSS) and area under the curve (AUC) test statistics for 8 SDMs predicted the distribution of 8 invasive bees to Hawai’i.
SpeciesInvasive Model TSSInvasive Model AUCNative Model TSSNative Model AUC
A. mellifera0.7310.9370.9210.995
C. dentipes0.9580.9930.8150.956
C. smaragdula0.9170.9870.7590.939
L. impavidum0.940.9750.8070.936
L. microlepoides0.9830.9970.6020.869
L. puteulanum0.9910.9980.7960.954
M. umbripennis0.9830.9980.9070.976
X. sonorina0.7530.9430.630.89
Table
Table 4. True skills statistic (TSS) thresholds to produce ensemble models for 8 invasive bees to Hawai’i in their native range and invasive Hawai’i range.
Table 4. True skills statistic (TSS) thresholds to produce ensemble models for 8 invasive bees to Hawai’i in their native range and invasive Hawai’i range.
SpeciesNative ModelInvasive Model
A. mellifera0.80.5
C. dentipes0.50.8
C. smaragdula0.80.8
L. impavidum0.50.5
L. microlepoides0.50.8
L. puteulanum0.50.8
M. umbripennis0.50.5
X. sonorina0.50.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tabor, J.A.; Koch, J.B. Ensemble Models Predict Invasive Bee Habitat Suitability Will Expand under Future Climate Scenarios in Hawai’i. Insects 2021, 12, 443. https://doi.org/10.3390/insects12050443

AMA Style

Tabor JA, Koch JB. Ensemble Models Predict Invasive Bee Habitat Suitability Will Expand under Future Climate Scenarios in Hawai’i. Insects. 2021; 12(5):443. https://doi.org/10.3390/insects12050443

Chicago/Turabian Style

Tabor, Jesse A., and Jonathan B. Koch. 2021. "Ensemble Models Predict Invasive Bee Habitat Suitability Will Expand under Future Climate Scenarios in Hawai’i" Insects 12, no. 5: 443. https://doi.org/10.3390/insects12050443

APA Style

Tabor, J. A., & Koch, J. B. (2021). Ensemble Models Predict Invasive Bee Habitat Suitability Will Expand under Future Climate Scenarios in Hawai’i. Insects, 12(5), 443. https://doi.org/10.3390/insects12050443

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop