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.

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.

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.

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.

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.

Wild insect diversity increases inter-annual stability in global crop pollinator communities.

While an increasing number of studies indicate that the range, diversity and abundance of many wild pollinators has declined, the global area of pollinator-dependent crops has significantly increased over the last few decades. Crop pollination studies to date have mainly focused on either identifying different guilds pollinating various crops, or on factors driving spatial changes and turnover observed in these communities. The mechanisms driving temporal stability for ecosystem functioning and services, however, remain poorly understood. Our study quantifies temporal variability observed in crop pollinators in 21 different crops across multiple years at a global scale. Using data from 43 studies from six continents, we show that (i) higher pollinator diversity confers greater inter-annual stability in pollinator communities, (ii) temporal variation observed in pollinator abundance is primarily driven by the three-most dominant species, and (iii) crops in tropical regions demonstrate higher inter-annual variability in pollinator species richness than crops in temperate regions. We highlight the importance of recognizing wild pollinator diversity in agricultural landscapes to stabilize pollinator persistence across years to protect both biodiversity and crop pollination services. Short-term agricultural management practices aimed at dominant species for stabilizing pollination services need to be considered alongside longer term conservation goals focussed on maintaining and facilitating biodiversity to confer ecological stability.

Opportunities to reduce pollination deficits and address production shortfalls in an important insect-pollinated crop.

Pollinators face multiple pressures and there is evidence of populations in decline. As demand for insect-pollinated crops increases, crop production is threatened by shortfalls in pollination services. Understanding the extent of current yield deficits due to pollination and identifying opportunities to protect or improve crop yield and quality through pollination management is therefore of international importance. To explore the extent of "pollination deficits," where maximum yield is not being achieved due to insufficient pollination, we used an extensive dataset on a globally important crop, apples. We quantified how these deficits vary between orchards and countries and we compared "pollinator dependence" across different apple varieties. We found evidence of pollination deficits and, in some cases, risks of overpollination were even apparent for which fruit quality could be reduced by too much pollination. In almost all regions studied we found some orchards performing significantly better than others in terms of avoiding a pollination deficit and crop yield shortfalls due to suboptimal pollination. This represents an opportunity to improve production through better pollinator and crop management. Our findings also demonstrated that pollinator dependence varies considerably between apple varieties in terms of fruit number and fruit quality. We propose that assessments of pollination service and deficits in crops can be used to quantify supply and demand for pollinators and help to target local management to address deficits although crop variety has a strong influence on the role of pollinators.

Crop pests and predators exhibit inconsistent responses to surrounding landscape composition.

The idea that noncrop habitat enhances pest control and represents a win-win opportunity to conserve biodiversity and bolster yields has emerged as an agroecological paradigm. However, while noncrop habitat in landscapes surrounding farms sometimes benefits pest predators, natural enemy responses remain heterogeneous across studies and effects on pests are inconclusive. The observed heterogeneity in species responses to noncrop habitat may be biological in origin or could result from variation in how habitat and biocontrol are measured. Here, we use a pest-control database encompassing 132 studies and 6,759 sites worldwide to model natural enemy and pest abundances, predation rates, and crop damage as a function of landscape composition. Our results showed that although landscape composition explained significant variation within studies, pest and enemy abundances, predation rates, crop damage, and yields each exhibited different responses across studies, sometimes increasing and sometimes decreasing in landscapes with more noncrop habitat but overall showing no consistent trend. Thus, models that used landscape-composition variables to predict pest-control dynamics demonstrated little potential to explain variation across studies, though prediction did improve when comparing studies with similar crop and landscape features. Overall, our work shows that surrounding noncrop habitat does not consistently improve pest management, meaning habitat conservation may bolster production in some systems and depress yields in others. Future efforts to develop tools that inform farmers when habitat conservation truly represents a win-win would benefit from increased understanding of how landscape effects are modulated by local farm management and the biology of pests and their enemies.

The effectiveness of flower strips and hedgerows on pest control, pollination services and crop yield: a quantitative synthesis.

Floral plantings are promoted to foster ecological intensification of agriculture through provisioning of ecosystem services. However, a comprehensive assessment of the effectiveness of different floral plantings, their characteristics and consequences for crop yield is lacking. Here we quantified the impacts of flower strips and hedgerows on pest control (18 studies) and pollination services (17 studies) in adjacent crops in North America, Europe and New Zealand. Flower strips, but not hedgerows, enhanced pest control services in adjacent fields by 16% on average. However, effects on crop pollination and yield were more variable. Our synthesis identifies several important drivers of variability in effectiveness of plantings: pollination services declined exponentially with distance from plantings, and perennial and older flower strips with higher flowering plant diversity enhanced pollination more effectively. These findings provide promising pathways to optimise floral plantings to more effectively contribute to ecosystem service delivery and ecological intensification of agriculture in the future.

