Wild pollinator activity negatively related to honey bee colony densities in urban context


Authors: Lise Ropars aff001;  Isabelle Dajoz aff002;  Colin Fontaine aff003;  Audrey Muratet aff004;  Benoît Geslin aff001
Authors place of work: IMBE, Aix Marseille Univ, Avignon Université, CNRS, IRD, Marseille, France aff001;  Institut d’Ecologie et des Sciences de l’Environnement de Paris (iEES Paris UMR 7618) Equipe Ecologie et Evolution des réseaux d’interactions, Université Paris Diderot, CNRS-SU, Paris, France aff002;  Centre d’Ecologie et des Sciences de la Conservation (CESCO UMR 7204), CNRS-Muséum National d’Histoire Naturelle-SU, Paris, France aff003;  Agence Régionale de la Biodiversité en Île-de-France (ARB ÎdF), Paris, France aff004;  Laboratoire Image, Ville, Environnement (LIVE UMR 7362), Université de Strasbourg, Strasbourg, France aff005
Published in the journal: PLoS ONE 14(9)
Category: Research Article
doi: 10.1371/journal.pone.0222316

Summary

As pollinator decline is increasingly reported in natural and agricultural environments, cities are perceived as shelters for pollinators because of low pesticide exposure and high floral diversity throughout the year. This has led to the development of environmental policies supporting pollinators in urban areas. However, policies are often restricted to the promotion of honey bee colony installations, which resulted in a strong increase in apiary numbers in cities. Recently, competition for floral resources between wild pollinators and honey bees has been highlighted in semi-natural contexts, but whether urban beekeeping could impact wild pollinators remains unknown. Here, we show that in the city of Paris (France), wild pollinator visitation rates are negatively correlated to honey bee colony densities present in the surrounding landscape (500m –slope = -0.614; p = 0.001 –and 1000m –slope = -0.489; p = 0.005). Regarding the morphological groups of wild pollinators, large solitary bee and beetle visitation rates were negatively affected by honey bee colony densities within a 500m buffer (slope = -0.425, p = 0.007 and slope = - 0.671, p = 0.002, respectively) and bumblebee visitation rates were negatively affected by honey bee colony density within a 1000m buffer (slope = - 0.451, p = 0.012). Further, lower interaction evenness in plant-pollinator networks was observed with high honey bee colony density within a 1000m buffer (slope = -0.487, p = 0.008). Finally, honey bees tended to focus their foraging activity on managed rather than wild plant species (student t-test, p = 0.001) whereas wild pollinators equally visited managed and wild species. We advocate responsible practices mitigating the introduction of high density of honey bee colonies in urban environments. Further studies are however needed to deepen our knowledge about the potential negative interactions between wild and domesticated pollinators.

Keywords:

Biology and life sciences – Organisms – Eukaryota – Animals – Invertebrates – Arthropoda – Insects – Hymenoptera – Bees – Honey bees – Bumblebees – Beetles – Plants – Flowering plants – Psychology – Behavior – Animal behavior – Foraging – Zoology – Animal behavior – Foraging – Plant science – Plant anatomy – Flowers – Social sciences – Psychology – Behavior – Animal behavior – Foraging – Research and analysis methods – Mathematical and statistical techniques – Mathematical models

Introduction

The recent decline of pollinating insect populations is driven by a conjunction of factors, including habitat fragmentation, use of pesticides, multiplication of pathogens, global warming and the decline of the wild flora [1]. Agricultural landscapes have changed, harbouring fewer floral resources and habitats to support diverse pollinating communities [2,3]. Consequently, many agricultural landscapes are becoming less conducive for pollinators and for beekeeping activities [4]. At the same time, areas that were previously rarely exploited by beekeepers are now under a strong pressure to receive apiaries; this is the case in natural habitats and cities [5,6]. Indeed, cities harbour diverse plant species flourishing all year long due to management practices [7] and heat island effect, thus providing resources throughout the year for pollinators [8]. The low pesticide policies applied in many conurbations may also create favourable conditions for the maintenance of diverse pollinator communities [9]. In the same time, honey bees are perceived as a symbol of biodiversity and ecosystem well-being by many city-dwellers and the media [10]. Many citizens have thus installed colonies as their own contribution to mitigate the pollinator decline [11,12] and urban introductions of honey bee colonies have been promoted by public authorities and decision makers. In many cities, this has translated into very recent and rapid increases in the number of honey bee colonies (e.g. 10 colonies per km2 in London–United Kingdom [13], 15 colonies per km2 in Brussels–Belgium [14]).

However, cities are not depauperate in wild pollinating insects and there is increasing evidence that they host diverse assemblages of wild bees [15,16]. This has led to rising concern about numerous introductions of honey bees in cities, that may negatively impact the wild pollinating fauna through competition for floral resources [11]. In other habitats, such as semi-natural (calcareous meadows [17] or scrubland [18,19]) or agricultural landscapes, several authors have detected exploitative competition between domesticated and wild pollinators through the monopolization of floral resources by honey bees [20,21]. However, we know of no studies that have assessed that honey bee introductions in cities could impact wild pollinator communities and their foraging activity on urban plant communities. Moreover, the effect of increasing honey bee densities has rarely been assessed using network approaches [11]. Massively introduced honey bees might impair the pollination function at community level by, for example, focusing their visits on managed (ornamental) plant species rather than wild ones [11]. Here, we explore those issues in the city of Paris (France), which has recently experienced a strong growth of its honey bee populations within a few years. In 2013, Paris hosted around 300 honey bee colonies, and in 2015 this figure had more than doubled, reaching 687 colonies, corresponding to 6.5 colonies.km-2 (data of the veterinary services of Paris; Fig 1), and has continued to increase since. In this context, our first objective was to analyze the effect of increasing honey bee colony densities on the visitation rates of wild pollinators at the community and morphological group levels. Secondly, we explored how the evenness of plant-pollinator networks was affected by increasing honey bee colony densities. The evenness index of plant-pollinator network reflects how balanced are the links realized by pollinators on plant communities. We expected that the interaction evenness of networks decreases along the gradient of increasing honey bee colony densities. Finally, we investigated the floral preferences of wild and domesticated pollinators for managed or wild plant species.

