Modeling the status, trends, and impacts of wild bee abundance in the United States

Edited by May R. Berenbaum, University of Illinois at Urbana–Champaign, Urbana, IL, and approved November 20, 2015 (received for review September 4, 2015)
December 22, 2015
113 (1) 140-145

Significance

In 2014, a presidential memorandum called for an assessment of the nation’s pollinators, in response to growing awareness of their economic importance and recent declines. We assess, for the first time to our knowledge, the status and trends of wild bee abundance and their potential impacts on pollination services across the United States. We develop national maps of wild bee abundance, report land-use–driven changes over time, and relate them to trends in agricultural demand for pollination. We estimate uncertainty in the findings, so future research can target the least-understood regions and topics. Our findings can also help focus conservation efforts where declines in bee abundance are most certain, especially where agricultural demand for pollination services is growing.

Abstract

Wild bees are highly valuable pollinators. Along with managed honey bees, they provide a critical ecosystem service by ensuring stable pollination to agriculture and wild plant communities. Increasing concern about the welfare of both wild and managed pollinators, however, has prompted recent calls for national evaluation and action. Here, for the first time to our knowledge, we assess the status and trends of wild bees and their potential impacts on pollination services across the coterminous United States. We use a spatial habitat model, national land-cover data, and carefully quantified expert knowledge to estimate wild bee abundance and associated uncertainty. Between 2008 and 2013, modeled bee abundance declined across 23% of US land area. This decline was generally associated with conversion of natural habitats to row crops. We identify 139 counties where low bee abundances correspond to large areas of pollinator-dependent crops. These areas of mismatch between supply (wild bee abundance) and demand (cultivated area) for pollination comprise 39% of the pollinator-dependent crop area in the United States. Further, we find that the crops most highly dependent on pollinators tend to experience more severe mismatches between declining supply and increasing demand. These trends, should they continue, may increase costs for US farmers and may even destabilize crop production over time. National assessments such as this can help focus both scientific and political efforts to understand and sustain wild bees. As new information becomes available, repeated assessments can update findings, revise priorities, and track progress toward sustainable management of our nation’s pollinators.
Bees and other flower-visiting animals provide essential pollination services to many US crops (1) and to wild plant species (2). Bees contributed an estimated 11% of the nation’s agricultural gross domestic product in 2009 (3), equal to $14.6 billion per year (4). Of this, at least 20% ($3.07 billion) is provided by wild pollinators that depend on suitable land for nesting and foraging (5). As the consumption of specialty fruit and vegetable crops has grown (6), the demand for pollination services has increased. However, the supply of managed honey bees (Apis mellifera L.) has not kept pace (7), due to management challenges and colony losses over the last decade (8). There is growing evidence that wild, unmanaged bees can provide effective pollination services where sufficient habitat exists to support their populations (9, 10). They can also contribute to the long-term stability of crop pollination, thereby reducing the risk of pollination deficits from variable supply or activity of honey bees (11, 12). As a result, wild pollinators should be integrated into crop pollination management plans as a supplement or alternative to managed bees (13).
Despite the agricultural importance of wild bees, there is increasing evidence that multiple species are declining in range or abundance. Some of the most important crop pollinators, such as bumble bees (Bombus spp.), have declined over past decades in the United States (1416). Among the numerous threats to wild bees, including pesticide use, climate change, and disease (17), habitat loss seems to contribute to most observed declines (18). Indeed, a National Research Council committee on the status of pollinators in North America reported that conserving and improving habitats for wild bees is important for ensuring continued pollination services and food security (19).
Recognizing both the growing need for pollination services and increasing threats to wild bees, a recent presidential memorandum called for a national assessment of the status of wild pollinators and available habitat in the United States (20). The resulting report sets a goal of 7 million acres of land for pollinators over the next 5 y (21). However, there has been no assessment at the national level of the current status of native pollinator habitat, where and at what rate this habitat is being degraded, and the impact of these changes on bee populations and the pollination services they provide.
A national assessment is challenging because plant–pollinator interactions and dynamics occur at relatively fine spatial scales. Wild bee populations are largely determined by the spatial distribution of habitat resources within their foraging range (2224), and this varies from ∼100–2,000 m (25, 26). Accordingly, most of our understanding of native bee populations is at the scale of landscapes and local sites. Several field-based assessments of habitat resources for native bee species have been developed at landscape scales (23, 2729). However, the required cost and time to scale this type of field assessment to cover all habitat types and bee species nationwide is logistically challenging and prohibitively expensive.
When field observations are lacking, careful use of expert-derived data has been shown to provide informative estimates that enable habitat assessments (30, 31), including studies on pollination (32, 33). Use of expert opinion may therefore be an efficient path to an initial nationwide assessment of pollinator habitat and abundance in the absence of consistent data across different land categories. Such an approach must include careful treatment of uncertainty that may arise from differences in expertise among regions, authorities, taxa, and so on (3436). Indeed, a robust analysis of uncertainty, and its implications for assessment findings, is a useful result in itself. It can help orient research toward addressing the most important gaps in our national knowledge of wild bees and their importance for crop pollination.
Here, we use a published model of bee abundance (32) and expert knowledge to assess the status, trends, and impact of wild bee abundance and associated uncertainties across the coterminous United States. The spatially explicit model predicts a relative index of wild bee abundance (hereafter, bee abundance) based on local nesting resources and the quality of surrounding forage (32). We parameterize the model with expert-derived estimates of nesting and forage quality for each of the main land-use types in each of the major ecoregions to construct a probability distribution for each parameter that captures estimates by multiple experts and their uncertainty. We first validate model predictions with bee collections and observations from a variety of landscape settings. We then map bee abundance, its uncertainty, and the agricultural demand for pollination across the United States to address the following questions: (i) What are the current status and trends of wild bee abundance across the coterminous United States? (ii) What land-use changes have driven these trends over a 5-y period (2008–2013)? (iii) Which regions and crops experience relatively low bee abundance compared with crop pollination demands? (iv) How does uncertainty in our knowledge affect these predictions? Responses to these questions will inform future research efforts and policy decisions to conserve native bees at the national level and can help guide a coordinated and ongoing nationwide assessment of wild bees.

Results

Bee Abundance.

Our model predicts generally high abundances of wild bees in areas rich in resources such as chaparral and desert shrublands, intermediate abundances in temperate forest and grassland/rangelands, and lower abundances in most agricultural areas (Fig. 1A). Patterns of wild bee abundance and expert uncertainty seem correlated (Fig. 1 A and B). In fact, whereas most areas with low bee abundance also present low uncertainty, only 5% of areas with high bee abundance have low uncertainty. This suggests that experts are more individually or collectively certain about uniformly poor bee habitats (e.g., corn fields) than they are about higher-quality habitats (e.g., shrublands), which can vary in quality over time and space (Discussion).
Fig. 1.
Maps of status, uncertainty, trend, and impacts of wild bees across the coterminous United States. (A) Status of wild bee abundance (relative index) for 2013. (B) Uncertainty (SD estimate) of wild bee abundance index for 2013. (C) Trends in wild bee abundance and its uncertainty (the likelihood of changes: pseudo-t values) between 2008 and 2013. (D) Status of supply of wild bees (model-predicted abundance from A) and demand for pollination services (summed area of animal-pollinated crops, weighted by their pollinator dependence) at a county scale for 2013. Counties with less than 1,000 ha of pollinator-dependency weighted crop area were left white. (E) Uncertainty in the supply of wild bees in 2013 for the counties identified as supply/demand mismatches in D. (F) Trend of supply and demand between 2008 and 2013 (zones I and II indicate high and low likelihood of decreases in supply, respectively).
Between 2008 and 2013, wild bee abundance was consistent in 67% of the US land area (−0.01 < index change < 0.01 in Fig. 1C). However, our model indicates decreases in 23% of the United States (index change < −0.01), and these decreases were highly likely in 9% of the United States (likelihood index ≤ −0.2 in Fig. 1C; Methods). Most of the areas of likely decrease occurred in agricultural regions of Midwestern and Great Plains states and in the Mississippi river valley. Eleven states [Minnesota, Texas (TX), Wisconsin (WI), South Dakota (SD), North Dakota (ND), Illinois, Missouri, Nebraska, Oklahoma, Kansas, and Louisiana] collectively accounted for 60% of the areas of predicted decrease in wild bee abundance. Over the 5-y period in these states, corn and grain cropland increased 200% and 100%, respectively, and mostly replaced grasslands and pasture (Fig. 2A and Fig. S1A). Bee abundance increased in 10% of the United States (index change > 0.01) and the increase was highly likely in 3% of the country (likelihood index ≥ 0.2 in Fig. 1C). Areas of likely increase in bee abundance were found in northern ND, eastern Washington (WA) and Pennsylvania (PA), southern Montana, parts of several states in the Great Plains, and in southeastern coastal areas (Fig. 1C). In these areas, grasslands, pastures, and corn/soy fields were converted to higher-quality habitat, such as shrublands or fallow crop fields (Fig. 2B and Fig. S1B).
Fig. 2.
Changes in land-use/cover corresponding to predicted changes in wild bee abundance. Bars represent land cover in pixels where decreases (A) and increases (B) of wild bee abundance are highly likely between 2008 and 2013 (i.e., bee abundance changes <−0.01 or >0.01 and the likelihood of changes ≤−0.2 or ≥0.2 in Fig 1C, respectively).
Fig. S1.
Land-cover changes that caused the largest predicted change in bee abundance from 2008 and 2013. (A) Top five land covers in 2013 where bee abundance index most likely decreased since 2008. (B) Top five land covers in 2013 where bee abundance index most likely increased since 2008.

