Abundance, diversity and body size: patterns from a grassland arthropod community

1. The empirical relationships among body size, species richness and number of individuals may give insight into the factors controlling species diversity and the relative abundances of species. To determine these relationships, we sampled the arthropods of grasslands and savannahs at Cedar Creek, MN using sweep nets (90 525 individuals of 1225 species) and pitfall traps (12 721 individuals of 92 species). Specimens were identified, enumerated and measured to determine body size.
2. Both overall and within abundant taxonomic orders, species richness and numbers of individuals peaked at body sizes intermediate for each group. Evolution could create unimodal diversity patterns by random diversification around an ancestral body size or from size-dependent fitness differences. Local processes such as competition or predation could also create unimodal diversity distributions.
3. The average body size of a species depended significantly on its taxonomic order, but on contemporary trophic role only within the context of taxonomic order.
4. Species richness ( S _i) within size classes was related to the number of individuals ( I _i) as S _i =  I i^0·5. This relationship held across a 100 000-fold range of body sizes. Within size classes, abundance distributions of size classes were all similar power functions. A general rule of resource division, together with similar minimum population sizes, is sufficient to generate the relationship between species richness and number of individuals.
5. Smaller bodied species had slightly shallower abundance distributions and may, in general, persist at lower densities than larger species.
6. Our results suggest there may be fewer undescribed small arthropod species than previously thought and that most undescribed species will be smaller than arthropods.     

The energetic equivalence rule rejected because of a potentially common sampling error: evidence from carabid beetles

Abundance data from pitfall traps are widely used to estimate the relationship between beetle body size and abundance. Such data probably overestimate densities of large bodied species and may overestimate slopes of size‐abundance relationships. Here, we test this idea by comparing size‐abundance patterns found using data from pitfall trapping with those found with data from a quantitative method of estimating abundance, quadrat sampling. We use data from a total of 47 communities. As expected, slopes of size‐abundance relationships are significantly more positive when estimated using data from pitfall traps compared to when using quadrat sampling data. This was seen when looking across different communities, within communities sampled by both methods and when focusing on the set of species found by both methods within a community. These results were also generally found regardless of method of analysis, which were done using regression with species values as independent data points and using the independent contrast method, and with slopes estimated using ordinary least square regression or the structural relation. Most important, slopes of size‐abundance relationships based on data from pitfall traps were on average significantly more positive than ?0.75 on log–log scales, and thus inconsistent with the energetic equivalence rule. Slopes based on quadrat sampling, on the other hand, were on average not significantly different from ?0.75. The rejection of the energetic equivalence rule based on data from pitfall traps here is therefore a sampling artefact. Similar problems may apply to abundance data from virtually all insect trapping methods, and should make us consider re‐examining many of the size‐abundance patterns documented so far. As a large proportion of all animal species are insects, and traps are widely used to estimate abundance, this is a potentially important problem for our general understanding of the relationship between species body size and abundance.     

Connecting species richness, abundance and body size in deep-sea gastropods

Aim This paper examines species richness, abundance, and body size in deep‐sea gastropods and how they vary over depth, which is a strong correlate of nutrient input. Previous studies have documented the empirical relationships among these properties in terrestrial and coastal ecosystems, but a full understanding of how these patterns arise has yet to be obtained. Examining the relationships among macroecological variables is a logical progression in deep‐sea ecology, where patterns of body size, diversity, and abundance have been quantified separately but not linked together. Location 196–5042 m depth in the western North Atlantic. Method Individuals analysed represent all Vetigastropoda and Caenogastropoda (Class Gastropoda) with intact shells, excluding Ptenoglossa, collected by the Woods Hole Benthic Sampling Program (3424 individuals representing 80 species). Biovolume was measured for every individual separately (i.e. allowing the same species to occupy multiple size classes) and divided into log₂ body size bins. Analyses were conducted for all gastropods together and separated into orders and depth regions (representing different nutrient inputs). A kernel smoothing technique, Kolmogorov‐Smirnov test of fit, and OLS and RMA were used to characterize the patterns. Results Overall, the relationship between the number of individuals and species is right skewed. There is also a positive linear relationship between the number of individuals and the number of species, which is independent of body size. Variation among these relationships is seen among the three depth regions. At depths inferred to correspond with intermediate nutrient input levels, species are accumulated faster given the number of individuals and shift from a right‐skewed to a log‐normal distribution. Conclusion A strong link between body size, abundance, and species richness appears to be ubiquitous over a variety of taxa and environments, including the deep sea. However, the nature of these relationships is affected by the productivity regime and scale at which they are examined.     

