Experimental dispersal reveals characteristic scales of biodiversity in a natural landscape
See allHide authors and affiliations
Edited by David Tilman, University of Minnesota, St. Paul, MN, and approved March 28, 2017 (received for review September 14, 2016)

Significance
Biological communities differ dramatically in numbers and identities of species, a pattern that could be explained by many mechanisms that each vary with spatial scale. Testing how ecological mechanisms transition among scales is key to understanding the maintenance of diversity but is infeasible in most systems. In a natural plant community, we experimentally enhanced seed dispersal over scales ranging from 1 m to 10 km using a seed vacuum. Our results indicate that pervasive dispersal limitation constrains local communities at and above 100 m, causing communities to contain half as many species as the environment is capable of supporting. Our results suggest that many species at risk for extirpation in this global biodiversity hotspot could be restored through managed dispersal.
Abstract
Ecological theory posits that dispersal among habitat patches links local communities and is a key “regional” process that maintains biological diversity. However, manipulations required to experimentally test regional processes are infeasible for most systems, and thus more work is needed to detect the scales at which regional processes manifest and their overall effect on diversity. In a Californian grassland, a hotspot for global biodiversity, we used a seed vacuum to increase dispersal at spatial scales varying from 1 m to 10 km while maintaining a realistic spatial structure of species pools and environmental conditions. We found that dispersal limitation has a profound influence on diversity; species richness increased with the spatial scale of seed mixing, doubling in plots that received seed from large (≥5 km) compared with small (≤5 m) scales. This increase in diversity corresponded to an increase in how well species distributions were explained by environmental conditions, from modest at small scales (R2 = 0.34) to strong at large scales (R2 = 0.52). Responses to the spatial scale of seed mixing were nonlinear, with no differences below 5 m or above 5 km. Nonlinearities were explained by homogeneity of environmental conditions below 5 m and by a lack of additional variation in the species pool above 5 km. Our approach of manipulating natural communities at different spatial scales reveals (i) nonlinear transitions in the importance of environmental sorting and dispersal, and (ii) the negative effects of dispersal limitation on local diversity, consistent with previous research suggesting that large numbers of species are headed toward regional extinction.
The problem of pattern and scale is the central problem in ecology, unifying population biology and ecosystem science, and marrying basic and applied ecology.
S. A. Levin (1992)
The processes that structure ecological populations, biodiversity, and ecosystem properties transition in importance across spatial scales (1, 2). As the spatial scale of observation increases, the range of environments sampled (1, 3) and the geographic distance separating localities (4) become increasingly important in shaping species distributions. As a consequence, the relative spatial scaling of different ecological processes is thought to underlie some of the most important patterns in ecology, such as species–area (5) and biodiversity–ecosystem function relationships (6). Because identifying the important scales is challenging, ecologists often compare ecological patterns among local and regional scales to simplify theoretical (7⇓–9) and empirical research (10, 11). However, how closely local and regional delineations match up with the actual scaling of ecological processes is rarely known (2, 12). Quantifying the spatial scaling of these processes promises to enrich our understanding of the mechanisms that maintain diversity, yet remains elusive even in biodiversity hotspots that require this information for conservation decisions (13).
A major challenge to testing how ecological processes transition among spatial scales in natural communities is that regional processes, unlike local processes, are not often amenable to experimental manipulation in the field due to the inability to move most communities of organisms. Manipulative tests of local and regional diversity are typically performed in mesocosm experiments on communities constructed using simplified environments, species pools, or dispersal patterns (14). Because it is often unclear how such simplifications affect experimental outcomes (15), mesocosm experiments allow essential tests of the range of potential outcomes under different sets of experimental conditions, but cannot capture the importance of processes that occur in nature. Finding new ways to combine the power of manipulative field experiments with the biological realism of natural landscape structure is necessary to address fundamental questions, such as: how widespread is dispersal limitation in a community and at what spatial scales does it manifest? How strongly do species sort along environmental gradients in the absence of dispersal limitation? And, even more basically, at what spatial scales are local and regional communities most appropriately defined?
The above questions can be tested by experimentally removing dispersal limitation through the homogenization of species pools, that is, by redistributing species equally across habitat patches while maintaining their abundance distributions. Theory makes distinct predictions for how such homogenization would affect species diversity, and whether it would increase or decrease the probability that species are found in environmentally suitable habitat patches (9, 16, 17). This latter concept, referred to here as “species–environment associations,” is commonly measured by the amount of variation in species distributions that is explained by environmental covariates (18). If dispersal chronically limits species movement and the landscape is patchy and heterogeneous, then homogenization of species pools would increase both species richness (19, 20) and the strength of species–environment associations (16) (SI Appendix, Fig. S1A). These changes occur as species access suitable localities that were previously inaccessible (9), causing the average number of species in a site to increase and the variation in species distributions explained by the environment to similarly increase.
Alternatively, if species are not dispersal-limited and already well-matched to environmental conditions (SI Appendix, Fig. S1B), the strength of species–environment associations should decrease even as local richness increases due to mass effects (16, 17), meaning that species are increasingly found in environments to which they are not suited. When species–environment associations are neutral with respect to species identity, species–environment associations should not change whether species have limited dispersal (SI Appendix, Fig. S1C; increased species richness) or are already dispersing evenly among communities (SI Appendix, Fig. S1D; no change in species richness) (8, 16). In the metacommunity ecology literature, these four potential outcomes of increasing dispersal would correspond to species sorting, mass effects, and neutral dynamics (with or without dispersal limitation), respectively (9, 16). The approach of homogenizing species pools and tracking diversity responses to local conditions is similar to that taken by many experiments that use mesocosm (e.g., ref. 11) or field (e.g., refs. 19⇓–21) manipulations with artificial species pools (14, 15, 22), but has yet to be implemented in the field with natural spatial structure in species pools, abundance distributions, and environments.
We tested the spatial scales at which dispersal and environmental heterogeneity impose constraints on species distributions and diversity using a “hay transfer” technique (23) in serpentine annual plant communities in California. This technique involved vacuuming the seed bank and other loose material from field plots, pooling the material among plots, and then redistributing it. Our “species-pooling” treatments captured five spatial scales (∼1 m, 5 m, 100 m, 5 km, and 10 km; SI Appendix, Fig. S2), with the 100-m treatment receiving a mix of seeds from all plots within 100 m, for example. In the annual plant communities studied, pooling the seed bank redistributes all plants during their dormant stage. In the following growing season, we collected data on species occupancy and environmental conditions to answer four questions: (i) How common is dispersal limitation in a metacommunity? (ii) How strongly do species sort along environmental gradients once potential dispersal limitation is experimentally removed? (iii) At what spatial scales do signals of dispersal limitation and environmental sorting manifest? (iv) Do these signals correspond to the natural structure of species pools and environmental conditions? Our surveys were conducted at peak flowering, after the filtering effects of competition and environment, which can be strongest at the seedling stage (24, 25), had time to take effect.
Our serpentine study system is ideal for testing questions of the spatial scaling of diversity for two reasons. First, the high occurrence of annual species and patchy distribution of serpentine habitat among a nonserpentine matrix is well-suited to testing spatial questions (26) with the hay transfer method (23). Second, our study took place in the California Floristic Province, a global biodiversity hotspot that is threatened by the pervasive invasion of European grasses (26, 27). Serpentine soils are hypothesized to act as spatial refugia for native species to escape the direct competitive effects of invasive grasses (28), but less often considered are the indirect negative effects of invaders on diversity through the isolation of habitat patches (13, 29).
Given that invasive grasses and other human impacts [road building, fire suppression (30)] have increased fragmentation of native plant communities in the past 200 y, we hypothesized that diversity is constrained by limited dispersal and that, if this dispersal constraint were removed, species would sort deterministically according to their environmental niche requirements (i.e., increased species–environment associations; SI Appendix, Fig. S1A). Environmental conditions were described by two composite variables—principal components analysis (PCA) axes—that summarized soil chemistry, fertility, moisture, and site topography. The strength of species–environment associations was the pseudo-R2 of logistic regressions that predicted the occurrences of all species with environmental conditions (PCA axes) used as predictors (see Data analysis in Materials and Methods and SI Appendix, Fig. S4 A–C for an example of data that underlie pseudo-R2 results). The degree of dispersal limitation, and the scales at which its effects manifest, are important to understanding biodiversity in this global hotspot and to projecting community stability over longer timescales (13).
Results
Increasing the spatial scale of species pooling caused a net increase of 12 species in plots (F4,112 = 69.4, P < 0.001; Fig. 1A). This increase was sigmoidal rather than linear, with no significant differences in species richness between the two smallest (1 m and 5 m; P = 0.72) or the two largest (5 km and 10 km; P = 0.34) spatial scales, but highly significant (all P < 0.001) differences among small (≤5 m), intermediate (100 m), and large (≥5 km) spatial scales. As a result, species richness in our plots doubled (F4,112 = 69.9, P < 0.001) from 10 species in the 1-m treatment to 21 species in the 10-km treatment (SI Appendix, Fig. S3A; all measures are for 0.75- × 0.75-m plots). Species richness postmanipulation was composed of species initially present at sites (not dispersal limited) plus those gained by species pooling (dispersal limited; SI Appendix, Fig. S3B).
Effect of spatial scale of species pooling on (A) net change in species richness (mean ± SE), relative to an unmanipulated control plot (C3 in SI Appendix, Fig. S2B) and (B) the strength of species–environment associations across a community of species (pseudo-R2 values); observed patterns correspond to prediction in SI Appendix, Fig. S1A. Points with the same letter were not significantly different in a multiple comparisons test. Net changes in species richness are driven by new species gained, which translates into a doubling in species richness when species are pooled at the 1-m versus at the 10-km scale (SI Appendix, Fig. S3).
To explore whether this increase in species richness corresponded with an increase or decrease in species–environment associations, we tested whether environmental conditions predicted the occurrences of all 73 annual species observed in our study. Environmental conditions were highly correlated and thus were summarized by the first and second axis scores of a principal components analysis (Materials and Methods). Axis 1 primarily summarized soil chemistry (i.e., Ca/Mg, Olsen-P, X-K, NO3-N, and organic matter), whereas axis 2 summarized elevation and soil moisture (SI Appendix, Fig. S5B). We found the proportion of variance explained (pseudo-R2) by our models increased from 0.34 to 0.52 with increasing spatial scale of species pooling (Fig. 1B); the sigmoidal response of species richness to the spatial scale of species pooling was closely mirrored by the sigmoidal increase in variation explained by the environment (Fig. 1). Supplementary analyses of species–environment associations with percentage of cover data, rather than presence/absence data, produced broadly similar increases with increasing spatial scale of species pooling (SI Appendix, Fig. S9 and discussion in SI Appendix, Supplementary Methods, Results, and Discussion).
In addition to plot-level species richness (alpha diversity), we tested the effects of our species-pooling treatments on other components of diversity. As with species richness, the composition of species in plots responded sigmoidally to the spatial scale of species pooling (Fig. 2A), as estimated using the axis 1 and 2 scores of a principal coordinates analysis (PCoA) (Materials and Methods). These axes were significantly associated with 42 plant species (P < 0.05), including many endemic serpentine-associated species (e.g., Clarkia gracilis, Collinsia sparsiflora, and Navarretia jepsonii). Beta diversity (i.e., among-plot turnover in species composition, measured with Jaccard’s dissimilarity index) decreased with the scale of species pooling from 0.81 to 0.56 between the 1-m to 10-km spatial scaling treatments (Fig. 2B). The size of the regional species pool (gamma diversity), estimated as the total number of unique species that emerged across plots, was not affected by the species-pooling treatments, with a mean pool size of 58.0 ± 0.5 SE across treatments (x-intercept; SI Appendix, Fig. S6). At a common gamma diversity, alpha and beta diversity are inversely related, as compositional dissimilarity among sites decreases as a consequence of species being found in more sites (31). Site occupancy (the number of sites occupied by each species) also increased with the spatial scale of species pooling (SI Appendix, Fig. S6), as did community evenness due to the increasing occurrences of regionally rare species (shallower slope in SI Appendix, Fig. S6). Not all 73 species were observed in all treatments; on average, each treatment had 1.2 species (±0.58 SE) that were not observed in any other treatment.
Effect of spatial scale of species pooling on (A) species composition and (B) the compositional dissimilarity of plots (beta diversity). Species composition shows the mean and SE of the first axis scores from a PCoA with Jaccard’s dissimilarity index (see biplot in SI Appendix, Fig. S5A); results with second-axis scores are qualitatively similar. Points with the same letter were not significantly different in a multiple comparisons test. Beta diversity patterns mirror changes in species richness, as increasing local (alpha) diversity with no change in regional (gamma) diversity reduces among-plot compositional dissimilarity.
Differences in the naturally occurring spatial structure of environmental conditions and species pools provide further evidence of how regional processes structure local patterns of diversity. We first tested for differences in environmental conditions among sites at spatial scales corresponding to our experimental setup (SI Appendix, Fig. S2) to characterize the spatial scaling of the environment. Environmental conditions differed significantly among groups of sites within halves (100-km scale; F = 14.05, P = 0.001), but differences between halves of the reserve (5-km scale; F = 247.68, P = 0.001) were larger, with limited overlap in environmental conditions (Fig. 3A). Differences among reserve halves were driven by axis 2 (% soil moisture, elevation) but not axis 1 (soil chemistry, fertility). For those environmental variables that were quantified with multiple within-site measurements (5-m scale), there were no statistically significant differences among plots that occurred at the same site (all P > 0.3; SI Appendix, Table S2). Despite these average changes in environmental conditions with distance, even distant sites (≥10 km apart) frequently shared similar environmental conditions (Fig. 3A and SI Appendix, Figs. S7 and S8), supporting the finding that dispersal limitation at these scales can have important consequences.
Natural differences among halves of the reserve in (A) environmental conditions (from a PCA) and (B) species composition (from a PCoA). The contours outline the bivariate normal distributions of each half of the reserve; dashed contours/upward triangles are sites that occur in the northern half, and solid contours/downward triangles are sites that occur in the southern half. Environmental conditions (P = 0.001) but not species (P = 0.175) significantly differed among halves (Results).
In contrast to environmental conditions, the composition of species significantly differed at the 100-m scale (F = 1.96, P = 0.001) but not between halves of the reserve (5-km scale; F = 1.26, P = 0.168). Despite the lack of a significant difference in species composition, the two halves of the reserve had only 62% (37/60) of species in common, with 20% (12/60) of species unique to the northern half and 18% (11/60) unique to the southern half; these numbers were calculated from the 1-m treatment and sum to less than 73 because not all species were observed in all treatments. The species unique to one-half of the reserve were significantly less associated with the axis 1 and 2 PCoA scores than species that were common to both halves (P = 0.023), which potentially reconciles the high compositional similarity among halves of the reserve despite the presence of unique species.
Discussion
Understanding the scale dependency of ecological processes is central to protecting biodiversity and a core interest of ecology research (1, 2, 32-34), yet experimental demonstrations have remained elusive for natural communities (2). By experimentally manipulating the spatial scale of seed dispersal in a biodiversity hotspot, our study resulted in two key findings: (i) compared with natural patterns, species richness and the strength of environmental matching doubled when seeds were distributed at large spatial scales, and (ii) nonlinear diversity responses to the spatial scale of species pooling suggests that there are distinct spatial transitions in the forces that structure diversity. These findings provide insight into how metacommunities and the distributions of their constituent species are structured, with implications for conservation.
In the absence of species pooling, species in our study showed modest associations with environmental conditions (blue point in Fig. 1B), consistent with the findings of several observational studies (e.g., refs. 34, 35). However, this type of observational evidence alone cannot identify the mechanisms underlying the pervasive lack of environmental associations, such as dispersal limitation, neutrality, and mass effects (SI Appendix, Fig. S1). The experimental enhancement of dispersal via seed pool manipulations allowed us to discriminate among potential mechanisms and supports a model in which dispersal limitation precludes the colonization of otherwise suitable sites (SI Appendix, Fig. S1A). Consistent with these results, a subsequent analysis showed that species with more effective dispersal modes (wind and vertebrate dispersed) tended to occupy more sites in the absence of species pooling (SI Appendix, Fig. S10 and Supplementary Methods, Results, and Discussion).
Similar evidence of the role of dispersal limitation has been obtained for individual species, such as the canary violet (4); we extend this perspective beyond individual species to an entire plant community and show a near doubling in the strength of species–environment associations that explains over half of the variation in how communities of species are distributed (Fig. 1B). An additional feature of our community-level seed additions is that they allowed natural spatial structure in both species composition and relative abundances, rather than fixed amounts of seed added for all sites and all species (19). Doing so most closely reflects natural variation in propagule pressure, an important determinant of dispersal limitation and community structure (20, 36). Although we cannot refute the possibility that the observed responses are transient, evidence from related seed addition research has shown that initial diversity effects are persistent, even in competitive environments (37), and increases in species–environment associations would not be expected in a transient scenario. Additionally, our results were insensitive to the exclusion of species with small populations, which are those most likely to be “sink” populations (SI Appendix, Table S3); we discuss these points at length in SI Appendix, Supplementary Methods, Results, and Discussion.
To what extent differences in species composition among local communities can be explained by distance or environment has been a major area of ecological research (38, 39), stimulated by the development of statistical methods to partition their relative influences from observational data (18, 40). However, variance-partitioning methods can produce biased estimates of spatial and environmental components of variation, as well as inflated model fits (41), and should be interpreted as rough estimates. Moreover, partitioning approaches cannot isolate the interactive effects of dispersal limitation and environmental components (18), yet it is these interactive effects that determine the persistence of species within landscapes (13, 29). In our study, we find that environmental conditions explain 32% of variation in species occupancy patterns in the absence of species pooling and an additional 20% when distance effects are experimentally removed (Fig. 1B). This increase when distance effects are removed quantifies the interactive effect of distance and environmental conditions, or the degree to which dispersal limitation prevents species from accessing suitable habitats.
Nearly all components of diversity that were examined showed consistent, nonlinear responses to the spatial scale of species pooling, including species richness (alpha diversity), species–environment associations, and compositional turnover (beta diversity). At small spatial scales (<5 m), conditions are homogenous, and most species are able to disperse these distances in one to few generations (42); our species-pooling manipulations had negligible effects at this scale (blue and green points, Figs. 1 and 2 and SI Appendix, Fig. S3). At intermediate (100 m) and larger (≥5 km) scales, our manipulations allowed seed to cross the nonserpentine grassland matrix to reach distant serpentine sites, resulting in increased species richness and species–environment associations consistent with dispersal limitation. Although environmental conditions also tended to be more different as distance increased, many sites were similar even when separated by ∼10 km (Fig. 3 and SI Appendix, Figs. S7 and S8), making dispersal across these distances important for accessing distant sites with similar environments.
Together, our evidence reveals three “characteristic scales” [sensu Levin (2)] or “domains of scale” [sensu Wiens (1)] among which strengths of ecological processes appear to differ. When seed pooling occurs at local scales ≤5 m, communities experience relatively homogenous environmental conditions, but the make-up of communities is severely limited by dispersal, such that species–environment associations frequently fail to emerge. As the scale of species pooling increases to 100 m and beyond, the increased environmental difference in new sites appears to be outweighed by the large numbers of similar sites (SI Appendix, Figs. S7 and S8) and the release from dispersal limitation. The net result is that species are better able to exploit the increased range of environmental heterogeneity, increasing species–environment associations. However, as the scale of species pooling increases beyond 5 km, dispersal limitation no longer alters community responses to continually increasing environmental heterogeneity.
The lack of differences among 5-km and 10-km scales was surprising given their differences in spatial extent and environmental conditions, particularly soil moisture and elevation (as summarized by PCA axis 2, Fig. 3A); there are two potential, nonmutually exclusive explanations for this trend. First, the high similarity in the composition of species pool among halves of the reserve (the 5-km scale, Fig. 3B) indicates that mixing seed among halves (10-km scale) would introduce little variation to the species pool, consistent with the lack of responses that we observed. Although we did observe a number of species that were unique to each half, these species were so rare that they did not contribute strongly to compositional variation among sites or the strengthening of species–environment associations in response to seed mixing. Second, it is possible that most species in the reserve are more specialized to soil fertility and biochemistry than to differences in elevation and soil moisture (Fig. 3, axis 1 and 2, respectively). Because the reserve halves did not significantly differ in axis 1 scores, specialization along this axis did not generate among-half differences in species composition.
Our finding that local communities are half as speciose as the environment is capable of supporting sheds light on the past and potential future of this global biodiversity hotspot. Californian landscapes no longer resemble those of 200 y ago; European invaders now form a matrix of unsuitable habitat (28) and have restricted native species to harsh refuge environments, such as serpentine (13, 26). Our findings of low species richness and high regional rarity (SI Appendix, Fig. S6) are consistent with those predicted for areas experiencing regional “extinction debts”—the delayed extinction of species from a region through chronic reductions in colonization rates (13). The prognosis for native diversity under this scenario is bleak because old invaders are adapting to tolerate harsh serpentine conditions (43) and new invaders that thrive in these harsh environments are spreading rapidly [e.g., barbed goatgrass (Aegilops triuncialis)]. The net effect is a shrinking and loss of native plant refuge patches that threaten the regional stability of native plant diversity.
In this scenario and others like it, a more regionally focused approach to conservation is needed to preserve many types of local communities and their constituent species, understanding that local persistence relies on colonization from adjacent patches (44). Management at the wrong scale may miss important covariation in species pools and environments, generating mismatches between species and local conditions (1, 4). In our system, management at a spatial extent of 5 km is likely to be most effective, particularly through assisted dispersal or creation of habitat patches that act as stepping stones between otherwise isolated patches. Although the exact scale of spatial delineations will differ among study systems depending on landscape structure, habitat matrix permeability, and dispersal distances of focal organisms, we have offered a clear example of how these critical scales can be identified. In other systems, approaches such as genetic analyses may be required to better understand characteristic scales of dispersal before assessing scales of environmental turnover.
The forces that dictate how species are distributed across landscapes have long fascinated ecologists. In the serpentine system, we have identified scale-specific processes that structure plant communities and how this spatial scaling is explained by the accumulation of species and environments across space. These results provide insights that are specific to an area that is considered a model system for biodiversity research (45) and more generally provide an experimental test of predictions from a well-developed body of theory (e.g., ref. 9). Broadly applying such approaches can better match conservation actions with ecological processes and promises to advance our understanding of one of the longest-standing challenges in ecology.
Materials and Methods
Study System.
The field experiment was conducted at the 2,800-ha McLaughlin Natural Reserve (nrs.ucdavis.edu/McL/) in northern California (38.8739° N, 122.4317° W). The region has a Mediterranean climate, featuring cool, wet winters (November–March) and hot, dry summers (April–October) with ∼750 mm of annual rainfall. The landscape is largely composed of chaparral, oak woodland, and grassland meadow habitat; common herbivores include mule deer, jack rabbits, and pocket gophers. The reserve lies on the San Andreas Fault and has unique soil chemistry owing to the emergence and erosion of the Earth’s mantle into the serpentine ultramafic soils that characterize the region. Serpentine soils have Ca/Mg ratios <1, as well as low levels of essential nutrients, high heavy metal content, and poor soil moisture retention.
We focused specifically on serpentine meadow habitat and observed 113 species in our plots from a potential pool of the 310 species that occur in all habitat types at the reserve; 73 of the 113 species were annual plants (SI Appendix, Table S1). The four most common species observed in unmanipulated plots were Vulpia microstachys (89% of sites), Hemizonia congesta (79%), Plantago erecta (69%), and Lasthenia californica (52%). Sixteen of the 113 total species could not be identified or classified as having an annual or perennial life history and were thus excluded from all analyses; however, these individuals occurred only at single sites and are thus unlikely to have large effects on diversity patterns.
Experimental Setup.
In May 2013, we surveyed the reserve for 30 serpentine meadow sites. Site locations were chosen such that they could be hierarchically grouped at five spatial scales, which resulted in three groups of five sites each at the northwestern and southeastern ends of the reserve (SI Appendix, Fig. S2). Within each site, we flagged and GPS-located eight 0.75- × 0.75-m plots arranged in a 2 × 4 block of plots, with plots separated by 1 m of bare ground. Each plot in a block was randomly assigned to receive a different experimental treatment (discussed below). In total, there were 240 plots (eight plots × 30 sites) and five spatial scales for comparison [at the level of the plot (1 m), a block of plots at a site (5 m), a group of sites (100 m), a reserve half (5 km), and a whole reserve (10 km)] (SI Appendix, Fig. S2). In late July 2013, after all winter annual species had senesced and the majority of summer annual species had set seed, we harvested all seed and standing vegetation from seven of the eight plots at each site using garden shears and a powerful gas-powered leaf vacuum (Stihl BG86); the eighth plot (C2; SI Appendix, Fig. S1B) was left unaltered to evaluate any unintended effects of the vacuuming procedure on plant diversity. All collected materials were stored outside in paper bags to allow natural heat stratification until they could be processed (<6 wk).
All collected materials were reintroduced to the seven vacuumed plots per site in two control and five species-pooling treatments (SI Appendix, Fig. S2B). The control treatments were as follows: a “vacuum without replacement” treatment (C1) to identify individuals that were left behind following vacuuming, and a “vacuum without movement” treatment (C3), where the collected material was homogenized at the plot level and redistributed back onto the source plot. The five spatial scaling treatments involved the pooling, homogenization, and redistribution of material collected from a single plot (1 m; blue plot in SI Appendix, Fig. S2B), multiple plots at a site (5 m; green plot), 5 sites of a single group (100 m; yellow plot), 15 sites from the same half of the reserve (5 km; orange plot), and all 30 sites across the entire reserve (10 km; red plot). Because treatments were nested within sites, and plots within sites are highly similar in environments and species composition, all treatments received a similar regional pool in terms of species richness and relative abundances. The redistributed material was secured to each recipient plot with twine.
We surveyed the plot-level community structure and corresponding environmental parameters during peak biomass in the following growing season. Plots were surveyed April 20 to May 2 in 2014, using percentage of cover estimates of each species because small-statured annuals can occur at densities up to 5,500 individuals per square meter (46); additional surveys were conducted later in the season to confirm the identities of late-flowering species. We surveyed the innermost 0.5 × 0.5 m2 of each plot to account for any edge effects in our analyses. We measured plot-level percentage of soil moisture content, understory photosynthetically active radiation (PAR) in full sun, and slope inclination, as well as site-level elevation, slope aspect, hillside slope steepness, and soil depth. We also performed site-level soil fertility analyses (NO3-N, Olsen-P, X-K, X-Na, X-Ca, X-Mg, pH, cation exchange capacity, organic matter; University of California at Davis Analytical Lab) on soil samples collected and pooled between four plots per site; X-Ca and X-Mg were converted to a ratio of Ca/Mg. Although site-level estimates of environmental conditions preclude finer-scale estimates, those variables for which we do have plot-level estimates showed no difference among plots that occur at the same site (all P > 0.35; SI Appendix, Table S2). This means that, although site-level estimates preclude estimates of within-site (among-plot) error in most of our environmental variables, all treatments and controls are subject to the same error because each is represented within all sites, and plots within sites do not significantly differ in environmental conditions (or species composition, Fig. 3).
Data Analysis.
A presence/absence matrix was created from percentage of cover estimates of all 73 annual species in 29 of the 30 sites that were sampled; data from one site were lost due to a corrupt data file. We used presence/absence data because it is most appropriate for tracking gains and losses of species in response to manipulations of species pools and is most comparable to other studies (14, 29).
We tested the responses of five components of diversity to our species pool manipulations: net changes in species richness, species richness, regional site occupancy, species composition (PCoA axis 1 and 2 scores), and beta diversity. The first four components were tested using separate linear mixed-effects models in the “lmerTest” R package, each with spatial species-pooling treatment as a fixed factor and plot nested within site as a random factor. The nested random effects account for the nonindependence of our nested experimental design in terms of both error structure and degrees of freedom (47). Because species-pooling treatment emerged as a significant predictor in all analyses, we used post hoc Tukey tests using the “multcomp” package to identify treatment levels that differed significantly from each other. Net changes in species richness were estimated by subtracting the control plot C3 (SI Appendix, Fig. S2B) presence/absence matrix from each species-pooling treatment matrix; doing so allowed us to decompose net changes into gains and losses of individual species (SI Appendix, Fig. S3B). Species composition was the first two axis scores of a PCoA using a Jaccard dissimilarity matrix (R package “vegan”; SI Appendix, Fig. S5A).
The last component, beta diversity, was calculated as a matrix of Jaccard dissimilarity coefficients among all pairwise site combinations and separately for each species-pooling treatment. We then performed a permutational multivariate analyses of variance (“adonis” function in R package vegan) to test differences in Jaccard dissimilarity among treatments, constrained to account for treatments being nested within sites. This analysis is analogous to a univariate analysis of variance, except expanded to handle multivariate response variables. Because the overall model was significant, we ran adonis with a Bonferroni correction on all pairwise treatment combinations to identify treatments that significantly differed.
Before testing the responses of species–environment associations to our species-pooling treatments, we used complementary multivariate methods in the R package vegan to reduce the dimensionality of the environmental variables that were measured. First, we used variance inflation factors (VIFs) to confirm that multicollinearity was low (VIFs < 10) and performed a constrained correspondence analysis (CCA), forward-selecting environmental variables using the R function “ordistep” to identify the subset of environmental variables that significantly influenced species occurrences (all but soil depth). Second, we ran a PCA on the subset of environmental variables identified as meaningful in the CCA, with the axis 1 and 2 scores of the PCA summarizing 51.4% of the variation in the among-site environmental dataset (SI Appendix, Fig. S5B).
We used the first (PCA1) and second (PCA2) PCA axis scores as composite environmental variables to examine how the species-pooling treatments affected the strength of species–environment associations. We performed a generalized linear model, separately for each species-pooling treatment, with species occurrences as the binomially distributed response variable and poly(PCA1,2)*species + poly(PCA2,2)*species as predictor variables. The poly(x,2) function in the R “stats” package calculates orthogonal first- and second-order polynomials of PCA1 and PCA2 to detect linear and quadratic relationships; species present in sites of intermediate environmental conditions but absent at the extremes would be best characterized by quadratic relationships. “Species” was a fixed factor in these analyses to facilitate model convergence and because we were directly interested in the amount of variation explained (pseudo-R2 [1-{residual deviance/null deviance}]) by species-specific differences in responses to environmental conditions (47). We then compared the amount of variation explained by composite environmental conditions among species-pooling treatments (18), with the prediction that pseudo-R2 values should increase with the spatial scale of species pooling. To further examine species-level responses (i.e., how strongly the occurrences of each species were associated with environmental conditions), we performed separate generalized linear models for each species and species-pooling treatment, testing the additive effects of poly(PCA1,2) and poly(PCA2,2). We visualized how the distribution of the pseudo-R2 values of the 73 species shifted across species-pooling treatments. We then visualized the distribution of pseudo-R2 values across species for each of the species-pooling treatments.
To quantify similarity among sites in environments and species that occurred naturally at distinct spatial scales, we performed permutational multivariate analyses of variance (adonis function in R package vegan). The response variables were matrices of Euclidian (for environment) and Jaccard (for species) distances among sites, and each matrix was tested for differences among groups of five sites (100-m scale) nested within each half of the reserve (5-km scale). The species distance matrix was created from plots that received the 1-m treatment only, as this treatment most closely reflected natural unmanipulated species distributions.
The control plots were used to assess the presence of any unintended effects of the hay transfer manipulation that were unrelated to the species-pooling treatment. We found no difference in species richness or composition between the unmanipulated (C2 in SI Appendix, Fig. S2) plots and our 1-m treatment plots (all P > 0.998), indicating that the hay transfer manipulation did not affect local diversity. Similarly, there was no difference among plots that received their own hay back (C3) and the 1-m treatment plots that received hay from a different adjacent plot (all P > 0.807); thus, the removal and replacement of hay did not affect plot diversity. As a result, we report only the results from the 1-m treatment plots. The removal without replacement control plot (C1) had significantly lower species richness than the unmanipulated control plot (C2) and all treatment plots (all P < 0.002), indicating that our seed vacuum was effective.
We incorporated percentage-of-cover data to assess whether our species-pooling manipulations caused small transient “sink” populations to establish, leading to biased diversity estimates. To do this, we revisited the original percentage-of-cover matrix from which the presence/absence matrix was generated and converted any percentage-of-cover entries less than or equal to four cutoff values (0.05, 0.10, 0.25, and 0.5%) to zero. For small-statured annual plants, 0.5% cover roughly translates into 10 individual plants; values greater than 0.5% risks removing small but stable populations. We then converted those four percentage-of-cover matrices to presence/absences and reanalyzed our species richness and species–environment association results (SI Appendix, Table S3).
Acknowledgments
We thank S. Harrison, B. Anacker, C. Kohler, C. Lee, E. Case, S. Copeland, P. Arenas, and L. Sergison for advice and/or assistance in the field; S.Y.S. and R.M.G. laboratory members for project feedback; and the editor, two anonymous reviewers, M. McPeek, and J. Stinchcombe for comments on an earlier version of the manuscript. This work was supported by the National Sciences and Engineering Research Council (NSERC)’s Canada Graduate Scholarship and Michael Smith Foreign Study Supplement, as well as Sigma Xi (R.M.G.), and by the California Agricultural Experiment Station (S.Y.S.).
Footnotes
- ↵1To whom correspondence should be addressed. Email: rgermain{at}zoology.ubc.ca.
Author contributions: R.M.G., S.Y.S., and B.G. designed research; R.M.G. performed research; R.M.G. and B.G. analyzed data; and R.M.G., S.Y.S., and B.G. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1615338114/-/DCSupplemental.
References
- ↵.
- Wiens JA
- ↵
- ↵.
- Bell G,
- Lechowicz MJ
- ↵
- ↵.
- Rosenzweig ML
- ↵.
- Tonkin JD,
- Death RG
- ↵.
- MacArthur RH,
- Wilson EO
- ↵.
- Hubbell SP
- ↵
- ↵
- ↵.
- Chase JM
- ↵
- ↵.
- Gilbert B,
- Levine JM
- ↵.
- Grainger TN,
- Gilbert B
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Germain RM, et al.
- ↵
- ↵.
- Coiffait-Gombault C,
- Buisson E,
- Dutoit T
- ↵
- ↵.
- Baldeck CA, et al.
- ↵
- ↵.
- Strauss SY,
- Webb CO,
- Salamin N
- ↵.
- Anacker BL
- ↵
- ↵.
- Szaro RC,
- Johnston DW
- Pearson S,
- Turner M,
- Gardner R,
- O’Neill R
- ↵.
- Chase JM,
- Kraft NJB,
- Smith KG,
- Velland M,
- Inouye B
- ↵
- .
- McGill BJ
- ↵
- ↵.
- Ramette A,
- Tiedje JM
- ↵.
- Levine JM
- ↵
- ↵.
- Tuomisto H,
- Ruokolainen K,
- Yli-Halla M
- ↵
- ↵
- ↵.
- Bennet JR,
- Gilbert B
- ↵
- ↵
- ↵
- ↵.
- Harrison S,
- Rajakaruna R
- ↵
- ↵
Citation Manager Formats
Article Classifications
- Biological Sciences
- Ecology