<h2>Location of honey bee colonies and study sites in the city of Paris.</h2>
Fig. 1.

Location of honey bee colonies and study sites in the city of Paris.


Vegetation height and land use maps were obtained from APUR database (http://opendata.apur.org/datasets/).

Methods

Study sites and plant-pollinator survey

The city of Paris (48°51′12″ N, 2°20′55″ E, Île-de-France, France) is a densely populated urban area (2 220 445 inhabitants in 2014, 105km2). In this city, for three consecutive years, we monitored plant-pollinator interactions in five (in 2014) to seven (in 2015 and 2016) green spaces. We chose these green spaces by their contrasted densities of honey bee colonies in their surroundings (Figs 1 and 2, S1 Data) and for their relative accessibility (access granted by the Bibliothèque nationale de France, campus of Paris Diderot University, Pierre et Marie Curie University, Descartes University, the Institut de Physique du Globe de Paris, and 2 gardens monitored by the Paris Direction des Espaces Verts et de l’Environnement). The distance between sites ranges from 410 to 6 264 meters (S2 Table). Honey bee colony densities were comprised between 0 and 28 colonies within 500m buffers around sites and between 7 and 53 colonies within 1000m buffers around sites (S1 Data). We chose to use the number of honey bee colonies around sites as it has been previously reported to be a good proxy to study potential competitive pressure exerted by honey bees on the wild pollinating fauna [11,18,21,22]. From May to July 2014 and from April to July 2015 and 2016, we carried 8, 11 and 13 observation rounds per green space respectively, spaced out at least by a week. For each round, in each site, we focused our observations on three one-meter square patches chosen to be the most flourished patches within flowerbeds. For each flower visited, we identified the visited plant to the lowest possible taxonomical level (from genus to horticultural variety) according to our knowledge and the taxonomic repository of France [23] and we classified it as managed or wild (S1 Table). Mean richness of visited plant species within patches could vary from 2.5 to 6.5 species depending on the flowering phenology of the plants present in the site. On each patch, we counted the number of visits realized by insect visitors during 5 minutes in 2014 and 2015 and 10 minutes in 2016. Each insect visitor was classified into one of these eight morphological groups: small and large solitary bees, honey bees–A. mellifera, bumblebees, beetles, butterflies, hoverflies and other flies [24]. Observation rounds were performed during warm sunny days (<15°C) with no wind and were carried out between 9 a.m. and 7 p.m. Because of daily variations in meteorological conditions, we alternated our samplings among sites between the morning and the afternoon from one week to the next.

<h2>Study sites within Paris (France) and related plant-pollinator networks.</h2>
Fig. 2.

Study sites within Paris (France) and related plant-pollinator networks.


Bipartite networks are retrieved from the compilation of all observed interactions between pollinators (top bar) and plants (bottom bar). Numbers under networks indicate the year. Each pollinator block represents a morphological group and each plant block represents a species. The width of links is proportional to the number of interactions (i.e. pollinating activity). Pollinator morphological groups are classified by order (left to right) colours and numbers: dark grey-1, small solitary bees; dark red-2, large solitary bees; black-3, syrphids; dark green-4, beetles; purple-5, butterflies; flies-6, light grey; orange-7, bumblebees; yellow-8, honey bees. Sites were represented by a capital letter to make the correspondance to S1 Data. Vegetation height and land use maps were obtained from APUR database (http://opendata.apur.org/datasets/).

Location of honey bee colonies

The law requires beekeepers to report their honey bee colonies to the veterinary services of the city. This is to our knowledge the most accurate data existing to date regarding the location of honey bee colonies within Paris—even if we are aware that some beekeepers may not report their colonies. We used these data to estimate honey bee colony densities within 500- and 1000-meter buffers centred on the study sites using the ArcGIS software (Version 10.2). We chose these buffer sizes to match the mean foraging distances of the majority of wild and domesticated bee species [25,26].

Statistical analyses

Spatial auto-correlation analysis

All statistical analyses were performed using R 3.2.2 (R Development Core Team, 2015). We checked the absence of spatial autocorrelation between our sites and honey bee colonies. We generated a matrix of distances between sites (S2 Table), and built matrices using the number of honey bee colonies in 500m and 1000m around our sites. Mantel tests were carried out between these matrices. No significant spatial autocorrelation was observed for both buffer sizes (500m –p = 0.749; 1000m –p = 0.204, respectively). We also tested the spatial autocorrelation of wild pollinator visitation rates. We found that the visitation rate of all wild pollinators together was spatially autocorrelated (p = 0.025) but not the morphological groups taken separately. Therefore, when analysing the visitation rate of all wild pollinators, we added an autoregressive process of order 1 correlation structure (addition of site coordinates and a random effect with sites nested within years) to deal with spatial autocorrelation.

Floral resources

Floral resources can affect the composition and activity of pollinator assemblages [27]. To account for the spatial heterogeneity of floral resources availability surrounding the study sites, we combined the area covered by the herbaceous, shrub and tree strata with an estimation of the average production of floral units per stratum along the observation period from February to July. We used a map of vegetation height, with a 50cm2/pixel resolution provided by the Parisian Urbanism Agency (APUR - http://opendata.apur.org/datasets/hauteur-vegetation-2015) to calculate the area covered by vegetation of less than 1 meter high (herbaceous stratum), between 1 and 10 meters (shrub stratum), and higher than 10 meters (tree stratum–Fig 1), this within buffers of 500 and 1000 meters centred on our study sites using Geographic Information Systems (GIS, ESRI ARC INFO v. 10.0). To estimate the floral resource production for shrubs and trees along the 6 months, we multiplied their area by their number of floral units/m2. The AgriLand database [28] (S3 Table) allowed us to estimate the number of floral units/m2 at a flowering peak. For these two strata, we considered that the flowering period lasted for 1 month [29]. For the herbaceous stratum, considering the flowering phenology, we modelled a normal distribution pattern (μ = 3; σ2 = 1.220) for 6 months, with the peak of floral production (2700 floral units/m2) occurring in the 3rd month. Using this method, we averaged the number of floral units at 1,371 per month for the 6 months of flowering. Although not targeting urban areas, AgriLand database is the most comprehensive database on floral unit production, thus allowing to account for differences in floral resource production among vegetation stratum.

To assess the local floral resources, we calculated a mean richness of visited plant species by pollinators corresponding to the cumulated visited plant species per site per year divided by the number of observation rounds carried per site during the considered year.

Foraging activity analysis

To standardize the observation effort among years, we calculated visitation rates as the number of visits per minute and per flower visitor group on each site and for each year. Visitation rates were analysed using linear mixed effects models (lme, package “nlme”, R Development Core Team, 2015) and log transformed to approach normality.

Thus, for each morphological group, fixed effects were a) the honey bee colony densities at 500 or 1000 meters around our sites, b) the estimation of the floral resources available in buffers of the same radius and c) the mean visited plant species richness of each site. We included the sites nested within years as a random effect to account for temporal repetitions. We performed model simplification based on the Akaike Information Criterion (AIC) and chose the best fit model with ΔAIC > 2 (dredge, package “MuMIn”, R Development Core Team 2015 –S6 Table) [30] and then obtained the associate P-values using the Maximum Likelihood method (ML). All variables were scaled to make their estimated effects comparable.