Richness, structure and functioning in metazoan parasite communities

Ecosystem functioning, characterized by components such as productivity and stability, has been extensively linked with diversity in recent years, mainly in plant ecology. The aim of our study was thus to quantify general relationships between diversity, community structure and ecosystem functions in metazoan parasite communities. We used data on parasite communities from 15 species of marine fish hosts from coastal Chile. The volumetric abundance (volume of all parasite species per individual host, in mm³) was used as a surrogate for productivity. Species diversity was measured using both species richness and evenness, while community structure was estimated using the co‐occurrence indices V‐ratio, C‐score and a new C‐score_s index standardized for the number of host replicates. After correcting for fish size, 47% of host species show no relationship, 13% show a hump shaped curve and 40% show positive monotonic relationships between productivity and parasite richness across all host individuals in a sample. We obtained a logarithmically decreasing relationship between evenness and productivity for all fish species, and propose a ‘dominance‐resistance’ hypothesis based on immunity to explain this pattern. The stability of the parasite community, measured as the coefficient of variation in productivity among individual hosts, was strongly and positively related to mean species richness across the 15 host species. The C‐score_s index, based on the number of checkerboard units in the host‐parasite presence/absence matrix, increases linearly with mean productivity across the 15 host species, suggesting that parasite communities tend to be more structured when they are more productive. This is the likely reason why linear relationships between richness and productivity were not observed consistently in all fish species. Parasite communities provide some clear patterns for the diversity–ecosystem functioning debate in ecology, although other factors, such as the history of community assembly, may also influence these patterns.     

Ecological scaling laws link individual body size variation to population abundance fluctuation

Scaling research has seen remarkable progress in the past several decades. Many scaling relationships were discovered within and across individual and population levels, such as species–abundance relationship, Taylor's law, and density mass allometry. However none of these established patterns incorporate individual variation in the formulation. Individual body size variation is a key evolutionary phenomenon and closely related to ecological diversity and species adaptation. Using a macroecological approach, I test 57 Long‐Term Ecological Research data sets and show that a power‐law and a generalized power‐law function describe well the mean‐variance scaling of individual body mass. This relationship connects Taylor's law and density mass allometry, and leads to a new scaling pattern between the individual body size variation and population abundance fluctuation, which is confirmed using freshwater fish and forest tree data. Underlying mechanisms and implications of the proposed scaling relationships are discussed. This synthesis shows that integration and extension of existing ecological laws can lead to the discovery of new scaling patterns and complete our understanding of the relation between individual trait and population abundance. Synthesis Scaling relationships are useful for community ecology as they reveal ubiquitous patterns across different levels of biological organizations. This work extends and integrates two existing scaling laws: Taylor's law and density‐mass allometry, and derives a new variance allometry between individual body mass and population abundance. The result shows that diverse individual body size is associated with stable population fluctuation, reflecting the effect of individual traits on population characteristics. Confirmed by several empirical data sets, these scaling relationships suggest new ways to study the underlying mechanisms of Taylor's law and have profound implications for fisheries and other applied sciences.     

Structure and resilience of a tidepool fish assemblage at Barbados

