Novel indices reveal that pollinator exposure to pesticides varies across biological compartments and crop surroundings.

Declines in insect pollinators have been linked to a range of causative factors such as disease, loss of habitats, the quality and availability of food, and exposure to pesticides. Here, we analysed an extensive dataset generated from pesticide screening of foraging insects, pollen-nectar stores/beebread, pollen and ingested nectar across three species of bees collected at 128 European sites set in two types of crop. In this paper, we aimed to (i) derive a new index to summarise key aspects of complex pesticide exposure data and (ii) understand the links between pesticide exposures depicted by the different matrices, bee species and apple orchards versus oilseed rape crops. We found that summary indices were highly correlated with the number of pesticides detected in the related matrix but not with which pesticides were present. Matrices collected from apple orchards generally contained a higher number of pesticides (7.6 pesticides per site) than matrices from sites collected from oilseed rape crops (3.5 pesticides), with fungicides being highly represented in apple crops. A greater number of pesticides were found in pollen-nectar stores/beebread and pollen matrices compared with nectar and bee body matrices. Our results show that for a complete assessment of pollinator pesticide exposure, it is necessary to consider several different exposure routes and multiple species of bees across different agricultural systems.

A blood test to monitor bee health across a European network of agricultural sites of different land-use by MALDI BeeTyping mass spectrometry.

There are substantial concerns about impaired honey bee health and colony losses due to several poorly understood factors. We used MALDI profiling (MALDI BeeTyping®) analysis to investigate how some environmental and management factors under field conditions across Europe affected the honey bee haemolymph peptidome (all peptides in the circulatory fluid), as a profile of molecular markers representing the immune status of Apis mellifera. Honey bees were exposed to a range of environmental stressors in 128 agricultural sites across eight European countries in four biogeographic zones, with each country contributing eight sites each for two different cropping systems: oilseed rape (OSR) and apple (APP). The full haemolymph peptide profiles, including the presence and levels of three key immunity markers, namely the antimicrobial peptides (AMPs) Apidaecin, Abaecin and Defensin-1, allowed the honey bee responses to environmental variables to be discriminated by country, crop type and site. When considering just the AMPs, it was not possible to distinguish between countries by the prevalence of each AMP in the samples. However, it was possible to discriminate between countries on the amounts of the AMPs, with the Swedish samples in particular expressing high amounts of all AMPs. A machine learning model was developed to discriminate the haemolymphs of bees from APP and OSR sites. The model was 90.6 % accurate in identifying the crop type from the samples used to build the model. Overall, MALDI BeeTyping® of bee haemolymph represents a promising and cost-effective "blood test" for simultaneously monitoring dozens of peptide markers affected by environmental stressors at the landscape scale, thus providing policymakers with new diagnostic and regulatory tools for monitoring bee health.

Land-use change in the past 40 years explains shifts in arthropod community traits.

Understanding how anthropogenic activities induce changes in the functional traits of arthropod communities is critical to assessing their ecological consequences. However, we largely lack comprehensive assessments of the long-term impact of global-change drivers on the trait composition of arthropod communities across a large number of species and sites. This knowledge gap critically hampers our ability to predict human-driven impacts on communities and ecosystems. Here, we use a dataset of 1.73 million individuals from 877 species to study how four functionally important traits of carabid beetles and spiders (i.e. body size, duration of activity period, tolerance to drought, and dispersal capacity) have changed at the community level across ~40 years in different types of land use and as a consequence of land use changes (that is, urbanisation and loss of woody vegetation) at the landscape scale in Switzerland. The results show that the mean body size in carabid communities declined in all types of land use, with particularly stronger declines in croplands compared to forests. Furthermore, the length of the activity period and the tolerance to drought of spider communities decreased in most land use types. The average body size of carabid communities in landscapes with increased urbanisation in the last ~40 years tended to decrease. However, the length of the activity period, the tolerance to drought, and the dispersal capacity did not change significantly. Furthermore, urbanisation promoted increases in the average dispersal capacities of spider communities. Additionally, urbanisation favoured spider communities with larger body sizes and longer activity periods. The loss of woody areas at the landscape level was associated with trait shifts to carabid communities with larger body sizes, shorter activity periods, higher drought tolerances and strongly decreased dispersal capacities. Decreases in activity periods and dispersal capacities were also found in spider communities. Our study demonstrates that human-induced changes in land use alter key functional traits of carabid and spider communities in the long term. The detected trait shifts in arthropod communities likely have important consequences for their functional roles in ecosystems.

