Microbial colonization dynamics in temporary aquatic systems

Most studies of community development in insular systems have investigated the colonization dynamics of only a portion of the biotic community using islands that have the benefit of prior biotic modification. Two experiments assessed the predictive ability of equilibrium island theory with regard to the development of the protistan community in temporary aquatic islands (100-1 plastic swimming pools) void of any prior biotic modification. The first experiment ran for 170 d (March–September 1985) and manipulated the access of islands to animals that represent potentially important dispersal pathways for microbes. A second experiment (May–June 1986) investigated in greater detail the early stages of species accrual in the absence of dispersal pathways other than via the atmosphere. Polyurethane foam substrates were used as sampling devices to increase the accuracy and replicability of sampling and provide habitat for colonists. Sampling was determined to be asymptotic. Species accrual was asymptotic in both experiments, although it initially lagged behind that predicted by theory. Autotrophs approached equilibrium faster than heterotrophs, but at a lower species richness. The predicted number of autotroph species at equilibrium was lower in islands that were excluded from contact with birds compared to less exclosed islands. The colonization dynamics of the entire community was not significantly different among islands having different degrees of exclosure. In both experiments, rates of species immigration were nonmonotonic with respect to time and species richness. This relationship appeared to reflect the importance of species interactions during the initial accrual phase before equilibrium. Rates of extinction were positively correlated with both these parameters, although they tended to decrease with time during the equilibrium period in less exclosed islands. Turnover at equilibrium was significant and resulted in directional changes in species composition over time. Assortative processes appeared to be important since later colonists exhibited greater persistence. Colonizing species generally fall into three autecological categories: 1) those that were ubiquitous and had temporally predictable patterns of immigration (successional species); 2) those possessing temporally predictable distributions but not spatially ubiquitous distributions (dispersal limited species); 3) those that showed little temporal or spatial predictability in immigration (transient or allochthonous species). Individual islands exhibited various degrees of fluctuation in species number during the predicted equilibrium which were poorly correlated with exogenous environmental variables and physicochemical habitat parameters. The presence of predacious mosquito larvae(Culex spp.) invariably resulted in a sharp decreases in microbial species richness, while documented contact with rodents was followed by an increase in species number in an island so contaminated. Several aspects of microbial colonization of temporary waters that contradict equilibrium predictions appear to be strongly influenced by microbe-microbe as well as macrobe-microbe interactions.     

Dispersal and species diversity: a meta-analysis

Species diversity in communities of interacting organisms is thought to be enhanced by dispersal, yet mechanisms predicting this have little to say about what effects differing rates of dispersal have on diversity and how dispersal affects diversity at larger spatial scales. I performed meta‐analyses on 23 studies comprising 50 experiments that manipulated species migration and measured community richness or diversity to test three hypotheses: that dispersal increases local diversity; that this effect depends on the rate of dispersal, specifically, that local diversity should be maximized at intermediate dispersal rates or else linearly related to dispersal rate; and that regional diversity may be either unaffected or negatively impacted by dispersal because dispersal tends to homogenize local communities. I found that immigration increased local diversity. Further, in animal studies, diversity appears maximized at intermediate dispersal rates but not with plant studies; however, more standardized studies are needed. Finally, results are ambiguous as to what happens at larger scales, with studies finding either declines or no change in regional diversity with dispersal. Taken together, these results reveal that dispersal has a complex, spatially contingent relationship with patterns of species diversity.     

High seed dispersal rates in faunally intact tropical rain forest: theoretical and conservation implications

Dispersal limitation may promote high tree-species diversity in rain forest by slowing local competitive exclusion, but evidence is scarce. By comparing the species lists of tree plots with those of nested seedling plots in rain forest at Gunung Palung, we found that the proportion of local seedlings arising from active dispersal events (via wind or animals, but not simply dropped from parent trees) was high: 68% of species and 46% of the individual seedlings. Local seedling species richness was not limited by the local richness of adults. Using these data, a spatially explicit simulation model indicated that dispersal limitation may not contribute substantially to the maintenance of tree diversity in this forest. We predict that the loss of animal seed dispersers would reduce local species richness of seedlings by 60%. While this reduction could possibly increase the influence of dispersal limitation, it would interfere with other mechanisms maintaining diversity.     