Synopsis Fish collections from 19 tidepools on a rock plateau at Martins Bay, on the east coast of Barbados, taken on three occasions (1981,1983 and 1987) contained 2078 individuals of 63 species. The number of species, individuals and total biomass increased with pool size. Partial residents, primarily juveniles of reef species, comprised 44% of species, 36% of numbers, and 26% of biomass. True and partial residents were of similar sizes. Most of the latter grow to larger sizes than those observed in the pools, indicating that the use of tidepools by fishes is size-dependent. Species richness, numbers of individuals and biomass in individual pools was positively associated with pool size. These relationships did not vary among sampling occasions. Species composition and relative abundance was also found to be similar among sampling occasions, leading to the conclusion that the tidepool assemblages are resilient and stable.     

The effect of fragment area on site‐level biodiversity

Habitat fragmentation accompanies habitat loss, and drives additional biodiversity change; but few global biodiversity models explicitly analyse the effects of both fragmentation and loss. Here we propose and test the hypothesis that, as fragment area increases, species density (the number of species in a standardised plot) will scale with an exponent given by the difference between the exponents of the species–area relationships for islands (z ~ 0.25) and in contiguous habitat (z ~ 0.15), and test whether scaling varies between land uses. We also investigate the scaling of overall abundance and rarefaction‐based richness, as some mechanisms make different predictions about how fragment area should affect them. The relevant data from the taxonomically and geographically broad PREDICTS database were used to model the three diversity measures, testing their scaling with fragment area and whether the scaling exponent varied among land uses (primary forest, secondary forest, plantation forest, cropland and pasture). In addition, the consistency of the response of species density to fragment area was tested across three well represented taxa (Magnoliopsida, Hymenoptera and ‘herptiles’). Species density and total abundance showed area‐scaling exponents of 0.07 and 0.16, respectively, and these exponents did not vary significantly among land uses; rarefaction‐based richness by contrast did not increase consistently with area. These results suggest that the area‐scaling of species density is driven by the area‐scaling of total abundance, with additive edge effects (species moving into the small fragments from the surroundings) opposing – but not fully overcoming – the effect of fragment area on overall density of individuals. The interaction between fragment area and higher taxon (plants, vertebrates and invertebrates), which remained in the rarefied richness model, indicates that mechanisms may vary among groups.     

Predicting abundance–body size relationships in functional and taxonomic subsets of food webs

Abundance–body size relationships are widely observed macroecological patterns in complete food webs and in taxonomically or functionally defined subsets of those webs. Observed abundance–body size relationships have frequently been compared with predictions based on the energetic equivalence hypothesis and, more recently, with predictions based on energy availability to different body size classes. Here, we consider the ways in which working with taxonomically or functionally defined subsets of food webs affected the relationship between the predicted and observed scaling of biomass and body mass in sediment dwelling benthic invertebrate communities at three sites in the North Sea. At each site, the energy available to body size classes in the “whole” community (community defined as all animals of 0.03125–32.0 g shell-free wet weight) and in three subsets was predicted from estimates of trophic level based on nitrogen stable isotope analysis. The observed and predicted scalings of biomass and body size were not significantly different for the whole community, and reflected an increase in energy availability with body size. However, the results for subsets showed that energy availability could increase or decrease with body size, and that individuals in the subsets were likely to be competing with individuals outside the subsets for energy. We conclude that the study of abundance–body mass relationships in functionally or taxonomically defined subsets of food webs is unlikely to provide an adequate test of the energetic equivalence hypothesis or other relationships between energy availability and scaling. To consistently and reliably interpret the results of these tests, it is necessary to know about energy availability as a function of body size both within and outside the subset considered.     

Size–abundance relationships in Florida ant communities reveal how ants break the energetic equivalence rule