Distribution of infectious and parasitic agents among three sentinel bee species across European agricultural landscapes.

Infectious and parasitic agents (IPAs) and their associated diseases are major environmental stressors that jeopardize bee health, both alone and in interaction with other stressors. Their impact on pollinator communities can be assessed by studying multiple sentinel bee species. Here, we analysed the field exposure of three sentinel managed bee species (Apis mellifera, Bombus terrestris and Osmia bicornis) to 11 IPAs (six RNA viruses, two bacteria, three microsporidia). The sentinel bees were deployed at 128 sites in eight European countries adjacent to either oilseed rape fields or apple orchards during crop bloom. Adult bees of each species were sampled before their placement and after crop bloom. The IPAs were detected and quantified using a harmonised, high-throughput and semi-automatized qPCR workflow. We describe differences among bee species in IPA profiles (richness, diversity, detection frequencies, loads and their change upon field exposure, and exposure risk), with no clear patterns related to the country or focal crop. Our results suggest that the most frequent IPAs in adult bees are more appropriate for assessing the bees' IPA exposure risk. We also report positive correlations of IPA loads supporting the potential IPA transmission among sentinels, suggesting careful consideration should be taken when introducing managed pollinators in ecologically sensitive environments.

Quantifying potential trade-offs and win-wins between arthropod diversity and yield on cropland under agri-environment schemes-A meta-analysis.

In Europe, agri-environment schemes (AES) are a key instrument to combat the ongoing decline of farmland biodiversity. AES aim is to support biodiversity and maintain ecosystem services, such as pollination or pest control. To what extent AES affect crop yield is still poorly understood. We performed a systematic review, including hierarchical meta-analyses, to investigate potential trade-offs and win-wins between the effectiveness of AES for arthropod diversity and agricultural yield on European croplands. Altogether, we found 26 studies with a total of 125 data points that fulfilled our study inclusion criteria. From each study, we extracted data on biodiversity (arthropod species richness and abundance) and yield for fields with AES management and control fields without AES. The majority of the studies reported significantly higher species richness and abundance of arthropods (especially wild pollinators) in fields with AES (31 % increase), but yields were at the same time significantly lower on fields with AES compared to control fields (21 % decrease). Aside from the opportunity costs, AES that promote out-of-production elements (e.g. wildflower strips), supported biodiversity (29-32 % increase) without significantly compromising yield (2-5 % increase). Farmers can get an even higher yield in these situations than in current conventional agricultural production systems without AES. Thus, our study is useful to identify AES demonstrating benefits for arthropod biodiversity with negligible or relatively low costs regarding yield losses. Further optimization of the design and management of AES is needed to improve their effectiveness in promoting both biodiversity and minimizing crop yield losses.

Distance functions of carabids in crop fields depend on functional traits, crop type and adjacent habitat: a synthesis.

Natural pest and weed regulation are essential for agricultural production, but the spatial distribution of natural enemies within crop fields and its drivers are mostly unknown. Using 28 datasets comprising 1204 study sites across eight Western and Central European countries, we performed a quantitative synthesis of carabid richness, activity densities and functional traits in relation to field edges (i.e. distance functions). We show that distance functions of carabids strongly depend on carabid functional traits, crop type and, to a lesser extent, adjacent non-crop habitats. Richness of both carnivores and granivores, and activity densities of small and granivorous species decreased towards field interiors, whereas the densities of large species increased. We found strong distance decays in maize and vegetables whereas richness and densities remained more stable in cereals, oilseed crops and legumes. We conclude that carabid assemblages in agricultural landscapes are driven by the complex interplay of crop types, adjacent non-crop habitats and further landscape parameters with great potential for targeted agroecological management. In particular, our synthesis indicates that a higher edge-interior ratio can counter the distance decay of carabid richness per field and thus likely benefits natural pest and weed regulation, hence contributing to agricultural sustainability.

Pesticide use negatively affects bumble bees across European landscapes.