The interplay of landscape composition and configuration: new pathways to manage functional biodiversity and agroecosystem services across Europe.

Managing agricultural landscapes to support biodiversity and ecosystem services is a key aim of a sustainable agriculture. However, how the spatial arrangement of crop fields and other habitats in landscapes impacts arthropods and their functions is poorly known. Synthesising data from 49 studies (1515 landscapes) across Europe, we examined effects of landscape composition (% habitats) and configuration (edge density) on arthropods in fields and their margins, pest control, pollination and yields. Configuration effects interacted with the proportions of crop and non-crop habitats, and species' dietary, dispersal and overwintering traits led to contrasting responses to landscape variables. Overall, however, in landscapes with high edge density, 70% of pollinator and 44% of natural enemy species reached highest abundances and pollination and pest control improved 1.7- and 1.4-fold respectively. Arable-dominated landscapes with high edge densities achieved high yields. This suggests that enhancing edge density in European agroecosystems can promote functional biodiversity and yield-enhancing ecosystem services.

The worldwide importance of honey bees as pollinators in natural habitats.

The western honey bee (Apis mellifera) is the most frequent floral visitor of crops worldwide, but quantitative knowledge of its role as a pollinator outside of managed habitats is largely lacking. Here we use a global dataset of 80 published plant-pollinator interaction networks as well as pollinator effectiveness measures from 34 plant species to assess the importance of A. mellifera in natural habitats. Apis mellifera is the most frequent floral visitor in natural habitats worldwide, averaging 13% of floral visits across all networks (range 0-85%), with 5% of plant species recorded as being exclusively visited by A. mellifera For 33% of the networks and 49% of plant species, however, A. mellifera visitation was never observed, illustrating that many flowering plant taxa and assemblages remain dependent on non-A. mellifera visitors for pollination. Apis mellifera visitation was higher in warmer, less variable climates and on mainland rather than island sites, but did not differ between its native and introduced ranges. With respect to single-visit pollination effectiveness, A. mellifera did not differ from the average non-A. mellifera floral visitor, though it was generally less effective than the most effective non-A. mellifera visitor. Our results argue for a deeper understanding of how A. mellifera, and potential future changes in its range and abundance, shape the ecology, evolution, and conservation of plants, pollinators, and their interactions in natural habitats.

Synergistic interactions of ecosystem services: florivorous pest control boosts crop yield increase through insect pollination.

Insect pollination and pest control are pivotal functions sustaining global food production. However, they have mostly been studied in isolation and how they interactively shape crop yield remains largely unexplored. Using controlled field experiments, we found strong synergistic effects of insect pollination and simulated pest control on yield quantity and quality. Their joint effect increased yield by 23%, with synergistic effects contributing 10%, while their single contributions were 7% and 6%, respectively. The potential economic benefit for a farmer from the synergistic effects (12%) was 1.8 times greater than their individual contributions (7% each). We show that the principal underlying mechanism was a pronounced pest-induced reduction in flower lifetime, resulting in a strong reduction in the number of pollinator visits a flower receives during its lifetime. Our findings highlight the importance of non-additive interactions among ecosystem services (ES) when valuating, mapping or predicting them and reveal fundamental implications for ecosystem management and policy aimed at maximizing ES for sustainable agriculture.

High effectiveness of tailored flower strips in reducing pests and crop plant damage.

Providing key resources to animals may enhance both their biodiversity and the ecosystem services they provide. We examined the performance of annual flower strips targeted at the promotion of natural pest control in winter wheat. Flower strips were experimentally sown along 10 winter wheat fields across a gradient of landscape complexity (i.e. proportion non-crop area within 750 m around focal fields) and compared with 15 fields with wheat control strips. We found strong reductions in cereal leaf beetle(CLB) density (larvae: 40%; adults of the second generation: 53%) and plant damage caused by CLB (61%) in fields with flower strips compared with control fields. Natural enemies of CLB were strongly increased in flower strips and in part also in adjacent wheat fields. Flower strip effects on natural enemies, pests and crop damage were largely independent of landscape complexity(8-75% non-crop area). Our study demonstrates a high effectiveness of annual flower strips in promoting pest control, reducing CLB pest levels below the economic threshold. Hence, the studied flower strip offers a viable alternative to insecticides. This highlights the high potential of tailored agri-environment schemes to contribute to ecological intensification and may encourage more farmers to adopt such schemes.