We also checked the correlations between honey bee visitation rate and the visitation rate of other morphological groups. Fixed effects were the visitation rate of honey bees and the mean visited plant species richness; sites nested within years was added as a random effect. We selected the best fit model based on the Akaike Information Criterion (AIC) with a ΔAIC > 2 and then obtained the associate P-values using the Maximum Likelihood method (ML). We only detected that visitation rates of wild pollinators and bumblebees were significantly correlated to the local visited plant species richness. Once the mean richness of visited plants was taken into account in our models, we did not detect any significant relation between honey bee visitation rates and wild pollinator visitation rates all together or separated by morphological groups (S4 Table).

Plant-pollinator network analysis

To determine the impact of honey bee colony densities on the structure of plant-pollinator networks, we constructed 19 quantified interaction networks linking flower visitor morphological groups excluding honey bees to plant species, one per site and year. Interaction frequencies were calculated as the number of visits per minute. The structure of the interaction networks was assessed by the interaction evenness using the “bipartite” package [31]. Interaction evenness is bounded between 0 and 1, and derived from the Shannon index, H = pijlog2pij/log2F, where F is the total number of plant–pollinator interactions in the matrix and pijis the proportion of those interactions involving plant i and pollinator j [32,33]. This index reflects how balanced are the interaction strengths between plants and pollinators. It decreases as the network is dominated by few highly frequent interactions and increases when the number of interactions is uniformly distributed [34]. We analyzed the interaction evenness using the same statistical models than for the visitation rate analysis, fixed effects were a) the honey bee colony densities at 500 or 1000 meters around sites, b) the estimation of the floral resources available in a buffer of the same radius and c) the mean plant species richness of each site. A model simplification based on the ΔAIC > 2 was used (S7 Table) [30] and the associated P-values were obtained using the Maximum Likelihood method (ML). We included the year nested within sites as a random effect to account for temporal repetitions.

Floral preferences analysis

To assess the pollinator floral preferences of both wild pollinators and honey bees, we summed their visitations on managed or wild plant species per site and per year (S1 Table). To consider the respective availability of both plant types, visitation rates of pollinator groups on managed or wild plants species were weighted by the percentage of managed and wild species recorded at each site and year. The number of each plant type sampled per year is available in the S5 Table. Floral preferences were tested using Student t-tests comparing the visitation rates of pollinator groups between managed and wild plants.

Results

Effect of honey bee colony densities on wild pollinator visitation rates

Pollinators were monitored for a total of 3,120 minutes during which we recorded 795 individual plant-pollinator links, totalling 32,694 visits on plants (16% of small solitary bees, 10% of large solitary bees, 12% of bumblebees, 1% of beetles, 6% of hoverflies, 4% of flies, 1% of butterflies and 50% of honey bees). 687 honey bee colonies were declared in Paris in 2015, which equates to an average density of 6.5 colonies/km2. Visitation rates of wild pollinators were negatively related to the density of honey bee colonies at both spatial scales (Figs 3 and 4, and Table 1, 500m –slope = -0.614; p = 0.001 –and 1000m– slope = -0.489; p = 0.005). Large solitary bees performed significantly fewer visits when the density of honey bee colonies increased within 500 meter buffers around our observation sites (Fig 3, Table 1, slope = -0.425; p = 0.007). This trend was also significant for beetles (Fig 3, Table 1, slope = - 0.671; p = 0.002). The visitation rate of bumblebees significantly decreased when the density of honey bee colonies increased within 1000 meter buffers (Fig 4, Table 1, slope = - 0.451; p = 0.012). The visitation rate of honey bees was positively correlated with the number of honey bee colonies within 500 meter buffers (Fig 3, Table 1 –slope = 0.501; p = 0.020). However, we did not record any significant increase in the visitation rate of honey bees with the increased density of colonies within 1000 meter buffers. Finally, we did not find any effects of honey bee colony densities on the visitation rate of other morphological groups of pollinators such as small solitary bees, flies, hoverflies and butterflies (ΔAIC < 2 between null models and convenience models or models containing only resources or richness variables–S6 and S7 Tables).