Sustainable agriculture requires balancing crop yields with the effects of pesticides on non-target organisms, such as bees and other crop pollinators. Field studies demonstrated that agricultural use of neonicotinoid insecticides can negatively affect wild bee species1,2, leading to restrictions on these compounds3. However, besides neonicotinoids, field-based evidence of the effects of landscape pesticide exposure on wild bees is lacking. Bees encounter many pesticides in agricultural landscapes4-9 and the effects of this landscape exposure on colony growth and development of any bee species remains unknown. Here we show that the many pesticides found in bumble bee-collected pollen are associated with reduced colony performance during crop bloom, especially in simplified landscapes with intensive agricultural practices. Our results from 316 Bombus terrestris colonies at 106 agricultural sites across eight European countries confirm that the regulatory system fails to sufficiently prevent pesticide-related impacts on non-target organisms, even for a eusocial pollinator species in which colony size may buffer against such impacts10,11. These findings support the need for postapproval monitoring of both pesticide exposure and effects to confirm that the regulatory process is sufficiently protective in limiting the collateral environmental damage of agricultural pesticide use.

Diverse pollen nutrition can improve the development of solitary bees but does not mitigate negative pesticide impacts.

Floral resource loss and pesticide exposure are major threats to bees in intensively managed agroecosystems, but interactions among these drivers remain poorly understood. Altered composition and lowered diversity of pollen nutrition may reinforce negative pesticide impacts on bees. Here we investigated the development and survival of the solitary bee Osmia bicornis provisioned with three different pollen types, as well as a mixture of these types representing a higher pollen diversity. We exposed bees of each nutritional treatment to five pesticides at different concentrations in the laboratory. Two field-realistic concentrations of three nicotinic acetylcholine receptor (nAChR) modulating insecticides (thiacloprid, sulfoxaflor and flupyradifurone), as well as of two fungicides (azoxystrobin and tebuconazole) were examined. We further measured the expression of two detoxification genes (CYP9BU1, CYP9BU2) under exposure to thiacloprid across different nutrition treatments as a potential mechanistic pathway driving pesticide-nutrition interactions. We found that more diverse pollen nutrition reduced development time, enhanced pollen efficacy (cocoon weight divided by consumed pollen weight) and pollen consumption, and increased weight of O. bicornis after larval development (cocoon weight). Contrary to fungicides, high field-realistic concentrations of all three insecticides negatively affected O. bicornis by extending development times. Moreover, sulfoxaflor and flupyradifurone also reduced pollen efficacy and cocoon weight, and sulfoxaflor reduced pollen consumption and increased mortality. The expression of detoxification genes differed across pollen nutrition types, but was not enhanced after exposure to thiacloprid. Our findings highlight that lowered diversity of pollen nutrition and high field-realistic exposure to nAChR modulating insecticides negatively affected the development of O. bicornis, but we found no mitigation of negative pesticide impacts through increased pollen diversity. These results have important implications for risk assessment for bee pollinators, indicating that negative effects of nAChR modulating insecticides to developing solitary bees are currently underestimated.

Estimation of the visual contribution to standing balance using virtual reality.

Sensory perturbations are a valuable tool to assess sensory integration mechanisms underlying balance. Implemented as systems-identification approaches, they can be used to quantitatively assess balance deficits and separate underlying causes. However, the experiments require controlled perturbations and sophisticated modeling and optimization techniques. Here we propose and validate a virtual reality implementation of moving visual scene experiments together with model-based interpretations of the results. The approach simplifies the experimental implementation and offers a platform to implement standardized analysis routines. Sway of 14 healthy young subjects wearing a virtual reality head-mounted display was measured. Subjects viewed a virtual room or a screen inside the room, which were both moved during a series of sinusoidal or pseudo-random room or screen tilt sequences recorded on two days. In a between-subject comparison of 10 [Formula: see text] 6 min long pseudo-random sequences, each applied at 5 amplitudes, our results showed no difference to a real-world moving screen experiment from the literature. We used the independent-channel model to interpret our data, which provides a direct estimate of the visual contribution to balance, together with parameters characterizing the dynamics of the feedback system. Reliability estimates of single subject parameters from six repetitions of a 6 [Formula: see text] 20-s pseudo-random sequence showed poor test-retest agreement. Estimated parameters show excellent reliability when averaging across three repetitions within each day and comparing across days (Intra-class correlation; ICC 0.7-0.9 for visual weight, time delay and feedback gain). Sway responses strongly depended on the visual scene, where the high-contrast, abstract screen evoked larger sway as compared to the photo-realistic room. In conclusion, our proposed virtual reality approach allows researchers to reliably assess balance control dynamics including the visual contribution to balance with minimal implementation effort.

Characterization Factors to Assess Land Use Impacts on Pollinator Abundance in Life Cycle Assessment.

