Introduction

Insects are declining worldwide (IPBES 2019) and the loss of insects acting as pollinators, such as wild bees, poses a threat to ecosystem services in natural and agricultural systems (Potts et al. 2010). In Germany alone, 53% of all native bee species are red listed (Westrich et al. 2011). Agricultural intensification through increased use of pesticides, the simplification of landscapes and the vanishing of nesting sites and food resources poses a big challenge to most pollinators (Potts et al. 2010; Vanbergen and Insect Pollinators Initiative 2013). The loss of suitable habitats can lead to a functional homogenization of insect communities through filtering for specific species traits but local effects can be mitigated by increasing landscape heterogeneity (Gámez-Virués et al. 2015; Perović et al. 2015).

The availability of floral resources is the main driver of wild bee abundances (Potts et al. 2003) and the loss of floral host species in the past century is a major reason for wild bee declines (Scheper et al. 2014). Heterogeneous landscapes with a high diversity of different habitats harbor diverse pollinator communities (Shackelford et al. 2013; Steckel et al. 2014) and increasing the amount of SNH in a landscape enhances wild bee abundance and species richness (e.g. Steffan-Dewenter et al. 2002; Nayak et al. 2015). In addition to the amount of SNH in a landscape, the quality of SNH, in terms of flower cover of important forage plants might be an essential factor shaping pollinator responses (Bartual et al. 2019) since both a high local flower cover and plant species richness are positively correlated with wild bee abundances and species richness (Potts et al. 2003; Kallioniemi et al. 2017).

Also, the cultivation of mass-flowering crops (MFC) can benefit wild bees by temporary providing high pollen and nectar amounts (Westphal et al. 2003; Diekötter et al. 2014; Rundlöf et al. 2014). The most commonly grown MFC in Europe is oilseed rape (OSR) (Brassica napus L.). But due to extreme drought during summer, cultivation area dropped by 28% from 2018 to 2019, accounting for 7.4% of the total German cropland (Destatis 2020). Other MFC, like grain legumes, are rarely cultivated (1.7% of cropland in Germany), although their cultivation area increased due to their promotion as ecological focus area (EFA) within the CAP greening reform of 2013 (Destatis 2020).

Successful crop pollination depends on the availability of suitable pollinators, whose morphological traits match those of the crops` flowers (Garibaldi et al. 2015). Vice versa, flower visiting pollinator communities of different MFC vary (Garratt et al. 2014). While OSR is visited by a diverse pollinator community (Garratt et al. 2014), only a few pollinator species, mainly honeybees and bumblebees, forage on faba beans (Vicia faba L.) (Garratt et al. 2014; Marzinzig et al. 2018). This is due to the complex flower morphology of faba beans and nectar located deep in the corolla tubes (Bond and Poulsen 1983). Force and a long tongue are necessary to access the flowers` reward (Bailes et al. 2018). Long-tongued bumblebees are therefore most effective pollinators of faba beans (Marzinzig et al. 2018), while solitary wild bees have been shown to be most effective pollinators for oilseed rape (Woodcock et al. 2013).

Regarding biodiversity benefits of MFC, the identity of the flowering crop might determine, which pollinator species are promoted with possible implications for functional trait compositions of pollinator communities (Diekötter et al. 2010). Beneficial effects of early flowering OSR on bumblebees and solitary wild bees are widely known (e.g. Westphal et al. 2003, 2009; Jauker et al. 2012; Diekötter et al. 2014), whereas there are only few studies, showing a positive effect of faba bean cultivation on bumblebees (Knight et al. 2009; Hanley et al. 2011). There is limited understanding of how MFC with differing flower morphologies, interactively affect wild bee communities in SNH of agricultural landscapes. Especially possible landscape-scale effects on the functional trait composition of wild bee communities are understudied.