Pollination Supply and Demand.

Bee abundance maps (Fig. 1A) can be interpreted as the potential “supply” of pollination services from wild bees. To compare this measure of supply to potential agricultural demand, we calculated the area of pollinator-dependent crops, weighted by each crop’s degree of pollinator dependence, for each US county in 2013 (Methods). By comparing the two maps, we identified counties with relatively high supply of wild bees and relatively low demand (Fig. 1D, light blue) and, conversely, where high demand occurs in counties with relatively low supply (Fig. 1D, purple). We identified 139 counties (which together comprise 39% of pollinator-dependent crop area) where high demand and low supply coincide (Fig. 1D, yellow outline) and 39 counties where this difference was particularly extreme (Fig. 1D, red outline). All of the 139 counties with a pollinator disparity had relatively low uncertainty for 2013 bee abundance (Fig. 1E), which indicates that there is high confidence in this mismatch. These counties tend to contain either a significant percentage of area that consists of highly pollinator-dependent crops [e.g., almonds, blueberries, and apples in California (CA), Oregon, and WA, respectively] or large amount of less-dependent crops (e.g., soybeans and canola in Midwestern states, cotton in northwest TX and the Mississippi Valley).
To examine changes in the relationship between wild bee supply and pollination demand, we combined the two trend maps (Methods). We found that 106 counties have simultaneously experienced increases in demand for pollination services and decreases in wild bee abundance (Fig. 1F, upper left quadrant). This represents 54% of the 195 counties that have experienced substantial changes in pollination demand (>500 ha of change). In 27 of these counties, declines in supply were highly likely (zone I in Fig. 1F legend), whereas in the remaining 79 counties declines were less certain (zone II in Fig. 1F legend). In counties of West Coast states and Michigan, increases in demand were mostly driven by increases in specialty crops such as almonds, cherries, blueberries, apples, watermelons, and squash. In contrast, demand increases in the Great Plains and Mississippi Valley were driven by increases in crops, such as sunflower, canola, soybeans, and cotton, with moderate to low pollinator dependency.
Trends in our measures of supply and demand vary widely among individual crops (Fig. 3). Most crops that require animal pollination have expanded in area (thus demand) between 2008 and 2013, whereas the predicted supply of wild bees in many of these cropped areas has declined. Specialty crops, such as pumpkins, blueberries, peaches, apples, and watermelons, are among the crops that present the strongest mismatch between changes in supply and demand. Others, such as canola, have experienced increases in both supply and demand. Of particular concern for future abilities to meet pollination demands, crops that are most dependent on pollinators (symbols in Fig. 3) tend to have experienced simultaneous declines in supply and increases in demand.
Fig. 3.
Nationwide changes in wild bee abundance and cultivated area for pollinator-dependent crops between 2008 and 2013. Symbols represent pollinator dependence for each crop reported by Klein et al. (49).

Discussion