While wild pollinators play a key role in global food production, their assessment is currently missing from the most commonly used environmental impact assessment method, Life Cycle Assessment (LCA). This is mainly due to constraints in data availability and compatibility with LCA inventories. To target this gap, relative pollinator abundance estimates were obtained with the use of a Delphi assessment, during which 25 experts, covering 16 nationalities and 45 countries of expertise, provided scores for low, typical, and high expected abundance associated with 24 land use categories. Based on these estimates, this study presents a set of globally generic characterization factors (CFs) that allows translating land use into relative impacts to wild pollinator abundance. The associated uncertainty of the CFs is presented along with an illustrative case to demonstrate the applicability in LCA studies. The CFs based on estimates that reached consensus during the Delphi assessment are recommended as readily applicable and allow key differences among land use types to be distinguished. The resulting CFs are proposed as the first step for incorporating pollinator impacts in LCA studies, exemplifying the use of expert elicitation methods as a useful tool to fill data gaps that constrain the characterization of key environmental impacts.

Nutritional stress exacerbates impact of a novel insecticide on solitary bees' behaviour, reproduction and survival.

Pesticide exposure and food stress are major threats to bees, but their potential synergistic impacts under field-realistic conditions remain poorly understood and are not considered in current pesticide risk assessments. We conducted a semi-field experiment to examine the single and interactive effects of the novel insecticide flupyradifurone (FPF) and nutritional stress on fitness proxies in the solitary bee Osmia bicornis. Individually marked bees were released into flight cages with monocultures of buckwheat, wild mustard or purple tansy, which were assigned to an insecticide treatment (FPF or control) in a crossed design. Nutritional stress, which was high in bees foraging on buckwheat, intermediate on wild mustard and low on purple tansy, modulated the impact of insecticide exposure. Within the first day after application of FPF, mortality of bees feeding on buckwheat was 29 times higher compared with control treatments, while mortality of FPF exposed and control bees was similar in the other two plant species. Moreover, we found negative synergistic impacts of FPF and nutritional stress on offspring production, flight activity, flight duration and flower visitation frequency. These results reveal that environmental policies and risk assessment schemes that ignore interactions among anthropogenic stressors will fail to adequately protect bees and the pollination services they provide.

Different roles of concurring climate and regional land-use changes in past 40 years' insect trends.

Climate and land-use changes are main drivers of insect declines, but their combined effects have not yet been quantified over large spatiotemporal scales. We analysed changes in the distribution (mean occupancy of squares) of 390 insect species (butterflies, grasshoppers, dragonflies), using 1.45 million records from across bioclimatic gradients of Switzerland between 1980 and 2020. We found no overall decline, but strong increases and decreases in the distributions of different species. For species that showed strongest increases (25% quantile), the average proportion of occupied squares increased in 40 years by 0.128 (95% credible interval: 0.123-0.132), which equals an average increase in mean occupancy of 71.3% (95% CI: 67.4-75.1%) relative to their 40-year mean occupancy. For species that showed strongest declines (25% quantile), the average proportion decreased by 0.0660 (95% CI: 0.0613-0.0709), equalling an average decrease in mean occupancy of 58.3% (95% CI: 52.2-64.4%). Decreases were strongest for narrow-ranged, specialised, and cold-adapted species. Short-term distribution changes were associated to both climate changes and regional land-use changes. Moreover, interactive effects between climate and regional land-use changes confirm that the various drivers of global change can have even greater impacts on biodiversity in combination than alone. In contrast, 40-year distribution changes were not clearly related to regional land-use changes, potentially reflecting mixed changes in local land use after 1980. Climate warming however was strongly linked to 40-year changes, indicating its key role in driving insect trends of temperate regions in recent decades.

Do pesticide and pathogen interactions drive wild bee declines?

There is clear evidence for wild insect declines globally. Habitat loss, climate change, pests, pathogens and environmental pollution have all been shown to cause detrimental effects on insects. However, interactive effects between these stressors may be the key to understanding reported declines. Here, we review the literature on pesticide and pathogen interactions for wild bees, identify knowledge gaps, and suggest avenues for future research fostering mitigation of the observed declines. The limited studies available suggest that effects of pesticides most likely override effects of pathogens. Bees feeding on flowers and building sheltered nests, are likely less adapted to toxins compared to other insects, which potential susceptibility is enhanced by the reduced number of genes encoding detoxifying enzymes compared with other insect species. However, to date all 10 studies using a fully-crossed design have been conducted in the laboratory on social bees using Crithidia spp. or Nosema spp., identifying an urgent need to test solitary bees and other pathogens. Similarly, since laboratory studies do not necessarily reflect field conditions, semi-field and field studies are essential if we are to understand these interactions and their potential effects in the real-world. In conclusion, there is a clear need for empirical (semi-)field studies on a range of pesticides, pathogens, and insect species to better understand the pathways and mechanisms underlying their potential interactions, in particular their relevance for insect fitness and population dynamics. Such data are indispensable to drive forward robust modelling of interactive effects in different environmental settings and foster predictive science. This will enable pesticide and pathogen interactions to be put into the context of other stressors more broadly, evaluating their relative importance in driving the observed declines of wild bees and other insects. Ultimately, this will enable the development of more effective mitigation measures to protect bees and the ecosystem services they supply.