Determinants of species richness, endemism and turnover in European longhorn beetles

This study assessed the diversity patterns of a large family of beetles, Cerambycidae, in Europe and tested the following hypotheses: 1) richness gradients of this hyperdiverse taxon are driven by water and energy variables; 2) endemism is explained by the same factors, but variation between areas also reflects post‐glacial re‐colonization processes; and 3) faunal composition is determined by the same climatic variables and, therefore, beta diversity (species turnover) is related to richness gradients. Species richness, endemism and beta diversity were modelled using inventories of 37 European territories, built from a database containing the distributions of 609 species. Area, spatial position, and nine topographical and climatic variables were used as predictors in regression and constrained analysis of principal coordinates modelling. Species richness was mostly explained by a temperature gradient, which produced a south‐to‐north decreasing richness gradient. Endemism followed the same pattern, but was also determined by longitudinal variation, peaking in the southwestern and southeastern corners of the continent. Faunal turnover was explained by an important purely spatial pattern and a spatially structured environmental gradient. Thus, contrary to other groups, cerambycid richness was mostly explained by environmental energy, but not by water availability. Endemism was concentrated in the Iberian and Greek peninsulas, but not in Italy. Thus, the latter area may have been the major source of post‐glacial re‐colonization for European longhorn beetles or, otherwise, a poor refuge during glaciations. Turnover patterns were independent of the richness gradient, because northern faunas are nested in southern ones. Turnover, in contrast to richness, was driven by both the independent effects of climate and geographic constraints that might reflect dispersal limitation or stochastic colonization events, suggesting that richness gradients are more environmentally deterministic phenomena than turnover patterns.     

Spatial heterogeneity, source-sink dynamics, and the local coexistence of competing species

Patch occupancy theory predicts that a trade-off between competition and dispersal should lead to regional coexistence of competing species. Empirical investigations, however, find local coexistence of superior and inferior competitors, an outcome that cannot be explained within the patch occupancy framework because of the decoupling of local and spatial dynamics. We develop two-patch metapopulation models that explicitly consider the interaction between competition and dispersal. We show that a dispersal-competition trade-off can lead to local coexistence provided the inferior competitor is superior at colonizing empty patches as well as immigrating among occupied patches. Immigration from patches that the superior competitor cannot colonize rescues the inferior competitor from extinction in patches that both species colonize. Too much immigration, however, can be detrimental to coexistence. When competitive asymmetry between species is high, local coexistence is possible only if the dispersal rate of the inferior competitor occurs below a critical threshold. If competing species have comparable colonization abilities and the environment is otherwise spatially homogeneous, a superior ability to immigrate among occupied patches cannot prevent exclusion of the inferior competitor. If, however, biotic or abiotic factors create spatial heterogeneity in competitive rankings across the landscape, local coexistence can occur even in the absence of a dispersal-competition trade-off. In fact, coexistence requires that the dispersal rate of the overall inferior competitor not exceed a critical threshold. Explicit consideration of how dispersal modifies local competitive interactions shifts the focus from the patch occupancy approach with its emphasis on extinction-colonization dynamics to the realm of source-sink dynamics. The key to coexistence in this framework is spatial variance in fitness. Unlike in the patch occupancy framework, high rates of dispersal can undermine coexistence, and hence diversity, by reducing spatial variance in fitness.     

Global patterns of amphibian phylogenetic diversity