1. Ants are among the most abundant terrestrial organisms, yet little is known of how ant communities divide resources because it is difficult to measure the number of individuals in colonies and the density of colonies. 2. The body size–abundance relationships of the ants of five upland ecosystems in Florida were examined. The study tested whether abundance, energy use, and total biomass were distributed among species and body sizes as predicted by Damuth's energetic equivalence rule. Estimates of average worker body size, colony size, colony mass, and field metabolic rates were used to examine the relationships among body sizes, energy use, and total biomass. 3. Analyses revealed significant variation in energy use and did not support the energetic equivalence hypothesis. Specifically, the energy use and total standing biomass of species with large workers and colonies was much greater than smaller species. 4. These results suggest that larger species with larger colonies account for a disproportionate fraction of the total abundance and biomass of ants. A general model of resource allocation in colonies provides a possible explanation for why ants do not conform to the predictions of the energetic equivalence rule and for why ants are so abundant.     

Plant species richness in continental southern Siberia: effects of pH and climate in the context of the species pool hypothesis

Aim Many high‐latitude floras contain more calcicole than calcifuge vascular plant species. The species pool hypothesis explains this pattern through an historical abundance of high‐pH soils in the Pleistocene and an associated opportunity for the evolutionary accumulation of calcicoles. To obtain insights into the history of calcicole/calcifuge patterns, we studied species richness–pH–climate relationships across a climatic gradient, which included cool and dry landscapes resembling the Pleistocene environments of northern Eurasia. Location Western Sayan Mountains, southern Siberia. Methods Vegetation and environmental variables were sampled at steppe, forest and tundra sites varying in climate and soil pH, which ranged from 3.7 to 8.6. Species richness was related to pH and other variables using linear models and regression trees. Results Species richness is higher in areas with warmer winters and at medium altitudes that are warmer than the mountains and wetter than the lowlands. In treeless vegetation, the species richness–pH relationship is unimodal. In tundra vegetation, which occurs on low‐pH soils, richness increases with pH, but it decreases in steppes, which have high‐pH soils. In forests, where soils are more acidic than in the open landscape, the species richness–pH relationship is monotonic positive. Most species occur on soils with a pH of 6–7. Main conclusions Soil pH in continental southern Siberia is strongly negatively correlated with precipitation, and species richness is determined by the opposite effects of these two variables. Species richness increases with pH until the soil is very dry. In dry soils, pH is high but species richness decreases due to drought stress. Thus, the species richness–pH relationship is unimodal in treeless vegetation. Trees do not grow on the driest soils, which results in a positive species richness–pH relationship in forests. If modern species richness resulted mainly from the species pool effects, it would suggest that historically common habitats had moderate precipitation and slightly acidic to neutral soils.     

How dispersal limitation shapes species-body size distributions in local communities

A critical but poorly understood pattern in macroecology is the often unimodal species-body size distribution (also known as body size-diversity relationship) in a local community (embedded in a much larger regional species pool). Purely neutral community models that assume functional equivalence among species are incapable of explaining this pattern because body size is the key determinant of functional differences between species. Several niche-based explanations have been offered, but none of them is completely satisfactory. Here we develop a simple model that unites a neutral community model with niche-based theory to explain the relationship. In the model, species of similar size are assumed to belong to the same size guild. Within a size guild, all individuals are equivalent in their competition for resources, sensu Hubbell's neutral community model; they have the same speciation rate and dispersal capacities. Between size guilds, however, the total number of individuals, the speciation rate, and the dispersal capacities differ, but using known allometric scaling laws for these properties, we can describe the differences between size guilds. Our model predicts that species richness reaches an optimum at an intermediate body size, in agreement with observations. The optimum at intermediate body size is basically the result of a trade-off between, on the one hand, allometric scaling laws for the number of individuals and the speciation rate that decrease with body size and, on the other hand, the scaling law for active dispersal that increases with body size.     

The effects of fertilization on the trait-abundance relationships in a Tibetan alpine meadow community