Our study is the first to our knowledge to map the status and trends of wild bees and their potential impacts on pollination services across the coterminous United States. By combining a spatial model with expert knowledge, we find highly heterogeneous patterns of both predicted abundance of wild bees and our uncertainty regarding those predictions. We also identify counties and crops of potential concern, where declines in wild bee abundance oppose increased need for crop pollination. These analyses form an important step toward a nationwide understanding of the status of wild pollinators. They can also help focus attention and future research toward regions of high uncertainty and to direct management efforts to areas of major concern.
Our mapped index of bee abundance (Fig. 1A) clearly shows that areas of intense agriculture (e.g., the Midwest Corn Belt and California’s Central Valley) are among the lowest in predicted wild bee abundance. Our predictions are also relatively certain in these areas (Fig. 1B). This reflects consensus among experts about the low suitability of intensively managed agricultural land for wild bees and is supported by an abundance of previous research on the negative effects of intensive agriculture on bee populations (37, 38). Recent trends (Fig. 1C) also correspond to increasing agricultural land use over time. Areas of bee abundance where declines are most certain tend to have experienced additional conversion of natural land covers to crops, especially corn (Fig. 2A). These results reinforce recent evidence that increased demand for corn in biofuel production has intensified threats to natural habitats in corn-growing regions (39). For example, a recent land-use simulation found that expansion of annual biofuel crops could reduce pollinator abundance and diversity at the state level (40). In areas where major land-use changes have gone in the opposite direction, however, bee abundance has tended to increase (Fig. 2B). These changes may represent detectable effects from the US Department of Agriculture Conservation Reserve Program, which compensates farmers for retiring marginal lands (41). Given the clear patterns in Fig. 1 AC, supported by other studies at finer spatial scales, this initial assessment can help set management priorities (e.g., habitat restoration or enhancement) to maintain populations of wild bees and other wildlife amid agricultural intensification (42, 43).
We put estimates of relative wild bee supply in the context of nationwide demand for pollination services, by comparing predicted bee abundance (Fig. 1A) to county-level information on crops. A total of 139 counties (Fig. 1D) contain almost half of pollinator-dependent crop area but support relatively low wild bee abundance. In these counties, there seems to be a significant mismatch between the supply of wild bees and demand for pollination services. Because our estimates are relative indices, they do not permit absolute comparisons of supply and demand that would determine whether pollinator abundances are adequate to pollinate crops fully. A more robust approach to locate regions of mismatch, therefore, is to identify counties in which supply and demand are changing in opposite directions (Fig. 1F). This comparison of trends pinpoints many of the same counties as Fig. 1D, and adds others. In these counties, regardless of whether demand for pollination services has already overtaken the ability of wild bees to supply them, recent trends indicate that the risk is growing over time (6). Growers of crops dependent on bees for pollination will need to depend more heavily on managed honey bees to supply pollination in the absence of abundant wild bee populations. We predict increasing demand (and rental fees) over time for honey bees in those regions highlighted in Fig. 1F, but a test of that prediction is beyond the scope of this paper. We also suggest that efforts to manage pollinator habitats, monitor bee populations, and evaluate pollen limitation in crops are most important in these regions.
The opposing trends of crop expansion and wild bee abundance may also be causally linked. Crop expansion probably contributed to the declining quality of bee habitats between 2008 and 2013; indeed, we find a negative correlation between changes in crop demand and bee abundance across all US counties (P < 0.01, Fig. S2). Studies from northern Europe have shown that mass-flowering crops can enhance wild bee abundance in surrounding landscapes (24, 44), but our analyses indicate the opposite relationship (perhaps because North America has larger-scale mass flowering crops) and emphasize the need for more careful assessment of North American systems.
Fig. S2.
Relationship between changes in dependency-weighted crop areas and wild bee abundance changes during 2008–2013.
Analysis of individual crops provides another perspective on potential mismatches between US wild bee supply and demand (Fig. 3). Crops that have decreasing wild bee abundance and increasing cultivated area (upper left quadrant of Fig. 3) tend to be those that are more dependent on bee-mediated pollination (symbols in Fig. 3). Pollination supply and demand are therefore mismatched for precisely the crops that most require pollination. Variability in US crop yields has been found to increase with greater dependence on pollinators (45), so these trends, if they continue, may destabilize crop production over time. To maintain stability in yields, farmers may need to maintain habitats for wild bees on and around their farms (46) or invest more heavily in managed pollinators.
We consider our estimates of uncertainty to be as informative as the bee abundance predictions themselves. All assessments involve uncertainty, but few report this crucial information with sufficient clarity and rigor (34). We are encouraged to note that our model validation supports the uncertainty estimates; expert-derived parameters improved model fit to a greater degree in areas where experts reported more certainty (Fig. S3). Quantifying uncertainty allows us to make initial predictions about the status and trends of pollinator abundance using uneven and incomplete information. It also helps identify regions where additional studies will most effectively improve our estimates and strengthen the national assessment over time. Highly uncertain regions are also those where the precautionary principle would be appropriate in land management strategies to prevent pollinator loss. In practice, uncertainty in our model can increase for three reasons: First, experts may not be certain about the resource quality of a particular land-cover type (e.g., idle cropland and woody wetland); next, individual experts are certain but disagree about the quality of resources available (e.g., developed open space or evergreen forest); and finally, experts acknowledge that a land-cover type is heterogeneous in its resource quality (e.g., grassland, deciduous and mixed forests, and developed open space). In our case, experts were less certain about the quality of nesting resources than of floral resources; this suggests a need to increase effort to understand the nesting biology of wild bees (29). Experts were also more certain about the quality of crops than of noncrop land covers (Fig. S4); this could reflect relative expertise among experts or greater spatial and temporal heterogeneity of natural land covers.
Fig. S3.
Fit between model predictions and field data for observed wild bee abundance and observed bumble bee abundance. Comparison of model fit for expert-informed and noninformative probability distributions for wild bee abundance data (A) and bumble bee abundance data (B). For bumble bees, we included expert-informed probability distribution of cavity nesting and summer floral resource. Bar plots represent the mean and SE of model fit (t value) resulting from 1,000 MC draws. Black circles depict fit when using the mean parameter value from experts (instead of random draws from expert-informed probability distributions). Additional circles represent fit with mean values, but separately for sites with low (red circles) and high (blue circles) model uncertainty. Sites were split into relatively low and high uncertainty group for wild bee abundance data (C) and bumble bee abundance data (D). The dashed line indicates split data as half for each uncertainty group. The relationship between bee abundance index (mean of 1,000 MC simulations) and observed abundance of wild bees (E) and bumble bees (F). Black line indicates the overall fixed effects. Two different colors (blue and red) indicate the model prediction for low and high uncertainty groups.
Fig. S4.
The relationship between the mean resource suitability estimates and uncertainty (SD) from expert-derived probability distributions of nesting (A) and floral (B) resources for wild bees for each of the land-cover categories.
Although our approach carefully captured expert uncertainty, three other sources of uncertainty arise from the data themselves. First, the Cropland Data Layer (CDL), like all land-cover and land-use data, contains classification errors (47), which contribute to the uncertainty in our estimates. For example, apparent land-use conversion from deciduous forest into woody wetlands contributed to predicted declines in bee abundance between 2008 and 2013, especially in Minnesota (Fig. 1C). Conversely, apparent conversion from grasslands into shrubs was the major driver in areas of increased pollinator abundances (Fig. 1C). Both changes, however, are partly a result of inconsistent classification, which led to apparent changes when none occurred. In addition, urban gardens could support a significant abundance of wild pollinators, but the CDL does not capture these specific features within “developed” categories (Table S1). Despite such inaccuracies, the CDL is the only available national coverage of land-uses/covers in agricultural as well as natural areas (48). Second, for our measures of pollination demand (Fig. 1D), for each crop we rely on Klein et al. (49) for estimates of pollinator dependence (Table S2). These estimates consist of simple percentages of yield and have been widely used in studies of pollination services (50, 51). They also contain some uncertainty, however, because each percentage represents the midpoint of a range reported originally in Klein et al. (49), whereas dependencies vary among crop varieties, climates, field settings, and cultivation practices. Because we focus on analyses of relative demand among crops and counties, our findings are likely robust to this uncertainty. Finally, we elicited expert parameters on nesting resources for different guilds and for floral resources at different seasons. However, we combined these estimates to produce a single probability distribution for each habitat type, which increased the uncertainty of our estimates (i.e., the SD of our probability distributions). In the future, more detailed assessments could integrate information on bee communities, nesting habits, and flight seasons to develop more refined probability distributions for each pixel. Indeed, our model predicted bumble bee abundances more accurately when we used parameters relevant to this genus (i.e., cavity-nesting species and summer floral resources), compared with averaged parameters (Fig. S3B). Although we have focused on bees, other taxa can be important crop pollinators (52). For simplicity in this initial nationwide assessment, we have also pooled all bee species into an overall abundance index, but bee taxa clearly vary in their importance as crop pollinators and their response to land use (53). Future work should distinguish pollinator taxa or guilds to model the trends and importance of each separately.
Table S1.
Reclassified crop categories and noncrop categories of NASS–CDL
No.Reclassified categoriesOriginal categories in CDL
Crop categories 
1AlfalfaAlfalfa
2CornCorn, Pop or Orn Corn, Sweet Corn, Sorghum
3BeanSoybeansL, Dry BeansL, PeanutL, Vetch
4BerriesBlueberriesG, CaneberriesG, CranberriesG
5StrawberriesStrawberriesM
6CitrusCitrusL, OrangesL
7CottonCottonL
8CucurbitsGourdsE, PumpkinsE, SquashE
9GrainBarley, Durum Wheat, Hops, Lentils, Millet, Oats, Other Small Grains, Rice, Rye, Speltz, Spring Wheat, Triticale, Winter Wheat, Sugarcane
10BuckwheatBuckwheatG
11GrassPasture/Grass, Sod/Grass Seed, Switchgrass
12GrapesGrapes
13HerbsHerbs, Mint
14MelonsCantaloupesE, Honeydew MelonsE, Dbl Crop Lettuce/CantaloupeM, CucumbersG
15WatermelonsWatermelonsE
16OilseedCamelinaL, CanolaL, FlaxseedL, Mustard, Rape SeedL
17FlowersSafflowerL, SunflowerM
18WildflowersClover/Wildflowers
19OrchardAlmondsG, ApplesG, ApricotsG, CherriesG, NectarinesG, PeachesG, PearsG, PlumsG, PomegranatesM, PrunesG
20Root VegetablesCarrots, Garlic, Onions, TurnipsG, Sugarbeets, Radishes
21SolanumsEggplantsM, Peppers, Potatoes, Sweet Potatoes, TomatoesL
22VegetablesCabbage, Cauliflower, Celery, Greens, Lettuce, Broccoli, Chick Peas, PeasL
23TobaccoTobacco
24Other CropsOther Crops*
25Vegetables and FruitsMiscellaneous Vegetables and FruitsL
26NutsPistachios, Walnuts, Pecans
27AsparagusAsparagus
28OlivesOlives
29Tree CropsOther Tree Crops
30Christmas TreesChristmas Trees
31Idle CroplandFallow/Idle Cropland
32Double (Dbl) CropDbl Crop Barley/Corn, Dbl Crop Barley/Sorghum, Dbl Crop Barley/SoybeansM, Dbl Crop Corn/SoybeansM, Dbl Crop Durum Wheat/Sorghum, Dbl Crop Lettuce/Barley, Dbl Crop Lettuce/Cotton, Dbl Crop Lettuce/Durum Wheat, Dbl Crop Oats/Corn, Dbl Crop Soybeans/CottonM, Dbl Crop Soybeans/OatsM, Dbl Crop Winter Wheat/Corn, Dbl Crop, Winter Wheat/Cotton, Dbl Crop Winter Wheat/Sorghum, Dbl Crop Winter Wheat/Soybeans
Noncrop categoriesCDL detailed categories description§
33Developed/Open SpaceDeveloped/Open Space: Areas with a mixture of some constructed materials, but mostly vegetation in the form of lawn grasses. Impervious surfaces account for less than 20% of total cover. These areas most commonly include large-lot single-family housing units, parks, golf courses, and vegetation planted in developed settings for recreation, erosion control, or aesthetic purposes.
34Developed/Low IntensityDeveloped/Low Intensity: Areas with a mixture of constructed materials and vegetation. Impervious surfaces account for 20–49% of total cover. These areas most commonly include single-family housing units.
35Developed/Med IntensityDeveloped/Med Intensity: Areas with a mixture of constructed materials and vegetation. Impervious surfaces account for 50–79% of the total cover. These areas most commonly include single-family housing units.
36Developed/High IntensityDeveloped/High Intensity: Highly developed areas where people reside or work in high numbers. Examples include apartment complexes, row houses, and commercial/industrial. Impervious surfaces account for 80–100% of the total cover.
37BarrenBarren: Areas of bedrock, desert pavement, scarps, talus, slides, volcanic material, glacial debris, sand dunes, strip mines, gravel pits and other accumulations of earthen material. Generally, vegetation accounts for less than 15% of total cover.
38Deciduous ForestDeciduous Forest: Areas dominated by trees generally greater than 5 m tall, and greater than 20% of total vegetation cover. More than 75% of the tree species shed foliage simultaneously in response to seasonal change
39Evergreen ForestEvergreen Forest: Areas dominated by trees generally greater than 5 m tall, and greater than 20% of total vegetation cover. More than 75% of the tree species maintain their leaves all year. Canopy is never without green foliage.
40Mixed ForestMixed Forest: Areas dominated by trees generally greater than 5 m tall, and greater than 20% of total vegetation cover. Neither deciduous nor evergreen species are greater than 75% of total tree cover.
41ShrublandShrubland: Areas dominated by shrubs less than 5 m tall with shrub canopy typically greater than 20% of total vegetation. This class includes true shrubs, young trees in an early successional stage or trees stunted from environmental conditions.
42Grass/PastureGrassland/Herbaceous: Areas dominated by gramanoid or herbaceous vegetation, generally greater than 80% of total vegetation. These areas are not subject to intensive management such as tilling but can be used for grazing.
  Pasture/Hay: Areas of grasses, legumes, or grass–legume mixtures planted for livestock grazing or the production of seed or hay crops, typically on a perennial cycle. Pasture/hay vegetation accounts for greater than 20% of total vegetation.
  Other Hay/Non Alfalfa: Mixed forage (Grass Mix below 25% Alfalfa, two or more Interseeded Coarse Grains, two or more Interseeded Grass Mix, two or more Interseeded Small Grains, two or more Legumes Interseeded), Grass/Small Grain Interseeding, Hay Oats and Peas, Legume/Coarse Grain, Legume/Grass Mixture, Legume/Small Grain, Legume/Small Grain/Grass, or Native Grass Interseeded.