Aim Phylogenetic diversity can provide insight into how evolutionary processes may have shaped contemporary patterns of species richness. Here, we aim to test for the influence of phylogenetic history on global patterns of amphibian species richness, and to identify areas where macroevolutionary processes such as diversification and dispersal have left strong signatures on contemporary species richness. Location Global; equal‐area grid cells of approximately 10,000 km². Methods We generated an amphibian global supertree (6111 species) and repeated analyses with the largest available molecular phylogeny (2792 species). We combined each tree with global species distributions to map four indices of phylogenetic diversity. To investigate congruence between global spatial patterns of amphibian species richness and phylogenetic diversity, we selected Faith’s phylogenetic diversity (PD) index and the total taxonomic distinctness (TTD) index, because we found that the variance of the other two indices we examined (average taxonomic distinctness and mean root distance) strongly depended on species richness. We then identified regions with unusually high or low phylogenetic diversity given the underlying level of species richness by using the residuals from the global relationship of species richness and phylogenetic diversity. Results Phylogenetic diversity as measured by either Faith’s PD or TTD was strongly correlated with species richness globally, while the other two indices showed very different patterns. When either Faith’s PD or TTD was tested against species richness, residuals were strongly spatially structured. Areas with unusually low phylogenetic diversity for their associated species richness were mostly on islands, indicating large radiations of few lineages that have successfully colonized these archipelagos. Areas with unusually high phylogenetic diversity were located around biogeographic contact zones in Central America and southern China, and seem to have experienced high immigration or in situ diversification rates, combined with local persistence of old lineages. Main conclusions We show spatial structure in the residuals of the relationship between species richness and phylogenetic diversity, which together with the positive relationship itself indicates strong signatures of evolutionary history on contemporary global patterns of amphibian species richness. Areas with unusually low and high phylogenetic diversity for their associated richness demonstrate the importance of biogeographic barriers to dispersal, colonization and diversification processes.     

Patch quality and context, but not patch number, drive multi-scale colonization dynamics in experimental aquatic landscapes

Colonization and extinction are primary drivers of local population dynamics, community structure, and spatial patterns of biological diversity. Existing paradigms of island biogeography, metapopulation biology, and metacommunity ecology, as well as habitat management and conservation biology based on those paradigms, emphasize patch size, number, and isolation as primary characteristics influencing colonization and extinction. Habitat selection theory suggests that patch quality could rival size, number, and isolation in determining rates of colonization and resulting community structure. We used naturally colonized experimental landscapes to address four issues: (a) how do colonizing aquatic beetles respond to variation in patch number, (b) how do they respond to variation in patch quality, (c) does patch context affect colonization dynamics, and (d) at what spatial scales do beetles respond to habitat variation? Increasing patch number had no effect on per patch colonization rates, while patch quality and context were critical in determining colonization rates and resulting patterns of abundance and species richness at multiple spatial scales. We graphically illustrate how variation in immigration rates driven by perceived predation risk (habitat quality) can further modify dynamics of the equilibrium theory of island biogeography beyond predator-driven effects on extinction rates. Our data support the importance of patch quality and context as primary determinants of colonization rate, occupancy, abundance, and resulting patterns of species richness, and reinforce the idea that management of metapopulations for species preservation, and metacommunities for local and regional diversity, should incorporate habitat quality into the predictive equation.     

Patterns of pollen dispersal in a small population of the Canarian endemic palm (Phoenix canariensis)

The genetic diversity of small populations is greatly influenced by local dispersal patterns and genetic connectivity among populations, with pollen dispersal being the major component of gene flow in many plants species. Patterns of pollen dispersal, mating system parameters and spatial genetic structure were investigated in a small isolated population of the emblematic palm Phoenix canariensis in Gran Canaria island (Canary Islands). All adult palms present in the study population (n=182), as well as 616 seeds collected from 22 female palms, were mapped and genotyped at 8 microsatellite loci. Mating system analysis revealed an average of 5.8 effective pollen donors (N_ep) per female. There was strong variation in correlated paternity rates across maternal progenies (ranging from null to 0.9) that could not be explained by the location and density of local males around focal females. Paternity analysis revealed a mean effective pollen dispersal distance of ∼71 m, with ∼70% of effective pollen originating from a distance of <75 m, and 90% from <200 m. A spatially explicit mating model indicated a leptokurtic pollen dispersal kernel, significant pollen immigration (12%) from external palm groves and a directional pollen dispersal pattern that seems consistent with local altitudinal air movement. No evidence of inbreeding or genetic diversity erosion was found, but spatial genetic structure was detected in the small palm population. Overall, the results suggest substantial pollen dispersal over the studied population, genetic connectivity among different palm groves and some resilience to neutral genetic erosion and subsequently to fragmentation.     

