Spatio-temporal dynamics of landscape use by the bumblebee Bombus pauloensis (Hymenoptera: Apidae) and its relationship with pollen provisioning

Caterpillar movement mediates spatially local interactions and determines the relationship between population density and contact

BACKGROUND

While interactions in nature are inherently local, ecological models often assume homogeneity across space, allowing for generalization across systems and greater mathematical tractability. Density-dependent disease models are a prominent example of models that assume homogeneous interactions, leading to the prediction that disease transmission will scale linearly with population density. In this study, we examined how the scale of larval butterfly movement interacts with the resource landscape to influence the relationship between larval contact and population density in the Baltimore checkerspot (Euphydryas phaeton). Our study was inspired by the recent discovery of a viral pathogen that is transmitted horizontally among Baltimore checkerspot larvae.

METHODS

We used multi-year larvae location data across six Baltimore checkerspot populations in the eastern U.S. to test whether larval nests are spatially clustered. We then integrated these spatial data with larval movement data in different resource contexts to investigate whether heterogeneity in spatially local interactions alters the assumed linear relationship between larval nest density and contact. We used Correlated Random Walk (CRW) models and field observations of larval movement behavior to construct Probability Distribution Functions (PDFs) of larval dispersal, and calculated the overlap in these PDFs to estimate conspecific contact within each population.

RESULTS

We found that all populations exhibited significant spatial clustering in their habitat use. Subsequent larval movement rates were influenced by encounters with host plants and larval age, and under many movement scenarios, the scale of predicted larval movement was not sufficient to allow for the "homogeneous mixing" assumed in density dependent disease models. Therefore, relationships between population density and larval contact were typically non-linear. We also found that observed use of available habitat patches led to significantly greater contact than would occur if habitat use were spatially random.

CONCLUSIONS

These findings strongly suggest that incorporating larval movement and spatial variation in larval interactions is critical to modeling disease outcomes in E. phaeton. Epidemiological models that assume a linear relationship between population density and larval contact have the potential to underestimate transmission rates, especially in small populations that are already vulnerable to extinction.

Utility of carbon and nitrogen stable isotopes for inferring wild bee (Hymenoptera: Apoidea) use of adjacent foraging habitats

Isotope analysis has proven useful for understanding diets of animals that are difficult to track for extended periods. Bees are small yet highly mobile and often forage from multiple habitats. However, current methods of assessing diet are limited in scope. Efficient methods of tracking bee diets that integrate across life stages, distinguish habitat use, and are sensitive to taxonomic differences will inform conservation strategies. We evaluated the utility of stable isotope analysis for estimating contributions of adjacent habitats to bees’ diets. We also investigated taxonomic variation in bee and flower isotope composition. We measured natural abundance of carbon and nitrogen stable isotopes in two body regions from three wild bee genera, as well as in 25 species of flowers that likely comprised their diets. Bee ∂13C and ∂15N varied with habitat and taxonomic groups (conflated with month), but did not match spatial or seasonal trends in their food plants. Flower ∂13C was lowest in the forest and in April–June, as expected if driven by water availability. However, bee ∂13C was elevated in the spring, likely from overwintering nutritional stress or unpredictable food availability. Bumble bees (Bombus) were enriched in ∂15N compared to others, possibly reflecting differences in larval feeding. Bee diet mixing models had high variation and should be interpreted with caution. Models estimated similar habitat contributions to diets of spring Andrena and overwintered Bombus queens. Summer Bombus queens and workers were indistinguishable. Sweat bees (Halictus) were estimated to use comparatively more field flowers than others. Overall, taxon more strongly influenced isotope composition than either foraging habitat or month, likely because of associated differences in sociality and timing of annual activity. Future studies seeking to reveal bee diets by isotope analysis may gain better resolution in more isotopically distinct habitats, in conjunction with controlled feeding or isotope labeling experiments.