<h2>Morphological group visitation rates along the gradient of honey bee colony number in 500m around our observation sites.</h2>
Fig. 3.

Morphological group visitation rates along the gradient of honey bee colony number in 500m around our observation sites.


Regressions of best-fit models were represented for each morphological group. When best-fit models included multiple explanatory variables, partial residual regressions were plotted.

<h2>Morphological group visitation rates along the gradient of honey bee colony number in 1000m around our observation sites.</h2>
Fig. 4.

Morphological group visitation rates along the gradient of honey bee colony number in 1000m around our observation sites.


Regressions of best-fit models were represented for each morphological group. When best-fit models included multiple explanatory variables, partial residual regressions were plotted.

Tab. 1.

Detailed effects of honey bee colony densities on wild pollinator visitation rates.

<h2>Detailed effects of honey bee colony densities on wild pollinator visitation rates.</h2>

Results of the best linear mixed-effects models containing the colonies number as response variable, and floral resources and richness, as covariables for each buffer scale. Model selection was performed according to the AIC criterion.

Effect of honey bee colony densities on network structure

Regarding the structure of the pollination networks, we found that the evenness of interactions between wild pollinators, excluding honey bees, and plants was negatively related to honey bee colony density within 1000 meter buffers (1000m –slope = -0.487; p = 0.008 –Fig 5, Table 2). We did not find any variation of the interaction evenness with the increased number of honey bee colonies within the 500 meter buffers. The model containing colony densities was equivalent to the null model (ΔAIC < 2—S7 Table).

<h2>Partial regression of the interaction evenness along the gradient of increased honey bee colony numbers in 1000m buffers around our observation sites, once the effect of resources was taken into account.</h2>
Fig. 5.

Partial regression of the interaction evenness along the gradient of increased honey bee colony numbers in 1000m buffers around our observation sites, once the effect of resources was taken into account.


Tab. 2.

Detailed effets of honey bee colony densities on the interaction evenness.

<h2>Detailed effets of honey bee colony densities on the interaction evenness.</h2>

Results of best linear mixed effects models on interaction eveness containing colonies number and floral resources as covariates. Model selection was performed according to the AIC criterion.

Floral preferences between wild and domesticated pollinators

We found that wild pollinators visited significantly more wild plant species than honey bees (t-test, p = 0.022). Furthermore, honey bees significantly preferred foraging on managed plant species than on wild ones (t-test, p = 0.001; Fig 6) whereas wild pollinators had no preference for a particular plant group, managed and wild plant species being equally visited (t-test, p = 0.745; Fig 6).

<h2>Floral preferences of honey bees and wild pollinators.</h2>
Fig. 6.

Floral preferences of honey bees and wild pollinators.


Mean visitation rate on wild and managed flora weighted by the percentage of managed and wild plant species sampled at each site. Mean and 95% confidence intervals are represented. P-values were obtained with Student’s t-test.

Discussion

We showed that in the city of Paris, the visitation rate of wild pollinators and especially the pollinating activity of large solitary bees, bumblebees and beetles, was negatively related to the density of honey bee colonies in the surrounding landscape. This first finding resonates with a growing body of literature highlighting a negative effect of high honey bee colony densities on the wild pollinating fauna [11,21]. Although our study is correlative and does not provide direct evidences, our results are consistent with the hypothesis that honey bee might outcompete the wild pollinating fauna by exploiting flowering rewards (nectar and pollen) [18,19,35].

The negative correlation between the visitation rates of the total wild fauna and the honey bee colony density was found for both scales, within 500 and 1000 meter buffers. When focusing on each pollinator morphological group, this effect was however scale dependent. The visitation rate of large solitary bees and beetles was negatively correlated to honey bee colony density within 500 meter buffers whereas the visitation rate performed by honey bees increased. The bumblebee visitation rates were negatively correlated to the honey bee colony density within 1000 meter buffers. Those differences might be partly due to the foraging abilities of these groups. The large solitary bees includes numerous species which can forage from few hundred meters to several kilometers from their nest, depending on the species considered and the landscape context [25,36]. Bumblebees on the other hand are known to forage at large scales, up to 2800 meters from their nest [37]. Large solitary bees, bumblebees and honey bees have similar dietary requirements, exploiting the same floral resources (pollen and nectar) [17,38]. As summarized in Wojcik et al. 2018 [39], previous studies have found that adding honey bee colonies may negatively affect wild bees and particularly bumblebees especially due to this overlap in resource use.