Effects of temporal floral resource availability and non-crop habitats on broad bean pollination.

Context: Flowering plants can enhance wild insect populations and their pollination services to crops in agricultural landscapes, especially when they flower before the focal crop. However, characterizing the temporal availability of specific floral resources is a challenge. Objectives: Developing an index for the availability of floral resources at the landscape scale according to the specific use by a pollinator. Investigating whether detailed and temporally-resolved floral resource maps predict pollination success of broad bean better than land cover maps. Methods: We mapped plant species used as pollen source by bumblebees in 24 agricultural landscapes and developed an index of floral resource availability for different times of the flowering season. To measure pollination success, patches of broad bean (Vicia faba), a plant typically pollinated by bumblebees, were exposed in the center of selected landscapes. Results: Higher floral resource availability before bean flowering led to enhanced seed set. Floral resource availability synchronous to broad bean flowering had no effect. Seed set was somewhat better explained by land cover maps than by floral resource availability, increasing with urban area and declining with the cover of arable land. Conclusions: The timing of alternative floral resource availability is important for crop pollination. The higher explanation of pollination success by land cover maps than by floral resource availability indicates that additional factors such as habitat disturbance and nesting sites play a role in pollination. Enhancing non-crop woody plants in agricultural landscapes as pollen sources may ensure higher levels of crop pollination by wild pollinators such as bumblebees. Supplementary information: The online version contains supplementary material available at 10.1007/s10980-022-01448-2.

No evidence for impaired solitary bee fitness following pre-flowering sulfoxaflor application alone or in combination with a common fungicide in a semi-field experiment.

Pesticide exposure is considered a major driver of pollinator decline and the use of neonicotinoid insecticides has been restricted by regulatory authorities due to their risks for pollinators. Impacts of new alternative sulfoximine-based compounds on solitary bees and their potential interactive effects with other commonly applied pesticides in agriculture remain unclear. Here, we conducted a highly replicated full-factorial semi-field experiment with the solitary bee Osmia bicornis, an important pollinator of crops and wild plants in Europe, and Phacelia tanacetifolia as a model crop. We show that spray applications of the insecticide sulfoxaflor (product Closer) and the fungicide azoxystrobin (product Amistar), both alone and combined, had no significant negative impacts on adult female survival or the production, mortality, sex ratio and body size of offspring when sulfoxaflor was applied five days before crop flowering. Our results indicate that for O. bicornis (1) the risk of adverse impacts of sulfoxaflor (Closer) on fitness is small when applied at least five days before crop flowering and (2) that azoxystrobin (Amistar) has a low potential of exacerbating sulfoxaflor effects under field-realistic conditions.

Flowering resources modulate the sensitivity of bumblebees to a common fungicide.

Bees are exposed to various stressors, including pesticides and lack of flowering resources. Despite potential interactions between these stressors, the impacts of pesticides on bees are generally assumed to be consistent across bee-attractive crops, and regulatory risk assessments of pesticides neglect interactions with flowering resources. Furthermore, impacts of fungicides on bees are rarely examined in peer-reviewed studies, although these are often the pesticides that bees are most exposed to. In a full-factorial semi-field experiment with 39 large flight cages, we assessed the single and combined impacts of the globally used azoxystrobin-based fungicide Amistar® and three types of flowering resources (Phacelia, buckwheat, and a floral mix) on Bombus terrestris colonies. Although Amistar is classified as bee-safe, Amistar exposure through Phacelia monocultures reduced adult worker body mass and colony growth (including a 55% decline in workers and an 88% decline in males), while the fungicide had no impact on colonies in buckwheat or the floral mix cages. Furthermore, buckwheat monocultures hampered survival and fecundity irrespective of fungicide exposure. This shows that bumblebees require access to complementary flowering species to gain both fitness and fungicide tolerance and that Amistar impacts are flowering resource-dependent. Our findings call for further research on how different flowering plants affect bees and their pesticide tolerance to improve guidelines for regulatory pesticide risk assessments and inform the choice of plants that are cultivated to safeguard pollinators.