43Woody WetlandsWoody Wetlands: Areas where forest or shrubland vegetation accounts for greater than 20% of vegetative cover and the soil or substrate is periodically saturated with or covered with water.
44Herbaceous WetlandsHerbaceous Wetlands: Areas where perennial herbaceous vegetation accounts for greater than 80% of vegetative cover and the soil or substrate is periodically saturated with or covered with water.
45Open WaterOpen Water: Lakes and ponds
Four levels of pollination demand are indicated by L (little = 0.05), M (modest = 0.25), G (great = 0.65), and E (essential = 0.95), which values follow pollinator dependency rates reported by Klein et al. (49). Note: pollinator-dependent crops for seed production such as alfalfa, onion, and asparagus were ignored in this study. Crop categories with bold font indicate that the pollinator dependency rate was modified from Klein et al. (49) to be more conservative for soybean, cotton, and canola/rape seed crops in the United States.
*,†,‡These crop categories have multi minor crops that are already clumped by NASS. Some crops may repeat in other categories. For example, there is “Pears” category in CDL but “Other Tree Crops” category also includes Pears. However, “Other Tree Crops” include many other tree crops, so that the Pears in “Other Tree Crops” is a minor crop. See the following detailed description.
*
Sorghum, Forage, Bamboo Shoots, Buckwheat, Guar, Ginger, Tea, Tannier, Lespedeza, Mushrooms, Indigo, Kenaf, Jojoba, Guayule, Carob, Jerusalem Artichokes/Sunchoke, Salsify (“Oyster Plant”), Taro, Rice (Wild), Yu Cha (“Tea Tree Oil”-oilseed plant), Crambe (Colewort), Psyllium, Quinoa, Meadowfoam, Lesquerella, Hesperaloe/Agave, Chia, Nursery, Teff, Kochia (Prostrata), Milkweed, Niger Seed, Cactus, Flowers (Horticulture), Sunn Hemp, Aloe Vera, or Canary Seed.
Elderberries, Avocados, Brussel Sprouts, Cucumbers, Pohole, Aronia (Chokeberry), Cassava, Pineapple, Okra, Currants, Rhubarb, Mulberries, Kohlrabi, Leeks, Gooseberries, Artichokes, Tangos, Dates, Shallots, Water Cress, Huckleberries, Sprite Melon, Broccoflower, Gailon/Gai Lein/Chinese Broccoli, Antidesma, Jujube, Pejibaye (Heart of Palm), Tomatillos, Scallions, Melongene, Israel Melons, Calaloo, Mayhaw Berries, Korean Golden Melon, Crenshaw Melon, Citron Melon, Chinese Bitter Melon, Casaba Melon, Canary Melon, or Calabaza Melon.
Other Tree Crops: Maple Sap, Coconuts, Chestnuts, Hazel Nuts, Macadamia Nuts, Ti, Cashew, Figs, Pears, Acerola (Barbados Cherry), Bananas, Coffee, Papaya, Plantain, Kiwifruit, Mangos, Persimmons, Plumcots, Quinces, Kumquats, Guava, Loquats, Passion Fruits, Atemoya (Custard Apple), Sapote, Carambola (Star Fruit), Caimito, Guanabana/Soursop, Breadfruit, Genip, Guavaberry, Jack Fruit, Rambutan, Mangosteen, Wampee, Longan, Lychee, Sapodilla, Cherimoya (Sugar Apple), or Canistel.
§
Detailed categories description referred to National Land Cover Database (www.mrlc.gov/nlcd06_leg.php).
According to NASS, Grassland/Pasture category collapses the following historical CDL categories, Grassland Herbaceous (code 171) and Pasture/Hay (code 181). In this study we also collapsed Other Hay/Non Alfalfa because it represents pasture areas although it is originally belonged to crop categories in CDL. However, we asked experts in all of the three categories and then averaged it for Grass/Pasture category.
Table S2.
Levels of significance for the differences of mean estimates among three ecoregions (Temperate Forests, Great Plains, and Mediterranean California)
CategoriesNesting resourcesFloral resourcesFloral bloom duration
 GroundCavityStemWoodSpringSummerFallSpringSummerFall
Developed/Open Space      *   
Developed/Low Intensity*     * ***
Developed/Medium Intensity *        
Developed/High Intensity *        
Barren          
Deciduous Forest          
Evergreen Forest          
Mixed Forest          
Shrubland          
Grass/Pasture          
Grassland/Herbaceous      ***  
Pasture/Hay*    ** * 
Other Hay/Non Alfalfa          
Woody Wetlands      **  
Herbaceous Wetlands          
Open Water          
Significant differences: *P < 0.05 and **P < 0.01.
Beyond these uncertainties, three additional caveats deserve mention. First, our assessment is based on a simple landscape model that predicts relative abundance of bees based on nesting resources, floral resources, and foraging distance. Although this model has proved to be informative in a variety of settings (32, 33, 54), it neither captures abundances of individual bee species nor reports visitation rates, pollination efficiency, or other variables important for realized pollination services. Second, although the model validation explained significant amounts of variance in field data, substantial variance remained unexplained. Clearly, other factors influence bee abundance in landscapes, but this study is intended as an initial national assessment of wild pollinators in general. Third, we evaluate trends over only 5 y; analysis of longer-term changes in both wild bee populations and land cover will provide a more robust assessment.
This first national assessment of status and trends of wild bee abundance will be valuable as a response to the recent federal mandates (20, 21) to direct additional research and management attention toward pollinators. A national program to detect future changes in bee populations has been estimated to cost $2,000,000 (55) and to require 5–10 y. Our national assessment can be used to focus such a costly effort, targeting bee and habitat surveys on regions that show high uncertainty, especially where agricultural demand for pollination services is high. Counties with mismatched levels of relative pollinator “supply” and “demand” warrant priority efforts to conserve and restore habitats for pollinators as well as other actions that can affect bees. As such efforts proceed, national assessments can be repeated with new information to update estimates, revise priorities, and track progress toward sustainable management of our nation’s wild pollinators.