On the contrary, flies, syrphids and butterflies do not exclusively rely on floral resources, especially during their larval life-stage, which might explain the absence of negative interactions with honey bees [40]. Small solitary bees do require pollen and nectar for both larva and adult stages but because body size and mouthparts length are correlated traits [41], small solitary bees might prefer to seek resources preferentially on shallow flowers [24]. Conversely, larger pollinators such as honey bees and bumblebees, could preferentially forage on plants best adapted to their morphology (preferentially deep flowers—see [42,43]). In that way, small solitary bees might be less sensitive to the increase in honey bee colony densities. The sharp decline of beetles’ foraging activity with the honey bee colony density within 500 meter buffers is more surprising. There is little literature on floral preferences of beetles. Also, their foraging range seems to be highly variable. As examples, Englund (1993) found that Cetonia aurata had a 18m foraging range and Juhel et al. (2017) estimated the foraging range of Brassicogethes aeneus up to 1.2km [44,45]. This underlines the difficulty to relate scale dependent ecological effects with ecological traits of species. For honey bees, we did not detect any increase in their visitation rate with honey bee colony density within 1000 meter buffers. Honey bee foraging range seems to be highly context dependent, from several hundred meters to several kilometers [46,47]. Additionally, Couvillon et al. 2015 demonstrated that honey bee foraging distances both depend on the type of rewards that honey bees seek (nectar or pollen) and on the month considered [48]. The scale to which organisms respond to landscape characteristics thus appear dependent of the context and sensitive to various components acting together. In dense urban habitats, pollinator’s foraging distance might also be sensitive to building height, width or to the spatial distribution of green spaces and floral resources [49].

We also recorded a decrease in the evenness of plant-pollinator interaction networks with the honey bee colony density within 1000 meter buffers. Interaction evenness decreases when the network is dominated by few and/or highly frequent interactions. A high evenness has been previously associated with a good network stability [50,51]. Being opposite, a low interaction evenness has been highlighted in degraded ecosystems [52] and in invaded networks [53]. In a previous meta-analysis [11], we showed that the honey bee position in interaction networks is comparable to that often found for invasive pollinators. Here, the lower evenness at high colonies densities within 1000 meter buffers could be due to the decrease of wild pollinators and particularly of bumblebee’s visitation rate. This questions the potential impact of urban beekeeping on the whole interaction network and urges once again the need for news studies regarding this topic.

In parallel, we showed that honey bees tended to significantly focus their visits on managed plant species, whereas wild pollinators did not show preferences between managed and wild plants. Honey bees often focus their visits on the most abundant resources to cover the colony needs [54] and ornamental flowerbeds might thus be attractive for them. Among the species most visited by honey bees, we indeed found Lavandula sp. and Geranium sanguineum which are common in ornamental flowerbeds. In the other hand, wild flowers received significantly more visits from wild pollinators and might rely more on the wild fauna for pollination. The observed decline of the wild fauna visitation rate associated with high colony densities may have negative consequences for the reproduction of this wild flora. Nevertheless, several other factors might explain insect’s flowers preferences and foraging choices such as the morphology, the color, the amount of resources, or the life span of flowers [55].

As this study took place in a city, urban environment may provide a large range of confounding factors such as pollutions, pesticides or floral resources quantity, which could also explain the observed decline in the foraging activity of some morphological groups. Few studies have explored the impact of pollutants on bees in urban areas, Moroń et al., 2012 & 2014 demonstrated that heavy metal pollution decreased the diversity and the abundance of solitary bees and can reduce the fitness of Osmia rufa, a common Megachilidae in urban environments [56,57]. Concerning pesticides, in Paris, public parks do not use any pesticides in their management since 2008, which limits the impact of this factor. Finally, we observed pollinators activity from April to July corresponding to the peak of floral resources, and consequently at this period there is a large amount of nectar and pollen available to flower visitors. However, along a year, the quantity of floral resources fluctuates and leads to peaks and gaps in floral resources abundance [58]. Here, we found a negative relation between wild pollinator foraging activities and honey bee colony densities but the intensity of this relation could be modulated by the amount of resources available that could be less abundant in early spring or late summer for example.

In this study, the proxy we did use to study the potential impact of honey bees on the wild fauna was the density of honey bee colonies around our sites, and this proxy has been used in several other studies [18,22]. We did not however find any significant correlations between the visitation rate of honey bees and the visitation rate of wild pollinators. We advocate here that the honey bee colony density variable represents a more continuous pressure on the wild fauna that simply the foraging activity of honey bees at a given time. However, at this stage, and as underlined by Mallinger et al. (2017), other critical parameters such as wild pollinator reproductive success (fitness), population or community dynamics are yet rarely explored [21]. This lack of knowledge impedes us to have a more comprehensive view of the potential impact of high honey bee colony densities on the wild pollinating fauna.

Nevertheless, numerous cities around the globe have experienced recent and fast increases in honey bee colony densities. The average colony density in Paris (6.5 colonies/km2) is higher than the national level (2.5 colonies/km2) but is far below other cities such as Brussels (15 colonies/km2) or London (10 colonies/km2) [13,14,59]) and cities also harbor a non-negligible diversity of wild pollinators [6062]. Altogether, our results not only question the fast development of urban beekeeping and the enthusiasm of citizens and mass media for the installation of hives in cities, but also some of urban management practices supposedly conducted to sustain biodiversity. This underlines the need for new studies exploring how domestic and wild pollinators coexist in urban habitats. In conclusion, we suggest that stakeholders should take into account the impacts that apiaries could have on the wild fauna [19,35]. If the capacity of urban ecosystems to provide the pollination function is to be preserved, land owners may focus their management practices on increasing floral resources and nesting habitats for pollinators in urban environments instead of adding honey bee colonies.

Supporting information

S1 Table [docx]
List of plant species with their status.

S2 Table [docx]
Distance between study sites in meters.

S3 Table [docx]
Open floral unit number per m of vegetative cover at the peak of flowering season from AgriLand Database.

S4 Table [docx]
Detailed effects of honey bee visitation rates on wild pollinator visitation rates.

S5 Table [docx]
Number of plants occurrences sampled per year, according to their status (wild or managed species).

S6 Table [docx]
Results of morphological groups’ model selection based on AIC criterion.

S7 Table [docx]
Results of interaction evenness model selection based on AIC criterion.

S1 Data [csv]
Dataset of the study.