Bee Tracker-an open-source machine learning-based video analysis software for the assessment of nesting and foraging performance of cavity-nesting solitary bees.

The foraging and nesting performance of bees can provide important information on bee health and is of interest for risk and impact assessment of environmental stressors. While radiofrequency identification (RFID) technology is an efficient tool increasingly used for the collection of behavioral data in social bee species such as honeybees, behavioral studies on solitary bees still largely depend on direct observations, which is very time-consuming. Here, we present a novel automated methodological approach of individually and simultaneously tracking and analyzing foraging and nesting behavior of numerous cavity-nesting solitary bees. The approach consists of monitoring nesting units by video recording and automated analysis of videos by machine learning-based software. This Bee Tracker software consists of four trained deep learning networks to detect bees that enter or leave their nest and to recognize individual IDs on the bees' thorax and the IDs of their nests according to their positions in the nesting unit. The software is able to identify each nest of each individual nesting bee, which permits to measure individual-based measures of reproductive success. Moreover, the software quantifies the number of cavities a female enters until it finds its nest as a proxy of nest recognition, and it provides information on the number and duration of foraging trips. By training the software on 8 videos recording 24 nesting females per video, the software achieved a precision of 96% correct measurements of these parameters. The software could be adapted to various experimental setups by training it according to a set of videos. The presented method allows to efficiently collect large amounts of data on cavity-nesting solitary bee species and represents a promising new tool for the monitoring and assessment of behavior and reproductive success under laboratory, semi-field, and field conditions.

Comparing floral resource maps and land cover maps to predict predators and aphid suppression on field bean.

CONTEXT: Predatory insects contribute to the natural control of agricultural pests, but also use plant pollen or nectar as supplementary food resources. Resource maps have been proposed as an alternative to land cover maps for prediction of beneficial insects. OBJECTIVES: We aimed at predicting the abundance of crop pest predating insects and the pest control service they provide with both, detailed flower resource maps and land cover maps. METHODS: We selected 19 landscapes of 500 m radius and mapped them with both approaches. In the centres of the landscapes, aphid predators - hoverflies (Diptera: Syrphidae), ladybeetles (Coleoptera: Coccinellidae) and lacewings (Neuroptera: Chrysopidae) - were surveyed in experimentally established faba bean phytometers (Vicia faba L. Var. Sutton Dwarf) and their control of introduced black bean aphids (Aphis fabae Scop.) was recorded. RESULTS: Landscapes with higher proportions of forest edge as derived from land cover maps supported higher abundance of aphid predators, and high densities of aphid predators reduced aphid infestation on faba bean. Floral resource maps did not significantly predict predator abundance or aphid control services. CONCLUSIONS: Land cover maps allowed to relate landscape composition with predator abundance, showing positive effects of forest edges. Floral resource maps may have failed to better predict predators because other resources such as overwintering sites or alternative prey potentially play a more important role than floral resources. More research is needed to further improve our understanding of resource requirements beyond floral resource estimations and our understanding of their role for aphid predators at the landscape scale. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10980-021-01361-0.

Fungicide and insecticide exposure adversely impacts bumblebees and pollination services under semi-field conditions.

Sulfoximines, the next generation systemic insecticides developed to replace neonicotinoids, have been shown to negatively impact pollinator development and reproduction. However, field-realistic studies on sulfoximines are few and consequences on pollination services unexplored. Moreover, the impacts of other agrochemicals such as fungicides, and their combined effects with insecticides remain poorly investigated. Here, we show in a full factorial semi-field experiment that spray applications of both the product Closer containing the insecticide sulfoxaflor and the product Amistar containing the fungicide azoxystrobin, negatively affected the individual foraging performance of bumblebees (Bombus terrestris). Insecticide exposure further reduced colony growth and size whereas fungicide exposure decreased pollen deposition. We found indications for resource limitation that might have exacerbated pesticide effects on bumblebee colonies. Our work demonstrates that field-realistic exposure to sulfoxaflor can adversely impact bumblebees and that applications before bloom may be insufficient as a mitigation measure to prevent its negative impacts on pollinators. Moreover, fungicide use during bloom could reduce bumblebee foraging performance and pollination services.