Methods

Pollination Model.

The spatially explicit model of wild bee abundance (ref. 32; hereafter, the Lonsdorf model) generates an index of relative bee abundance at each spatial unit (e.g., map pixel). The model assumes that bees forage from a nest site to acquire floral resources in the surrounding landscape and the probability of acquiring resources declines exponentially with increasing distance between the nest site and floral resources. The model also assumes that nesting and floral resources vary among land-cover types in the landscape. To apply these model assumptions to the United States and evaluate their accuracy, we needed to identify a standard land-cover map, estimate the nesting and floral resources of each land cover, and validate the predictions with observations.

Data Sources.

We used the CDL (30-m resolution) to provide land-use and -cover types. This is the only such dataset produced annually at the national scale by the National Agricultural Statistics Service (NASS) since 2008. We reduced the number of crop cover types from over 100 to 32 representative categories based on shared crop characteristics and we retained 13 noncrop categories that are derived from the National Land Cover Database (Table S1). Based on a synthesis study (26), we applied an average foraging distance (670 m) of temperate wild bees as an input parameter for the forage distance function in the model.

Expert Opinion of Nesting and Floral Resources.

For each of the reclassified 45 land-use categories, a panel of 14 experts evaluated nesting suitability for four bee nesting guilds (ground, cavity, stem, and wood) and floral resource availability for three foraging seasons (spring, summer, and fall). Experts selected one of five options to represent nesting suitability or floral resource production (0.05, 0.25, 0.5, 0.75, or 0.95). For floral resources they selected the proportion of each 12-wk season in which the cover produced such resources (1–12 wk). For each estimate, panel members also specified one of four levels of certainty (none, low, medium, or high; SI Methods, Expert Survey and Table S2). We represented experts’ estimates and uncertainties as a continuous beta probability distribution (hereafter “pd”; SI Methods, Determining Final Probability Distribution of Resource Suitability). Ultimately we generated a single nesting suitability pd by summarizing across all experts and nesting guilds, and a floral resource pd in the same manner using floral seasons (SI Methods, Determining Final Probability Distribution of Resource Suitability and Fig. S5).
Fig. S5.
Example of determining floral and nesting resource probability distributions from expert elicitation. (A) Distributions of four types of nesting recourses and their average for orchard cover. (B) Distributions of seasonal floral resources and their average for watermelon cover. BD, bloom duration.

Modeling and Uncertainty.

The expert-informed probability distributions (pds) of nesting and floral resources for all land-use categories of the CDL were used as input parameters of the Lonsdorf model to predict a relative index (0–1) of wild bee abundance at each parcel of land (120 m × 120 m, one pixel). Because these input parameters are probability distributions, we can also express the bee abundance index as a probability distribution. We used Monte Carlo (MC) simulation to estimate the mean and SD for bee abundance at each parcel. These may be interpreted as the best estimate and the uncertainty of the index. Modeling uncertainty with probability distributions, however, bounds the uncertainty (measured as SD) possible for low and high estimates. This tends to result in greater estimates of uncertainty for moderate parameter values (Fig. S4), where bounding effects are not as important (30).

Model Validation.

We validated the model prediction and its uncertainty with field data of wild bee abundance. We used several data sets (SI Methods, Validation Data). All wild bees were observed at 180 sites on crop fields and seminatural and natural areas in six states between 2008 and 2013 (12, 5660). We also used a separate data set of bumble bees at 343 sites along roadsides in 40 states between 2008 and 2009 (15). We compared the model predictions based on expert-derived parameters and CDL corresponding to the year in which data were collected with the field data. Through the extensive model validation process, we verified that predicted bee abundance and its uncertainty respect current knowledge on wild bees (SI Methods, Model Validation Process and Fig. S3).

Mapping Status and Trends.

Status.

We used the expert-informed probability distributions (pds) and 2013 CDL as inputs to the Lonsdorf model to generate maps of the mean and uncertainty of bee abundance at 120-m resolution across the coterminous United States. For each pixel, we approximated the mean abundance index using the means of expert-informed pds and we represented uncertainty by estimating the SD of bee abundance indices again by using the expert-informed pds (SI Methods, Estimation of Mean and SD and Fig. S6). We recognize that model uncertainty may also have other sources, including the accuracy of classification for land-cover maps, but an examination of these effects on model uncertainty was beyond the scope of this study.
Fig. S6.
Estimation of the mean and SD of bee abundance index. (A) The relationship between bee abundance index using means of expert-informed probability distributions and mean from a 1,000 MC simulation for 10,000 randomly selected locations across the United States. (B) The relationship between SDI and SD of bee abundance index from a 1,000 MC simulation.

Trends.

We assessed trends in wild bee abundance as the differences in the mean bee abundance index between 2008 and 2013. To assess the uncertainty of trends, we calculated a pseudo-t value of the difference, by dividing the mean difference between the two years by the variation of the difference using the SD estimate for the two years (SI Methods, Likelihood of Index Change). High positive or negative values in the likelihood of change indicate a high likelihood of increase or decrease in the mean wild bee abundance index, respectively. Finally, we examined which land-use changes occurred in the counties whose predicted bee abundance changed the most, whether the abundance increased or decreased.

Supply and Demand Analysis.

We summarized the supply as the relative abundance of wild bees for each US county by averaging the bee abundance index and its uncertainty for all pixels within that county (SI Methods, Supply Assessment). We analyzed supply separately for 2008 and 2013. To assess the demand for pollination in each US county in 2008 and 2013, we summed the dependency-weighted area of all pollinator-dependent crops (49) for that county (SI Methods, Demand Assessment and Table S1). To assess the current status of supply and demand and to identify those counties with relatively low supply and high demand, we compared the average bee abundance with the dependency-weighted crop area. We also identified counties with relatively high uncertainty in the supply. To assess the trends in supply and demand between 2008 and 2013, we compared the likelihood of changes in bee abundance and the dependency-weighted crop area (SI Methods, Likelihood of Changes in Supply). Finally, we analyzed the trend of supply and demand for individual crops by comparing the likelihood of changes from 2008 to 2013 in wild bee abundance and dependency-weighted crop area across the entire coterminous United States.

SI Methods

Expert Survey.

We used a beta probability distribution to represent the experts’ estimates of the relative quality of each land-use type in terms of their provision of nesting and floral resources for wild bees. Their initial estimates of nesting (0.05, 0.25, 0.5, 0.75, or 0.95) quality were used as the mean values. Each expert provided two components of floral resource quality for every land-cover type: first, floral resource quality during the season (0.05, 0.25, 0.5, 0.75 or 0.95), and then the proportion of each of three seasons in which flowering occurred (spring, March–May; summer, June–August; and autumn, September– November). The experts’ mean value was simply the average of the two components (0–1). We translated their qualitative level of certainty into a SD, which equaled 0.15, 0.1, and 0.05 for low, medium, and high certainty, respectively. In all cases, we used a uniform probability distribution between 0 and 1; when experts were completely uncertain, this was 0 (none).
The panel of experts provided their estimates and uncertainty about nesting and floral resources for each of the noncrop categories in the five ecoregions (seven in Temperate Forest, two in Great Plains, one in Desert and Semiarid, one in Northwestern Forested Mountain, and five in Marine West Coast and Mediterranean California regions). To determine whether there was a difference in mean estimates among the ecoregions that would prevent us from pooling them, we used a Kruskal–Wallis rank sum test statistic to determine significance among the means resulting from expert estimates (i.e., the five levels of suitability values) for nesting and floral resources and used a one-way ANOVA to determine significance of the differences among means of the expert estimates for seasonal floral duration among the three ecoregions that had more than two responses. We found that only a few noncrop categories had significant differences in mean estimates of nesting and floral resources and floral duration among ecoregions (Table S2).

Determining Final Probability Distribution of Resource Suitability.

