Large parks and city-wide tree cover boost butterfly diversity across 22 major U.S. cities
Ulrich, J.; Cheung, Y. Y. J.; Cosma, C. T.; Kharouba, H.; Guzman, L. M.
Show abstract
Accelerating global urbanization necessitates a better understanding of how to manage cities that promote biodiversity. However, we currently lack multi-year, multi-city studies, which limits a generalizable understanding of how both within and between city differences impact the spatial and temporal dynamics of urban biodiversity. Here, we tested hypotheses about the drivers of butterfly diversity within and across urban parks by applying Bayesian occupancy models to five years of iNaturalist community science data from 2,550 parks in 22 major U.S. cities. We found that cities with bigger parks supported more species per park, including more disturbance- and edge-avoidant species. This was driven by a positive effect of park size on butterfly species colonization rates. We also found that attributes of habitat quality (plant diversity within parks and tree cover surrounding parks) contributed to butterfly species occupancy. Park connectivity increased species persistence, but the overall effects on butterfly species occupancy varied across cities. Finally, we found that the total area of tree cover throughout a city, rather than the size or connectivity of individual parks, was the primary determinant of city-wide diversity: Increasing total tree canopy cover from below-average (~6%) to above-average (~22%) increased city-wide species richness by ~10%. These findings highlight the need for cities to maintain large parks while also increasing city-wide tree cover to support biodiversity across local to regional scales. By integrating high-resolution community science data across the continental U.S., this study provides mechanistic insight into how cross-scale processes shape urban biodiversity dynamics and identifies generalizable recommendations for improving urban conservation management.
Matching journals
The top 8 journals account for 50% of the predicted probability mass.