Aims Comparisons of the trait–abundance relationships from various habitat types are critical for community ecology, which can offer us insights about the mechanisms underlying the local community assembly, such as the relative role of neutral vs. niche processes in shaping community structure. Here, we explored the responses of trait–abundance relationships to nitrogen (N), phosphorus (P) and potassium (K) fertilization in an alpine meadow.Methods Five fertilization treatments (an unfertilized control and additions of N, P, K and NPK respectively) were implemented using randomized block design in an alpine Tibetan meadow. Species relative abundance (SRA), plant above-ground biomass and species richness were measured in each plot. For 24 common species, we measured species functional traits: saturated height, specific leaf area (SLA) and leaf dry matter content (LDMC) in each treatment but seed size only in the unfertilized control. Standard major axis (SMA) regression and phylogenetically independent contrasts (PICs) analysis were used to analyse species trait–abundance relationships in response to different fertilization treatments.Important findings Positive correlations between SRA and saturated height were raised following N, P and NPK fertilizations, which indicated an increase in light competition in these plots. In P fertilized plots, SRA was also positively correlated with LDMC because tall grasses with a nutrients conservation strategy often have a relative competitive advantage in capturing limited light and soil nutrients. In K fertilized plots, neither the trait–abundance relationships nor above-ground biomass or species richness significantly differed from that in the control, which suggests that K was not a limiting resource in our study site. These significant correlations between species traits and relative abundance in fertilized treatment suggest that trait-based selection plays an important role in determining species abundance within local communities in alpine meadows.     

Body size distributions in North American freshwater fish: large‐scale factors

Aim To document continental‐ and regional‐scale variation in the size distributions of freshwater fish and examine some energetic, evolutionary and biogeographic explanations for these patterns. Location North America. Methods Regional species lists, coupled with habitat and body size information, were used to document the spatial patterns. Results At the continental scale, riverine specialist fishes show a unimodal, right‐skewed, body size distribution whereas habitat generalist and lacustrine specialist species exhibit bimodal size distributions, with only a slight preponderance of small‐mode species. Most large‐mode species are migratory. Resident species, unlike migratory ones, show a latitudinal increase in mean size, but the size increase across all species is steeper because the importance of large migratory species increases with latitude. Size distributions change from right‐ to left‐skewed with increasing latitude. Maximum body size does not change with increasing family richness but minimum size declines and skewness increases, consistent with diversification of small species. Skewness does not vary with mean family body size. Main conclusions Post‐glacial recolonization by large, habitat generalist, migratory species is the main determinant of latitudinal size distribution trends. There is little support for the energetic hypothesis, but the data are consistent with a negative Cope's rule.     

Energetic inequivalence in eusocial insect colonies

The energetic equivalence rule states that population-level metabolic rate is independent of average body size. This rule has been both supported and refuted by allometric studies of abundance and individual metabolic rate, but no study, to my knowledge, has tested the rule with direct measurements of whole-population metabolic rate. Here, I find a positive scaling of whole-colony metabolic rate with body size for eusocial insects. Individual metabolic rates in these colonies scaled with body size more steeply than expected from laboratory studies on insects, while population size was independent of body size. Using consumer-resource models, I suggest that the colony-level metabolic rate scaling observed here may arise from a change in the scaling of individual metabolic rate resulting from a change in the body size dependence of mortality rates.     

Environmental correlates of body size distributions of European springtails (Hexapoda: Collembola)