For each land-cover type, expert input provided 56 estimates for nesting resources (14 experts × 4 nesting guilds) and 42 estimates for floral resources (14 experts × 3 floral seasons). We ultimately used a single beta distribution based on a composite of estimates from all experts across all guilds or seasons for each land-cover type. We first determined the experts’ mean estimate for each guild or season before determining the average across all guilds and seasons.
To generate a composite beta distribution for each nesting guild or floral season, we fit a beta distribution to the collated values resulting from 1,000 iterations of a randomly selected expert and drew a random suitability value based on the selected expert’s beta distribution for each nesting guild. This resulted in four distributions for nesting and three distributions for floral resource quality. To obtain the final probability distribution (pd) across the four nesting resources (Fig. S5A), we repeated the process again except that the probability of drawing a particular nesting guild reflected the proportion of native bee species that are considered members of that guild: 70%, 20%, 5%, and 5% for ground, cavity, stem, and wood nesters, respectively (61). To obtain a final pd across the three floral resources (Fig. S5B), we fit a beta distribution to the collated values resulting from 1,000 iterations of a randomly drawn suitability value based on each of the season-specific beta distributions. We used the fitdistr function from the MASS R package (62) to fit a beta distribution to the probability density of the collated samples for the nesting and floral resources.

Characteristics of Final Resource Suitability Distributions.

As is characteristic of the relationship between the mean and variance of a binomial distribution, we found a parabolic relationship between mean suitability and the uncertainty level (Fig. S4). Experts evaluated some noncrop land uses such as deciduous and coniferous forests, shrublands, and grasslands, as well as some crop land uses such as pollinator-dependent croplands (orchard and oilseed crops) as relatively high suitability, but with high uncertainty (Fig. S4). The nesting and floral resources in forests and shrubs are associated with the canopy cover and tree density that determine understory herb cover and diversity (63, 64). For ground-nest resources, soil type is an important factor (65). Habitat resources in croplands may be more associated with farming practices such as pesticide use and tillage that can directly influence populations of native bees (33, 66).

Validation Data.

We used wild bee abundance data from six previously published studies across three ecoregions, two each from the Mediterranean ecoregion, the Great Plains, and Temperate Forests. In the Mediterranean regions, we used the observation data from 50 sites in watermelon fields (12) and 42 sites near hedgerow fields in CA (56, 57). In the Great Plains, we used the observation data from 13 sites in orchards and prairies in Iowa (58, 59) and 12 sites in corn, soy, and pasture fields in SD (available at figshare.com/s/e865f26c9a9e11e5b86e06ec4b8d1f61). In Temperate Forests, we used the data from 45 sites in cranberry fields in WI (60) and 19 sites in watermelon fields in PA and New Jersey (12). We also used a nationwide study of bumble bee abundance (15). All 343 sites across 40 states fell into three regions: 153 sites in the West (Mediterranean and Northwestern Forested Mountains), 70 sites in the Middle (Desert and Semiarid regions and the Great Plains), and 120 sites in the East (Temperate Forests). All validation data sets were collected between 2008 and 2013 and at locations where the CDLs are available.

Model Validation Process.

We used a generalized mixed effects model to fit the model prediction, wild bee abundance index, to field collected bee count data. We coded study, site-within-study, year, and study-within-year as random effects. This accounted for different numbers of field samples and studies among multiple years and variability in survey methods and sampling units among studies. In fitting the model for bumble bee data, we considered the year to be a random effect and we considered the bee abundance index to be a fixed effect on the field observation of bee abundance. We represented the strength of the model’s fit using Student’s t value of the coefficient for the bee abundance index. To examine and compare the strength of model fit between noninformative and expert-informed distributions for wild bee data, we generated input parameter sets that were randomly drawn from noninformative distributions (i.e., flat probability distributions) and from expert-informed beta distributions for nesting and floral resources. Then we fitted the model index to field data and investigated the range of the t values for each of the 1,000 combinations of parameter sets. In the model fitting, we used log transformed-count data for wild bees and square root transformed-count data for bumble bees to meet assumptions of normality. The generalized mixed effects model was performed using the R package lmer4 and Student t values were measured with the R package lmerTEST.
We were interested in evaluating three questions in the validation: (i) How well do randomly drawn expert-derived estimates correlate with observed data compared with noninformative estimates? (ii) How much does the mean (i.e., best) estimate improve the fit of the model? (iii) Do our expert-derived assessments of uncertainty predict the fit of the model with data? Using an MC simulation (n = 1,000) we found that model estimates using expert-informed pd were more positively correlated with observed bee abundance than those using noninformative pd (i.e., random values between 0 and 1) for both wild bees (Fig. S3A) and for bumble bees (Fig. S3B). The mean estimate of bee abundance for the MC simulation improved model prediction compared with field-observed bee data. Model fits were also better when using sites with low uncertainty (i.e., relatively low SD of the wild bee abundance index) than when using sites with high uncertainty (Fig. S3 CF). In addition, model estimates better predicted observed bumble bee abundances when they used parameters specific to the guild of cavity nesting bees, which are active throughout the summer (Fig. S3B).

Estimation of Mean and SD.

We developed an analytical estimate of variation in predicted bee abundance caused by parametric uncertainty because simulating the variation across the entire coterminous United States is simply not feasible. Instead, we randomly sampled 10,000 locations across the United States and determined a “true” estimate of mean and SD of model prediction and compared them to our analytical proxies. For each location, we ran an MC simulation (n = 1,000) to obtain the mean and SD of model prediction with input parameter values (i.e., nesting and floral suitability values for each land use). These were randomly drawn from the beta probability distributions described earlier. Using these N samples for the estimation, we then calculated the mean (M) and the SD (SD) of the bee abundance index. These represent, respectively, the mean estimate of the model prediction given the expert-derived distributions and the amount of uncertainty for the estimation:
M=1NiNBi
[S1]
SD=1N1iN(BiM)2,
[S2]
where Bi is the model prediction (i.e., the bee abundance index) from i of N parameter sets that were randomly drawn from prior distributions and N is the number of simulations.
We found that M was similar to the single model prediction based on means of expert-informed pds (Fig. S6A). To generate an analytical estimate of SD, we used the three estimates (mean and 25% and 75% bounds) from the beta distribution to generate three bee abundance indices (Bm, B25, and B75). Then we measured the variance of the two samples around the mean (SD index, SDI):
SDI=(B75Bm)2+(B25Bm)2.
[S3]
We found a significant nonlinear positive relationship between the estimated SD index and SD from the MC analysis (Fig. S6B). Finally, we used the fitted curve between SDI and SD to estimate SD (hereafter, SDE) of the mean wild bee abundance index.

Likelihood of Index Change.

We used our analytical estimate of uncertainty to develop a pseudo-t value to measure the likelihood of bee abundance index changes. We used the SD estimates for the 2 y of land-cover data used in the analysis (SDE08 and SDE13) to calculate SE of the difference of the mean index between 2008 and 2013 (SEdiff):
SEdiff=SDE0822+SDE1322.
[S4]
Then, we calculated the pseudo-t value by dividing the mean difference of the index by the SE:
pseudo-t=M13M08SEdiff,
[S5]
where M13 and M08 are the mean bee abundance indexes in 2013 and 2008.

Supply Assessment.

We used the mean wild bee abundance index (P) and SD estimate of the mean index (SDE) to assess the county-level status and trend of pollination service provision by wild bees to pollinator-dependent croplands across the United States. In each county, we calculated M and SDE in all NC pixels of each pollinator-dependent crop, C. Using these two indices, we measured the combined mean (MC) and combined SD estimate (SDEC) of pollination service-provision for each crop, c, as follows:
MC=n=1NcPnNC
[S6]
SDEC=(n=1NC[(SDEn)2+(PnMC)2])NC.
[S7]
We used these two measurements to represent mean and uncertainty of county-level pollination service provision by native bees for each crop. Then we measured the mean (MCO) and SD estimate (SDECO) of the bee abundance index for all number of types of pollinator-dependent crops (TN) at a county level:
MCO=n=1TN(NCn×MCn)n=1TNNCn
[S8]
SDECO=n=1TN[NCn(SDECn2+(MCnMCO)2)]n=1TNNCn.
[S9]

Likelihood of Changes in Supply.