Consequences of plant invasions on compartmentalization and species' roles in plant-pollinator networks.

Compartmentalization-the organization of ecological interaction networks into subsets of species that do not interact with other subsets (true compartments) or interact more frequently among themselves than with other species (modules)-has been identified as a key property for the functioning, stability and evolution of ecological communities. Invasions by entomophilous invasive plants may profoundly alter the way interaction networks are compartmentalized. We analysed a comprehensive dataset of 40 paired plant-pollinator networks (invaded versus uninvaded) to test this hypothesis. We show that invasive plants have higher generalization levels with respect to their pollinators than natives. The consequences for network topology are that-rather than displacing native species from the network-plant invaders attracting pollinators into invaded modules tend to play new important topological roles (i.e. network hubs, module hubs and connectors) and cause role shifts in native species, creating larger modules that are more connected among each other. While the number of true compartments was lower in invaded compared with uninvaded networks, the effect of invasion on modularity was contingent on the study system. Interestingly, the generalization level of the invasive plants partially explains this pattern, with more generalized invaders contributing to a lower modularity. Our findings indicate that the altered interaction structure of invaded networks makes them more robust against simulated random secondary species extinctions, but more vulnerable when the typically highly connected invasive plants go extinct first. The consequences and pathways by which biological invasions alter the interaction structure of plant-pollinator communities highlighted in this study may have important dynamical and functional implications, for example, by influencing multi-species reciprocal selection regimes and coevolutionary processes.

Does the aboveground herbivore assemblage influence soil bacterial community composition and richness in subalpine grasslands?

Grassland ecosystems support large communities of aboveground herbivores that are known to directly and indirectly affect belowground properties such as the microbial community composition, richness, or biomass. Even though multiple species of functionally different herbivores coexist in grassland ecosystems, most studies have only considered the impact of a single group, i.e., large ungulates (mostly domestic livestock) on microbial communities. Thus, we investigated how the exclusion of four groups of functionally different herbivores affects bacterial community composition, richness, and biomass in two vegetation types with different grazing histories. We progressively excluded large, medium, and small mammals as well as invertebrate herbivores using exclosures at 18 subalpine grassland sites (9 per vegetation type). We assessed the bacterial community composition using terminal restriction fragment length polymorphism (T-RFLP) at each site and exclosure type during three consecutive growing seasons (2009-2011) for rhizosphere and mineral soil separately. In addition, we determined microbial biomass carbon (MBC), root biomass, plant carbon:nitrogen ratio, soil temperature, and soil moisture. Even though several of these variables were affected by herbivore exclusion and vegetation type, against our expectations, bacterial community composition, richness, or MBC were not. Yet, bacterial communities strongly differed between the three growing seasons as well as to some extent between our study sites. Thus, our study indicates that the spatiotemporal variability in soil microclimate has much stronger effects on the soil bacterial communities than the grazing regime or the composition of the vegetation in this high-elevation ecosystem.

Species traits, landscape quality and floral resource overlap with honeybees determine virus transmission in plant-pollinator networks.

Emerging infectious diseases pose a threat to pollinators. Virus transmission among pollinators via flowers may be reinforced by anthropogenic land-use change and concomitant alteration of plant-pollinator interactions. Here, we examine how species' traits and roles in flower-visitation networks and landscape-scale factors drive key honeybee viruses-black queen cell virus (BQCV) and deformed wing virus-in 19 wild bee and hoverfly species, across 12 landscapes varying in pollinator-friendly (flower-rich) habitat. Viral loads were on average more than ten times higher in managed honeybees than in wild pollinators. Viral loads in wild pollinators were higher when floral resource use overlapped with honeybees, suggesting these as reservoir hosts, and increased with pollinator abundance and viral loads in honeybees. Viral prevalence decreased with the amount of pollinator-friendly habitat in a landscape, which was partly driven by reduced floral resource overlap with honeybees. Black queen cell virus loads decreased with a wild pollinator's centrality in the network and the proportion of visited dish-shaped flowers. Our findings highlight the complex interplay of resource overlap with honeybees, species traits and roles in flower-visitation networks and flower-rich pollinator habitat shaping virus transmission.

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.

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.