Predicting local–regional richness relationships using island biogeography models

Local species richness frequently is linearly related to the richness of the regional species pool from which the local community was presumably assembled. What, if anything, does this pattern imply about the relative importance of species interactions and dispersal as determinants of local species richness? Two recent papers by Hugueny and Cornell and He et al. propose that the classical island biogeography model of MacArthur and Wilson can help answer this question, by serving as a null model of the relationship between local (island) and regional (mainland) species richness in the absence of local species interactions. The two models make very different predictions, despite being derived from apparently‐similar assumptions. Here we reinterpret these two models and show that their contrasting predictions can be regarded as arising from different, implicit assumptions about how species abundances vary with species richness on the mainland. We derive a more general island biogeography model of local–regional richness relationships that explicitly incorporates mainland species abundance and subsumes the two previous models as limiting cases. The new model predicts that the local–regional richness relationship can range from nearly linear to strongly curvilinear, depending on how species abundances on the mainland vary with mainland richness, as well as on rates of immigration to and extinction from islands. Local species interactions are not necessary for producing curvilinear local–regional richness relationships. We discuss the implications of our new model for the interpretation of local–regional richness relationships.     

Regional zooplankton dispersal provides spatial insurance for ecosystem function

Changing environmental conditions are affecting diversity and ecosystem function globally. Theory suggests that dispersal from a regional species pool may buffer against changes in local community diversity and ecosystem function after a disturbance through the establishment of functionally redundant tolerant species. The spatial insurance provided by dispersal may decrease through time after environmental change as the local community monopolizes resources and reduces community invasibility. To test for evidence of the spatial insurance hypothesis and to determine the role dispersal timing plays in this response we conducted a field experiment using crustacean zooplankton communities in a subarctic region that is expected to be highly impacted by climate change – Churchill, Canada. Three experiments were conducted where nutrients, salt, and dispersal were manipulated. The three experiments differed in time‐since‐disturbance that the dispersers were added. We found that coarse measures of diversity (i.e. species richness, evenness, and Shannon–Weiner diversity) were generally resistant to large magnitude disturbances, and that dispersal had the most impact on diversity when dispersers were added shortly after disturbance. Ecosystem functioning (chl‐a) was degraded in disturbed communities, but dispersal recovered ecosystem function to undisturbed levels. This spatial insurance for ecosystem function was mediated through changes in community composition and the relative abundance of functional groups. Results suggest that regional diversity and habitat connectivity will be important in the future to maintain ecosystem function by introducing functionally redundant species to promote compensatory dynamics.     

Dispersal and establishment: spatial patterns and species–area relationships

According to the equilibrium theory of island biogeography, high colonization ability of species is associated with low exponents (z) of the species–area relationship (SAR) and weak spatial patterns in species number and dissimilarity. However, the relationship between z and the strength of these spatial patterns has not been investigated systematically. We used a multispecies metapopulation model to investigate these relationships in an archipelago of islands. We conclude that this relationship can only be predicted if either the dispersal ability or the power of establishment of species is known. With species richness limited by establishment, we generated high z‐values associated with weak spatial patterns in species number and dissimilarity. If species richness was constrained by the dispersal ability of species, we observed low to medium z‐values but strong spatial patterns. If the dispersal ability and the abilities of species to establish were both high, z‐values and spatial pattern tend to be low and species numbers were limited by the size of the regional species pool.     

