A conceptual overview of our methods.

We aggregated individual eBird checklists across the world (shown on the map), represented by the three colored insets which show the relationship between global range size (x-axis) and local abundance (y-axis) and the associated correlation value. We then aggregated these checklist level measures for 16,562,995 eBird checklists into the largest-ever meta-analysis to find the global level relationship between global range size and local abundance.

Funnel plots.

(A) the relationship between 16,562,995 effect sizes (Fisher’s; x-axis) and their precision (the square root of the inverse of the sampling variance; y-axis). (B) the relationship between 16,562,995 correlations based on 3,005,668,285 observations of 7,635 species (Pearson’s correlation coefficients; x-axis) and the number of species – 3, which is the inverse of the sampling variances for Zr (y-axis). Both plots consist of data points with the red dashed line indicating zero effect.

The relationship between species average (mean) density and species range size.

We calculated the mean density of a species in 5-degree grids where species occurred (y-axis) while the species range size (x-axis) was estimated by the sum of the percentage occurrence of the species multiplied by the grid size (km2) across all the 575 grids (7,464 species). The blue line indicates an average slope line from phylogenetic comparative models with 100 different posterior phylogenetic trees.

Results of the intercept meta-analytic model using the asreml function

Results of the meta-regression model with ‘checklist duration’ as a moderator, using the asreml function

Results of the meta-regression model with ‘sampling variance’ as a moderator, using the asreml function

Results of the meta-regression model with ‘checklist duration’ and ‘sampling variance’ as moderators, using the asreml function

Results of the phylogenetic regression using the phylolm and miInference functions