Extreme species density of bees (Apiformes, Hymenoptera) in the warm deserts of North America

Despite the long intertwined evolutionary histories of bees and plants, bee diversity peaks in the xeric areas of the eastern and western hemispheres and not the tropics, where plant diversity is greatest. Intensive sampling in the northeast Chihuahuan Desert of Mexico and the United States provide the first quantitative estimate of bee species richness where high diversity had been predicted in North America from museum records. We find that the density of bee species in a limited area of 16 km2 far exceeds any other site in the world and amounts to approximately 14% of the bee species described from the United States. Long-term studies of bees and other pollinators from areas that are minimally impacted by humans provide much-needed baseline data for studies of bees where human impacts are more severe and as climate change accelerates.


Introduction
The keystone role of bees in ecosystems has led to greater appreciation and a deeper understanding of some previously understudied aspects of their biology. These include the contribution of non-managed bee species to crop pollination (Kremen et al. 2002;Winfree et al. 2008;Garibaldi et al. 2013), ecological and life history variables that predict vulnerability to environmental change among species (Winfree et al. 2009;Williams et al. 2010), and the environmental conditions that favor increases in bee populations and their species richness (Cane et al. 2006;Williams et al. 2010;Ollerton et al. 2014). Studies that have compared historical records from museums to more recent samples have concluded that bees along with other pollinators are in decline due to human activity (Biesmeijer et al. 2006;Lebuhn et al. 2012;Burkle et al. 2013;Ollerton et al. 2014; Bartomeus et al. 2018;Mathiasson and Rehan 2019). Comparisons of bee faunas at sites where land-use histories differ have reached similar conclusions (reviewed in Winfree et al. 2011: but see Winfree et al. 2007;Minckley 2014). However, with some notable recent exceptions (Herrera 2019;Meiners et al. 2019 and references therein), most of this recent work has not extended beyond one or several years and has been focused on faunas in or near human-disturbed habitats. Needed are more long-term studies from habitats that are little disturbed by humans to provide baseline data to appropriately gauge how bee dynamics are influenced where there is anthropogenic disturbance (Winfree 2010; Archer et al. 2014;De Palma et al. 2016). Herrera (2019) has gone further to emphasize that studies are needed particularly where diversity is highest.
This study was done to thoroughly document the number of bee species (=species richness) in a bee biodiversity hotspot and serves to address the concern of Herrera (2019) mentioned above. We use repeat sampling to estimate the total number of bee species (observed as well as those not captured) from a defined area so both the number of species (species richness, an area free estimate) and number of species from a known area (species density) can be compared to data from other studies that are focused on elucidating biogeographic patterns of bees. The amount of area sampled strongly predicts species richness in all groups (Rozenzweig 1995) and allows species density to be evaluated separately from species richness (Gotelli and Colwell 2001). Global patterns of species richness remain one of the central questions to evolutionary ecology (Colwell and Lees 2000;Lomolino et al. 2017) and anomalous diversity patterns sometimes provide unusual insights into the mechanisms that generate diversity (Kindlmann et al. 2012). Bees are a species-rich group (N = 20,000 species) and have very different and counterintuitive biogeographic patterns from their hosts, the flowering plants. Using species lists and museum records, Michener (1979) showed that the warm desert areas of North America and xeric regions around the Mediterranean Sea of the eastern hemisphere were unusually rich in bee species. Why bee diversity does not follow the temperate to tropic increase in species richness of their floral hosts remains poorly understood (Orr et al. 2020). One hypothesis is the persistently wet soils of tropical regions favor fungi and other pathogens of bees that nest in the soil (Michener 1979). This is consistent with the observation that many bee groups that occur in the tropics have their nests above-ground in trees and broken stems where moisture is reduced (Michener 2007). Bee diversity in the tropics may also be depressed because many tropical groups are highly eusocial with colonies that can be large. If floral resources are limiting, at least sometimes, social bees may successfully outcompete solitary bees. Also possible is that some feature of the ecosystems where warm xeric climates occur favor bee diversification. In parts of the Chihuahuan and Sonoran Deserts of North America, and around the Mediterranean Sea in the eastern hemisphere rainfall and temperature patterns result in separate spring and summer blooms. Where this occurs, bee activity is also biseasonal and many of these bee species are active in only one season. The temporally separated but geographically co-occurring bee fauna may allow more species to occur in one area than could coexist otherwise (Michener 1979;Minckley 2008). Finally, topographical complexity and biogeography may also contribute to global patterns of bee diversity. Where bee species richness peaks are also areas with ongoing mountain building and associated tectonics (Casas-Sainz and de Vicente 2009; González-León et al 2011) that have been shown to favor diversification in other groups (Rahbek et al. 2019), and are where biomes converge (Spector 2002). Although it is possible that one factor is responsible for most variation in bee diversity, it is more likely that some combination of the factors listed above underlie where and how many bees occur worldwide.
We examined the hypothesis that bee richness peaks in the warm deserts of North America using intensive, standardized methods in one area of the Chihuahuan Desert along the international border between the United States and Mexico. Where this study was done had been cattle ranchland since the early 1800's, but has not been grazed for the past 20 years in Mexico and the past 50 years in the United States (Minckley 2013). This region is a biogeographical crossroads where the subtropical dry forest, Great Plains, Chihuahuan and Sonoran Desert intersect (Brown 1994). The pollinator fauna is largely native and strongly dominated by bees. The honey bee, Apis mellifera, is the sole introduced bee species representing 3.5% of all bee captures. The vast majority of collections in this study were within an area of approximately 16 km 2 on the floor of the San Bernardino Valley. Repeat sampling using the same protocol allowed us to use abundance data to make a quantitative estimate of bee species density including species not captured. We then compared estimates of bee species density from this 1-year intensive sampling to nine years of collecting using sampling by several methods from this same area. Combining sampling methods that vary in effort and efficacy violates some assumptions of the models used to estimate the richness of observed plus unsampled species (Gotelli and Colwell 2001). However, that bee communities, particularly those in deserts, vary dramatically across all spatial and temporal scales (Minckley et al. 1999;Williams et al. 2001;Russo et al. 2015;Herrera 2019) and in response to sampling effort, protocol, and technique (Toler et al. 2005;Roulston et al. 2007;Westphal et al. 2008;Droege et al. 2010;Cane et al. 2013;Prendergast et al. 2020), suggests that multiple approaches will yield more bee species, and more accurate estimates of the number of species. For an unusually well-studied, species-rich community of ants in Costa Rica, Longino et al. (2002) showed that the actual species richness was greatly underestimated with datasets from single sampling protocols; and only by combining different datasets was it possible to obtain an accurate estimate. In this study, bee species density estimates based on repeat samples with one protocol is compared to estimates based on a larger dataset from the same area that combined sampling protocols from more sites over nine years.
Our prediction was that intensive sampling using one technique (pan traps) over one year would underestimate the observed species richness obtained over more years using multiple sampling techniques from the same area. However, we also predicted this intensive pan trapping would provide a reasonable estimate of species richness using statistical methods that extrapolate from species observed to include species not captured but present (Gotelli and Colwell 2001).
We also gauged how well our samples from the valley floor estimate the number of species in a broader area, by compiling a dataset from collections made in two nearby areas and making an assumption that this list of species from all collections combined was the best estimate for bee richness in this region. Extrapolation to the expected species richness from the multiyear samples made on the valley floor closely approximated the observed species richness of the larger regional dataset. In comparison to 21 other long-term studies of bees, the density of bee species in this part of the Chihuahuan Desert is greater than that reported from any other area.