We used these two indices as the mean and SD at a county level to measure the uncertainty in changes of bee supply (i.e., likelihood of changes) by dividing the mean difference between 2013 and 2008 with a variance as follows:
MCO13MCO08(SDECO13N13)2+(SDECO08N08)2,
[S10]
where N13 and N08 mean relative sample size of all parcels of pollinator-dependent croplands in 2013 and 2008, respectively (e.g., if there were 110 and 100 parcels in 2013 and 2008, then N13 and N08 are 1.1 and 1.0, respectively).

Demand Assessment.

We used the median value of the range of pollination dependency rate for CDL crop categories (Table S1) reported by Klein et al. (49) to calculate the demand-weighted crop area (DCO) at a county scale across the United States. We calculated the demand-weighted crop areas in 2008 and 2013 at a county level as follows:
DCO=i=1TNDRiAi,
[S11]
where DRi is the pollinator dependency rate of crop i and Ai is the area of crop i. The county-level change in pollination demand was measured by calculating the differences in demand-weighted crop areas between 2013 and 2008 (i.e., DCO13DCO08).

Acknowledgments

We thank the nonauthor expert survey participants: J. Cane, J. Cruz, E. Evans, K. Gill, J. Hemberger, T. Harrison, J. Hopwood, H. Sardinas, C. Stanley-Stahr, M. Vaughan, and M. Veit. We also thank C. Kremen, S. Hendrix, R. Winfree, C. Mogren, H. Gaines-Day, C. Gratton, H. Sardinas, A. Sciligo, and K. Nemec, who provided their field-observation datasets. We thank P. Willis at the USDA National Agricultural Statistics Service for his assistance regarding our questions about NASS-Cropland Data Layer and M. O’Neal and D. Cohen for contributing their knowledge of pollination of soybean and canola crops. We thank L. Richardson and C. Nicholson for comments that improved our manuscript. This research was supported by the USDA-NIFA Specialty Crop Research Initiative, from Project 2012-51181-20105: Developing Sustainable Pollination Strategies for US Specialty Crops.

Supporting Information

Supporting Information (PDF)
Supporting Information
pnas.1517685113.sd01.xlsx

References

1
R Winfree, NM Williams, H Gaines, JS Ascher, C Kremen, Wild bee pollinators provide the majority of crop visitation across land-use gradients in New Jersey and Pennsylvania, USA. J Appl Ecol 45, 793–802 (2008).
2
J Ollerton, R Winfree, S Tarrant, How many flowering plants are pollinated by animals? Oikos 120, 321–326 (2011).
3
S Lautenbach, R Seppelt, J Liebscher, CF Dormann, Spatial and temporal trends of global pollination benefit. PLoS One 7, e35954 (2012).
4
RA Morse, NW Calderone, The value of honey bees as pollinators of U.S. crops in 2000. Bee Culture 128, 1–15 (2000).
5
JE Losey, M Vaughan, The economic value of ecological services provided by insects. Bioscience 56, 311–323 (2006).
6
NW Calderone, Insect pollinated crops, insect pollinators and US agriculture: Trend analysis of aggregate data for the period 1992-2009. PLoS One 7, e37235 (2012).
7
MA Aizen, LD Harder, The global stock of domesticated honey bees is growing slower than agricultural demand for pollination. Curr Biol 19, 915–918 (2009).
8
NA Steinhauer, et al., A national survey of managed honey bee 2012-2013 annual colony losses in the USA: Results from the Bee Informed Partnership. J Apic Res 53, 1–18 (2014).
9
R Winfree, NM Williams, J Dushoff, C Kremen, Native bees provide insurance against ongoing honey bee losses. Ecol Lett 10, 1105–1113 (2007).
10
LA Garibaldi, et al., Wild pollinators enhance fruit set of crops regardless of honey bee abundance. Science 339, 1608–1611 (2013).
11
C Brittain, C Kremen, AM Klein, Biodiversity buffers pollination from changes in environmental conditions. Glob Change Biol 19, 540–547 (2013).
12
R Rader, J Reilly, I Bartomeus, R Winfree, Native bees buffer the negative impact of climate warming on honey bee pollination of watermelon crops. Glob Change Biol 19, 3103–3110 (2013).
13
C Kremen, Managing ecosystem services: What do we need to know about their ecology? Ecol Lett 8, 468–479 (2005).
14
JC Grixti, LT Wong, SA Cameron, C Favret, Decline of bumble bees (Bombus) in the North American Midwest. Biol Conserv 142, 75–84 (2009).
15
SA Cameron, et al., Patterns of widespread decline in North American bumble bees. Proc Natl Acad Sci USA 108, 662–667 (2011).
16
I Bartomeus, et al., Historical changes in northeastern US bee pollinators related to shared ecological traits. Proc Natl Acad Sci USA 110, 4656–4660 (2013).
17
D Goulson, E Nicholls, C Botías, EL Rotheray, Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. Science 347, 1255957 (2015).
18
SG Potts, et al., Global pollinator declines: Trends, impacts and drivers. Trends Ecol Evol 25, 345–353 (2010).
19
; Committee on the Status of Pollinators in North America, National Research Council Status of Pollinators in North America (National Academy, Washington, DC, 2007).
20
; The White House, Presidential memorandum – Creating a federal strategy to promote the health of honey bees and other pollinators (The White House, Washington, DC). (2014).
21
; Pollinator Health Task Force National strategy to promote the health of honey bees and other pollinators (The White House, Washington, DC, 2014).
22
TH Roulston, K Goodell, The role of resources and risks in regulating wild bee populations. Annu Rev Entomol 56, 293–312 (2011).
23
NM Williams, J Regetz, C Kremen, Landscape-scale resources promote colony growth but not reproductive performance of bumble bees. Ecology 93, 1049–1058 (2012).
24
M Rundlof, AS Persson, HG Smith, R Bommarco, Late-season mass-flowering red clover increases bumble bee queen and male densities. Biol Conserv 172, 138–145 (2014).
25
SS Greenleaf, NM Williams, R Winfree, C Kremen, Bee foraging ranges and their relationship to body size. Oecologia 153, 589–596 (2007).
26
TH Ricketts, et al., Landscape effects on crop pollination services: Are there general patterns? Ecol Lett 11, 499–515 (2008).
27
HM Hines, SD Hendrix, Bumble bee (Hymenoptera: Apidae) diversity and abundance in tallgrass prairie patches: Effects of local and landscape floral resources. Environ Entomol 34, 1477–1484 (2005).
28
C Carvell, et al., Bumble bee species’ responses to a targeted conservation measure depend on landscape context and habitat quality. Ecol Appl 21, 1760–1771 (2011).
29
HS Sardiñas, C Kremen, Evaluating nesting microhabitat for ground-nesting bees using emergence traps. Basic Appl Ecol 15, 161–168 (2014).
30
GM Lee, D Molen, HP Boogaard, H Klis, Uncertainty analysis of a spatial habitat suitability model and implications for ecological management of water bodies. Landscape Ecol 21, 1019–1032 (2006).
31
SJ O’Neill, TJ Osborn, M Hulme, I Lorenzoni, AR Watkinson, Using expert knowledge to assess uncertainties in future polar bear populations under climate change. J Appl Ecol 45, 1649–1659 (2008).
32
E Lonsdorf, et al., Modelling pollination services across agricultural landscapes. Ann Bot (Lond) 103, 1589–1600 (2009).
33
CM Kennedy, et al., A global quantitative synthesis of local and landscape effects on wild bee pollinators in agroecosystems. Ecol Lett 16, 584–599 (2013).
34
CJ Johnson, MP Gillingham, Mapping uncertainty: Sensitivity of wildlife habitat ratings to expert opinion. J Appl Ecol 41, 1032–1041 (2004).
35
JV Murray, et al., How useful is expert opinion for predicting the distribution of a species within and beyond the region of expertise? A case study using brush-tailed rock-wallabies Petrogale penicillata. J Appl Ecol 46, 842–851 (2009).
36
CA Czembor, WK Morris, BA Wintle, PA Vesk, Quantifying variance components in ecological models based on expert opinion. J Appl Ecol 48, 736–745 (2011).
37
F Hendrickx, et al., How landscape structure, land-use intensity and habitat diversity affect components of total arthropod diversity in agricultural landscapes. J Appl Ecol 44, 340–351 (2007).
38
V Le Feon, et al., Intensification of agriculture, landscape composition and wild bee communities: A large scale study in four European countries. Agric Ecosyst Environ 137, 143–150 (2010).
39
CK Wright, MC Wimberly, Recent land use change in the Western Corn Belt threatens grasslands and wetlands. Proc Natl Acad Sci USA 110, 4134–4139 (2013).
40
AB Bennett, TD Meehan, C Gratton, R Isaacs, Modeling pollinator community response to contrasting bioenergy scenarios. PLoS One 9, e110676 (2014).
41
; USDA Economic Research Service, The Conservation Reserve Program is regionally concentrated. Available at www.ers.usda.gov/data-products/chart-gallery/detail.aspx?chartId=40030&ref=collection&embed=True. Accessed August 1, 2015. (2013).
42
JL Hopwood, The contribution of roadside grassland restorations to native bee conservation. Biol Conserv 141, 2632–2640 (2008).
43
BR Blaauw, R Isaacs, Flower plantings increase wild bee abundance and the pollination services provided to a pollination-dependent crop. J Appl Ecol 51, 890–898 (2014).
44
C Westphal, I Steffan-Dewenter, T Tscharntke, Mass flowering crops enhance pollinator densities at a landscape scale. Ecol Lett 6, 961–965 (2003).
45
S Sinnathamby, et al., Pollinator decline: US agro-socio-economic impacts and responses. J Nat Environ Sci 4, 1–13 (2013).
46
C Kremen, NM Williams, RW Thorp, Crop pollination from native bees at risk from agricultural intensification. Proc Natl Acad Sci USA 99, 16812–16816 (2002).
47
GM Foody, Status of land cover classification accuracy assessment. Remote Sens Environ 80, 185–201 (2002).
48
CK Wright, MC Wimberly, Reply to Kline et al.: Cropland data layer provides a valid assessment of recent grassland conversion in the Western Corn Belt. Proc Natl Acad Sci USA 110, E2864 (2013).
49
AM Klein, et al., Importance of pollinators in changing landscapes for world crops. Proc Biol Sci 274, 303–313 (2007).
50
N Gallai, JM Salles, J Settele, BE Vaissiere, Economic valuation of the vulnerability of world agriculture confronted with pollinator decline. Ecol Econ 68, 810–821 (2009).
51
TC Giannini, GD Cordeiro, BM Freitas, AM Saraiva, VL Imperatriz-Fonseca, The dependence of crops for pollinators and the economic value of pollination in Brazil. J Econ Entomol 108, 849–857 (2015).
52
KA Orford, IP Vaughan, J Memmott, The forgotten flies: The importance of non-syrphid Diptera as pollinators. Proc Biol Sci 282, 20142934 (2015).
53
JH Cane, RL Minckley, LJ Kervin, TH Roulston, NM Williams, Complex responses within a desert bee guild (Hymenoptera: Apiformes) to urban habitat fragmentation. Ecol Appl 16, 632–644 (2006).
54
TH Ricketts, E Lonsdorf, Mapping the margin: Comparing marginal values of tropical forest remnants for pollination services. Ecol Appl 23, 1113–1123 (2013).
55
G Lebuhn, et al., Detecting insect pollinator declines on regional and global scales. Conserv Biol 27, 113–120 (2013).
56
C Kremen, LK M’Gonigle, Small-scale restoration in intensive agricultural landscapes supports more specialized and less mobile pollinator species. J Appl Ecol 52, 602–610 (2015).
57
LA Morandin, C Kremen, Hedgerow restoration promotes pollinator populations and exports native bees to adjacent fields. Ecol Appl 23, 829–839 (2013).
58
SD Hendrix, KS Kwaiser, SB Heard, Bee communities (Hymenoptera: Apoidea) of small Iowa hill prairies are as diverse and rich as those of large prairie preserves. Biodivers Conserv 19, 1699–1709 (2010).
59
AN Harmon-Threatt, SD Hendrix, Prairie restorations and bees: The potential ability of seed mixes to foster native bee communities. Basic Appl Ecol 16, 64–72 (2015).
60
H Gaines-Day, Do bees matter to cranberry? The effect of bees, landscape, and local management on cranberry yield. PhD dissertation (Univ of Wisconsin, Madison, WI). (2013).
61
CD Michener The Bees of the World (Johns Hopkins Univ Press, Baltimore, 2007).
62
; R Development Core Team, R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna), Version 2.15.0. (2012).
63
I Steffan-Dewenter, et al., Tradeoffs between income, biodiversity, and ecosystem functioning during tropical rainforest conversion and agroforestry intensification. Proc Natl Acad Sci USA 104, 4973–4978 (2007).
64
R Winfree, T Griswold, C Kremen, Effect of human disturbance on bee communities in a forested ecosystem. Conserv Biol 21, 213–223 (2007).
65
JH Cane, Soils of ground-nesting bees (Hymenoptera, Apoidea) - texture, moisture, cell depth and climate. J Kans Entomol Soc 64, 406–413 (1991).
66
RE Shuler, TH Roulston, GE Farris, Farming practices influence wild pollinator populations on squash and pumpkin. J Econ Entomol 98, 790–795 (2005).

