Research Article |
Corresponding author: Stephanie Stiegel ( stiegel@uni-hildesheim.de ) Academic editor: Christopher K. Starr
© 2023 Laura Rojas-Arias, Daniel Gómez-Morales, Stephanie Stiegel, Rodulfo Ospina-Torres.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Rojas-Arias L, Gómez-Morales D, Stiegel S, Ospina-Torres R (2023) Niche modeling of bumble bee species (Hymenoptera, Apidae, Bombus) in Colombia reveals highly fragmented potential distribution for some species. Journal of Hymenoptera Research 95: 231-244. https://doi.org/10.3897/jhr.95.87752
|
Insect population decline has been reported worldwide, including those of pollinators important for ecosystem services. Therefore, conservation actions which rely on available rigorous species distribution data are necessary to protect biodiversity. Niche modeling is an appropriate approach to distribution maps, but when it comes to bumble bees, few studies have been performed in South America. We modeled ecological niches of nine Colombian Bombus species with MAXENT 3.4 software using bioclimatic variables available from WorldClim. This resulted in maps for each species that show the potential distribution area at the present time. Modeled species maps accurately represent potential niches according to the description of bioclimatic conditions in the species’ habitat. We grouped the species into three clusters based on our results, as well as on distributional information from literature on the topic: High Mountain, Mid- Mountain and inter-Andean, and the Amazon and Eastern Plains Basin. Niche modeling depicted bumble bee species’ distribution in Colombia, the results of which can serve as a useful tool for conservation policies in the country.
biogeography, distribution, maxent, native bees, pollinators
Bumble bee species (Bombus spp.) are one of the widest studied taxa of bees in the world and Colombia regarding their ecological traits, the genetic composition of populations, pathogens, and their role in pollination services (
However, a drastic decline in insects has been documented worldwide (Hallmann et al. 2017; Seibold et al. 2019; Klink et al. 2020). A massive decline in the species distribution of most bumble bees has been reported, especially in developed regions like Europe and North America (
Biogeographic information about species is vital for conservation planning. The study of species distribution patterns has two complementary approaches: historical biogeography, which elucidates the causal mechanisms of current distribution, and ecological biogeography, which evaluates the ecogeographic factors that are currently shaping the distribution of species, looking for patterns in characteristics that are required for a species’ long-term survival. It looks at abiotic characteristics like humidity, temperature, or salinity and ecological characteristics like interactions with other organisms or genetical features. Research in this last approach is usually performed at a local scale a short time frame (
Niche modeling is one of the methods used in ecological biogeography. It represents the fundamental niche without considering the realized niche, which requires detailed field research and confirmed species samples. Niche modeling aims to find the potential distribution area of a species (
Thus, niche modeling can provide potential distribution maps for bumble bees in Colombia, updating and complementing the previous available maps, made by
The geographical extent used for each species in the modeling process covers the entire continental Colombian territory, from the north (13°23.73'N) to the south (4°13.75'N) and from the west (81°44.13'W) to the east (66°50.63'W). No buffer was used.
There are nine species from the genus Bombus reported for Colombia: B. pauloensis Friese, 1913 (formerly B. atratus Franklin, 1913), B. excellens Smith, 1879, B. funebris Smith, 1854, B. hortulanus Friese, 1904, B. melaleucus Handlirschi, 1888, B. pullatus Franklin, 1913, B. robustus Smith, 1854, B. rubicundus Smith, 1854, and B. transversalis Oliver, 1789. Across the paramo ecosystem, four species can be found with different altitudinal distributions: B. funebris between 2500 and 4750 m, B. hortulanus between 2100 and 3600 m, B. robustus between 2100 and 3800 m, and B. rubicundus between 2500 and 3900 m (
We obtained occurrence points from the Wild Bee Research Lab of the Universidad Nacional de Colombia, Bogotá (
We used 19 environmental data layers for modeling (Table
Bioclimatic variables included in modeling were obtained from WorldClim.
Temperature | Precipitation |
---|---|
1 Annual Mean Temperature | 12 Annual Precipitation |
2 Mean Diurnal Range (Mean of monthly (max. temp.-min. temp.)) | 13 Precipitation of Wettest Month |
3 Isothermality (BIO2/BIO7) (×100) | 14 Precipitation of Driest Month |
4 Temperature Seasonality (standard deviation ×100) | 15 Precipitation Seasonality (Coefficient of Variation) |
5 Max. temperature of Warmest Month | 16 Precipitation of Wettest Quarter |
6 Min. Temperature of Coldest Month | 17 Precipitation of Driest Quarter |
7 Temperature Annual Range (5–6) | 18 Precipitation of Warmest Quarter |
8 Mean Temperature of Wettest Quarter | 19 Precipitation of Coldest Quarter |
9 Mean Temperature of Driest Quarter | |
10 Mean Temperature of Warmest Quarter | |
11 Mean Temperature of Coldest Quarter |
QGIS 2.8 (
Two different statistical analyses were performed per species: the Jackknife test to evaluate the weight or importance of each variable (
The models showed good results for all nine bumble bee species, with AUC values above 0.9, although the number of occurrences for Bombus excellens, Bombus melaleucus, and Bombus transversalis was low. Bombus pauloensis and Bombus hortulanus had the highest number of occurrences, being common in collections (Table
Bumble bee species data for the area under the curve (AUC obtained by niche modeling), number of occurrences used for modeling, and the distributional group selected for each species according to our results and previous results by
Species | AUC 1970-2000 | Distributional group |
---|---|---|
Bombus excellens | 0.951 | Mid-mountain and inter-Andean |
Bombus funebris | 0.976 | High Mountain |
Bombus hortulanus | 0.977 | High Mountain |
Bombus melaleucus | 0.965 | Mid-mountain and Inter-Andean |
Bombus pauloensis | 0.911 | Mid-mountain and Inter-Andean |
Bombus pullatus | 0.917 | Mid-mountain and Inter-Andean |
Bombus robustus | 0.983 | High Mountain |
Bombus rubicundus | 0.973 | High Mountain |
Bombus transversalis | 0.999 | Amazon and Eastern Plains Basin |
The potential distribution areas for B. funebris, B. hortulanus, B. robustus, and B. rubicundus occurred only in high mountain departments, along a range of different small and fragmented areas, with a high probability of occurrence in the central part of the Eastern Andes Range (Fig.
This group of species showed a wider altitudinal and potential distribution than the high mountain species, along and between the three mountain ranges (Fig.
For B. melaleucus, its potential distribution was indicated along the three mountain ranges. B. pauloensis exhibited a continuous distribution over the Andean region. B. pullatus was concentrated along low altitudinal areas of the inter-Andean valleys.
Bombus transversalis is the only species with a potential distribution along the Amazon Forest, the Andean foothills, and the Orinoquia region. The eastern part of Meta presented the highest altitudinal points for the species at 1150 m. The northern area of the biogeographic region of Chocó showed optimal bioclimatic conditions for B. transversalis, but no occurrence point was found or used there for this model (Fig.
The variables with the highest weight for each species are presented in Table
Bioclimatic variables had the highest weight in the model for each species, according to the results obtained by the Jackknife test.
Species (Distributional group) | Variables with the highest weight (Percent of contribution) | ||
---|---|---|---|
1st | 2nd | 3rd | |
Bombus excellens (Mid-Mountain and inter-Andean) | 8. Mean Temperature of Wettest Quarter 94.35% | 1. Annual Mean Temperature 94.27% | 6. Min. Temperature of Coldest Month 94.04% |
Bombus funebris (High Mountain) | 4. Temperature Seasonality 95.50% | 6. Min. Temperature of Coldest Month 94.54% | 1. Annual Mean Temperature 94.03% |
Bombus hortulanus (High Mountain) | 11. Mean Temperature of Coldest Quarter 96.72% | 1. Annual Mean Temperature 96.63% | 10. Mean Temperature of Warmest Quarter 96.61% |
Bombus melaleucus (Mid-Mountain and inter-Andean) | 5. Max. Temperature of Warmest Month 95.20% | 11. Mean Temperature of Coldest Quarter 95.16% | 1. Annual Mean Temperature 95.14% |
Bombus pauloensis (High Mountain) | 6. Min. Temperature of Coldest Month 91.30% | 11. Mean Temperature of Coldest Quarter 91.18% | 8. Mean Temperature of Wettest Quarter 91.09% |
Bombus pullatus (Mid-Mountain and inter-Andean) | 11. Mean Temperature of Coldest Quarter 83.91% | 8. Mean Temperature of Wettest Quarter 83.53% | 4. Temperature Seasonality 76.67% |
Bombus robustus (High Mountain) | 5. Max. Temperature of Warmest Month 93.74% | 6. Min. Temperature of Coldest Month 93.72% | 8. Mean Temperature of Wettest Quarter 93.51% |
Bombus rubicundus (High Mountain) | 6. Min. Temperature of Coldest Month 96.70% | 1. Annual Mean Temperature 96.70% | 8. Mean Temperature of Wettest Quarter 96.48% |
Bombus transversalis (Amazon and Eastern Plains Basin) | 13. Precipitation of Wettest Month 98.88% | 11. Mean Temperature of Coldest Quarter 97.87% | 6. Min. Temperature of Coldest Month 96.94% |
These potential distribution maps for bumble bees in Colombia improve the previous maps available (
The four high-mountain species are associated with the paramo and high Andean ecosystems (
The most important bioclimatic variables for bumble bees in Colombia are related to temperature (Table
The potential distribution for Colombian bumble bee species is reported here. Seven of them show a restricted distribution, shaped mainly by temperature restrictions. These results are a direct contribution to knowledge about Colombian bumble bee species, constituting useful knowledge for conservation, territorial planning, protection plans, and environmental management. Thus, by obtaining the most suitable areas for a species, this study can provide the ideal location to breed and reproduce the native species. These species contribute to improving the agricultural productivity of several regions by pollination. Likewise, this information can help avoid the introduction of foreign species for this purpose (
We would like to thank M. Sc. Andrea Lorena García and M. Sc. Daniela Hoyos-Benjumea from the Insect collection at Universidad del Quindío, Ph.D. Diego Riaño from the bee collection in Universidad Militar Nueva Granada, Ph.D. Francisco Serna from the entomological museum Universidad Nacional, School of Agronomy, Bogotá, and Ph.D. Sergio Orduz and biol. Carlos Londoño from the entomological museum Francisco Luis Gallego for their contributions to occurrence data. We are indebted to Ph.D. Sydney A. Cameron for comments and suggestions that improved this manuscript.
We acknowledge financial support by Stiftung Universität Hildesheim.
Points used for the Niche modeling of Bumblebee species (Hymenoptera, Apidae, Bombus) in Colombia
Data type: COL (excel document).
Explanation note: Ocurrence data used for the modeling.