Aim Species–body size distributions (SBDs) are plots of species richness across body size classes. They have been linked to energetic constraints, speciation–extinction dynamics and to evolutionary trends. However, little is known about the spatial variation of size distributions. Here we study SBDs of European springtails (Collembola) at a continental scale and test whether minimum, average and maximum body size and the shapes of size distributions change across latitudinal and longitudinal gradients and whether SBDs of islands and mainlands differ. We also test whether the island rule and the positive body size–range size relationship of vertebrates also holds for Collembola. Location Europe. Methods We use a unique data set on the spatial distributions of 2102 species of European springtails across 52 countries and larger islands together with associated data on body size, area, climate variables, longitude and latitude. Differences in the central moments of SBDs are inferred from simultaneous spatial autoregression models. Results The SBD of the European Collembola and its largest suborder Entomobryomorpha is unimodal and symmetrical. Average, minimum and maximum body weight and the skewness of the mainland/island SBDs peaked at intermediate latitudes. We could not find simple latitudinal gradients in minimum and maximum body weight. Average and maximum body size increased with country/island area in accordance with the island rule in vertebrates, while minimum body size did not significantly differ between islands and mainlands. Finally, we found a weak but statistically significant positive correlation of range size and body size. Main conclusions We provide evidence for differences in body size distributions between islands and mainlands that are in part in line with the island rule in invertebrates. We also find evidence for an interspecific body size–range size relationship similar to that of vertebrates although the vertebrate pattern is much stronger than the springtail pattern. Our results on latitudinal gradients of maximum and average body size imply the need to account for species richness and area effects in the study of latitudinal gradients in body size. We recommend implementing sample size and area effects in the study of body size distributions on islands and mainlands.     

Individual size variation reduces spatial variation in abundance of tree community assemblage,not of tree populations