Zdroje

1. Potts SG, Imperatriz-fonseca V, Ngo HT, Aizen MA, Biesmeijer JC, Breeze TD. Safeguarding pollinators and their values to human well-being. Nature. Nature Publishing Group; 2016; 1–10. doi: 10.1038/nature20588 27894123

2. Vanbergen AJ, Garratt MP, Vanbergen AJ, Baude M, Biesmeijer JC, Britton NF, et al. Threats to an ecosystem service: Pressures on pollinators. Front Ecol Environ. 2013;11: 251–259. doi: 10.1890/120126

3. Senapathi D, Carvalheiro G, Biesmeijer JC, Dodson C, Evans RL, Mckerchar M, et al. The impact of over 80 years of land cover changes on bee and wasp pollinator communities in England. Proc R Soc Biol Sci. 2015;282: 20150294–20150294. doi: 10.1098/rspb.2015.0294 25833861

4. Otto CR V, Roth CL, Carlson BL, Smart MD. Land-use change reduces habitat suitability for supporting managed honey bee colonies in the Northern Great Plains. Proc Natl Acad Sci U S A. 2016;113: 10430–5. doi: 10.1073/pnas.1603481113 27573824

5. Geldmann J, González-Varo JP. Conserving honey bees does not help wildlife. Science (80-). 2018;359: 392–393. doi: 10.1126/science.aar2269 29371456

6. Colla SR, MacIvor JS. Questioning public perception, conservation policy, and recovery actions for honeybees in North America. Conserv Biol. 2017;31: 1202–1204. doi: 10.1111/cobi.12839 27624856

7. Pardee GL, Philpott SM. Native plants are the bee ‘ s knees: local and landscape predictors of bee richness and abundance in backyard gardens. Urban Ecosyst. 2014;17: 641–659. doi: 10.1007/s11252-014-0349-0

8. Harrison T, Winfree R. Urban drivers of plant-pollinator interactions. Evans K, editor. Funct Ecol. 2015;29: 879–888. doi: 10.1111/1365-2435.12486

9. Fortel L, Henry M, Guilbaud L, Guirao AL, Kuhlmann M, Mouret H, et al. Decreasing Abundance, Increasing Diversity and Changing Structure of the Wild Bee Community (Hymenoptera: Anthophila) along an Urbanization Gradient. Smith MA, editor. PLoS One. 2014;9: e104679. doi: 10.1371/journal.pone.0104679 25118722

10. Smith TJ, Saunders ME. Honey bees: the queens of mass media, despite minority rule among insect pollinators. Leather SR, Packer L, editors. Insect Conserv Divers. 2016;9: 384–390. doi: 10.1111/icad.12178

11. Geslin B, Gauzens B, Baude M, Dajoz I, Fontaine C, Henry M, et al. Massively Introduced Managed Species and Their Consequences for Plant–Pollinator Interactions. Adv Ecol Res. 2017;57: 1–53. doi: 10.1016/j.talanta.2011.11.033 22265567

12. Stange E, Zulian G, Rusch G, Barton D, Nowell M. Ecosystem services mapping for municipal policy: ESTIMAP and zoning for urban beekeeping. One Ecosyst. 2017;2: e14014. doi: 10.3897/oneeco.2.e14014

13. Alton K, Ratnieks F. To bee or not to bee. Biol. 2016;60: 12–15.

14. L ‘ apiculture. Brusselles Environ. 2017; 1–2. Available: http://www.environnement.brussels/thematiques/alimentation/produire-mes-aliments/que-produire-en-ville/lapiculture#

15. Baldock KCR, Goddard MA, Hicks DM, Kunin E, Mitschunas N, Osgathorpe LM, et al. Where is the UK ‘ s pollinator biodiversity? The importance of urban areas for flower- visiting insects. Proc R Soc Biol Sci. 2015;282: 20142849. doi: 10.1098

16. Sirohi MH, Jackson J, Edwards M, Ollerton J. Diversity and abundance of solitary and primitively eusocial bees in an urban centre: a case study from Northampton (England). J Insect Conserv. 2015;19: 487–500. doi: 10.1007/s10841-015-9769-2

17. Steffan-Dewenter I, Tscharntke T. Resource overlap and possible competition between honey bees and wild bees in central Europe. Oecologia. 2000;122: 288–296. doi: 10.1007/s004420050034 28308384

18. Torné-Noguera A, Rodrigo A, Osorio S, Bosch J. Collateral effects of beekeeping: Impacts on pollen-nectar resources and wild bee communities. Basic Appl Ecol. Elsevier GmbH; 2016;17: 199–209. doi: 10.1016/j.baae.2015.11.004

19. Henry M, Rodet G. Controlling the impact of the managed honeybee on wild bees in protected areas. Sci Rep. 2018;8: 9308. doi: 10.1038/s41598-018-27591-y 29915290

20. Magrach A, González-Varo JP, Boiffier M, Vilà M, Bartomeus I. Honeybee spillover reshuffles pollinator diets and affects plant reproductive success. Nat Ecol Evol. 2017;1: 1299–1307. doi: 10.1038/s41559-017-0249-9 29046536

21. Mallinger RE, Gaines-Day HR, Gratton C. Do managed bees have negative effects on wild bees?: A systematic review of the literature. Raine NE, editor. PLoS One. 2017;12: e0189268. doi: 10.1371/journal.pone.0189268 29220412

22. Herbertsson L, Lindström SAM, Rundlöf M, Bommarco R, Smith HG. Competition between managed honeybees and wild bumblebees depends on landscape context. Basic Appl Ecol. Elsevier GmbH; 2016;17: 609–616. doi: 10.1016/j.baae.2016.05.001

