Unique Bee Communities within Vacant Lots and Urban Farms Result from Variation in Surrounding Urbanization Intensity

We investigated the relative importance of vacant lot and urban farm habitat features and their surrounding landscape context on bee community richness, abundance, composition, and resource use patterns. Three years of pan trap collections from 16 sites yielded a rich assemblage of bees from vacant lots and urban farms, with 98 species documented. We collected a greater bee abundance from vacant lots, and the two forms of greenspace supported significantly different bee communities. Plant–pollinator networks constructed from floral visitation observations revealed that, while the average number of bees utilizing available resources, niche breadth, and niche overlap were similar, the composition of floral resources and common foragers varied by habitat type. Finally, we found that the proportion of impervious surface and number of greenspace patches in the surrounding landscape strongly influenced bee assemblages. At a local scale (100 m radius), patch isolation appeared to limit colonization of vacant lots and urban farms. However, at a larger landscape scale (1000 m radius), increasing urbanization resulted in a greater concentration of bees utilizing vacant lots and urban farms, illustrating that maintaining greenspaces provides important habitat, even within highly developed landscapes.


Introduction
Urban greenspaces, once disregarded as homogeneous habitat dominated by a few, often invasive, species, are increasingly considered as potential habitat for species conservation [1,2].Part of their potential stems from the fact that urban greenspaces often support a high species richness that can include rare plants and animals [3][4][5].Urban agroecosystems, which contain a diversity of plants reliant on insect-mediate ecosystem services like pollination, are especially promising habitats for pollinator conservation [6,7].These agroecosystems are often transformed from vacant lots, which can themselves be valuable habitats for bees [2] and biodiversity more broadly [6,[8][9][10][11].Urban agroecosystems and vacant lots differ in their management and floral composition, both of which are important drivers of species abundance and richness [11,12].At the same time, both of these habitat types are located within the urban landscape that can act as a strong environmental filter for communities [13].Understanding how both the landscape context and the local design and management decisions ongoing within vacant lots and urban agroecosystems act to filter a city's bee community is key to incorporating conservation into urban greenspace planning [14,15].
The landscape structure surrounding an urban greenspace can strongly influence conservation-driven management outcomes.Habitat fragmentation and loss that result from urbanization are known contributors to patterns of species richness and abundance across rural to urban gradients [16,17].At the same time, landscape features at the city-wide scale can act as strong environmental filters, influencing if taxa from a given species pool are able to colonize distinct greenspace patches [12].For instance, isolation from other greenspaces can lead to reduced bee species richness within urban forest fragments [18] and connectivity with other greenspaces positively predicts bee abundance on green roofs [19].In some cases, landscape variables also influence the functional trait distribution found within groups [17,20].For instance, highly fragmented urban landscapes tend to favor smaller species and cavity nesting species [21,22].
At a local scale, species interactions among bees and floral resources can drive patterns of pollinator community assembly [12].Increasing bee species richness is often linked to floral species richness [23], which frequently increases with the degree of urbanization [4,24,25].Within a city, urban greenspaces offer a varied breadth and quality of dietary resources [12,[26][27][28][29]. Exotic plant species can make up a substantial portion of the floral resources found in urban areas [30].Although some exotic species common within vacant lots have been found to be highly attractive [31,32], these resources may also be most suitable for polylectic [33] and/or exotic bees [31].Urban agroecosystems are likely to have a wider breadth of floral resources, including food crops, native and exotic ornamentals, and many of the same weedy species found in vacant lots.Adding flowering resources, principally native species aimed at supporting pollinators and other beneficial insects, increases available nectar and pollen resources [34] and results in a greater richness and abundance of pollinators within a patch for some [35,36], but not all [37] cases.
Herein, we examine the role of landscape structure on the taxonomic and functional trait distributions of bees among urban greenspaces, and evaluate how the design and management of these habitats influences the richness, abundance, dietary breadth and niche overlap within bee assemblages.Our study was conducted within the city of Cleveland, Ohio, a shrinking city that experienced significant population decline resulting in an overabundance of infrastructure, which was demolished over time resulting in vacant land.Worldwide, over 450 cities have experienced significant population decline, including 59 within the United States [38].Many of these cities have large holdings of reclaimed greenspace, for instance, Cleveland currently encompasses over 27,000 vacant lots that totals more than 1500 ha.[39].While the spatial extent of vacant land within shrinking cities represents a potentially important conservation resource, the relative importance of local habitat management and landscape context for structuring bee assemblages is largely unknown.Therefore, our objectives were to: (1) compare the richness, abundance and composition of bees found within vacant lots and urban agroecosystems, (2) determine if the intensity of urbanization surrounding these sites influenced bee assemblages and (3) examine how the conversion of vacant land to urban agriculture influenced foraging patterns, including the dietary breadth and overlap of bee species.Our hypotheses were that (i) urban agroecosystems would support a greater richness of pollinator species as a result of enhanced dietary niche breadth and reduced overlap, and (ii) that densely urbanized landscapes would supply a reduced taxonomic richness, higher proportion of exotic species, reduced average body size, and reduced proportion of ground nesting species relative to cavity nesting species to vacant lots and urban farms.

Methods
Study sites: We selected eight vacant lots and eight urban community gardens, educational and for-profit farms (hereafter referred to collectively as urban farms) within the city of Cleveland, OH, USA for this study in 2010.Site selection was constrained by the availability of vacant lots and the willingness of urban farm managers to participate.Vacant lots were parcels seeded with a fescue grass seed mixture following the demolition of a single-family structure.These properties were managed by periodic mowing, approximately once per month, by the City of Cleveland Land Bank.Urban farms consisted of six sites managed as learning farms by the Cleveland Botanical Garden, one for-profit farm, and one non-profit community garden.All urban farm sites had been converted from vacant lots and managed for food production for at least two years prior to the beginning of this study.No pesticides were used within any vacant lot or urban farm for the duration of the study.Property areas of vacant lots ranged from 843 m 2 to 6951 m 2 and from 271 m 2 to 4614 m 2 for urban farms.The distance between neighboring sites varied from 145 m to 1856 m.Two vacant lots were replaced after the 2010 sampling season: one due to construction of a park on the site, and the other due to repeated vandalism of sampling equipment.One urban farm site was replaced following 2011 data collection due to repeated vandalism of sampling equipment.
Pollinator richness and abundance: We collected bees three times per year during 2010-2012 using pan traps on calm, sunny days.Collections occurred in early, mid, and late summer.Following Droege [40], pan traps consisted of 96 mL white soufflé cups (Solo cup company, Lake Forest, IL, USA), which we kept white or painted either fluorescent blue or fluorescent yellow (Guerra Paint and Pigment, New York, NY, USA) [40].We divided sites into four, equally-sized plots and placed one of each color pan trap at the center of each plot 1 m apart and filled with a 1% dish soap solution (Dawn© Proctor and Gamble, Cincinnati, OH, USA).After 48 h, we collected the traps and returned them to the laboratory, where we washed individuals in water and dried them with a blow dryer.We pinned and identified bees to the species level using the Discover Life website [41].All specimen identifications were verified by either Dr. Jason J. Gibbs  Bee functional traits: To each individual collected in a pan trap that could be identified to species, we assigned three functional traits: origin, nesting guild, and body length.The origin of each individual was classified as either native or exotic to the United States.Nesting guild assignment was modified from [42] to include: CARPENTER (excavator in woody substrate), MINER (excavator in the ground), NONEST (cleptoparasitic), PITH (nests in pithy stems of plants), and RENT (nests in pre-existing cavities).We used the Discover Life website [41] to determine a mean body length (we calculated the mean of the reported male and female body length) for each species.Measurements were not available for 15 species collected as part of our study, so we measured collected specimens (five individuals per species, unless fewer were collected) and calculated a mean body length.
Resource use: To determine patterns of resource use, we measured floral visitation in 2011 and 2014 within the urban farms and vacant lot sites.Depending on weather conditions, we visited site on a weekly or bi-weekly basis from 20th May to 14th August in 2011, and from 3rd June to 5th August in 2014.As urban farms and vacant lots are heterogeneous habitats, we felt that timed plant species-specific sampling was more appropriate than using transects [43].Thus, we identified all flowering plant species in bloom within a site and observed blooms of each species for 10 min.For each plant species, the observer either monitored blooms on one or more plants within a 1 m 2 area for 10 min, or if multiple patches were present, the observation time was split and two or more 1 m 2 patches were observed.During each observation period, we collected any bee that landed on an open flower within the observation area using "insect vacuums" (BioQuip Products, Compton, CA, USA).We then pinned and identified each individual in the same manner as described for specimen collected in pan traps.
Landscape classification: To examine the effect of landscape structure on bee assemblages and patterns of floral resource use, we quantified the landscape surrounding each site using remotely sensed images at 1-2 m resolution [44][45][46].We used eCognition Developer 8.7 (Trimble, Westminster, CO, USA) and ArcGIS 10.1 (ESRI, Redlands, CA, USA) for land cover identification and landscape composition analysis, respectively [47,48].We classified land cover into four categories: buildings, pavement, fine vegetation (grasses and low-growing forbs), and coarse vegetation (trees and woody shrubs).From these data, we calculated the following predictor variables: proportion impervious surface (area of buildings and pavement), number of buildings, and number of distinct greenspace patches.We measured these variables within circular buffers with radii of 100 m and 1000 m surrounding each site.We selected these scales to measure how local land use (100 m) as well as larger-scale landscape heterogeneity and urbanization (1000 m) influence bee species assemblages.

Statistical analysis:
We used generalized linear models to examine the relationship between bee communities and local and landscape factors in R 3.4.3[49].All models initially included the fixed effect of habitat type (urban farm or vacant lot), and the following predictor variables at both the local (100 m) and landscape (1000 m) scale: the proportion of impervious surface, the number of buildings, and the number of distinct greenspace patches.Collinearity analyses revealed a strong correlation between the number of buildings at 1000 m and other landscape variables.As a result, we removed this variable from the analysis.All models also included interactions between habitat type and each of the landscape variables.We pooled data across years and collection time periods.Additionally, because the number of pan traps collected at each site varied due to loss and vandalism (median = 78, min = 45, max = 94 collected), we included the number of pan traps per site as an offset in all models.Models evaluating count data (e.g., abundance and richness) were first fit with a Poisson distribution.If we detected overdispersion, we refit models using a negative binomial generalized linear model with the glm.nb function.Models analyzing mean body length were assessed for normality and analyzed using the glm function.We then used backwards model selection with the step function to determine the optimal fixed effects structure for each model.Finally, if we detected a significant interaction between habitat type and a continuous predictor variable, we evaluated whether the slope for either habitat type was significantly different from zero by manually computing the slopes, standard errors, and resulting p-values.
To determine if resource use differed between habitat types and across landscape structure, we used generalized linear mixed models in the lme4 package [50] to analyze the abundance and richness of the bee communities visiting flowers at each site.Models included the same fixed effects as those used to analyze pan trap data, and we included the time observing a particular floral resource as an offset and site as a random effect.Similar to the pan trap analysis, we first fit models with a Poisson distribution using the glmer function and then with a negative binomial distribution using the glmer.nbfunction if we detected overdispersion.We determined the optimal fixed effects structure for each model using the drop1 function for backwards model selection and assessed the significance of interaction slopes by manually computing the slopes, standard errors, and resulting p-values for both habitat types.We analyzed differences in community composition using non-metric multidimensional scaling and analysis of similarities using the metaMDS and adonis functions in the vegan package [51].
Network analysis: To describe the plant-pollinator networks at each site, we used the networklevel function in the bipartite package [52,53] to calculate the following indices: (i) generality: the mean number of plant species per pollinator, and (ii) niche overlap: the mean similarity in interaction patterns between pollinators.Because network indices are sensitive to web dimension and sampling intensity [53], we calculated 1000 null models for each web using the nullmodel function with the 'r2dtable' method [52] and corrected our index values as z-scores.For a given index, the z-score was calculated as the raw score minus the mean index score from the null models divided by the standard deviation of the null models.As in Dormann et al. [54], we then determined whether a given score was significantly different from the mean null model z-score using the number of cases where a null model's score was greater than or equal to the observed value.Finally, we used t-tests to compare corrected indies from vacant lots and urban farms to determine if patterns in these networks differed between habitat types.

Results
We collected a total of 2733 bees from pan traps, identifying 2671 of those individuals to species.Overall, we collected bees from 98 species across five families (Table 1).The five most abundant species collected were Agapostemon virescens, Halictus ligatus, Ceratina calcarata, Halictus confusus, and Melissodes bimaculatus.The majority of individuals found occupying urban farms and vacant lots were in the MINER nesting guild, representing 82% of all individuals (Table 1).Insufficient numbers and representation by the CARPENTER (<1% of individuals), NONEST (<1% of individuals), and PITH (6% of all individuals) nesting guilds across sites prompted us to remove these guilds from further analysis.Exotic species represented 8% of the collected bee fauna, with 208 individuals belonging to 10 species (Table 1).These species included the first known North American collection of Hylaeus pictipes; one specimen was collected from a vacant lot site annually in 2010, 2011, and 2012.Nest parasites represented 14% of collected species with 14 unique species found.
Table 1.Species list of the bee fauna collected from pan traps within urban farms and vacant lots.We list the total number of specimens of each species collected over the 2010-2012 sampling period as well as the assigned nesting guild (modified from [41]), status in North America, and mean body size.Bee composition, richness, and abundance: The composition of bee species collected from pan traps within vacant lots versus urban farms was significantly different (F = 2.06, p = 0.001, Figure 1).We found that total abundance (z = −2.16,p = 0.03) and native abundance (z = −2.72,p = 0.02) of bees collected per pan trap was greater in vacant lots versus urban farms.Bee abundance (z = 3.32, p < 0.001) and native bee abundance (z = 3.11, p = 0.002) were also negatively correlated with the proportion of impervious surface present within a 100 m radius of our study sites.The opposite pattern was found at 1000 m, with bee abundance (z = 2.41, p = 0.02) and native bee abundance (z = 2.12, p = 0.03) increasing with the proportion of impervious surface in the surrounding landscape.

Family
We found significant interactions between habitat type and landscape variables for exotic bee abundance and bee species richness.For exotic bee abundance, we found a significant interaction between habitat type and the number of greenspaces (z = −2.53,p = 0.01) at a 1000 m radius.Exotic bee abundance in vacant lots was negatively correlated with the number of greenspace patches present at 1000 m (t = −3.06,p = 0.01), but no such relationship existed for urban farms (t = 0.85, p = 0.42) (Figure 2A).We also found a significant interaction between habitat type and the proportion of impervious surface (z = −3.48,p = 0.001), where exotic bee abundance in urban farms (t = 3.60, p = 0.01), but not vacant lots (t = −1.49,p = 0.17), was positively correlated to the proportion of impervious surface present at 1000 m (Figure 2B).Exotic bee abundance within both habitats was also positively correlated with the number of buildings at 100 m (z = 2.04, p = 0.04).For bee species richness, we found a significant interaction between habitat and the proportion of impervious surface at 1000 m (z = −3.33,p = 0.001), with richness increasing significantly in urban farms (t = 2.45, p = 0.03), but not in vacant lots (t = −2.17,p = 0.06), with greater impervious surface in the surrounding landscape at a radius of 1000 m (Figure 2C).Within both vacant lots and urban farms, the number of greenspace patches in the surrounding landscape was positively correlated with bee richness at a 100 m radius (z = 2.09, p = 0.04) and negatively correlated at 1000 m (z = −2.02,p = 0.04).
in vacant lots (t = −2.17,p = 0.06), with greater impervious surface in the surrounding landscape at a radius of 1000 m (Figure 2C).Within both vacant lots and urban farms, the number of greenspace patches in the surrounding landscape was positively correlated with bee richness at a 100 m radius (z = 2.09, p = 0.04) and negatively correlated at 1000 m (z = −2.02,p = 0.04).Bee functional trait distributions: We found a marginally significant difference in the abundance of bees in the MINER nesting guild among habitats (z = −1.80,p = 0.07), with these ground nesting species being less frequently collected from urban farms.Within both vacant lots and urban farms, the abundance of MINER bees decreased with the proportion of impervious surface within 100 m and increased with the proportion of impervious surface within 1000 m of research sites.For the RENT nesting guild, we found a significant interaction between habitat type and the proportion of impervious surface at 1000 m (z = −2.66,p = 0.01), with the abundance of the RENT guild within urban farms increased significantly with the proportion of impervious surface (t = 2.40, p = 0.05).We did not detect a relationship between RENT abundance within vacant lots and the proportion of impervious surface at 1000 m (t = −1.27,p = 0.25) (Figure 2D).For mean body size, we found a significant interaction between habitat and the number of buildings within 100 m of a focal site (z = −2.65,p = 0.04), with an increasing relationship found in urban farms (t = 2.98, p = 0.02) and no relationship found in vacant lots (t = 0.66, p = 0.54) (Figure 2E).Bee functional trait distributions: We found a marginally significant difference in the abundance of bees in the MINER nesting guild among habitats (z = −1.80,p = 0.07), with these ground nesting species being less frequently collected from urban farms.Within both vacant lots and urban farms, the abundance of MINER bees decreased with the proportion of impervious surface within 100 m and increased with the proportion of impervious surface within 1000 m of research sites.For the RENT nesting guild, we found a significant interaction between habitat type and the proportion of impervious surface at 1000 m (z = −2.66,p = 0.01), with the abundance of the RENT guild within urban farms increased significantly with the proportion of impervious surface (t = 2.40, p = 0.05).We did not detect a relationship between RENT abundance within vacant lots and the proportion of impervious surface at 1000 m (t = −1.27,p = 0.25) (Figure 2D).For mean body size, we found a significant interaction between habitat and the number of buildings within 100 m of a focal site (z = −2.65,p = 0.04), with an increasing relationship found in urban farms (t = 2.98, p = 0.02) and no relationship found in vacant lots (t = 0.66, p = 0.54) (Figure 2E).
There was no difference in the breadth of floral visitation among habitats (generality, t = −1.27,p = 0.23) or the extent of overlap in resource use (niche overlap, t = 0.59, p = 0.57).The number of bees observed visiting floral resources per 10-min observation period did not differ among urban farms and vacant lots and was not influenced by landscape variables.In contrast, both overall species richness and native species richness were influenced by landscape variables.As these responses were virtually identical, we report only the results from native species richness.Within both vacant lots and urban farms, species richness of native bee floral visitors per 10 min was negatively correlated with the proportion of impervious surface in the surrounding landscape at 100 m (z = −2.07,p = 0.04).Each species is represented by a black rectangle, and the width of that rectangle illustrates that relative frequency that a particular species interacts with a member of the opposite trophic level.A list of the full species names is provided in Appendix A: Tables A1 and A2.There was no difference in the breadth of floral visitation among habitats (generality, t = −1.27,p = 0.23) or the extent of overlap in resource use (niche overlap, t = 0.59, p = 0.57).The number of bees observed visiting floral resources per 10-min observation period did not differ among urban farms and vacant lots and was not influenced by landscape variables.In contrast, both overall species richness and native species richness were influenced by landscape variables.As these responses were virtually identical, we report only the results from native species richness.Within both vacant lots and urban farms, species richness of native bee floral visitors per 10 min was negatively correlated with the proportion of impervious surface in the surrounding landscape at 100 m (z = −2.07,p = 0.04).We found a significant interaction between habitat type and the number of greenspaces in the surrounding landscape at 100 m (z = 2.32, p = 0.02), with a significant negative relationship between the number of greenspaces and native visitor species richness in urban farms (t = −2.27,p = 0.04) but no such relationship existed for vacant lots (t = 0.49, p = 0.63) (Figure 4).There was also a significant interaction between habitat and the number of greenspaces at 1000 m (z = −2.65,p = 0.01), with the native species richness of bees collected per 10 min in urban farms increasing with the number of greenspaces (t = 2.98, p = 0.01) but not in vacant lots (t = 0.10, p = 0.92) (Figure 4).

Discussion
A meaningful shift to value historically overlooked urban habitats as resources for conservation has occurred [2,55].This is especially true within shrinking cities, where thousands of minimallymanaged vacant lots present a vast opportunity to conserve species and deliver needed ecosystem services such as local cooling, carbon sequestration, storm water infiltration and local food production [2,6,10,56].Ideally, managing urban greenspaces for species conservation promotes the provisioning of ecosystem services, but in reality this is hardly a causal relationship.Understanding how to design urban greenspaces to achieve these distinct goals presents a challenge to conservation practitioners [57].In this study, we examined how transforming vacant land into urban agroecosystems influenced bee assemblages and their foraging patterns.Furthermore, we measured if surrounding urbanization both locally (100 m radius) and at a larger landscape scale (1000 m radius) altered the conservation potential of these two types of urban greenspaces.
Urban greenspaces have repeatedly demonstrated an ability to support a high bee richness [12,24,35,58,59], and arthropod assemblages often vary among distinct forms of greenspace within an urban landscape [11,[59][60][61].In total, we documented 98 bee species, with the MINER nesting guild dominating both habitat assemblages.We also found 14 unique species within the NO NEST guild (14% of bee species); these nest parasites can be an indicator of habitat quality and their presence has been noted to promote community stability by influencing competition among hosts and reducing dominance of common species [62].We found that vacant lots and urban farms supported significantly different bee assemblages, with 14 species unique to urban farms and 29 species collected solely from vacant lots.We also found a higher abundance of bees within vacant lots versus urban farms.This was due in part to a marginally significant reduction in the MINER guild within urban farm sites, 66.7% fewer individuals from 34 species were collected relative to vacant lots.Interestingly, communal and eusocial species [63][64][65] such as A. virescens, H. ligatus, H. confusus, and L. leucocomum were among the most common MINER species collected within vacant lots, and a 36-65% reduction in the collection of these species was found within urban farms.Frequent soil disturbances resulting from tillage, weeding, irrigation, and mulching practices within urban farms may discourage the establishment and impact the survival of mining species in general and could have a particularly negative influence on communal and eusocial species.Given the nesting behavior of these species, small-scale disturbances such as tillage of a single bed could remove several females

Discussion
A meaningful shift to value historically overlooked urban habitats as resources for conservation has occurred [2,55].This is especially true within shrinking cities, where thousands of minimally-managed vacant lots present a vast opportunity to conserve species and deliver needed ecosystem services such as local cooling, carbon sequestration, storm water infiltration and local food production [2,6,10,56].Ideally, managing urban greenspaces for species conservation promotes the provisioning of ecosystem services, but in reality this is hardly a causal relationship.Understanding how to design urban greenspaces to achieve these distinct goals presents a challenge to conservation practitioners [57].In this study, we examined how transforming vacant land into urban agroecosystems influenced bee assemblages and their foraging patterns.Furthermore, we measured if surrounding urbanization both locally (100 m radius) and at a larger landscape scale (1000 m radius) altered the conservation potential of these two types of urban greenspaces.
Urban greenspaces have repeatedly demonstrated an ability to support a high bee richness [12,24,35,58,59], and arthropod assemblages often vary among distinct forms of greenspace within an urban landscape [11,[59][60][61].In total, we documented 98 bee species, with the MINER nesting guild dominating both habitat assemblages.We also found 14 unique species within the NO NEST guild (14% of bee species); these nest parasites can be an indicator of habitat quality and their presence has been noted to promote community stability by influencing competition among hosts and reducing dominance of common species [62].We found that vacant lots and urban farms supported significantly different bee assemblages, with 14 species unique to urban farms and 29 species collected solely from vacant lots.We also found a higher abundance of bees within vacant lots versus urban farms.This was due in part to a marginally significant reduction in the MINER guild within urban farm sites, 66.7% fewer individuals from 34 species were collected relative to vacant lots.Interestingly, communal and eusocial species [63][64][65] such as A. virescens, H. ligatus, H. confusus, and L. leucocomum were among the most common MINER species collected within vacant lots, and a 36-65% reduction in the collection of these species was found within urban farms.Frequent soil disturbances resulting from tillage, weeding, irrigation, and mulching practices within urban farms may discourage the establishment and impact the survival of mining species in general and could have a particularly negative influence on communal and eusocial species.Given the nesting behavior of these species, small-scale disturbances such as tillage of a single bed could remove several females and their offspring from a local population.Furthermore, the extended gap in time between foundress and worker foraging periods that exists for eusocial species [63] may make these species more likely to be active and negatively impacted when crop management practices occur at their nest site.
We observed floral visitation within both habitats and did not find evidence to support our prediction that differences in bee richness and abundance observed among vacant lots and urban farms was due to variation in dietary niche breadth or overlap.We found no difference in the breadth of floral visitation per bee species (as measured by "generality") or the extent of overlap in resource use among species (described by "niche overlap") within urban farms versus vacant lots.The average number of bees visiting plants during our observation periods was also consistent across habitats.However, vacant lots and urban farms offered different floral resources, possibly contributing to the variation found in bee community composition we observed.For instance, the composition of common foragers utilizing plant resources varied, with only Apis mellifera shared among the three most common species observed visiting flowers in vacant lots and urban farms.Our findings also highlight that a diversity of plant species can contribute to bee conservation in cities. Red and white clover (T.pretense and T. repens) were the most attract plants within vacant lots, demonstrating the role that weedy exotic plants can play in supporting an urban bee species pool [30,32,66] in self-assembled, minimally-managed greenspaces.In urban farms, the most frequently visited plants per observation period included both harvestable crops and ornamental flowers.These findings support recommendations directed at large-scale rural agroecosystems, suggesting that bee conservation in agricultural habitats could be aided by practices such as the establishment of additional floral resources and relaxed weed control measures [67].Furthermore, variation in visitation frequency among planted ornamental species found in our study and others [26,27,68,69] highlights the importance of plant selection toward the conservation value of gardens and urban agroecosystems.In our study, black-eyed Susan (Rudbeckia hirta), oxeye daisy (Leucanthemum vulgare), and purple coneflower (Echinacea purpurea) received the most visitors per 10 min observation.Measuring bee visitation rates to inform conservation focused plantings provides important information, for instance, Garbuzov et al. [27] found a 1000-fold difference in attractiveness among plants surveyed in public urban gardens.Similarly, among the 168 plants surveyed in allotment gardens, 50% of all bumble bee visits were to just 14 species [68].
We predicted that urbanization intensity would negatively impact the abundance and taxonomic richness of bees found within urban greenspaces, which was based on evidence that isolation from natural habitat sources can reduce pollinator abundance, richness, and the provision of pollination services [70,71].We found support for this prediction within the localized 100 m radius surrounding a site, with bee abundance and richness generally declining with the proportion of impervious surface and increasing with the number of greenspaces patches surrounding urban farms and vacant lots.The richness of bees observed visiting plants within urban farms and vacant lots also decreased with more localized impervious surface present.Together, these findings suggest that patch isolation due to a high concentration of localized impervious surface may limit greenspace colonization by bees.
However, many of our findings, particularly at the 1000 m radius scale, are in contrast to our hypothesis that densely urbanized landscapes supply a reduced taxonomic richness, higher proportion of exotic species, reduced average body size, and reduced proportion of ground nesting species.Increased impervious surface within 1000 m of a focal site was positively associated with total bee, native bee, and MINER nesting guild abundance.Total species richness, exotic bee abundance, and RENT nesting guild abundance within urban farms were also positively related to urbanization at a 1000 m radius scale and average bee body length increased with a greater concentration of buildings surrounding urban farms at 100 m.These findings are contrary to most previous studies who have either found no strong influence of landscape composition [7,37,72] or that amount of impervious surface surrounding greenspaces is negatively correlated with the abundance, and richness of bee species, with ground nesting bees being particularly sensitive [12,73,74].Furthermore, landscapes with increasing proportions of urban land typically use support assemblages with a smaller average body size [21].Cavity nesting species are one exception, which have been found to increase with urbanization [74,75] potentially due to human-derived nesting resources such as fences and buildings.Our results are consistent with these findings as RENT nesting guild abundance was positively correlated with the proportion of impervious surface at 1000 m within urban farm sites.
For most bee species, however, it is unlikely that urbanization intensity actually aids bee conservation.Instead, the large-scale landscape patterns found in our study are more likely due to a concentration effect, wherein a greater proportion of the landscape-scale species pool relies on a specific greenspace patch as a foraging habitat when fewer sites are available [76].This phenomenon may be most likely to occur in shrinking cities, where the prevalence of minimally-managed vacant lots can maintain a species pool able to concentrate within available greenspace even when the proportion of the landscape composed of impervious surface is high.For instance, at 1000 m, our most densely developed landscape was comprised of 71% impervious surface, and vacant land made up a considerable portion of the remaining 29% of available greenspace, which may provide enough habitat to support many urban-adapted bee species.This pattern warrants further study but suggests that, even within highly urbanized landscapes, maintaining vacant land provides a functional habitat with conservation value for bees in shrinking cities.

Conclusions
In conclusion, our study adds to the growing body of literature advocating for the maintenance of minimally-managed vacant lot habitats as a conservation resource [6,[8][9][10][11].These sites support a high pollinator richness and abundance and require relatively little financial investment.Furthermore, a unique bee community is found within vacant lots versus urban farms, illustrating the value in diversifying the design and management of greenspaces in cities to maximize conservation gains and the delivery of ecosystem services.Maintaining vacant lot greenspaces may be most important for bees within highly urbanized landscapes where few foraging patches exist; however, patches surrounded locally by high concentrations of impervious surface may limit colonization and thus not be ideal spaces for focused conservation investment.Finally, the strength and relationship of landscape effects can be habitat dependent; in several cases, we found a significant influence of landscape variables on assemblages within urban gardens, but not vacant lots.This supports the idea that landscape context influences the impact of local-scale management on insect community dynamics [77], a hypothesis that has not been supported in traditional agricultural systems [78].

Figure 1 .
Figure 1.Non-metric multidimensional scaling ordination of the community composition of bees collected from urban farm and vacant lot greenspaces.Stress of the final configuration was acceptable (stress = 0.19) and the ordination was plotted in two dimensions.

Figure 1 .
Figure 1.Non-metric multidimensional scaling ordination of the community composition of bees collected from urban farm and vacant lot greenspaces.Stress of the final configuration was acceptable (stress = 0.19) and the ordination was plotted in two dimensions.

Figure 2 .
Figure 2. Habitat specific effects of landscape variables.Relationship between exotic bee abundance and (A) the number of greenspace patches, or (B) the proportion of impervious surface within 1000 m of a focal site.Relationship between (C) bee species richness and the proportion of impervious surface within 1000 m of a focal site, (D) the abundance of the RENT nesting guild and the proportion of impervious surface within 1000 m of a focal site, and (E) the mean body length and the number of buildings within a 100 m radius.Best fit lines (±95% confidence interval) for urban farms are represented in red and vacant lots are represented in blue.An asterisk next to the best fit line indicates that the slope of that line is significantly (* p < 0.05; ** p < 0.01) different from zero.Bee foraging patterns: Across 2011 and 2014, the number of bee species observed visiting available floral resources ranged from 5-22 per site.Plants with the most visits per 10-min observation were milkweed (Asclepias syriaca; mean = 2.0 visits/10 min), chives (Allium schoenoprasum; mean = 1.8 visits/10 min), and squash (Cucurbita pepo; mean = 1.8 visits/10 min) in urban farms and red clover (Trifolium pretense; mean = 2.4 visits/10 min), white clover (Trifolium repens; mean = 2.2 visits/10 min), and Queen Anne's lace (Daucus carota; mean = 2.1 visits/10 min) in vacant lots (Figure3).The most commonly observed bees were squash bees (Peponapis pruinosa; mean = 1.9 observations/10 min), honey bees (Apis mellifera; mean = 1.4 observations/10 min), and European wool carder bees (Anthidium.manicatum;mean = 1.2 observations/10 min) in urban farms and common eastern bumble bees (Bombus impatiens; mean = 2.3 observations/10 min), a species of mining bee (Calliopsis

Figure 2 .
Figure 2. Habitat specific effects of landscape variables.Relationship between exotic bee abundance and (A) the number of greenspace patches, or (B) the proportion of impervious surface within 1000 m of a focal site.Relationship between (C) bee species richness and the proportion of impervious surface within 1000 m of a focal site, (D) the abundance of the RENT nesting guild and the proportion of impervious surface within 1000 m of a focal site, and (E) the mean body length and the number of buildings within a 100 m radius.Best fit lines (±95% confidence interval) for urban farms are represented in red and vacant lots are represented in blue.An asterisk next to the best fit line indicates that the slope of that line is significantly (* p < 0.05; ** p < 0.01) different from zero.Sustainability 2018, 10, x FOR PEER REVIEW 9 of 17andreniformis; mean = 2.0 observations/10 min), and honey bees (Apis mellifera; mean = 1.5 observations/10 min) in vacant lots (Figure3).

Figure 3 .
Figure 3. Quantitative plant-pollinator networks in (A) vacant lots and (B) urban farms.Pollinators are represented along the top, plants along the bottom.Each species is represented by a black rectangle, and the width of that rectangle illustrates that relative frequency that a particular species interacts with a member of the opposite trophic level.A list of the full species names is provided in Appendix A: TablesA1 and A2.

Figure 3 .
Figure 3. Quantitative plant-pollinator networks in (A) vacant lots and (B) urban farms.Pollinators are represented along the top, plants along the bottom.Each species is represented by a black rectangle, and the width of that rectangle illustrates that relative frequency that a particular species interacts with a member of the opposite trophic level.A list of the full species names is provided in Appendix A: TablesA1 and A2.

Figure 4 .
Figure 4. Habitat specific relationship between native species richness and the number of greenspaces within (A) 100 m and (B) 1000 m radii.Predictor variables were scaled to ensure model convergence.Best fit lines (±95% confidence interval) for urban farms are represented in red and vacant lots are represented in blue.An asterisk next to the best fit line indicates that the slope of that line is significantly (* p < 0.05; ** p < 0.01) different from zero.

Figure 4 .
Figure 4. Habitat specific relationship between native species richness and the number of greenspaces within (A) 100 m and (B) 1000 m radii.Predictor variables were scaled to ensure model convergence.Best fit lines (±95% confidence interval) for urban farms are represented in red and vacant lots are represented in blue.An asterisk next to the best fit line indicates that the slope of that line is significantly (* p < 0.05; ** p < 0.01) different from zero.
(University of Manitoba) or Dr. Sam Droege (USGS Patuxent Wildlife Research Center).We deposited voucher specimens in the C.A. Triplehorn Insect Collection within the Museum of Biological Diversity at The Ohio State University.

Table A2 .
List of species codes and corresponding scientific names for plant species.