Research on individual trait variation has gained much attention because of its implication for ecosystem functions and community ecology. The effect of individual variation on population and community abundance (number of individuals) variation remains scarcely tested. Using two established ecological scaling laws (Taylor's law and abundance–size relationship), we derived a new scaling relationship between the individual size variation and spatial variation of abundance. Tested against multi‐plot tree data from Diaoluo Mountain tropical forest in Hainan, China, the new scaling relationship showed that individual size variation reduced the spatial variation of community assemblage abundance, but not of taxon‐specific population abundance. The different responses of community and population to individual variation were reflected by the validity of the abundance–size relationship. We tested and confirmed this scaling framework using two measures of individual tree size: aboveground biomass and diameter at breast height. Using delta method and height‐diameter allometry, we derived the analytic relation of scaling exponents estimated under different individual size measures. In addition, we used multiple regression models to analyze the effect of taxon richness on the relationship between individual size variation and spatial variation of population or community abundance, for taxon‐specific and taxon‐mixed data, respectively. This work offers empirical evidence and a scaling framework for the negative effect of individual trait variation on spatial variation of plant community. It has implications for forest ecosystem and management where the role of individual variation in regulating population or community spatial variation is important but understudied.     

Size-Energy Relationships in Ecological Communities

Hypotheses that relate body size to energy use are of particular interest in community ecology and macroecology because of their potential to facilitate quantitative predictions about species interactions and to clarify complex ecological patterns. One prominent size-energy hypothesis, the energetic equivalence hypothesis, proposes that energy use from shared, limiting resources by populations or size classes of foragers will be independent of body size. Alternative hypotheses propose that energy use will increase with body size, decrease with body size, or peak at an intermediate body size. Despite extensive study, however, size-energy hypotheses remain controversial, due to a lack of directly-measured data on energy use, a tendency to confound distinct scaling relationships, and insufficient attention to the ecological contexts in which predicted relationships are likely to occur. Our goal, therefore, was to directly evaluate size-energy hypotheses while clarifying how results would differ with alternate methods and assumptions. We comprehensively tested size-energy hypotheses in a vertebrate frugivore guild in a tropical forest in Madagascar. Our test of size-energy hypotheses, which is the first to examine energy intake directly, was consistent with the energetic equivalence hypothesis. This finding corresponds with predictions of metabolic theory and models of energy distribution in ecological communities, which imply that body size does not confer an advantage in competition for energy among populations or size classes of foragers. This result was robust to different assumptions about energy regulation. Our results from direct energy measurement, however, contrasted with those obtained with conventional methods of indirect inference from size-density relationships, suggesting that size-density relationships do not provide an appropriate proxy for size-energy relationships as has commonly been assumed. Our research also provides insights into mechanisms underlying local size-energy relationships and has important implications for predicting species interactions and for understanding the structure and dynamics of ecological communities.     

Body size and measurement of species diversity in large grazing mammals

Species are by definition different from each other. This fact favours ranking rather than additive indices. However, ecologists have measured species diversity in terms of species richness, or by combining species richness with the relative abundance of species within an area. Both methods raise problems: species richness treats all species equally, while relative abundance is not a fixed property of species but varies widely temporally and spatially, and requires a massive sampling effort. The functional aspect of species diversity measurement may be strengthened by incorporating differences between species such as body size as a component of diversity. An index of diversity derived from a measure of variation in body size among species is proposed for large grazing mammals. The proposed diversity index related positively to species abundance, indicating that the use of body size as a surrogate for diversity is adequate. Because the proposed index is based on presence or absence data, the expensive and time consuming counting of individuals per species in each sampling unit is not necessary.     

Spatial variation of species diversity across scales in an old-growth temperate forest of China

Species richness and abundance are the two most important diversity variables. Species abundance is additive when aggregated across spatial scale, whereas species richness is non-additive. This study analyzes the effect of spatial scale and site on species abundance and richness in a 25-ha temperate forest plot in the Changbai Mountains, northeastern China. The result shows that species abundance and richness are not only dependent on spatial scales, but also dependent on site. Species abundance responds linearly to changes of spatial scale with no intersection in different sites of the study area. However, although species richness also increases with the increase of spatial scale, there are some intersections for the different sites, suggesting that a species-rich site does not always have a high value if the spatial scale is changed. In all, with respect to additive variables, it is relatively easy to extrapolate them from one spatial scale to another spatial scale, as they and the spatial scale usually form a linear relationship. In contrast, non-additive variables are difficult to extrapolate across spatial scales, because they often respond nonlinearly to spatial scale changes. In order to extrapolate these non-additive variables across spatial scales, it is necessary to estimate the relationships between them and spatial scales. As a result, extrapolation of information among spatial scales may be possible, but very difficult, especially for non-additive variables. Because the 25-ha Changbai plot is very small compared to the extent of the world temperate forests, and the vegetation is a relatively uniform type, more such studies in other ecosystems are needed before theories and generalization about scaling effects can be formulated.     

Inter‐assemblage facilitation: the functional diversity of cavity‐producing beetles drives the size diversity of cavity‐nesting bees

Inter‐specific interactions are important drivers and maintainers of biodiversity. Compared to trophic and competitive interactions, the role of non‐trophic facilitation among species has received less attention. Cavity‐nesting bees nest in old beetle borings in dead wood, with restricted diameters corresponding to the body size of the bee species. The aim of this study was to test the hypothesis that the functional diversity of cavity‐producing wood boring beetles ‐ in terms of cavity diameters ‐ drives the size diversity of cavity‐nesting bees. The invertebrate communities were sampled in 30 sites, located in forested landscapes along an elevational gradient. We regressed the species richness and abundance of cavity nesting bees against the species richness and abundance of wood boring beetles, non‐wood boring beetles and elevation. The proportion of cavity nesting bees in bee species assemblage was regressed against the species richness and abundance of wood boring beetles. We also tested the relationships between the size diversity of cavity nesting bees and wood boring beetles. The species richness and abundance of cavity nesting bees increased with the species richness and abundance of wood boring beetles. No such relationship was found for non‐wood boring beetles. The abundance of wood boring beetles was also related to an increased proportion of cavity nesting bee individuals. Moreover, the size diversity of cavity‐nesting bees increased with the functional diversity of wood boring beetles. Specifically, the mean and dispersion of bee body sizes increased with the functional dispersion of large wood boring beetles. The positive relationships between cavity producing bees and cavity nesting bees suggest that non‐trophic facilitative interactions between species assemblages play important roles in organizing bee species assemblages. Considering a community‐wide approach may therefore be required if we are to successfully understand and conserve wild bee species assemblages in forested landscapes.