Methods
The study was done in the San Bernardino Valley in the northeastern Chihuahuan Desert of western North America near the confluence of Arizona and New Mexico, USA and Chihuahua and Sonora, Mexico. The San Bernardino Valley runs north-south across the Mexico-United States border in northeastern Sonora, Mexico and southeastern Arizona, USA. Elevation is approximately 1070 m and climate is xeric temperate with an annual average precipitation of 360 mm/ year. Approximately 50% of the annual precipitation occurs in July and August and the two driest months are April and May based on long-term weather records 30 kilometers west of the San Bernardino Valley at Douglas, Arizona. Minckley (2008Minckley ( , 2013 provides a detailed description of the climate, area and vegetation. The main study area contains a number of habitats ranging from desert scrub to permanent springs that are the headwaters of the San Bernardino River. In 2000, seven to eleven permanently marked 1-ha sites were established in five different vegetation types (Fig. 1). Two vegetation types were in proximity to surface (riparian) or close subsurface (mesquite forest) water and were more mesic than other vegetation types. Two vegetation types were upland and were much drier (grassland, desert scrub) with cacti and shrubs characteristic of Chihuahuan Desert, such as Acacia constricta Benth. (cat claw), Larrea tridentata (DC.) Coville (creosote bush), and Flourensia cernua DC. (tarbush). One habitat was abandoned agricultural fields. All sites were within a 16 km 2 area and ranged from 1100 to 1220 m in elevation.
Collections of bees were made at all of the permanent sites by pan trap in 2001 and at approximately half of these sites from 2002 to 2008. Pan trap sampling was done using three yellow, three white and three blue 177 ml. pan traps placed out in alternating colors and spaced 3-4 m apart along a transect following the protocol of LeBuhn et al. (2007). All bee specimens were identified to species and sex by taxonomic experts and published keys (Suppl. material 1: Table S1). Most were identified to recognized species, except for 73 species that were recognized as morphospecies. These morphospecies are 1.43% of all individuals, and most are members of Lasioglossum (Dialictus), Dufourea and Megachile. All species and their abundance are listed in Table 1 and are held in the Minckley collection at the University of Rochester. Data are available from the Global Biodiversity Information Facility (GBIF) under the RL Minckley Insect and Plant Collection (https://www.gbif.org/dataset/84b69033-f9de-40bc-8d5f-8ee41987d6cc).

Species density
Bee species density in the San Bernardino Valley was first estimated from the bee collections captured in pan traps at permanent sites in 2001 (Fig. 1). These samples were made at all 45 permanent sites on 10 to 14-day intervals from mid-April to mid-September, except for three weeks in June, the driest time of the year (Tables 2, 3).
A second estimate of bee species density was done from all samples in the same area circumscribed by the permanent sites described above (hereafter referred to as the  Bossert et al. (2020).    (Cockerell, 1924) 118   (Michener, 1954) 13 13 1 3 A. (Isosmia) rubrella (Michener, 1943) 0 2 0 5 Atoposmia (Atoposmia) anthodyta (Michener, 1943) 0 (Michener, 1954) 2 2 1 0 A. (Eremosmia) beameri (Michener, 1951) 1 4 1 0 A. (Eremosmia) hypostomalis (Michener, 1954)  "core area") made from 2000 to 2008 (Fig. 1). Not included were the samples made by pan traps at permanent sites in 2001. In this dataset, samples were included regardless of sampling technique (pan trap, aerial netting at flowering plants) or duration, and included collections made at the permanent sites and intervening sites (Tables 2, 3). We tested if results from the 2001 intensive pan trap study differed from the results from the same "core area" sampled from 2000-2008 by comparing the accumulation curves from both datasets after they were standardized for equal sample size and coverage following the approach by Chao et al. (2014). Predicted species richness was estimated after the two samples were extrapolated to 42,798 individuals, twice the abundance of the smaller sample (pan traps from 2001 at the permanent sites, N observed = 21,399 individuals) as recommended by Chao et al. (2014) and implemented in the program iNext (http://chao.stat.nthu.edu.tw/wordpress/software_download/inext-online/). Confidence intervals and standard errors were calculated from 50 bootstrap replicates. Curves were significantly different if confidence intervals did not overlap. The contribution of vegetation types to bee species richness is not considered in this study. An earlier study found that there was little difference among vegetation types in bee species richness and composition, most likely because of the close proximity of the habitats and sites sampled (Minckley 2014).

Species richness
To examine if the samples from the valley floor approximated the species richness of the fauna in a larger regional area, we compared species composition and estimated species richness to our complete dataset of bees from the San Bernardino Valley that includes all bee samples including two from nearby areas. One area extended 18 km east from the valley floor area up to desert grassland habitat. We refer to the collections from this area as from "peripheral" sites (Fig. 1). These samples were made during the months when samples were made in the valley floor and used pan traps and aerial netting. A second set of samples was made from 2009-2014, as part of a study of cactus reproductive biology by W. Radke, at a limestone outcrop dominated by desert vegetation approximately 10.5 km north-northeast of the valley floor area (Fig. 1). Most, but not all, of these samples were made in the spring months and included sampling by pan traps and nets. Tables 2, 3 show details on the number of collections and sampling protocols used for these two areas.
Here, and as described above, the standardized curves of species interpolation and extrapolation and sample coverage were generated by the program iNEXT (http:// chao.stat.nthu.edu.tw/wordpress/software_download/inext-online/), and comparisons of datasets and estimates of species richness were done as recommended by Chao et al. (2014). Predicted species richness from samples made in the core area, core area plus  2001 pan traps, and samples from all sites made in the region was estimated after the samples were extrapolated to 92,240 individuals, or twice the abundance of the smallest sample (core area collections, N observed = 47,120 individuals). Confidence intervals and standard errors were calculated from 50 bootstrap replicates, and accumulation curves were considered significantly different if confidence intervals did not overlap.
Comparison to other studies Meiners et al. (2019) recognized that comparisons among long-term bee studies would be more informative if differences in the area sampled were considered. To standardize species-area relationships, they generated the species-area curve from 22 long-term bee studies in the United States and then ranked each study based on the percent difference between observed and predicted species density. Meiners et al. (2019) included the San Bernardino Valley study in their analysis but calculated species density using the area of the San Bernardino National Wildlife Refuge (1088 km 2 ), instead of the area among permanent sampling sites (16 km 2 ), and an earlier estimate of bee species richness (N = 383) reported in Minckley (2008) than known at present (N = 473). We used the same data and study sites included in the Meiners et al. (2019) study to compare new results on species density from the San Bernardino Valley to the other studies. Table 4 in Meiners et al. (2019) with changes updated to reflect the findings in this study are in Suppl. materials 1: Table S2.

Results
Samples limited only to pan traps at permanent sites in 2001 yielded a total of 21,399 specimens representing 277 bee species (Table 1) from 848 samples (day x site). Samples over nine years and more collections (N = 1370) with both pan traps and netting at flowers (Table 2), and from the same area defined by the permanent sites, yielded more than twice the specimens and 447 species (Table 3), including 58 species not collected in the 2001 pan trap dataset. Of the total specimens represented when the two datasets are combined, 31% were from the intensive pan trap sampling in 2001 and 69% were from sampling using nets and pan traps across multiple years. The species accumulation curves from these two sets of samples differ significantly after 200 individuals, as judged by non-overlapping confidence intervals (Fig. 2), despite having been collected from the same area (Fig. 1) and representing robust samples based on estimators that include sample coverage (Suppl. material 1: Figures S1-S4). The core area samples rarefied to the number of individuals captured in the 2001 pan trap samples (N = 42,798) yielded significantly more species (N = 437.6) than the 321.8 species predicted from the 2001 pan trap samples (Fig. 2). Overall, the observed species richness from the intensive 2001 pan trap sampling yielded 62% of the species collected from the core area over all 9 years. Far fewer collections of bees were made from sites outside the core area and these amounted to fewer individuals and species of bees (Tables 2, 3). Nevertheless, the 221 species from the peripheral and limestone samples are a large fraction (79%) of the 277 species captured by pan traps in 2001 and about half the number of species represented in the core area dataset (Table 3). Twelve species were unique to the peripheral site samples. The limestone outcrop samples were the smallest in terms of specimens collected and species (Tables 2, 3), however, 12 of the species were unique even though sampling was focused primarily on pan traps and bees were recorded from only two floral hosts (Table 2).
Including all samples at all sites in the area defined by the permanent sites (2001 pan traps and core area) yielded 473 bee species (Table 1). With the additional sampling made at the peripheral sites and the limestone outcrop there are 497 bee species known from the San Bernardino Valley (Table 1). Overall, the two datasets from the area defined by the permanent sites (2001 pan traps and core area) captured most of the bee species richness in the desert areas of the San Bernardino Valley; samples for sites peripheral to the core area increased the total number of bee species by only 12 species.
Pooled across all sites in the region the collections total to 74,147 individuals and 497 bee species (Table 1). Although the core plus 2001 pan trap dataset has many more individuals (N = 68,519) than the core dataset (N = 42,120), it results in a species accumulation curve that is significantly lower than the curve based on all datasets combined (Fig. 3). Up to approximately 42,000 individuals, the species accumulation curve from the core dataset is significantly lower than the curve based on all datasets combined. However, at greater numbers of individuals the curves are statistically indistinguishable based on overlapping confidence intervals (Fig. 3). Extrapolations of species richness to 94,240 individuals predicted 508.4 species for the core area dataset and 498.3 species for the 2001 pan trap combined with the core area dataset and closely approximate the estimate of 516.1 bee species based on all samples made in the San Bernardino Valley (Fig. 3). The increasingly flattened accumulation curves as more individuals are sampled indicate these estimates of bee species richness are unusually robust (see Williams et al. 2001), and biologically reasonable, given that they closely approximate the number of bee species (N = 497) from all samples combined made by us in the San Bernardino Valley.
Among the 22 long-term bee studies that Meiners et al. (2019) listed, three (Yosemite National Park, Clark County, Nevada and Grand Staircase Escalante National Monument) reported more bee species and only one (Black Rock Forest Reserve, New York) was from a smaller area than the core area sampled in the San Bernardino Valley (Suppl. material 1: Table S2). The updated information from the San Bernardino Valley changed the species area curve and the rank order of species density for the 22 longterm bee studies as judged by the percentage change between observed and predicted species per area (Fig. 4a, b). Meiners et al (2019, see their fig. 5a, b) concluded that the six sites with the greatest percentage increase over that predicted from the species-area curve was Pinnacles National Park, California followed in descending order by Grand Staircase Escalante National Monument, Utah, Yosemite National Park, California, Carlsbad Caverns National Park, New Mexico, Clark County, Nevada and the San Bernardino Valley, Arizona/ Sonora. In our study with the new information reported herein, the San Bernardino Valley, Arizona/ Sonora is the highest ranked study followed in descending order by Grand Staircase Escalante National Monument, Utah, Pinnacles National Park, California, Yosemite National Park, California, Clark County, Nevada, and Carlsbad Caverns National Park, New Mexico. With these changes, bee species density from the valley floor samples in the San Bernardino Valley is greater than any other area in North America surveyed to date (Fig. 4b).

Discussion
The San Bernardino Valley of the Chihuahuan Desert is extremely rich for bees. An area of 16 km 2 that varies in elevation by less than 120 m yielded 277 species in    Table S2, Meiners et al. 2019). When the effect of area alone is standardized, the density of bee species in the San Bernardino Valley, is the highest reported from any long-term study of bees in North America. For questions that rely on accurate estimates of species composition and richness, this study highlights the importance of using a variety of sampling protocols and longterm sampling. Our dataset based on samples from one passive sampling protocol (pan traps) in 2001 that recovered more than 21,000 specimens yielded 93 fewer species than we found from samples made over nine years that included pan traps and aerial netting at flowers (Table 3). If we assume that the collections from 2001 pan traps and core area are complete measures of species richness, pan traps underestimated species richness by 38%. The significantly lower species accumulation curve from the 2001 pan trap samples than the long-term dataset from the same area (Fig. 2) is further evidence that surveys based only on pan traps underestimate bee species richness. Biases of pan trap sampling for bees have been well documented (reviewed in Packer and Darla-West 2021) but not from a long-term study in a bee diversity hotspot.
There are numbers of reasons why species richness estimates for desert bees are sensitive to sampling duration that broadly apply to biodiversity studies of most desert insects and bees in general. For one, many species, especially those bees that are specialists of one or several host plants, emerge as adults for a brief period per year, are small-bodied with limited dispersal capabilities and occur patchily (Minckley et al. 1999;Kazenel et al. 2020). Many species are rare: singleton (represented by one specimen) and doubleton (represented by two specimens) species represent 83 (16.7%) and 52 (10.5%) of all species in our full dataset, respectively. Finally, bee populations vary dramatically between years in response to changes in floral resources (Crone 2013) and propensity to remain in diapause . Variation of species abundance taken at the same location across years is often marked (Stubblefield et al. 1993;Meiners et al. 2019). An extreme case occurred in spring 2006 when a drought resulted in the complete absence of 11 of 46 bee species that had been abundant in the San Bernardino Valley in the near-normal rainfall years before and after . Of the absent species, seven were pollen specialists of plants that require rainfall to trigger growth and flower (annual species and creosotebush [Larrea tridentata]), and four were pollen generalists. In contrast, the five pollen specialist bee species of mesquite (Prosopis velutina), a deep-rooted perennial that blooms regardless of winter rainfall, were active during the drought.
The ability to undergo facultative diapause is common among desert bee species (Hurd, 1957;Rust et al. 1988;Danforth et al. 2019) and suggests a bet-hedging life history strategy is a common adaptation for persisting in these areas . Danforth (1999) showed in a laboratory experiment that some Macrotera portalis (Timberlake, 1954), a specialist of mallows (Sphaeralcea spp), emerged in response to moisture in the late-summer (when it is active) but others remained in diapause for one or two years. Desert plant species vary tremendously in their propensity to emerge and bloom (Bowers and Dimmitt 1994;Bowers 2005), even in common gardens (Gremer and Venable 2014;Gremer et al. 2016), which suggests that activity of the specialist bees that depend on different hosts will differ in their activity between years and across habitats. Sampling that is only occasional through the season or widely-spaced across an area may not detect many resident bee species. In the San Bernardino Valley, there is 80% turnover of pollen specialist species every two weeks in the spring (Minckley 2008), suggesting that samples for these and other short-lived solitary bees should be more than bimonthly. The samples in this study made over 9 years and using pan traps and nets minimize these biases and account for the higher estimates of bee species richness than from very intensive pan trapping in one year even when both datasets are from the same location.
The species richness estimates using extrapolation methods from the samples on the valley floor suggest approximately 500 bee species occur in the San Bernardino Valley, of which we captured 497 when all samples in the region were combined (Fig. 3, Table 1). Extrapolations of all samples suggests species richness is 516.1, not substantially more than the number of species observed, and approximately 14% of the 3562 species known for all of North America north of the United States-Mexico border . The number of bee species from the San Bernardino Valley and their relative abundance sets a reasonable expectation for communities in other undisturbed desertscrub/ desert grassland habitats of the Chihuahuan Desert and can be usefully compared to areas with more severe human impacts. For example, a compilation limited to described bee species from southeast Arizona in the American Museum entomology collection and from the San Bernardino Valley yielded 540 bee species when desert scrub, mid-elevation grassland and high montane ecosystems were included (Minckley and Ascher 2013). Considered individually, the three desert areas north of the San Bernardino Valley had between 150-180 species but did not include some genera that had yet to be curated in the American Museum of Natural History collection. Thus, further taxonomic and curatorial work should substantially increase the bee species richness from these well-collected sites. Similarly, a 12-year study of bees on the extreme northern edge of the Chihuahuan Desert in New Mexico (approximately 510 km NNW of the San Bernardino Valley) captured 308 bee species using a modified pan trap design (Kazenel et al. 2020). Because bee species richness is often underestimated by pan traps, the reported number of species suggest this bee fauna is much richer and may be comparable to that in the San Bernardino Valley.
Museum records, taxonomic revisions, and surveys have long suggested that bee diversity runs counter to that of the plants they pollinate in that bees have peak diversity in the western hemisphere in the Chihuahuan, Mojave and Sonoran deserts and in the eastern hemisphere around the Mediterranean Sea (Michener 1979(Michener , 2007. The repeat sampling from a defined area in the Chihuahuan Desert provides empirical data that supports this counterintuitive biogeographical pattern. However, more sampling of bees is needed in tropical and subtropical areas, as well as a more complete picture of herbaceous plant biogeography. Plant biogeography still rests largely on what we know about tree species, and when short-lived plant species are better documented (Massante et al. 2019) their distribution may more closely mirror the distribution of their primary pollinators, the bees. sites over 9 years (the "core" area). 95% confidence intervals calculated from 50 bootstrap replications are not prominent because of their small size. Both samples have greater that 95% coverage and broadly overlap up to their reference sizes (21, 403 individuals for 2001 pan trap samples, 47,120 for core area samples) and when extrapolated to twice the reference sample size. Figure S3. Comparison of the coverage-based rarefaction (solid line) and extrapolation (dashed line), for the pan trap samples at permanent sites in 2001 and the samples made by netting and pan traps throughout the same area encompassed by the permanent sites over 9 years (the "core" area) for all three orders of q. In each of the three graphs the lower line is the 2001 pan trap samples and the upper line is the samples from the core area. Note that the lines do not intersect at any sample coverage. Figure S4. Comparison of the coverage-based rarefaction (solid line) and extrapolation (dashed line), for the samples made in the core area (middle line), core area combined with the pan trap sampling at permanent sites (blue line on bottom), and all samples made in the San Bernardino Valley (top line) at q =0. The confidence intervals are very small they do not show on the graph. Sample coverage is above 95% for diversity values below 200 species for the datasets and the extrapolations overlap broadly at higher sampling coverage indicating all three datasets predict very similar bee species richness in the San Bernardino Valley. Note that the analyses at other orders of q are not shown because the results match that shown here. Copyright notice: This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited. Link: https://doi.org/10.3897/jhr.82.60895.suppl1