In this study, we aim to investigate how crop identity of two MFC species affects wild bee densities, species richness and their functional traits in SNH at landscape-scale. Moreover, we focus on the relative importance of local flower cover of SNH, the amount of SNH and landscape diversity for wild bees. We established 30 paired 1 km2 study landscapes, i.e. 15 study landscapes with V. faba fields and 15 without grain legumes. We expect that faba bean cultivation increases mainly the number of bumblebees, the main pollinators of V. faba, by offering mid-seasonal resources. Further, we hypothesize that all wild bees benefit from an enhanced coverage of OSR, from an enhanced amount and local flower cover of SNH in the study landscapes and high landscape diversity. We expect local and landscape resource availability to modify the functional trait composition of wild bee communities, regarding bee size, foraging preferences, sociality and bumblebee tongue length. With respect to faba bean flower morphology, bigger pollinators, which are heavier and have proportional longer tongues and bees known to forage on Fabaceae should be favored by its cultivation. In contrast, OSR with nectar being easily accessible from flowers might promote predominantly small bees with proportional shorter tongues and bee species with no preference for Fabaceae. We expect social bees to be attracted by and to benefit from both MFC stronger than solitary bees. In addition, we hypothesize, that especially long tongued and big bumblebees, might be favored by faba bean cultivation because their morphology facilitates access to the reward of the flowers.

Methods

Study design

We selected 30 paired 1 km × 1 km study landscapes in three different regions in Germany (see Beyer et al. 2020) (Fig. 1). Twelve landscapes were situated close to Göttingen in Central Germany (51°32´N, 9°56´E), ten landscapes close to the North Sea in vicinity of the city Itzehoe (53°55´N, 9°30´E) and eight in eastern Germany close to Leipzig (51°20´N, 12°22´E). 15 landscapes were surrounding a conventionally managed faba bean field (faba bean landscapes: FB +). We ensured that the respective farmers applied diversified crop rotations, i.e. conventional crop sequences that included grain legumes as an additional component for at least three years. Respecting control landscapes (FB−) did not include any grain legumes. Landscape selection was based on a similar landscape composition of the study landscapes within a pair. For further details see Beyer et al. (2020).

Fig. 1
figure 1

Map showing the location of the study landscapes in the three regions of Germany. Landscapes with faba bean cultivation (FB +) are depicted in green and control landscapes (FB−) in orange. On the right, an example of a study landscape is shown with six transects located in grassy and woody semi-natural habitats

All landscape elements with a minimum size of 10 m2 were mapped within the study landscapes by site inspections. As landscape composition measures, we calculated the percentage cover of semi-natural and flower-rich habitats (SNH, Table S1), landscape diversity (Shannon index of 14 land-use types; see appendix for further information) and percentage cover of mass-flowering OSR with QGIS version 2.18 (QGIS Development Team 2016). Our study landscapes covered an independent gradient of SNH and OSR (Fig. S1).

Wild bee sampling

The study was conducted during the summer of 2017 and wild bee sampling was repeated three times between May and August (early summer: 22.05.–19.06., mid-summer: 22.06.–20.07., late summer: 27.07.–22.08.). Oilseed rape full bloom was over when sampling started. Faba beans were in full bloom during the first sampling run in early summer and did not flower during the last two sampling runs for most of the study landscapes. During every sampling run six variable transect walks (50 m × 5 m, 5 min. excluding handling time) were conducted within the most attractive flower rich SNH in each study landscape to cover the entire spectrum of suitable foraging patches for wild bees at landscape scale (Westphal et al. 2008) (Fig. 1). Due to the high spatial and temporal floral variability in SNH, we think that varying the transect location according to floral availability is a good approach to capture the entire pollinator community in each study landscape across the season. To cover different structural SNH types, we distributed the transect walks, if possible, equally within both, woody and grassy SNH types (for details see Table S1 and S2).

Transect walks were conducted between 9 a.m. and 7 p.m. on days with suitable weather conditions (no rain, low wind speed, temperatures > 15 °C). Walking slowly through the transect in one direction, the number of flower-visiting bees was counted. Species, that could not be identified in the field, were caught and taken to the lab for further identification. To avoid excessive killing, we collected only the first five bumblebee individuals of every species/morphospecies per transect for confirmatory identification. Based on similar coloration, we defined three morphospecies: (1) Bombus terrestris, B. lucorum, B. hortorum, B. bohemicus and B. sylvestris, (2) B. lapidarius, B. ruderarius and B. rupestris, (3) B. pascuorum and B. muscorum. All other bumblebee species could be distinguished as separate species. We identified all bumblebees (Mauss 1994) and all non-Bombus wild bees (Amiet et al. 2001, 2007, 2014) in the lab. Because Bombus terrestris, B. lucorum, B. cryptarum and B. magnus are difficult to distinguish, we aggregated them. The percentage of local flower cover within the transect area was estimated during each transect walk (mean: 28.7%, range: 0–95%). Flowering plants, on which we observed pollinator visits were recorded on family level (Fig. S2).