Genetic influences on social dominance: cow wars

Analyzing population dynamics in changing habitats is a prerequisite for population dynamics forecasting. The recent development of metapopulation modeling allows the estimation of dispersal kernels based on the colonization pattern but the accuracy of these estimates compared with direct estimates of the seed dispersal kernel has rarely been assessed. In this study, we used recent genetic methods based on parentage analysis (spatially explicit mating models) to estimate seed and pollen dispersal kernels as well as seed and pollen immigration in fragmented urban populations of the plant species Crepis sancta with contrasting patch dynamics. Using two independent networks, we documented substantial seed immigration and a highly restricted dispersal kernel. Moreover, immigration heterogeneity among networks was consistent with previously reported metapopulation dynamics, showing that colonization was mainly due to external colonization in the first network (propagule rain) and local colonization in the second network. We concluded that the differences in urban patch dynamics are mainly due to seed immigration heterogeneity, highlighting the importance of external population source in the spatio-temporal dynamics of plants in a fragmented landscape. The results show that indirect and direct methods were qualitatively consistent, providing a proper interpretation of indirect estimates. This study provides attempts to link genetic and demographic methods and show that patch occupancy models may provide simple methods for analyzing population dynamics in heterogeneous landscapes in the context of global change.     

Congruence,but no cascade—Pelagic biodiversity across three trophic levels in Nordic lakes

Covariation in species richness and community structure across taxonomical groups (cross‐taxon congruence) has practical consequences for the identification of biodiversity surrogates and proxies, as well as theoretical ramifications for understanding the mechanisms maintaining and sustaining biodiversity. We found there to exist a high cross‐taxon congruence between phytoplankton, zooplankton, and fish in 73 large Scandinavian lakes across a 750 km longitudinal transect. The fraction of the total diversity variation explained by local environment alone was small for all trophic levels while a substantial fraction could be explained by spatial gradient variables. Almost half of the explained variation could not be resolved between local and spatial factors, possibly due to confounding issues between longitude and landscape productivity. There is strong consensus that the longitudinal gradient found in the regional fish community results from postglacial dispersal limitations, while there is much less evidence for the species richness and community structure gradients at lower trophic levels being directly affected by dispersal limitation over the same time scale. We found strong support for bidirectional interactions between fish and zooplankton species richness, while corresponding interactions between phytoplankton and zooplankton richness were much weaker. Both the weakening of the linkage at lower trophic levels and the bidirectional nature of the interaction indicates that the underlying mechanism must be qualitatively different from a trophic cascade.     

Species richness patterns and metapopulation processes – evidence from epiphyte communities in boreo-nemoral forests

For several epiphyte species, dispersal limitation and metapopulation dynamics have been suggested. We studied the relative importance of local environmental conditions and spatial aggregation of species richness of facultative and obligate epiphytic bryophytes and lichens within two old‐growth forests in eastern Sweden. The effect of the local environment was analyzed using generalized linear models (GLM). We tested whether species richness was spatially structured by fitting variogram models to the residuals of the GLM. In addition, we analyzed the species‐area relationship (area=tree diameter). Different environmental variables explained the richness of different species groups (bryophytes vs lichens, specialists vs generalists, sexual vs asexual dispersal). In most groups, the total variation explained by environmental variables was higher than the variation explained by the spatial model. Spatial aggregation was more pronounced in asexually than in sexually dispersed species. Bryophyte species richness was only poorly predicted by area, and lichen species richness was not explained by area at all. Spatial aggregation may indicate effects of dispersal limitation and metapopulation dynamics on community species richness. Our results suggest that species groups differ in habitat requirements and dispersal abilities; there were indications that presence of species with different dispersal strategies is linked to the age of the host tree. Separate analyses of the species richness of species groups that differ in the degree of habitat specialization and dispersal ability give insights into the processes determining community species richness. The poor species‐area relationship, especially in lichens, may indicate species turnover rather than accumulation during the lifetime of the host tree. Epiphyte species extinctions may be mainly caused by deterministic processes, e.g. changes in habitat conditions as the host tree grows, ages and dies, rather than by stochastic population processes.     