Information & Authors

Information

Published in

Go to Proceedings of the National Academy of Sciences
Proceedings of the National Academy of Sciences
Vol. 113 | No. 1
January 5, 2016
PubMed: 26699460

Classifications

Submission history

Published online: December 22, 2015
Published in issue: January 5, 2016

Keywords

  1. crop pollination
  2. ecosystem services
  3. habitat suitability
  4. land-use change
  5. uncertainty

Acknowledgments

We thank the nonauthor expert survey participants: J. Cane, J. Cruz, E. Evans, K. Gill, J. Hemberger, T. Harrison, J. Hopwood, H. Sardinas, C. Stanley-Stahr, M. Vaughan, and M. Veit. We also thank C. Kremen, S. Hendrix, R. Winfree, C. Mogren, H. Gaines-Day, C. Gratton, H. Sardinas, A. Sciligo, and K. Nemec, who provided their field-observation datasets. We thank P. Willis at the USDA National Agricultural Statistics Service for his assistance regarding our questions about NASS-Cropland Data Layer and M. O’Neal and D. Cohen for contributing their knowledge of pollination of soybean and canola crops. We thank L. Richardson and C. Nicholson for comments that improved our manuscript. This research was supported by the USDA-NIFA Specialty Crop Research Initiative, from Project 2012-51181-20105: Developing Sustainable Pollination Strategies for US Specialty Crops.

Notes

This article is a PNAS Direct Submission.
Dr. Mogren's unpublished bee observation data has been deposited online at figshare.com/s/e865f26c9a9e11e5b86e06ec4b8d1f61.

Authors

Affiliations

Gund Institute for Ecological Economics, University of Vermont, Burlington, VT 05405;
Eric V. Lonsdorf
Gund Institute for Ecological Economics, University of Vermont, Burlington, VT 05405;
Biology Department, Franklin and Marshall College, Lancaster, PA 17604;
Neal M. Williams
Department of Entomology and Nematology, University of California, Davis, CA 95616;
Claire Brittain
Department of Entomology and Nematology, University of California, Davis, CA 95616;
Rufus Isaacs
Department of Entomology, Michigan State University, East Lansing, MI 48824;
Jason Gibbs
Department of Entomology, Michigan State University, East Lansing, MI 48824;
Taylor H. Ricketts
Gund Institute for Ecological Economics, University of Vermont, Burlington, VT 05405;
Rubenstein School of Environment and Natural Resources, University of Vermont, Burlington, VT 05405

Notes

1
To whom correspondence should be addressed. Email: [email protected].
Author contributions: I.K., E.V.L., N.M.W., C.B., R.I., J.G., and T.H.R. designed research; I.K., E.V.L., N.M.W., C.B., and T.H.R. performed research; I.K., E.V.L., and T.H.R. analyzed data; and I.K., E.V.L., N.M.W., C.B., R.I., J.G., and T.H.R. wrote the paper.

Competing Interests

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


Citation statements




Altmetrics

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

    Loading...

    View Options

    View options

    PDF format

    Download this article as a PDF file

    DOWNLOAD PDF

    Get Access

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to access the full text.

    Single Article Purchase

    Modeling the status, trends, and impacts of wild bee abundance in the United States
    Proceedings of the National Academy of Sciences
    • Vol. 113
    • No. 1
    • pp. 1-E104

    Media

    Figures

    Tables

    Other

    Share

    Share

    Share article link

    Share on social media