Using Temporally Resolved Floral Resource Maps to Explain Bumblebee Colony Performance in Agricultural Landscapes

: Wild bumblebees are key pollinators of crops and wild plants that rely on the continuous availability of ﬂoral resources. A better understanding of the spatio-temporal availability and use of ﬂoral food resources may help to promote bumblebees and their pollination services in agricultural landscapes. We placed colonies of Bombus terrestris L. in 24 agricultural landscapes with various degrees of ﬂoral resource availability and assessed di ﬀ erent parameters of colony growth and ﬁtness. We estimated pollen availability during di ﬀ erent periods of colony development based on detailed information of the bumblebee pollen diet and the spatial distribution of the visited plant species. Total pollen availability did not signiﬁcantly explain colony growth or ﬁtness. However, when using habitat maps, the weight gain of colonies, the number of queen cells, and colony survival decreased with increasing distance from the forest. The better explanation of bumblebee performance by forest proximity than by (plant-inferred) pollen availability indicates that other functions of forests than pollen provision were important. The conservation of forests next to agricultural land might help to sustain high populations of these important wild pollinators and enhance their crop pollination services. Combining di ﬀ erent mapping approaches might help to further disentangle complex relationships between B. terrestris and their environment in agricultural landscapes.


Introduction
Animal pollination is crucial for around one third of worldwide food production, with 85% of leading global crop types relying to varying degrees on pollination [1].In addition to domestic honeybees (Apis mellifera L.), wild bees greatly enhance and stabilize crop pollination, and they are often the most effective pollinators [2][3][4][5][6].However, habitat loss and intensive agricultural practices contribute to pollinator declines in different regions of the world [7], while the dependency of global agriculture on pollinators is increasing [8].Hence, supporting wild bees in agricultural landscapes is crucial to future economic and environmental stability [9,10].
Bumblebees are important wild pollinators that increase the yield of many crops, e.g., fruit trees, pepper, pumpkin, strawberries, and tomatoes [11,12].The buff-tailed bumblebee (Bombus terrestris L.) is one of the dominant crop pollinators in Europe [13].Like all wild bees, wild bumblebees rely on foraging and nesting resources provided by the surrounding landscape [14,15].Thus, favorable foraging habitats can enhance pollinator populations and crop pollination at the landscape scale [10,[16][17][18].As habitat types differ in resource availability, habitat maps using broad land use categories are commonly used to predict wild bee development and the service they provide [5].For example, habitat maps explained faster growing and heavier colonies of B. terrestris in suburban gardens than elsewhere [19].In different studies, the pollen deposition of B. terrestris was higher with a decreasing proportion of cropland in the surrounding landscape [20], and B. terrestris colonies had higher reproductive success and survived longer in urban areas than agricultural areas [21].In intensively managed agricultural landscapes, floral resources are mainly offered by small fragments of semi-natural habitats [4,22,23] or single mass flowering crops.Consequently, in addition to habitat maps, maps of single mass flowering species, e.g., apple, oilseed rape, and strawberry, have been found to explain colony development of B. terrestris [24][25][26].Bombus terrestris mainly collects pollen on woody plants [27,28].No single plant accounts for >15% of the total pollen diet of B. terrestris [28].Furthermore, the flowering of single plant species is temporally restricted, and most pollinators use a sequence of specific plant species during their flight season [28].Resources may be limited, especially during early stages of colony development [24,[29][30][31].Hence, comprehensive floral resource maps at the landscape scale that give information on temporal resource availability may further improve our ability to predict the growth and reproductive success of bumblebee colonies [32].
Here, we placed colonies of B. terrestris in 24 agricultural landscapes in southwest Germany, quantified their pollen use, mapped the most commonly used plants (71 species derived from 30 pollen types, offering 95% of the total pollen diet), and calculated a weighed pollen availability index for the early and late flight period of the species in each landscape.We tested the effects of detailed plant-inferred pollen availability versus classical land use maps (i.e., distances from and proportions of land use categories in the landscapes) on colony development (i.e., colony growth and longevity).We tested the following hypotheses: (1) pollen resources for bumblebees are mostly provided by woody semi-natural habitats; (2) a high availability of pollen resources in the landscape enhances colony development; (3) early pollen resources have stronger effects on the weight gain of colonies than late pollen resources; and (4) floral resource maps predict colony development better than habitat maps.

Study Landscapes and Creation of Floral Resource Maps
The study was conducted in the surroundings of Landau in the Upper Rhine Valley, Rhineland-Palatinate, Germany (Electronic Supplementary Materials, Figure S1).A total of 24 agricultural landscapes of 500 m radius were selected to represent gradients of the amounts of early and late pollen resources for B. terrestris [33].Landscape centers were at least 800 m away from each other and were placed in grassy field margins.The cover of plants offering selected pollen types and the major land use types were mapped according to field inspections between June and November 2017 (woody plants) and between April and June 2018 (land use and annual herbaceous plants; for further details see Supplement S1).We used Copernicus (Sentinel-2; L2A-L2C) optical satellite imagery (2016; 10 m resolution), processed by the Federal Agency for Cartography and Geodesy (BKG), as orientation in the field during mapping.We ground-truthed each landscape element via field inspection during our mapping process in 2018 and noted every change if necessary.Obtained data were later digitized as polygon (land use classes) and point layers (single resources), and they were analyzed with QGIS 3.6.2[34] using the above-mentioned satellite imagery as a base map.To create land use/land cover (LULC) maps (hereafter: "habitat maps"), the following habitat types were distinguished: arable land, permanent crops, forest edges (i.e., the first ten meter into forests), forest (i.e., forest interiors without its edges), open woody semi-natural habitat (i.e., hedgerows and single standing trees including those of semi-natural orchards), herbaceous semi-natural habitat (i.e., intensively and extensively managed meadows and pastures), and built-up area (i.e., rural settlements).The proportions of habitat types across the landscapes are shown in Table S2.We also measured the Euclidean distances of colonies to forests because these were important predictors of wild bee development or performance in previous studies [23,[35][36][37].To analyze single resources, the cover of and distance to dominant pollen types were used: Brassica napus L. (hereafter: oilseed rape), because it has been an important predictor of bumblebee densities and colony development in earlier studies [24,38,39], and Cornus, Lonicera, Prunus, Rubus, and Tilia because each of these plant taxa contributed more than 10% of pollen availability to B. terrestris in at least one season during our study (Electronic Supplementary Materials, Table S1).To create floral resource maps, the area covered by plant species offering pollen types that constitute at least 5% of the pollen diet of B. terrestris in our study region in at least one season, and all remaining woody plant species were mapped [28].These maps accounted for the area covered by 71 plant species derived from 30 pollen types.Plants in forest interiors (i.e., deeper than ten meters into a forest) were not mapped because most bees [40], and B. terrestris in particular, prefer open habitat for foraging [41][42][43].Floral resource indices were used as described in [33].Indices considered the relative cover of plants offering pollen types multiplied by their utilization (i.e., percentage of the total collected pollen volume) by B. terrestris during a specific period, summed up over all pollen types.An index value of 1 denotes the average pollen availability of all studied landscapes, while larger or smaller values describe the proportional difference of resource availability in a given landscape relative to the average.Pollen availability was calculated over the whole duration of B. terrestris colony field placement (from mid-April to mid-June) and separately for the early and late phase of colony development.The early season started with the placement of colonies in the field (mid-April) and ended when they reached their maximum weight (end of May).The late season was from when the colonies had their maximum weight until colony termination (i.e., when no evidence of living bumblebees could be detected at the nest boxes during measurements; mid-June).Floral resource maps accounted for 96.3% of early, 83.2% of late, and 94.9% of the total pollen diet of B. terrestris according to the analysis of pollen diets from the 48 colonies when considering 45,900 pollen grains out of 306 samples (Table S1).The composition of pollen diet excluded from index calculations is given in Table S3.

Placement of B. terrestris Colonies
Two commercially bred B. terrestris colonies (STB Control, Aarbergen, Germany) were established in the center of each of the 24 landscapes (48 colonies in total) on a grassy field margin in mid-April 2018.Colonies were even-aged, consisted of one queen and approximately 50-60 workers, and were embedded in a nest box.For acclimatization, a tank filled with sugar water was provided for each colony.Nesting box and tank were enclosed by a cardboard box to protect colonies against unfavorable weather conditions.To protect colonies from soil moisture, boxes were placed on wooden plates that were 10 cm above ground level.Colony entrances were facing south-east.The day after placement, colonies were weighted and colony entrances were opened.After one week of acclimatization, sugar tanks were closed.The weight of nest boxes containing colonies was measured every second week until colony termination.After that, colonies were harvested and frozen at −18 • C for later dissection in the laboratory.Weight gain was obtained by subtracting the maximum colony weight by its initial weight.After colony dissection, cells were counted and assigned to two classes: queen cells (diameter ≥ 11 mm or length ≥ 19 mm if diameter was not measurable because cells were fragmented) and other cells (i.e., male/worker cells with diameter < 11 mm and length < 19 mm; [44]).Male and worker cells could not be consistently differentiated and were therefore grouped together [19,29,31].Days of survival were measured from the day of colony placement until colony termination.Weight gain, the number of cells, and the survival of the two colonies in each landscape were averaged.

Pollen Collection, Preparation and Determination
The pollen diet of B. terrestris was recorded at up to four sampling dates between 24 April and 16 June 2018.At each nest, pollen loads were collected from up to four workers returning from foraging trips per sampling date.The start of the sampling period coincided with the peak flowering of oilseed rape and Crataegus spec., while at the end, Rubus fruticosus L. and Tilia spec.were flowering, which are dominant flowering resources of B. terrestris in the respective seasons [28].Pollen was stored in water, frozen at −18 • C, and then acetolyzed [45].Acetolyzed pollen was mounted on permanent slides in glycerin, and 150 pollen grains per sample were counted starting at a random position of the slide and identified to the highest possible taxonomic resolution using a light microscope (400× magnification), a palynological key [46], a photo atlas [47], and the reference pollen collection of the Institute of Plant Sciences of the University of Bern.Pollen types that contributed more than 3% of the diet of B. terrestris colonies in any period (i.e., before colonies reached their maximum weight or from this moment until colony termination) were classified as key pollen types.Deformed pollen grains that could not be assigned to pollen types (2.2% of detected pollen) and pollen fragments were excluded from the analysis.The pollen type Brassicaceae could not be identified to the species level.However, oilseed rape accounted for 98.3% of the cover of Brassicaceae in the landscapes.

Statistical Analysis
Correlations between colony response (i.e., weight gain, queen cells, and survival) and explanatory variables (i.e., variables used to describe floral resource maps, habitat maps, or single resources) were tested with linear regression.To visualize the interaction among variables, a correlation matrix among all variables was drawn (Figure S3).All variables were standardized before the analysis using the standardize function of the arm package [48] in order to allow for a comparison between effect sizes between models.To compare the predictive power of the mapping approaches, three different sets of linear candidate models were set up for each colony response variable.The explanatory variables of the candidate models were either the pollen availability indices, the parameters derived from habitat maps, or the single pollen resources.The correlations of variables in models were below |r| ≤ 0.7 (Figure S3).Models were compared via the Akaike second-order Information Criterion (AICc) [49][50][51][52] using the dredge function from the MuMin package [53].We compared the overall best model(s) using ∆ i < 2 as a cutoff rule [50,52].In addition to this statistically based model comparison, models with single explanatory variables that contain either one of the focal predictors of this study (resource availability and dominant single pollen resources) or key variables reported in the literature (distance to forests, cover of arable land and built-up area, see introduction) are displayed.To compare the predictive power of the different mapping approaches, the best models of each approach were compared to each other.Data analyses were conducted in R 4.0 [54].Model diagnostic plots were visually checked (residuals vs. fitted values and normal Quantile-Quantile plots).A pollen network graph was created using the package bipartite [55].The plotting of bar plots and linear models was done using the package ggplot2 [56], and the correlation matrix was drawn using the corrplot package [57].
After giving each pollen type a weight proportional to its use by bumblebees, open woody semi-natural habitats (woody semi-natural habitat excluding forests and forest edges) provided 75.2% of the total pollen availability (expressed with indices) in the 24 studied landscapes (Figure 1).In open, woody semi-natural habitats, plants offering Cornus sanguinea, Lonicera xylosteum type, and Prunus type pollen were the most important contributors to pollen availability in the early season (39.8%, 15.4%, and 6.8%, respectively).In the late season, Tilia, L. xylosteum type and Rubus were the most important contributors (25.2%, 23.2%, and 17.4%, respectively).Forest edges contributed to 11.3% of total pollen availability.Additionally in the early season, C. sanguinea and L. xylosteum type were most important contributors (5.0% and 1.6%, respectively); and in the late season Rubus, Tilia, and L. xylosteum type (9.7%, 4.2%, and 2.4%, respectively).The herbaceous semi-natural habitat contributed 3.5% to early pollen availability and 8.5% to late pollen availability.Here, the main contributors in the early season were Papaver rhoeas type, Lamium album type, and Trifolium pratense type (1.4%, 1.3%, and 0.6%, respectively), and the main contributors in the late season were T. pratense type and Phacelia tanacetifolia (5.3% and 2.8%, respectively).In permanent crops (average cover: 4.0%), the most important contributors were fruit trees in the early season (5.2%) and Vitis in the late season (1.8%).Though arable land covered the major part of the landscapes, its contribution to pollen availability over the seasons did not exceed 3.2%.The most important crop pollen resources were oilseed rape in the early season (2.5%) and Asparagus officinalis type in the late season (1.1%).The contributions of built-up area to pollen availability were low (<2.1% in any season).Figure 1 shows the average area of habitat categories in the studied landscapes, as well as their average contributions to pollen availability (indices) during different periods.For a list of the plants detected in the landscapes offering collected pollen types and their contributions to pollen availability during different time periods, see Table S1.

Pollen Availability and Colonies
The total pollen availability in the landscapes during any time did not significantly affect colony weight gain, survival, or the number of queen cells (p > 0.35; Table 1 and Figure S3).Using habitat maps, we found that the weight gain of colonies, queen cells, and survival decreased with distance to forests (Figure 2).In addition to the distance from forests, the best models contained either increases of survival with distance to built-up area or decreased survival with an increasing proportion of it in the landscapes (Table 1).Using important predictors from previous studies, we found that the number of queen cells tended to decrease with the proportion of arable land (Table 1).However, this model was distinctly worse than models containing distance to forests (Δi ≥ 2; Table 1).

Pollen Availability and Colonies
The total pollen availability in the landscapes during any time did not significantly affect colony weight gain, survival, or the number of queen cells (p > 0.35; Table 1 and Figure S3).Using habitat maps, we found that the weight gain of colonies, queen cells, and survival decreased with distance to forests (Figure 2).In addition to the distance from forests, the best models contained either increases of survival with distance to built-up area or decreased survival with an increasing proportion of it in the landscapes (Table 1).Using important predictors from previous studies, we found that the number of queen cells tended to decrease with the proportion of arable land (Table 1).However, this model was distinctly worse than models containing distance to forests (Δi ≥ 2; Table 1).

Pollen Availability and Colonies
The total pollen availability in the landscapes during any time did not significantly affect colony weight gain, survival, or the number of queen cells (p > 0.35; Table 1 and Figure S3).Using habitat maps, we found that the weight gain of colonies, queen cells, and survival decreased with distance to forests (Figure 2).In addition to the distance from forests, the best models contained either increases of survival with distance to built-up area or decreased survival with an increasing proportion of it in the landscapes (Table 1).Using important predictors from previous studies, we found that the number of queen cells tended to decrease with the proportion of arable land (Table 1).However, this model was distinctly worse than models containing distance to forests (∆ i ≥ 2; Table 1).Regarding single pollen resources, weight gain decreased with increasing distance to oilseed rape, and no significant relationship was found with any other variable (Table 1).The number of queen cells showed no significant relationship with any of the single pollen resources.Colony survival decreased with increasing distance to Cornus sanguinea in all of the best-ranking models (Table 1).Alternatively, a negative correlation also existed between colony survival and the distance to Rubus (Δi = 2.11).For correlations among all observed variables, see Figure S3.Regarding single pollen resources, weight gain decreased with increasing distance to oilseed rape, and no significant relationship was found with any other variable (Table 1).The number of queen cells showed no significant relationship with any of the single pollen resources.Colony survival decreased with increasing distance to Cornus sanguinea in all of the best-ranking models (Table 1).Alternatively, a negative correlation also existed between colony survival and the distance to Rubus (∆ i = 2.11).For correlations among all observed variables, see Figure S3.
Table 1.Comparison of models explaining bumblebee colony response variables with pollen availability indices, classical habitat maps, and single resources (see main text for description of mapping approaches).The abbreviation "SNH" is used for "Semi-natural habitat".The Akaike second-order Information Criterion (AICc) and the dredge function from the MuMin package [53] were used to select the best models (∆ i < 2).Delta weight (∆ i ) is the difference between the AICc of the model and the best model.Best models containing pollen availability and important predictors using habitat maps or single pollen resources are displayed regardless of AICc values (see main text).Models listed below a dashed line are not included in the best model set (∆ i < 2).Variables were standardized [48].

Discussion
As expected, B. terrestris mostly used pollen from woody semi-natural habitats.Surprisingly, however, species-specific floral resource maps at the landscape level did not explain the colony development of B. terrestris, although they accounted for the vast part of their pollen diet.Instead, all three parameters of colony performance declined with the distance to forests.This suggests that forest edges had other important functions for bumblebees than pollen provision, such as protection from adverse weather or nectar provision.Alternatively, our results could indicate that the proximity to floral resources is more important than their amount in the landscape, which was expressed by the pollen availability index.In line with a high importance of proximity rather than amount, colony survival decreased with increasing distance to major pollen sources like Cornus sanguinea and Rubus.In turn, the cover of Rubus and Tilia declined with distance from the forest (see Supplementary Materials, Figure S3).Our findings therefore suggest that, in addition to pollination [37], pollinator species richness and flower visitation rate of pollinators [23], bumble bee colony growth and fitness might be positively affected by proximity to forests.Positive effects of forests on B. terrestris may be surprising, because the species is considered an inhabitant of open landscapes [41][42][43].Following the literature on B. terrestris habitat use, we only mapped pollen resources in open habitats and in the first 10 m of forests ("forest edge").Still, we cannot fully exclude that bumblebees were exploiting floral resources, e.g., Acer spec., Castanea sativa Mill., Tilia spec., and Rubus fruticosus, within forests and that our floral resource index is thus incomplete.Flower use by pollinators in the forest canopy is difficult to quantify, and we are unaware of studies that have comprehensively described pollinator communities of trees growing in European forest interiors.Furthermore, other pollen resource plants commonly used by B. terrestris such as Asparagus officinalis L., Brassica napus, Cornus sanguinea L., Papaver rhoeas L., Phacelia tanacetifolia Benth., and Trifolium pratense L. are absent or rare in forest interiors.Of the most important pollen resource plants collected by bumblebees in our study, only Lonicera xylosteum L. and Rubus fruticosus can be commonly found in forest interiors, but they are equally found along forest edges, in hedgerows, and in gardens.Typically, forests are semi-natural habitats that often have positive effects on pollinator richness, visitation rate, or pollination service [5,23,58,59].Apart from food availability, possible benefits of forests for pollinators were summarized in [60]: already established bumblebee colonies may benefit from the reduced daytime temperature in forest interiors in comparison to open habitats during summers (e.g., [61]), and B. terrestris might have benefited from microclimatic conditions along and inside forests during its foraging flights.Other benefits of forests include the reduction of air movement, which leads to reduced energetic costs of foraging flights compared to open habitats (e.g., [60,61]).In addition, in the same study year, B. terrestris was found collecting honeydew from a colony of the giant willow aphid in England, probably due to an increase of nectar sugar concentration while floral nectar resources were simultaneously restricted, followed by exceptional hot and dry weather [62].Our study year was extremely hot, with the highest average temperature ever recorded in Germany along with drought due to low summer rainfall in combination with a high sunshine duration [63].These extreme weather conditions could have caused a shift in limitations from pollen towards nectar resources and potentially enhanced the collection of honeydew in forest interiors.Weather conditions that exceed the thermal tolerance limits of species are likely to increase with climate change [64,65].
Bombus terrestris tend to forage on close-by patches with high resource densities [27,42,66], and their average flight distances lie below or close to 500 m if rewarding resources are available [67,68].However, they were also found foraging up to several kilometers from their nests [42,69].Hence, although bumblebees prefer patches of abundant floral resources close to their colonies, the unexpected low effect of local floral resource availability on colony development might partly be explained by the potentially long foraging distances of B. terrestris.On the other hand, during early phases of colony development, the number of workers is still low and every lack of resources is detrimental to colony development [30,70].In addition, long distance flights are more energy-consuming than short distance flights, making close resources more valuable than resources further away from the nest [70].
Further, a mean foraging range of 275.3 ± 18.5 m with a range of 70-631 m was observed for B. terrestris, indicating that the major part of their foraging may occur at that scale [71].Thus, we believe that our landscape radius of 500 m is still appropriate.
The negative effect of built-up area on colony termination contrasts with findings of increased weight gain, higher numbers of males and queens, higher queen survival, and more food stores in colonies of B. terrestris in suburban areas and sites with varying degrees of urbanization [19,21].Bumblebees can profit from large proportions of beneficial habitats at the outskirts of urban areas related to urban sprawl [72].In our landscapes, built-up areas had a minor role in providing pollen resources compared to some previous studies [73], probably because landscapes were selected to be dominated by agricultural land use.
The increasing amount of arable land in the landscape provides an alternative explanation for the decreased colony performance with increasing distance to forests.The negative effects of arable land on pollinators have commonly been observed (e.g., [20]) and can be explained by, e.g., the negative effects of pesticides or the scarcity of floral resources in intensive agriculture.Pollen availability in our study was lowest in arable land.In addition, longer foraging trips and foraging on scarce, widely distributed, and distant pollen resources (thus less efficient foraging flights) might have had a negative influence on colony development.Large fields in arable land largely lack the beneficial microclimate offered by rural settlements or woody semi-natural habitats [61,72].The decreasing number of queen cells with arable land was in accordance with the decrease of seed set in Vicia faba L. with arable land observed in [33].When we excluded forest distance from our models, weight gain increased with the proximity to oilseed rape fields in the landscape.This resembled the results of [24,38], the authors of which found increasing abundance and weight gain in colonies of B. terrestris with larger amounts of oilseed rape in the landscape.The less strong effect on weight gain on B. terrestris observed in our study may be explained by the lower amounts of oilseed rape in our study (the mean proportion of oilseed rape was 1.8%, compared to 7.0% in [24]).In addition, weight gain during the early season does not account for differences in colony growth during the different phases of early colony development.In our study, pollen from Cornus sanguinea and from oilseed rape were collected in similar amounts, despite oilseed rape, as a mass flowering crop, theoretically being highly attractive for B. terrestris.In addition, Brassicaceae pollen has a lower mean pollen grain volume than Cornus sanguinea pollen.Thus, despite its benefits for colony development, oilseed rape has a rather low contribution to early pollen availability indices.The high effect of oilseed rape on colony growth despite low pollen use indicates that B. terrestris might visit oilseed rape mostly for nectar rather than for its pollen [27].Thus, the true effect of oilseed rape and of plants visited for nectar rather than for pollen might also be underestimated when using the pollen diet as base for floral resource indices.
The high importance of pollen from woody plants in the diet of B. terrestris, especially in the early season (mid-March to end of May) is in line with other studies [27,28].A positive effect of woody floral resources on the development on wild pollinators might be stronger in the beginning of the early season, with the full flowering of Salix, Acer, and Prunus (especially Prunus spinosa, Prunus domestica, and Prunus avium), which are important floral resource plants of B. terrestris in the early season [28].

Conclusions
Classical habitat maps predicted the colony development of B. terrestris better than detailed landscape-scale floral resource maps based on pollen use.This indicates that high amounts of attractive pollen food resources in the landscape alone are not sufficient to ensure a high fitness of bumblebee colonies.Still, the floral resource maps and diet analyses provided information that was not accessible only through classical habitat maps, e.g., that hedgerows play an overriding role in pollen availability to bumblebees in our study region despite their very small cover.More knowledge of floral resource use by B. terrestris inside forest areas is needed to better understand its effects on colony development.Microclimatic conditions in or along forests may help bumblebees to better survive hot and dry weather periods and counteract possible stressors like pesticide exposure in arable land.Overall, our study demonstrates how predictors created by different mapping approaches are needed to complement each other and help to explain their complex relationships between B. terrestris colonies and their development in agricultural landscapes.Using a combination of predictors created by different mapping approaches might help to clearly identify dominant drivers of wild pollinator development and their service in crop pollination in agricultural landscapes.

Supplementary Materials:
The following materials are openly available in figshare (https://figshare.com/) at DOI:10.6084/m9.figshare.13233893,Figure S1: Location of landscapes, Supplement S1: Floral resource mapping, Figure S2: Pollen collection network, Figure S3: Pearson correlation matrix between colony parameter and predictors using different mapping approaches, Table S1: Use and availability of key pollen types included in the study and used for calculating the resource availability index for B. terrestris, Table S2: Proportions of habitat types across landscapes, Table S3: Composition of remaining pollen diet of returning foragers of B. terrestris excluded from index calculation, data presented in this study, and R script used for the analysis.

Figure 1 .
Figure 1.Average (±1 standard error) area of habitat categories (arable, permanent crops, forest edges, herbaceous semi-natural habitats, and woody semi-natural habitats but excluding forests and builtup area) in the studied landscapes (a), their average contribution to pollen availability per landscape (b), in the early season (c), and in the late season (d).Pollen availability is based on the relative area covered by plant species offering a pollen type in the respective habitat type across all landscapes multiplied by the relative collected pollen volume of each type over the whole season.

Figure 1 .
Figure 1.Average (±1 standard error) area of habitat categories (arable, permanent crops, forest edges, herbaceous semi-natural habitats, and woody semi-natural habitats but excluding forests and builtup area) in the studied landscapes (a), their average contribution to pollen availability per landscape (b), in the early season (c), and in the late season (d).Pollen availability is based on the relative area covered by plant species offering a pollen type in the respective habitat type across all landscapes multiplied by the relative collected pollen volume of each type over the whole season.

Figure 1 .
Figure 1.Average (±1 standard error) area of habitat categories (arable, permanent crops, forest edges, herbaceous semi-natural habitats, and woody semi-natural habitats but excluding forests and built-up area) in the studied landscapes (a), their average contribution to pollen availability per landscape (b), in the early season (c), and in the late season (d).Pollen availability is based on the relative area covered by plant species offering a pollen type in the respective habitat type across all landscapes multiplied by the relative collected pollen volume of each type over the whole season.