Habitat diversity associated to island size and environmental filtering control the species richness of rock‐savanna plants in neotropical inselbergs

Disentangling the multiple factors controlling species diversity is a major challenge in ecology. Island biogeography and environmental filtering are two influential theories emphasizing respectively island size and isolation, and the abiotic environment, as key drivers of species richness. However, few attempts have been made to quantify their relative importance and investigate their mechanistic basis. Here, we applied structural equation modelling, a powerful method allowing test of complex hypotheses involving multiple and indirect effects, on an island‐like system of 22 French Guianan neotropical inselbergs covered with rock‐savanna. We separated the effects of size (rock‐savanna area), isolation (density of surrounding inselbergs), environmental filtering (rainfall, altitude) and dispersal filtering (forest‐matrix openness) on the species richness of all plants and of various ecological groups (terrestrial versus epiphytic, small‐scale versus large‐scale dispersal species). We showed that the species richness of all plants and terrestrial species was mainly explained by the size of rock‐savanna vegetation patches, with increasing richness associated with higher rock‐savanna area, while inselberg isolation and forest‐matrix openness had no measurable effect. This size effect was mediated by an increase in terrestrial‐habitat diversity, even after accounting for increased sampling effort. The richness of epiphytic species was mainly explained by environmental filtering, with a positive effect of rainfall and altitude, but also by a positive size effect mediated by enhanced woody‐plant species richness. Inselberg size and environmental filtering both explained the richness of small‐scale and large‐scale dispersal species, but these ecological groups responded in opposite directions to altitude and rainfall, that is positively for large‐scale and negatively for small‐scale dispersal species. Our study revealed both habitat diversity associated with island size and environmental filtering as major drivers of neotropical inselberg plant diversity and showed the importance of plant species growth form and dispersal ability to explain the relative importance of each driver.     

When do localized natural enemies increase species richness?

The Janzen–Connell hypothesis states that local species‐specific density dependence, mediated through specialist enemies of offspring such as fungal pathogens and insect seed predators, can facilitate coexistence of species by preventing recruitment near conspecific adults. We use spatially explicit simulation models and analytical approximations to evaluate how spatial scales of offspring and enemy dispersal affect species richness. In comparison with model communities in which both offspring and enemies disperse long distances, species richness is substantially decreased when offspring disperse long distances and enemies disperse short distances. In contrast, when both offspring and enemies disperse short distances species richness more than doubles and adults of each species are highly spatially clumped. For the range of conditions typical of tropical forests, locally dispersing specialist enemies may decrease species richness relative to enemies that disperse long distances. In communities where dispersal distances of both offspring and enemies are short, local effects may enhance species richness.     

Clumped versus scattered: how does the spatial correlation of disturbance events affect biodiversity?

In this study, we systematically explore the effects of rate and spatial correlation (level of clumping) of disturbance events on a community of sessile species differing in their life history traits. A spatially explicit individual-based model shows that long-term coexistence is very sensitive to spatial correlation when the trade-off in life history traits includes differences in dispersal distances. Highest biodiversity emerges at highly correlated disturbances of intermediate rates. Diversity peaks shift to larger rates when clumping decreases. Scattered disturbances lead to competitive exclusion. Interestingly, we observed additional peaks in the diversity–disturbance curves at certain levels of clumping. Thus, subject to the differences in life history traits, particular combinations of disturbance rate and spatial correlation may enable subsets of species to coexist, which opens new possibilities for explaining diversity. Our results suggest that observation of high biodiversity under spatially correlated disturbances points to a competition–colonisation trade-off, which includes dispersal distances.     

Spatial averaging and disturbance lead to high productivity in aquatic metacommunities