Functional trait composition of pollinator communities

To study the effects of faba bean cultivation and landscape metrics on the functional trait composition of wild bee communities, we chose three functional traits, likely to be affected by local and landscape resource availability: body size, foraging preferences and sociality. Trait analyses were conducted for all sampled wild bees. Since we collected only the first five individuals of every Bombus color group per transect, the sampled bumblebee individuals do not mirror the actual structure of the bumblebee community in a certain study landscape. To account for the community structure, we only considered the first five collected bumblebee individuals per transect for the trait analyses. We consider this sub-sample to be representative as the mean density of recorded bumblebees per transect was 3.8 and in approx. 80% of the cases we only caught five or less bumblebee individuals per transect. We excluded cleptoparasites from the trait analyses (six species) because of their different life-history strategy and undefined foraging preferences. We measured the intertegular distance (ITD) of at least one and maximum five worker individuals per species (57 species in total) as a measure of body size (Table S3). For four species we caught only one single male individual and thus were not able to measure ITD. For the foraging preferences, wild bees were categorized into species collecting pollen from Fabaceae (any species, known to forage on Fabaceae, not specifically Fabaceae specialists) and species not known to exploit Fabaceae (specialists of a different plant family or generalists, not known to visit Fabaceae). For sociality, wild bees were grouped into social and solitary bees. Foraging preference and sociality information were extracted from the literature (Westrich 2018) (Table S3).