23. Gargominy O, Tercerie S, Régnier C, Ramage T, Dupont P, Daszkiewicz P, et al. TAXREF v11, référentiel taxonomique pour la France: méthodologie, mise en oeuvre et diffusion. Museum national d’Histoire naturelle, Paris. Rapport Patrinat 2017–116. 2017. doi: 10.1109/TDEI.2009.5211872

24. Geslin B, Gauzens B, Thébault E, Dajoz I. Plant Pollinator Networks along a Gradient of Urbanisation. Ollerton J, editor. PLoS One. 2013;8: e63421. doi: 10.1371/journal.pone.0063421 23717421

25. Zurbuchen A, Landert L, Klaiber J, Müller A, Hein S, Dorn S. Maximum foraging ranges in solitary bees: only few individuals have the capability to cover long foraging distances. Biol Conserv. Elsevier Ltd; 2010;143: 669–676. doi: 10.1016/j.biocon.2009.12.003

26. Wright IVR, Roberts SPM, Collins BE. Evidence of forage distance limitations for small bees (Hymenoptera: Apidae). Eur J Entomol. 2015;112: 303–310. doi: 10.14411/eje.2015.028

27. Hegland SJ, Boeke L. Relationships between the density and diversity of floral resources and flower visitor activity in a temperate grassland community. Ecol Entomol. 2006;31: 532–538. doi: 10.1111/j.1365-2311.2006.00812.x

28. Baude M, Kunin WE, Boatman ND, Conyers S, Davies N, Gillespie MAK, et al. Historical nectar assessment reveals the fall and rise of floral resources in Britain. Nature. Nature Publishing Group; 2016;530: 85–88. doi: 10.1038/nature16532 26842058

29. Badeau V, Bonhomme M, Bonne F, Carre J, Cecchini S, Chuine I, et al. Les plantes au rythme des saisons. Biotope. 2017.

30. Logan M. Biostatistical Design and Analysis Using R [Internet]. Oxford, UK: Wiley-Blackwell; 2010. doi: 10.1002/9781444319620

31. Dormann CF, Fruend J, Gruber B. Visualising bipartite networks and calculating some (ecological) indices. Bipartite Ref Man. 2009; 1–73. doi: 10.1002/sim.4177 21484849

32. Vasquez DP, Chacoff NP, Cagnolo L. Evaluating multiple determinants of plant–animal mutualistic networks. Ecology. 2009;90: 2039–2046. doi: 10.1890/08-1837.1 19739366

33. Tylianakis JM, Laliberté E, Nielsen A, Bascompte J. Conservation of species interaction networks. Biol Conserv. Elsevier Ltd; 2010;143: 2270–2279. doi: 10.1016/j.biocon.2009.12.004

34. Bluthgen N, Fründ J, Vasquez D, Menzel F. What do interaction network metrics tell us about specialization and biological traits? Ecology. 2008;89: 3387–3399. doi: 10.1890/07-2121.1 19137945

35. Cane JH, Tepedino VJ. Gauging the Effect of Honey Bee Pollen Collection on Native Bee Communities. A J Soc Conserv Biol. 2016;0: 1–6. doi: 10.1111/conl.12263

36. Gathmann A, Tscharntke T. Foraging ranges of solitary bees. J Anim Ecol. 2002;71: 757–764.

37. Knight ME, Martin AP, Bishop S, Osborne JL, Hale RJ, Sanderson RA, et al. An interspecific comparison of foraging range and nest density of four bumblebee (Bombus) species. Mol Ecol. 2005;14: 1811–1820. doi: 10.1111/j.1365-294X.2005.02540.x 15836652

38. Goulson D, Stout JC, Kells AR. Do alien bumblebees compete with native flower-visiting insects in Tasmania? J Insect Conserv. 2002;6: 179–189.

39. Wojcik VA, Morandin LA, Davies Adams L, Rourke KE. Floral Resource Competition Between Honey Bees and Wild Bees: Is There Clear Evidence and Can We Guide Management and Conservation? Environ Entomol. 2018; 1–12. doi: 10.1093/ee/nvx174 29145607

40. Balfour NJ, Garbuzov M, Ratnieks FLW. Longer tongues and swifter handling: Why do more bumble bees (Bombus spp.) than honey bees (Apis mellifera) forage on lavender (Lavandula spp.)? Ecol Entomol. 2013;38: 323–329. doi: 10.1111/een.12019

41. Stang M, Klinkhamer PGL, van der Meijden E. Size constraints and flower abundance determine the number of interactions in a plant-flower visitor web. Oikos. 2006;112: 111–121. doi: 10.1111/j.0030-1299.2006.14199.x

42. Fontaine C, Dajoz I, Meriguet J, Loreau M. Functional diversity of plant-pollinator interaction webs enhances the persistence of plant communities. PLoS Biol. 2006;4: 0129–0135. doi: 10.1371/journal.pbio.0040001 16332160

43. Rodríguez-Gironés MA, Santamaría L. Models of optimal foraging and resource partitioning: Deep corollas for long tongues. Behav Ecol. 2006;17: 905–910. doi: 10.1093/beheco/arl024

44. Juhel AS, Barbu CM, Franck P, Roger-Estrade J, Butier A, Bazot M, et al. Characterization of the pollen beetle, Brassicogethes aeneus, dispersal from woodlands to winter oilseed rape fields. PLoS One. 2017;12. doi: 10.1371/journal.pone.0183878 28841712

45. Englund R. Movement patterns of Cetonia beetles (Scarabaeidae) among flowering Viburnum opulus (Caprifoliaceae)—Option for long-distance pollen dispersal in a temperate shrub. Oecologia. 1993;94: 295–302. doi: 10.1007/BF00341330 28314045