Dispersal in heterogeneous ecosystems, such as coastal metacommunities, is a major driver of diversity and productivity. According to theory, both species richness and spatial averaging shape a unimodal relationship of productivity with dispersal. We experimentally tested the hypothesis that disturbances acting on local patches would buffer the loss of productivity at high dispersal by preventing synchronized species oscillations. To simulate these disturbances, our experimental assemblages involved species that self‐organized in isolation under three inflow pulsing frequencies, where hydraulic displacement and nutrient loading affected assemblage diversity and composition. At steady‐state, the emerging isolated assemblages were connected at three levels of dispersal creating three metacommunities of different connectivity. Consistent with theory, as dispersal increased, species richness in the metacommunity declined; productivity however remained high. This occurred because the most productive species in our study (which dominated the isolated patch of intermediate inflow pulsing frequency) dominated all three patches (low, intermediate and high inflow pulsing frequencies) after dispersal commenced in our metacommunities. This experimental result provides empirical support for the mechanism of spatial averaging. Furthermore, disturbances, in the form of localized pulsed inflows, prevented population oscillation synchrony caused by homogenization. Overall, our observations suggest that localized environmental fluctuations and the identity of species seem to be more influential than dispersal in shaping the diversity and composition of phytoplankton assemblages and stabilizing productivity.     

Positive interactions and the emergence of community structure in metacommunities

The significant role of space in maintaining species coexistence and determining community structure and function is well established. However, community ecology studies have mainly focused on simple competition and predation systems, and the relative impact of positive interspecific interactions in shaping communities in a spatial context is not well understood. Here we employ a spatially explicit metacommunity model to investigate the effect of local dispersal on the structure and function of communities in which species are linked through an interaction web comprising mutualism, competition and exploitation. Our results show that function, diversity and interspecific interactions of locally linked communities undergo a phase transition with changes in the rate of species dispersal. We find that low spatial interconnectedness favors the spontaneous emergence of strongly mutualistic communities which are more stable but less productive and diverse. On the other hand, high spatial interconnectedness promotes local biodiversity at the expense of local stability and supports communities with a wide range of interspecific interactions. We argue that investigations of the relationship between spatial processes and the self-organization of complex interaction webs are critical to understanding the geographic structure of interactions in real landscapes.     

Geographical,Temporal and Environmental Determinants of Bryophyte Species Richness in the Macaronesian Islands

Species richness on oceanic islands has been related to a series of ecological factors including island size and isolation (i.e. the Equilibrium Model of Island Biogeography, EMIB), habitat diversity, climate (i.e., temperature and precipitation) and more recently island ontogeny (i.e. the General Dynamic Model of oceanic island biogeography, GDM). Here we evaluate the relationship of these factors with the diversity of bryophytes in the Macaronesian region (Azores, Madeira, Canary Islands and Cape Verde). The predictive power of EMIB, habitat diversity, climate and the GDM on total bryophyte richness, as well as moss and liverwort richness (the two dominant bryophyte groups), was evaluated through ordinary least squares regressions. After choosing the best subset of variables using inference statistics, we used partial regression analyses to identify the independent and shared effects of each model. The variables included within each model were similar for mosses and liverworts, with orographic mist layer being one of the most important predictors of richness. Models combining climate with either the GDM or habitat diversity explained most of richness variation (up to 91%). There was a high portion of shared variance between all pairwise combinations of factors in mosses, while in liverworts around half of the variability in species richness was accounted for exclusively by climate. Our results suggest that the effects of climate and habitat are strong and prevalent in this region, while geographical factors have limited influence on Macaronesian bryophyte diversity. Although climate is of great importance for liverwort richness, in mosses its effect is similar to or, at least, indiscernible from the effect of habitat diversity and, strikingly, the effect of island ontogeny. These results indicate that for highly vagile taxa on oceanic islands, the dispersal process may be less important for successful colonization than the availability of suitable ecological conditions during the establishment phase.