Since bumblebees are the main pollinators of faba beans, we additionally tested whether MFC affect the functional trait composition of bumblebee communities concerning bumblebee body size and tongue length. The tongue length of maximum 12 bumblebee workers per species was measured and the tongue length relative to the bees` body sizes were calculated. According to the mean relative tongue length, bumblebee species were assigned to the long- or short-tongued category (Table S4).

Statistical analyses

All analyses were conducted in R version 3.6.1 (R Core Team 2019) using generalized linear mixed models with the glmmTMB package (Brooks et al. 2017) and linear mixed effect models with the nlme package (Pinheiro et al. 2019). We followed a multimodel inference approach (Burnham and Anderson 2002). Based on global models (see details below), we fitted candidate models containing all possible combinations of explanatory variables and their interactions with the help of the dredge function of the MuMIn package (Barton 2019). The number of variables in the single candidate models was restricted to a maximum of seven. Candidate models were ranked by Bayesian Information Criterion (BIC), which is an information criterion similar to the more commonly used Akaike information criterion (AIC), but with a penalty term that also considers the number of samples. BIC thus tends to select less complex models in comparison to AIC (Burnham and Anderson 2002; Link and Barker 2006). We considered all models with a delta BIC < 2 compared to the best fitting model to have important empirical support (Raftery 1995). We used BIC model weights (wi), to assess the posterior probability of each model (Link and Barker 2006) and derived the sum of BIC model weights (Σwi) across all models to rank the importance of explanatory variables.

We used residual plots to validate model assumptions. As a measure of the model´s goodness of fit we calculated marginal (Rm2) and conditional (Rc2) R2 values with the MuMIn package. To test for collinearity, we used the variance inflation factor (VIF), which was below 3 for all cases (Zuur et al. 2010). To test for differences between faba bean and control landscapes or between different runs we performed post hoc tests with the emmeans package (Lenth 2019) at alpha level of 0.05. All graphics were generated with ggplot2 package (Wickham 2016) and predictions and confidence intervals in scatter plots were obtained from the effects package (Fox and Weisberg 2019) and in the boxplots from the emmeans package.

Wild bee densities and species richness

To analyze the effect of local and landscape resource availability on pollinator densities, we performed two separate models for bumblebees and all other wild bees (non-Bombus wild bees). While all observed wild bees (excluding honeybees) were included in the wild bee densities, only identified individuals were used for the wild bee species richness. Overall species richness (see Fig. S3 for species accumulation curve) was analyzed at landscape scale (due to very low species numbers per transect), i.e. richness was calculated across all six transects per study landscape and run and SNH flower cover was averaged across all transects per study landscape and run. Our response variables were: density of bumblebees, density of non-Bombus wild bees (per 250 m2 transect area and 5 min), species richness of all wild bees per study landscape (per 1500 m2 transect area and 30 min), and additionally, density of the long-tongued B. pascuorum/ B. muscorum color group (hereinafter referred to as B. pascuorum densities, since B. muscorum made up only a small proportion of that group). We included faba bean cultivation (FB + / FB−), local flower cover of the SNH transect area (FC), sampling run (early/mid/late summer), percentage of SNH and of OSR per study landscape and landscape diversity as explanatory variables in the global models. Moreover, we added all possible two-way interactions, as we consider each of them as ecologically relevant. All continuous explanatory variables were standardized to a mean of zero and a standard deviation of 1 to improve model convergence. As nested random effects, we included study region, landscape pair, study landscape and sampling run (region/landscape pair/landscape/run) into the pollinator density models and region, landscape pair and study landscape (region/landscape pair/landscape) into the species richness model.

Global generalized linear mixed models with all independent variables and two-way interactions were fitted with Poisson and in case of overdispersion with negative binomial distribution. For pollinator densities, models with and without zero inflation terms were compared and the model with lowest BIC was selected for multimodel inference. The zero inflation was either assumed to be constant across the data or to be related to the local flower cover (Brooks et al. 2017). The model with negative binomial distribution and without zero inflation term was selected for the bumblebee densities and for non-Bombus densities we chose the model with negative binomial distribution and zero inflation related to the local flower cover. For the species richness, we chose the model with Poisson distribution. To test for competition effects of honeybees on wild bee densities, we repeated density analyses as described above, with honeybee densities as additional predictor.

Functional trait composition of wild bee communities

To analyze the effects of faba bean cultivation and landscape metrics on the functional traits of wild bees at landscape level, we aggregated data across transects for each run. Using the FD package (Laliberté and Legendre 2010), we calculated the community weighted means (CWM) for ITD for the entire pollinator community and for bumblebees seperately. We applied linear mixed effect models with Gamma distribution and log link for the entire wild bee community and with normal distribution for the analysis of CWM ITD of bumblebees. Moreover, we used generalized linear mixed models with binomial familiy and logit link to test the effect of local and landscape wide floral resources on the following response variables: the proportion of wild bees collecting pollen from Fabaceae (Fabaceae foragers) and the proportion of social wild bees relative to all bees as well as the proportion of long-tongued bumblebees relative to all bumblebees. The explanatory variables, random effects and multimodel inference procedure were the same as described above for the species richness.

Results

In total we observed 2620 wild bees of which 78.0% were bumblebees. 1342 wild bees were caught and used for species identification, resulting in 66 recorded species (Table S3). The most frequently caught wild bee was Bombus terrestris/lucorum accounting for about 29.7% of all sampled bees, making up 64.6% of all short-tongued bumblebees. B. pascuorum accounted for 24.4% of all wild bees and represented 82.8% of all long-tongued bumblebees. 12.2% of all collected bees were Bombus lapidarius. Lasioglossum pauxillum was the most commonly found non-Bombus bee and accounted for 5.2% of all wild bees. Most visited families of all visited plants in early summer were Fabaceae, Boranginaceae and Rosaceae, while in mid and late summer Fabaceae, Boranginaceae and Asteraceae were most visited (Fig. S2). 44.3% of all non-cuckoo species were categorized as Fabaceae foragers, which corresponded to 87.8% of all caught individuals.

Wild bee densities and species richness

The multimodel inference revealed two models within dBIC < 2 explaining the effect of faba bean cultivation (Σwi = 0.96), local flower cover (Σwi = 1.00) and run (Σwi = 1.00) on bumblebee densities (Tables 1, S5, S6). We found on average 1.6 times more bumblebees in faba bean compared to control landscapes (Fig. 2a). The highest bumblebee densities were observed in run 2 and lowest in run 3 (Fig. S4). Moreover, bumblebee densities increased with increasing local flower cover and this effect was stronger in control compared to bean landscapes (Σwi = 0.66) (Fig. 2b).

Table 1 Summary of the best fitting candidate models (dBIC < 2) and null models for the effects of faba bean cultivation and landscape metrics on bee densities, species richness and functional traits of all wild bees and bumblebees
Fig. 2
figure 2

Effect of faba bean cultivation on bumblebee densities in landscapes with (FB +) and without (FB−) faba bean cultivation (p < 0.001) (a). Effect of local flower cover (%) on bumblebee densities in FB + and FB− landscapes (b). Effect of local flower cover (%) on non-Bombus wild bee densities (c) and on species richness of all wild bees (d). Effect of landscape diversity on wild bee species richness (e). Predicted mean values and 95% confidence intervals are displayed (in red for (a)). Different letters above boxplots indicate significant differences between FB + and FB− landscapes. Predictions are based on the respective best fitting models m.bb1, m.nb1 and m.spr1 (Table 1). Y-axes are sqrt-transformed in (ac) for graphical reasons

There was a single model within dBIC < 2, showing positive effects of local SNH flower cover (Σwi = 1.00), on the density of non-Bombus wild bees (Tables 1, S5, S6, Fig. 2c). There were two best-fitting models, showing a positive effect of local SNH flower cover (Σwi = 1.00) and of landscape diversity (Σwi = 0.67) on the species richness of all wild bees (Tables 1, S5, S6, Fig. 2d, e). Honeybees did not affect wild bee densities (see Tables S7 and S8).

Functional trait composition of wild bee communities

All wild bees

Functional composition of bee communities with respect to the bees’ body sizes (CWM of ITD) were not affected by local or landscape resource availability. No fitted model was superior to the null model (Table 1).

Three models with dBIC < 2 explained faba bean cultivation (Σwi = 0.99), run (Σwi = 1.00), OSR (Σwi = 0.32) and SNH cover (Σwi = 0.40) effects on the proportion of social wild bees within the communities (Tables 1, S5, S6). We found a higher proportion of social bees in control compared to faba bean landscapes in run 1, whereas a higher proportion of social bees was recorded in bean landscapes in comparison with controls in the second and third run (Fig. 3a). The proportion of social bees decreased with increasing amounts of OSR in the surrounding study landscape (Fig. 3b). We found proportions of social bees to increase with increasing percentage of SNH in run 1 and to decrease in run 2 and 3 (Fig. S5).

Fig. 3
figure 3

Effect of (a) faba bean cultivation and of (b) oilseed rape cover on the proportion of social bees (for the three different sampling runs in (a)). Different letters above boxplots indicate significant differences between faba bean (FB +) and control (FB−) landscapes within each run. Predicted mean values and 95% confidence intervals are displayed (in red for (a)). Predictions are based on the second-best fitting model m.s2 (Table 1)

The multimodel inference resulted in one best model explaining faba bean cultivation (Σwi = 0.83) and run (Σwi = 1.00) effects on the proportion of Fabaceae foragers (Tables 1, S5, S6). Patterns were similar as for the proportion of social bees, which might be due to the correlation of the two response variables (spearman r = 0.74, p < 0.001). Proportions of Fabaceae foragers were higher in faba bean landscapes compared to controls in run 2 and 3, whereas the opposite was true for run 1 (Fig. S6). Besides, composition of the bee communities was unaffected by faba bean cultivation but changed by OSR cover, mainly due to common bumblebee species (for further details see supplementary material Table S9 and Fig. S7).

Bumblebees

There were two models within dBIC < 2 explaining the effect of faba bean cultivation (Σwi = 0.40), sampling run (Σwi = 1.00) and local flower cover (Σwi = 1.00) on the proportion of long-tongued bumblebees (Tables 1, S5, S6). We found 51% of the bumblebee community to be long-tongued in faba bean landscapes and 42% in control landscapes (p = 0.07) (Fig. 4a). Highest proportions of long-tongued bumblebees were present in the last sampling run (Fig. S8) and the proportions decreased with increasing local flower cover (Fig. 4b). Focusing on the most dominant long-tongued species, we found higher B. pascuorum densities in faba bean compared to control landscapes (p < 0.01) (Fig. S9). Densities were highest in run 2 and lowest in 1.

Fig. 4
figure 4

Effect of faba bean cultivation (p = 0.07) (a) and of local flower cover (%) (b) on the proportions of long-tongued bumblebees. Effect of faba bean cultivation (p = 0.12) (c) and of local flower cover (%) (d) for faba bean (FB +) and control (FB−) landscapes on the community weighted means (CWM) for intertegular distance (ITD) of bumblebees. Same letters above boxplots indicate non-significant differences between FB + and FB− landscapes. Predicted mean values and 95% confidence intervals are displayed (in red for (a and c)). Predictions are based on the second-best model m.ltb2 for the proportion of long-tongued bumblebees and on the best fitting model m.itdb1 for CWM ITD of Bombus (Table 1)

The CWM for ITD of bumblebees was fitted by a single best model including faba bean cultivation (Σwi = 0.90), local flower cover (Σwi = 0.88) and their interaction (Σwi = 0.86) as predictors (Tables 1, S5, S6). Mean size of bumblebees did not differ between faba bean and control landscapes (Fig. 4c). CWM ITD increased with increasing flower cover in control landscapes, whereas it decreased in faba bean landscapes (Fig. 4d). CWM ITD was higher in faba bean compared to control landscapes when local flower cover of SNH was low (p < 0.01), while at high local flower cover of SNH CWM ITD did not differ significantly (p = 0.07).

Discussion

Wild bee densities and species richness

We showed that faba bean cultivation enhances bumblebee densities in SNH, which is in accordance with our hypothesis. Accordingly, bumblebee densities were enhanced in faba bean landscapes, when sampling within field margins (Beyer et al. 2020). Bumblebees are known to preferentially forage on flowers, that offer high quality, protein-rich pollen, like Fabaceae (Somme et al. 2015) and are, next to honeybees, the most common flower visitors to faba beans in temperate regions (Bond and Poulsen 1983; Garratt et al. 2014; Marzinzig et al. 2018). Due to their tongue length and body size, bumblebees` morphology fits the morphology of faba beans with its very long corolla-tubes. Force and a long tongue is necessary to easily access nectar of faba bean flowers (Bailes et al. 2018). An increased colony density of long-tongued Bombus pascuorum in the vicinity of faba bean fields (Knight et al. 2009) and a higher density of bumblebees in field margins adjacent to bean compared to wheat fields (Hanley et al. 2011) has already been shown. However, the effect of faba bean cultivation on biodiversity at landscape-scale is widely unknown. Our study provides evidence that conventionally managed faba beans increase bumblebee densities at landscape-scale by providing resources during mid-summer after the flowering of OSR. Resource continuity is of major importance for pollinators to persist in simplified agricultural landscapes (Westphal et al. 2009; Schellhorn et al. 2015) and March, June and August/September have been identified as critical times of nectar scarcity for bumblebees (Timberlake et al. 2019). Thus, faba beans might provide essential resources for wild bees in times of food scarcity during June. Resource continuity throughout the season might be especially important for colonial species, while most solitary bees have shorter life cycles and depend on resource availability at the time of their activity period. This might be one reason for why non-Bombus wild bees did not benefit from faba bean cultivation. Another reason is the beans` flower morphology with rewards not easily accessible for small bees with short tongues. Non-Bombus wild bees and bumblebees were closely linked to the quantity of floral resources within the SNH of our study landscapes. Bumblebee densities increased more strongly with increasing local flower cover in control compared to bean landscapes. Faba bean cultivation might have mitigated the concentration of bumblebees on highly attractive SNH through the provisioning of an alternative resource. In addition to common grassland species, woody plant species were frequently used by our wild bees (Rosaceae in run 1), emphasizing their importance as wild bee forage plants (Kämper et al. 2016; Leidenfrost et al. 2020). While SNH flower cover was the most important determinant for non-Bombus wild bee densities recorded in the respective SNH, the cover of SNH had no effect on wild bees. However, when investigating wild bees outside SNH, such as field margins or crop fields, SNH cover has been shown to be an important driver of wild bee densities (Nayak et al. 2015; Beyer et al. 2020). Wild bees have been shown to shift between different crop and non-crop habitats in agricultural landscapes according to their flowering phenology and changing resource availability (Mandelik et al. 2012). Similarly, local floral resource availability drove wild bees in our study. However, only bumblebees but not non-Bombus wild bees were additionally affected by landscape metrics, which might be explained by the larger foraging ranges of bumblebees compared to smaller and less mobile solitary bees (Greenleaf et al. 2007).

In addition to the local flower cover within SNH, landscape diversity enhanced wild bee species richness. More complex landscapes might offer more diverse nesting and foraging habitats supporting diverse pollinator communities with different resource requirements (Shackelford et al. 2013; Perović et al. 2015). Our study provides important evidence that diverse and flower-rich landscapes are prerequisites to promote and conserve species-rich pollinator communities in often species-poor agricultural landscapes.

Contrary to our expectations, the cover of OSR did neither affect bumblebee, non-Bombus wild bee densities, nor wild bee species richness. OSR has been shown to be beneficial for wild bees (Westphal et al. 2003; Diekötter et al. 2014). Nevertheless, Westphal et al. (2009) showed that early positive effects on bumblebee colonies do not translate into a higher colony reproductive success later in the season. Similarly, the early resources provided by OSR in our study did not affect wild bees in SNH. Our sampling started after most OSR fields had ceased blooming in mid-May. Thus, a possible boost through an early resource pulse in April did not translate into persistent effects in our wild bee communities.

Functional trait composition of wild bee communities

All wild bees

Faba bean cultivation affected the proportion of social bees and of faba bean foragers in SNH, depending on the sampling time. While we observed higher proportions in control than in bean landscapes when faba beans were in full bloom, the effect was reversed afterwards as also shown in Beyer et al. (2020). While a higher proportion of social bees might result from an increase of social bees or a decrease of solitary bees or both, those effects were probably driven by bumblebees, the most dominant group of social bees in our study landscapes. They are known to frequently forage on Fabaceae (Somme et al. 2015) and are attracted by mass-flowers (Walther-Hellwig and Frankl 2000). Most of them were likely to forage in faba bean fields instead of SNH, where we sampled, when beans were blooming. Such dilution effects driven by MFC have been already shown for OSR and might lead to diminishing pollination services for less attractive wild plant populations (Holzschuh et al. 2016). Similar effects are conceivable concerning faba bean flowering, pulling social bees away from SNH during their blooming period with possible negative implications for wild Fabaceae populations that are adapted and dependent on bumblebee pollinators. After the flowering of V. faba, we observed a higher proportion of social bees and Fabaceae foragers in bean than in control landscapes. This indicates that faba bean cultivation generally promotes social bees and bees foraging on Fabaceae. Again, this effect might be driven by the bumblebees through an enhanced colony growth facilitated by the additional resources provided by beans in June, as already shown for later flowering crops (Rundlöf et al. 2014).

High amounts of OSR in the surrounding landscape had a slightly negative effect on the proportion of social bees in wild bee communities found in SNH. OSR is visited by a variety of different pollinator groups, with solitary wild bees as most effective pollinators (Woodcock et al. 2013). Solitary bees, which are known to be promoted by OSR (Diekötter et al. 2014; Jauker et al. 2012) might have benefited more than social bees, like bumblebees, which are associated with Fabaceae (Goulson et al. 2005). Our results indicate that high covers of OSR can modify the functional trait composition of wild bee communities in terms of sociality with resulting effects on the community composition.

In contrast to faba beans that enhanced the proportion of social bees and Fabaceae foragers, oilseed rape promoted solitary wild bees. Consequently, a single species of MFC cannot conserve functionally diverse wild bee communities in agricultural landscapes. Biodiversity conservation measures should therefore aim at increasing crop diversity, which is already one of the three CAP greening practices, farmers receive payments for. The introduction of the crop diversification measure in fact might have led to a decreased cultivation of common cereals in favor of an increased cultivation area of rape, sunflower and leguminous crops in 2015 (European Commission 2017). Nevertheless, agricultural management should also target the cultivation of morphologically different MFC to increase functional farmland plant diversity, which enhances pollinator diversity (Sutter et al. 2017).

Bumblebees

In addition to higher bumblebee densities, there was a slight tendency of enhanced proportions of long-tongued bumblebees in faba bean compared to control landscapes. Additionally, we found enhanced densities of B. pascuorum in bean landscapes, which was the most dominant species of the long-tongued bumblebees (more than 80%) in our study sites. Thus, faba bean cultivation promotes particularly long-tongued bumblebee species. Long-tongued bumblebees have been shown to be most effective faba bean pollinators, while short-tongued species often rob nectar by biting holes into the corolla base (Bond and Poulsen 1983; Marzinzig et al. 2018). Many long-tongued bumblebee species are threatened and the decrease of bumblebee forage plants, especially legumes, during the last decades has been postulated as a reason for their decline (Goulson et al. 2005). Most long-tongued bumblebee species emerge late (May/June) in contrast to the early emerging short-tongued bumblebees (Goulson et al. 2005). Likewise, we found highest proportions of long-tongued bumblebees in the latest run from mid-July to mid-August. Dietary specialization, a long tongue and a late emergence time have been related to rarity in bumblebees (Goulson et al. 2005). Hence, the cultivation of faba beans could help to promote rare long-tongued bumblebee species, for which the faba bean provides valuable resources for colony establishment.

An increasing local flower cover led to higher proportions of short-tongued bumblebees. Short-tongued bumblebee species are known to be particularly attracted by mass-flowers (Walther-Hellwig and Frankl 2000) which explains their dominance in bumblebee communities when flower cover of SNH was high. The cover of OSR did not affect the functional composition of the bumblebee communities. Contrary, Diekötter et al. (2010) found decreased densities of long-tongued bumblebees with increasing amounts of OSR in the landscape. This pattern probably arose from competition effects through short-tongued bumblebees, which strongly profited from OSR, leading to higher rates of nectar robbing on red clover (Diekötter et al. 2010).

We found bumblebees to have higher mean sizes in faba bean compared to control landscapes when local flower cover of SNH was low. Faba beans provided additional resources in simplified agricultural landscapes, from which bumblebee colonies were likely to benefit (Knight et al. 2009). Bumblebees in complex landscapes with more floral resources have been shown to be larger than in simple landscapes (Persson and Smith 2011) since fewer resources result in smaller offspring (Oliveira et al. 2016; Radmacher and Strohm 2010). Another explanation might be that bigger bumblebees have proportional longer tongues and can apply more force, which is required to access the V. faba flower rewards (Bailes et al. 2018). The improved resource exploitation of big bumblebees is likely to have driven their dominance of the community. However, when local flower cover of SNH increased, the difference between bean and control landscapes regarding mean bee size disappeared. If bumblebees were resource limited in our study landscapes, more flower-rich SNH were likely to enable a bigger worker size through increased nectar and pollen supply, explaining increasing ITD of bumblebees in control landscapes. Accordingly, Grab et al. (2019) showed that the provisioning of floral resources can buffer against a decrease of bee body sizes due to the simplification of landscapes. Contrary, mean bumblebee size in faba bean landscapes decreased with increasing SNH flower cover. Increased alternative resource availability probably compensated for the dominance of big bumblebees, benefiting from faba beans and exerting competitive pressure on smaller bees. Potential wild bee functional trait composition changes driven by faba beans are therefore likely to be mitigated by a high floral resource supply of SNH, i.e. heterogenous, resource-rich landscapes.

Conclusions

Our study shows that the identity of MFC moderates functional trait compositions of pollinator communities. Different MFC promote different functional groups of wild bees and the dominance of one MFC cannot conserve diverse pollinator communities. High OSR covers shifted bee communities towards a higher proportion of solitary wild bees. Faba bean cultivation enhanced bumblebee densities at landscape-scale and filtered for functional traits associated with Fabaceae. Social wild bees, bees foraging on Fabaceae and long-tongued bumblebee species were promoted. Increased floral resource availability can buffer against changes of the functional trait composition regarding bumblebee size, caused by a dominance of faba beans in resource-poor landscapes. It is already known that the loss of suitable foraging and nesting habitats in simplified agricultural landscapes results in a functional homogenization of insect communities (Gámez-Virués et al. 2015; Perović et al. 2015). Since crop yields are positively correlated with functional pollinator diversity (Garibaldi et al 2015; Woodcock et al. 2019), it is desirable to conserve diverse pollinator communities, which requires a high landscape heterogeneity (Steckel et al. 2014; Perović et al. 2015) and a high functional diversity of farmland plants (Sutter et al. 2017). Farmland floral resources can be provided by either promoting SNH or by the cultivation of MFC and both approaches have already been shown to benefit pollinators (e.g. Westphal et al. 2003; Diekötter et al. 2014). Biodiversity conservation measures should therefore include both, off-field practices targeting SNH as well as on-field practices aiming at crop diversification, especially with respect to flowering crops with various flower morphologies. A high crop and landscape diversity, combined with a high quantity and quality of SNH can help to promote common crop pollinators, to conserve functional diverse pollinator communities and to ultimately secure crop yields.