46. Steffan-Dewenter I, Münzenberg U, Bürger C, Thies C, Tscharntke T. Scale-dependant effects of landscape context on three pollinator guilds. Ecology. 2002;83: 1421–1432. doi: 10.1890/0012-9658(2002)083[1421:SDEOLC]2.0.CO;2

47. Cunningham SA, Fournier A, Neave MJ, Le Feuvre D. Improving spatial arrangement of honeybee colonies to avoid pollination shortfall and depressed fruit set. J Appl Ecol. 2016;53: 350–359. doi: 10.1111/1365-2664.12573

48. Couvillon MJ, Riddell Pearce FC, Accleton C, Fensome KA, Quah SKL, Taylor EL, et al. Honey bee foraging distance depends on month and forage type. Apidologie. 2015;46: 61–70. doi: 10.1007/s13592-014-0302-5

49. Levé M, Baudry E, Bessa-Gomes C. Domestic gardens as favorable pollinator habitats in impervious landscapes. Sci Total Environ. Elsevier B.V.; 2019;647: 420–430. doi: 10.1016/j.scitotenv.2018.07.310 30086494

50. Orford KA, Vaughan IP, Memmott J. The forgotten flies: the importance of non-syrphid Diptera as pollinators. Proc R Soc B Biol Sci. 2015;282: 20142934. doi: 10.1098/rspb.2014.2934 25808886

51. Ballantyne G, Baldock KCR, Willmer PG. Constructing more informative plant–pollinator networks: visitation and pollen deposition networks in a heathland plant community. Proc R Soc B Biol Sci. 2015;282: 20151130. doi: 10.1098/rspb.2015.1130 26336181

52. Tylianakis JM, Tscharntke T, Lewis OT. Habitat modification alters the structure of tropical host–parasitoid food webs. Nature. 2007;445: 202–205. doi: 10.1038/nature05429 17215842

53. Kaiser-Bunbury CN, Valentin T, Mougal J, Matatiken D, Ghazoul J. The tolerance of island plant-pollinator networks to alien plants. J Ecol. 2011;99: 202–213. doi: 10.1111/j.1365-2745.2010.01732.x

54. Holzschuh A, Dainese M, González-Varo JP, Mudri-Stojnić S, Riedinger V, Rundlöf M, et al. Mass-flowering crops dilute pollinator abundance in agricultural landscapes across Europe. Novotny V, editor. Ecol Lett. 2016;19: 1228–1236. doi: 10.1111/ele.12657 27531385

55. Fenster CB, Armbruster WS, Wilson P, Dudash MR, Thomson JD. Pollination Syndromes and Floral Specialization. Annu Rev Ecol Evol Syst. 2004;35: 375–403. doi: 10.1146/annurev.ecolsys.34.011802.132347

56. Moroń D, Szentgyörgyi H, Skórka P, Potts SG, Woyciechowski M. Survival, reproduction and population growth of the bee pollinator, Osmia rufa (Hymenoptera: Megachilidae), along gradients of heavy metal pollution. Insect Conserv Divers. 2014;7: 113–121. doi: 10.1111/icad.12040

57. Moroń D, Grześ IM, Skórka P, Szentgyörgyi H, Laskowski R, Potts SG, et al. Abundance and diversity of wild bees along gradients of heavy metal pollution. J Appl Ecol. 2012;49: 118–125. doi: 10.1111/j.1365-2664.2011.02079.x

58. Timberlake TP, Vaughan IP, Memmott J. Phenology of farmland floral resources reveals seasonal gaps in nectar availability for bumblebees. J Appl Ecol. 2019; 1–12. doi: 10.1111/1365-2664.13403

59. Chauzat M-P, Cauquil L, Roy L, Franco S, Hendrikx P, Ribière-Chabert M. Demographics of the European Apicultural Industry. vanEngelsdorp D, editor. PLoS One. 2013;8: e79018. doi: 10.1371/journal.pone.0079018 24236084

60. Ropars L, Dajoz I, Geslin B. La diversité des abeilles parisiennes. Osmia. 2018;7: 14–19. doi: 10.1051/lhb/2010051

61. Ropars L, Dajoz I, Geslin B, Ropars L, Dajoz I, Geslin B. La ville un désert pour les abeilles sauvages? J Bot. 2017;79: 29–35.

62. Geslin B, Le Féon V, Kuhlmann M, Vaissière BE, Dajoz I. The bee fauna of large parks in downtown Paris, France. Ann la Société Entomol Fr. 2015;51: 487–493. doi: 10.1080/00379271.2016.1146632


Článek vyšel v časopise

PLOS One


2019 Číslo 9

Nejčtenější v tomto čísle

Tomuto tématu se dále věnují…


Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Léčba bolesti v ordinaci praktického lékaře
nový kurz
Autoři: MUDr. PhDr. Zdeňka Nováková, Ph.D.

Revmatoidní artritida: včas a k cíli
Autoři: MUDr. Heřman Mann

Jistoty a nástrahy antikoagulační léčby aneb kardiolog - neurolog - farmakolog - nefrolog - právník diskutují
Autoři: doc. MUDr. Štěpán Havránek, Ph.D., prof. MUDr. Roman Herzig, Ph.D., doc. MUDr. Karel Urbánek, Ph.D., prim. MUDr. Jan Vachek, MUDr. et Mgr. Jolana Těšínová, Ph.D.

Léčba akutní pooperační bolesti
Autoři: doc. MUDr. Jiří Málek, CSc.

Nové antipsychotikum kariprazin v léčbě schizofrenie
Autoři: prof. MUDr. Cyril Höschl, DrSc., FRCPsych.

Všechny kurzy
Kurzy