What drives elevational patterns of diversity? A test of geometric constraints,climate and species pool effects for pteridophytes on an elevational gradient in Costa Rica

Aim We studied pteridophyte species richness between 100 m and 3400 m along a Neotropical elevational gradient and tested competing hypotheses for patterns of species richness. Location Elevational transects were situated at Volcán Barva in the Braulio Carrillo National Park and La Selva Biological Station (100–2800 m) and Cerro de la Muerte (2700–3400 m), both on the Atlantic slope of Costa Rica, Central America. Method We analysed species richness on 156 plots of 20 × 20 m and measured temperature and humidity at four elevations (40, 650, 1800 and 2800 m). Species richness patterns were regressed against climatic variables (temperature, humidity, precipitation and actual evapotranspiration), regional species pool, area and predicted species number of a geometric null model (the mid‐domain effect, MDE). Results The species richness of the 484 recorded species showed a hump‐shaped pattern with elevation with a richness peak at mid‐elevations (c. 1700 m). The MDE was the single most powerful explanatory variable in linear regression models, but species richness was also associated strongly with climatic variables, especially humidity and temperature. Area and species pool were associated less strongly with observed richness patterns. Main conclusions Geometric models and climatic models exclusive of geometric constraints explained comparable amounts of the elevational variation in species richness. Discrimination between these two factor complexes is not possible based on model fits. While overall fits of geometric models were high, large‐ and small‐ranged species were explained by geometric models to different extents. Species with narrow elevational ranges clustered at both ends of the gradient to a greater extent than predicted by the MDE null models used here. While geometric models explained much of the pattern in species richness, we cannot rule out the role of climatic factors (or vice versa) because the predicted peak in richness from geometric models, the empirical peak in richness and the overlap in favourable environmental conditions all coincide at middle elevations. Mid‐elevations offer highest humidity and moderate temperatures, whereas at high elevations richness is reduced due to low temperatures, and at low elevations by reduced water availability due to high temperatures.     

Continental and regional ranges of North American mammals: Rapoport's rule in real and null worlds

Aim To assess the relationship between species richness and distribution within regions arranged along a latitudinal gradient we use the North American mammalian fauna as a study case for testing theoretical models. Location North America. Methods We propose a conceptual framework based on a fully stochastic mid‐domain model to explore geographical patterns of range size and species richness that emerge when the size and position of species ranges along a one‐dimensional latitudinal gradient are randomly generated. We also analyse patterns for the mammal fauna of North America by comparing empirical results from a biogeographical data base with predictions based on randomization null models. Results We confirmed the validity of Rapoport's rule for the mammals of North America by documenting gradients in the size of the continental ranges of species. Additionally, we demonstrated gradients of mean regional range size that parallel those of continental range. Our data also demonstrated that mean range size, measured both as a continental or a regional variable, is significantly correlated with the geographical pattern in species richness. All these patterns deviated sharply from null models. Main conclusions Rapoport's statement of an areographic relationship between species distribution and richness is highly relevant in modern discussions about ecological patterns at the geographical scale.     

Midpoint attractors and species richness: Modelling the interaction between environmental drivers and geometric constraints

We introduce a novel framework for conceptualising, quantifying and unifying discordant patterns of species richness along geographical gradients. While not itself explicitly mechanistic, this approach offers a path towards understanding mechanisms. In this study, we focused on the diverse patterns of species richness on mountainsides. We conjectured that elevational range midpoints of species may be drawn towards a single midpoint attractor – a unimodal gradient of environmental favourability. The midpoint attractor interacts with geometric constraints imposed by sea level and the mountaintop to produce taxon‐specific patterns of species richness. We developed a Bayesian simulation model to estimate the location and strength of the midpoint attractor from species occurrence data sampled along mountainsides. We also constructed midpoint predictor models to test whether environmental variables could directly account for the observed patterns of species range midpoints. We challenged these models with 16 elevational data sets, comprising 4500 species of insects, vertebrates and plants. The midpoint predictor models generally failed to predict the pattern of species midpoints. In contrast, the midpoint attractor model closely reproduced empirical spatial patterns of species richness and range midpoints. Gradients of environmental favourability, subject to geometric constraints, may parsimoniously account for elevational and other patterns of species richness.     

Biogeographical patterns and Rapoport's rule in southeastern Pacific benthic polychaetes of the Chilean coast

Recently three biogeographical units were identified along the Chilean coast (the Magellanic Province, an Intermediate Area, and the Peruvian Province), however few studies have focused on the factors and dynamic processes that formed these spatial units (e.g. Rapoport's rule and its causal mechanisms). In this study we used benthic polychaetes of the Chilean coast to evaluate patterns of latitudinal distribution and species richness, and the existence of the three main biogeographical provinces described for the Chilean coast. Additionally, we evaluated the latitudinal Rapoport effects and geometric constraint as a null hypothesis explaining the species richness distribution.
We found that benthic polychaete diversity increased towards southern latitudes. The cluster and ordination (non-metric MultiDimensional Scaling, nMDS) analyses of the distribution data, presented only two statistically significant (bootstrapping techniques) biogeographic provinces along the Chilean coast, with a break occurring between 41° and 42°S. While, our results did not support a latitudinal Rapoport effect, they did support the view that latitudinal Rapoport effects are a local phenomenon, occurring only for the Northeastern Pacific marine taxa. The relationship between latitudinal range extent and mean latitude indicated the existence of two hard boundaries at either extreme of the Chilean coast, limiting the geographical ranges of the species. However, geometric constraints tested using a Monte Carlo simulation approach showed a weak level of mid-domain effect on species richness. Finally, we propose that geometric constraint together with the geomorphology and historical characteristics of the Chilean coast explain the biogeographical patterns of benthic polychaete taxa in Chile.     

The mid-domain effect applied to elevational gradients: species richness of small mammals in Costa Rica

Aim The objective of this study was to comprehensively document and examine the alpha and gamma patterns of species richness in non-volant, small mammals (rodents, shrews and mouse opossums) along a tropical elevational gradient. These data were used to determine the support for existing hypotheses of species richness encompassing mid-domain null models, as well as climatic, and community overlap hypotheses. Location Field studies were conducted along a Caribbean slope of the Río Peñas Blancas watershed in the north-eastern region of Costa Rica between 750 and 1850 m at 10 sampling sites. Methods Species richness and abundances of small mammals were surveyed for four seasons including three temporal replicates at each of five elevational sites: late wet season (2000), early wet season (2001), and dry season (2002), and one spatial replicate at five different sites within the same elevations during the late wet season (2001). Species richness at elevations below 700 m was compiled from specimen records from 23 US national and international collections. Predictions of a null model based solely on geometric constraints were examined using a Monte Carlo simulation program, Mid-Domain Null. Results In 16,900 trap nights, 1561 individuals from 16 species were captured. Both alpha and gamma species richness peaked at mid-elevation between 1000 and 1300 m, with richness declining both at higher and lower elevations. Most of the empirical curves of species richness occur within 95% prediction curves of the mid-domain model, although deviations from the null model exist. Regression of the empirical richness on the null model predictions explained nearly half of the variation observed (r² = 0.45, P = 0.002). Main conclusions The geometric constraints of montane topography appear to influence the diversity pattern of small mammals, although climatic conditions including an intermediate rainfall and temperature regime, and distance from the persistent cloud cap also are correlated with the pattern of species richness. The predictions of productivity, and community overlap hypotheses are not supported with the empirical data.     

Comment on ‘Island biogeography: patterns of marine shallow‐water organisms’ by Hachich et al., Journal of Biogeography (2015)

In a recent article, Hachich et al. (2015, Journal of Biogeography, 42 , 1871–1882) studied the large‐scale biogeographical patterns of the species–area, species–island age and species–isolation relationships associated with marine shallow‐water groups (reef fish, gastropods and seaweeds) from 11 Atlantic archipelagos. We here express our concerns regarding the data accuracy used to compute the different models that tested the null hypothesis of species richness being independent of the selected variables. In our commentary, we focus mainly on the use of out‐of‐date checklists of gastropod and seaweed species from different archipelagos, but we also point out inaccuracies in some island age estimates and explain our disagreement with the use of the 200 m depth limit for the shallow‐water gastropods and seaweeds.     

Biogeographical patterns of marine larval trematode parasites in two intermediate snail hosts in Europe

Aim We used published inventories of trematodes in Littorina littorea (L.) and Hydrobia ulvae (Pennant) in European seas to search for two basic biogeographical patterns in the spatial occurrence of various trematode species: (1) do parasite distribution and richness patterns in the two host snails overlap with known ecoregions of free‐living organisms; and (2) does trematode species richness in the snails follow latitudinal or longitudinal gradients? Location North East Atlantic. Methods We used multidimensional scaling (MDS), analysis of similarity (ANOSIM) and analysis of variance (ANOVA) to test whether there were overlaps of parasite distribution and richness with known ecoregions of free‐living organisms. In addition, we used linear regression analyses to test whether trematode richness in snails (corrected for sampling effort) was correlated with the latitude or longitude of the sampling sites. Results When corrected for sampling effort, mean trematode species richness per site did not differ among the different ecoregions in L. littorea. In contrast, in H. ulvae, mean species richness was much lower for sites from the Celtic Sea compared with sites from the Baltic Sea and the North Sea. Based on the results of MDS analyses, trematode species composition was distinct among ecoregions; in particular, communities from the Baltic Sea differed markedly from communities in the Celtic Sea, for both snail species. Latitude and longitude were not significantly correlated with parasite species richness in either snail species. Most trematode species had restricted distributions, and only three species in L. littorea and five species in H. ulvae occurred at more than 50% of the sites. Main conclusions There is more structure in the large‐scale distribution of trematodes in gastropods than one would expect from the large‐scale dispersal capabilities of their bird and fish final hosts. We propose mechanisms based both on limited dispersal via fish and bird final hosts and on gradients in environmental factors to explain the observed patterns.     

Atlas versus range maps: robustness of chorological relationships to distribution data types in European mammals

Aim Chorological relationships describe the patterns of distributional overlap among species. In addition to revealing biogeographical structure, the resulting clusters of species with similar geographical distributions can serve as natural units in conservation planning. Here, we assess the extent to which temporal, methodological and taxonomical differences in the source of species’ distribution data can affect the relationships that are found. Location Western Europe. Methods We used two data sets – the Atlas of European mammals and polygon range maps from the IUCN Global Mammal Assessment – both as presence–absence data for UTM 50 km × 50 km squares. We performed pairwise comparisons among 156 species for each data set to build matrices of the similarity in distribution across species, using both Jaccard’s and Baroni‐Urbani & Buser’s indices. We then compared these similarity matrices (chorological relationships), as well as the species richness and occurrence patterns from the two data sets. Results As expected, range maps increased both the mean prevalence per species and mean species richness per grid cell in comparison to atlas data, reflecting the general view that these data types respectively over‐ and underestimate species occurrence. However, species richness and occurrence patterns in atlas and range map data were positively associated and, most importantly, the chorological relationships underlying the two data sets were highly similar. Main conclusions Despite many methodological, temporal and taxonomical differences between atlas data and range maps, the chorological relationships encountered between species were similar for both data sets. Chorological analyses can thus be robust to the data source used and provide a solid basis for analytical biogeographical studies, even over broad spatial scales.     

Type and spatial structure of distribution data and the perceived determinants of geographical gradients in ecology: the species richness of African birds

Aim Studies exploring the determinants of geographical gradients in the occurrence of species or their traits obtain data by: (1) overlaying species range maps; (2) mapping survey‐based species counts; or (3) superimposing models of individual species’ distributions. These data types have different spatial characteristics. We investigated whether these differences influence conclusions regarding postulated determinants of species richness patterns. Location Our study examined terrestrial bird diversity patterns in 13 nations of southern and eastern Africa, spanning temperate to tropical climates. Methods Four species richness maps were compiled based on range maps, field‐derived bird atlas data, logistic and autologistic distribution models. Ordinary and spatial regression models served to examine how well each of five hypotheses predicted patterns in each map. These hypotheses propose productivity, temperature, the heat–water balance, habitat heterogeneity and climatic stability as the predominant determinants of species richness. Results The four richness maps portrayed broadly similar geographical patterns but, due to the nature of underlying data types, exhibited marked differences in spatial autocorrelation structure. These differences in spatial structure emerged as important in determining which hypothesis appeared most capable of explaining each map's patterns. This was true even when regressions accounted for spurious effects of spatial autocorrelation. Each richness map, therefore, identified a different hypothesis as the most likely cause of broad‐scale gradients in species diversity. Main conclusions Because the ‘true’ spatial structure of species richness patterns remains elusive, firm conclusions regarding their underlying environmental drivers remain difficult. More broadly, our findings suggest that care should be taken to interpret putative determinants of large‐scale ecological gradients in light of the type and spatial characteristics of the underlying data. Indeed, closer scrutiny of these underlying data — here the distributions of individual species — and their environmental associations may offer important insights into the ultimate causes of observed broad‐scale patterns.     

The need for richness‐independent measures of turnover when delineating biogeographical regions

Delineating biogeographical regions is one of the primary steps when analysing biogeographical patterns. In their proposed quantitative framework, Kreft & Jetz (2010, Journal of Biogeography, 37 , 2029–2053) recommended the use of the β_sim index to delineate biogeographical regions because this turnover measure is weakly affected by differences in species richness between localities. A recent study by Carvalho et al. (2012, Global Ecology and Biogeography, 21 , 760–771) critiziced the use of β_sim in ecological and biogeographical studies, and proposed the β_‐3 index. Here we used simple numerical examples and an empirical case study (European freshwater fishes) to highlight potential pitfalls associated with the use of β_‐3 for bioregionalization. We show that β_‐3 is not a richness‐independent measure of species turnover. We also show that this index violates the ‘complementarity’ property, namely that localities without species in common have the largest dissimilarity, which is an essential prerequisite for beta diversity studies.     

The river domain: why are there more species halfway up the river?

Biologists have long noted higher levels of species diversity in the longitudinal middle‐courses of river systems and have proposed many explanations. As a new explanation for this widespread pattern, we suggest that many middle‐course peaks in richness may be, at least in part, a consequence of geometric constraints on the location of species’ ranges along river courses, considering river headwaters and mouths as boundaries for the taxa considered. We demonstrate this extension of the mid‐domain effect (MDE) to river systems for riparian plants along two rivers in Sweden, where a previous study found a middle‐course peak in richness of natural (non‐ruderal) species. We compare patterns of empirical richness of these species to null model predictions of species richness along the two river systems and to spatial patterns for six environmental variables (channel width, substrate fineness, substrate heterogeneity, ice scour, bank height, and bank area). In addition, we examine the independent prediction of mid‐domain effects models that species with large ranges, because the location of their ranges is more constrained, are more likely to produce a mid‐domain peak in richness than are species with small ranges. Species richness patterns of riparian plants were best predicted by models including both null model predictions and environmental variables. When species were divided into large‐ranged and small‐ranged groups, the mid‐domain effect was more prominent and the null model predictions were a better fit to the empirical richness patterns of large‐ranged species than those of small‐ranged species. Our results suggest that the peak in riparian plant species richness in the middle courses of the rivers studied can be explained by an underlying mid‐domain effect (driven by geometric constraints on large‐ranged species), together with environmental effects on richness patterns (particularly on small‐ranged species). We suggest that the mid‐domain effect may help to explain similar middle‐course richness peaks along other rivers.     

Richness patterns, species distributions and the principle of extreme deconstruction

Aim To analyse the global patterns in species richness of Viperidae snakes through the deconstruction of richness into sets of species according to their distribution models, range size, body size and phylogenetic structure, and to test if environmental drivers explaining the geographical ranges of species are similar to those explaining richness patterns, something we called the extreme deconstruction principle. Location Global. Methods We generated a global dataset of 228 terrestrial viperid snakes, which included geographical ranges (mapped at 1° resolution, for a grid with 7331 cells world‐wide), body sizes and phylogenetic relationships among species. We used logistic regression (generalized linear model; GLM) to model species geographical ranges with five environmental predictors. Sets of species richness were also generated for large and small‐bodied species, for basal and derived species and for four classes of geographical range sizes. Richness patterns were also modelled against the five environmental variables through standard ordinary least squares (OLS) multiple regressions. These subsets are replications to test if environmental factors driving species geographical ranges can be directly associated with those explaining richness patterns. Results Around 48% of the total variance in viperid richness was explained by the environmental model, but richness sets revealed different patterns across the world. The similarity between OLS coefficients and the primacy of variables across species geographical range GLMs was equal to 0.645 when analysing all viperid snakes. Thus, in general, when an environmental predictor it is important to model species geographical ranges, this predictor is also important when modelling richness, so that the extreme deconstruction principle holds. However, replicating this correlation using subsets of species within different categories in body size, range size and phylogenetic structure gave more variable results, with correlations between GLM and OLS coefficients varying from –0.46 up to 0.83. Despite this, there is a relatively high correspondence (r = 0.73) between the similarity of GLM‐OLS coefficients and R² values of richness models, indicating that when richness is well explained by the environment, the relative importance of environmental drivers is similar in the richness OLS and its corresponding set of GLMs. Main conclusions The deconstruction of species richness based on macroecological traits revealed that, at least for range size and phylogenetic level, the causes underlying patterns in viperid richness differ for the various sets of species. On the other hand, our analyses of extreme deconstruction using GLM for species geographical range support the idea that, if environmental drivers determine the geographical distribution of species by establishing niche boundaries, it is expected, at least in theory, that the overlap among ranges (i.e. richness) will reveal similar effects of these environmental drivers. Richness patterns may be indeed viewed as macroecological consequences of population‐level processes acting on species geographical ranges.     

The mid-domain effect in ectomycorrhizal fungi: range overlap along an elevation gradient on Mount Fuji,Japan

Mid-domain effect (MDE) models predict that the random placement of species'' ranges within a bounded geographical area leads to increased range overlap and species richness in the center of the bounded area. These models are frequently applied to study species-richness patterns of macroorganisms, but the MDE in relation to microorganisms is poorly understood. In this study, we examined the characteristics of the MDE in richness patterns of ectomycorrhizal (EM) fungi, an ecologically important group of soil symbionts. We conducted intensive soil sampling to investigate overlap among species ranges and the applicability of the MDE to EM fungi in four temperate forest stands along an elevation gradient on Mount Fuji, Japan. Molecular analyses using direct sequencing revealed 302 EM fungal species. Of 73 EM fungal species found in multiple stands, 72 inhabited a continuous range along the elevation gradient. The maximum overlap in species range and the highest species richness occurred at elevations in the middle of the gradient. The observed richness pattern also fit within the 95% confidence interval of the mid-domain null model, supporting the role of the MDE in EM fungal richness. Deviation in observed richness from the mean of the mid-domain null estimation was negatively correlated with some environmental factors, including precipitation and soil C/N, indicating that unexplained richness patterns could be driven by these environmental factors. Our results clearly support the existence of microbial species'' ranges along environmental gradients and the potential applicability of the MDE to better understand microbial diversity patterns.     

Biogeographical patterns of algal communities in the Mediterranean Sea: Cystoseira crinita‐dominated assemblages as a case study

Aim The aim of this study was to describe the composition, community structure and biogeographical variation of subtidal algal assemblages dominated by the brown alga Cystoseira crinita across the Mediterranean Sea. Location The Mediterranean coast, from Spain (1°25′ E) to Turkey (30°26′ E). Methods Data on the species composition and structure of assemblages dominated by the species C. crinita were collected from 101 sites in nine regions across the Mediterranean Sea. Multivariate and univariate statistical tools were used to investigate patterns of variation in the composition of the assemblages among sites and regions, and to compare these with previously defined biogeographical regions. Linear regressions of species richness versus longitude and versus latitude were also carried out to test previously formulated hypotheses of biodiversity gradients in the Mediterranean Sea. Results The main features characterizing C. crinita‐dominated assemblages across the Mediterranean included a similar total cover of species, a similar cover of C. crinita, and consistency in the presence of the epiphyte Haliptilon virgatum. Biogeographical variation was detected as shifts in relative abundances of species among regions, partly coinciding with previously described biogeographical sectors. A significant positive correlation was found between species richness and latitude, while no significant correlation was detected between species richness and longitude. Main conclusions The patterns of variation in community structure detected among the studied regions reflected their geographical positions quite well. However, latitude seemed to contribute more to the explanation of biological patterns of diversity than did geographical distances or boundaries, which classically have been used to delimit biogeographical sectors. Moreover, the positive correlation between species richness and latitude reinforced the idea that latitude, and possibly temperature as a related environmental factor, plays a primary role in structuring biogeographical patterns in the Mediterranean Sea. The lack of correlation between species richness and longitude contradicts the notion that there is a decrease in species richness from west to east in the Mediterranean, following the direction of species colonization from the Atlantic.     

Interpreting biogeographical boundaries among Afrotropical birds: spatial patterns in richness gradients and species replacement

Aim We analyse the geographical distribution of 1911 Afrotropical bird species using indices of three simple biogeographic patterns. The first index, the frequency of species with range edges (T_e), is formulated to map directly the density of species distribution limits, for comparison with the results of traditional biogeographical classification and ordination procedures, in order to show variations in the strength and breadth of transition zones. The other two indices are formulated to seek to distinguish as directly as possible between two components within these transition-zone patterns: contributions from gradients in species richness (T_g); and contributions from replacements among species (T_r). We test the ability of these indices to discover the same boundaries among Afrotropical bird faunas as one popular procedure for classifying areas (TWINSPAN) and then use them to look for geographical trends in the different kinds of transition zones. Location The analysis is restricted to the sub-Saharan or Afrotropical region, excluding the Arabian Peninsula, Madagascar and all offshore islands. Methods We record the presence of each species in 1961 1°×1° grid cells of the map. To apply the three indices, each (core) grid cell in turn is compared with its neighbouring eight cells in the grid. The range edges index (T_e) counts the number of species with range edges between the core cell and the surrounding cells. The richness gradients index (T_g) counts the largest difference in species richness measured diametrically across the core cell in any direction when there is a consistent trend in richness along this line of three cells. The species replacements index (T_r) counts the number of species pairs recorded within a nine-cell neighbourhood that are not corecorded within any of the cells. Values for each of the 1961 grid cells are calculated and used to produce colour-scale maps of transition zones. Results Large-scale spatial patterns of variation in density of range edges (T_e) are consistent with classifications of the same data and with most previous biogeographical classifications proposed for the region. Variation in richness gradients (T_g) and species replacements (T_r) explain different parts of this pattern, with transition zones around humid forests in the equatorial region being dominated by species replacement, and transition zones around deserts (most extensive in the north and south) being dominated by richness gradients. Main conclusions The three indices distinguish the spatial arrangement and intensity of different kinds of transition zones, thereby providing a first step towards a more rigorous mechanistic understanding of the different processes by which they may have arisen and are maintained. As an example of one such pattern shown by our analyses of Afrotropical birds, there is evidence for a broad latitudinal trend in the nature of transition zones in faunal composition (following the latitudinal distribution of the different kinds of habitat transitions), from being dominated by species replacements near the equator to being dominated by richness gradients further from the equator.     

Reconciling topographic and climatic effects on widespread and range-restricted species richness

Aim To identify the reasons behind differing geographical species richness patterns of range‐restricted and widespread species. Location The Western Hemisphere. Methods We used regression to determine the strongest environmental predictors of richness for widespread and range‐restricted mammal species in 10,000 km² quadrats in the continental Americas. We then used range‐placement models to predict the expected correlation between range‐restricted and widespread species richness were they to be determined by identical, random, or contrasting environmental factors. Finally, to determine the reasons underlying deviations from these predictions, we divided the Americas into 5% quantiles based on temperature and topographic heterogeneity and correlated richness of these two assemblages across quantiles – an approach that avoids constraints on statistical testing imposed by low potential for range overlap among range‐restricted species. Results Minimum annual temperature was the strongest predictor of widespread species richness while topographic heterogeneity was the best, although weak, predictor of range‐restricted species richness in conventional regression analysis. Our models revealed that the observed correlation between range‐restricted and widespread species richness was similar to what would be observed if both range‐restricted and widespread species richness were determined by temperature. Patterns of range‐restricted and widespread species richness were highly correlated across temperature quantiles, but range‐restricted species uniquely showed an increasing pattern across heterogeneity quantiles. Main conclusions Species richness gradients among range‐restricted species differ from those of widespread species, but not as extensively or for the reasons reported previously. Instead, these assemblages appear to share some but not all underlying environmental determinants of species richness. Our new approach to examining species richness patterns reveals that range‐restricted and widespread species richnesses share a common response to temperature that conventional analyses have not previously revealed. However, topographic heterogeneity has assemblage‐specific effects on range‐restricted species.     

The coincidence of rarity and richness and the potential signature of history in centres of endemism

We investigate the relative importance of stochastic and environmental/topographic effects on the occurrence of avian centres of endemism, evaluating their potential historical importance for broad‐scale patterns in species richness across Sub‐Saharan Africa. Because species‐rich areas are more likely to be centres of endemism by chance alone, we test two null models: Model 1 calculates expected patterns of endemism using a random draw from the occurrence records of the continental assemblage, whereas Model 2 additionally implements the potential role of geometric constraints. Since Model 1 yields better quantitative predictions we use it to identify centres of endemism controlled for richness. Altitudinal range and low seasonality emerge as core environmental predictors for these areas, which contain unusually high species richness compared to other parts of sub‐Saharan Africa, even when controlled for environmental differences. This result supports the idea that centres of endemism may represent areas of special evolutionary history, probably as centres of diversification.     

Constraining null models with environmental gradients: a new method for evaluating the effects of environmental factors and geometric constraints on geographic diversity patterns

The relative effects of climate and geometric constraints on geographic diversity patterns have long been controversial. We developed a new method to assess the role of geometric constraints in shaping altitudinal richness patterns. We showed how plant species richness on four mountains in southwest China are shaped by geometric constraints and environmental gradients together. Contrary to previous studies, our results suggested that: 1) small‐ and large‐ranged species richness were largely controlled by the same environmental gradients, and differed mainly in the effect of geometric constraints. 2) The contribution of geometric constraints (in addition to environmental gradients) to explaining species richness was greater when species richness peaked at low altitudes than at mid‐altitudes, suggesting that geometric constraints may be very important when richness peaks are far away from mid‐domains. 3) Relating species richness directly to environmental factors (the most widely used method in biodiversity studies) may be misleading when geometric constraints may be affecting the richness pattern, because this method may overestimate the effect of environmental factors by failing to distinguish the confounding effect of geometric constraints. Instead, the effect of environmental factors can be evaluated with an underlying gradient derived from small‐ranged species. 4) The geometric constraints effect cannot be fully evaluated by pure geometric constraints models, and is better evaluated with range‐based models constrained with environmental gradients. 5) If the generality of our findings is supported for other taxa on other gradients, then many previous studies on the effects of climate and of geometric constraints on geographic diversity patterns may need to be re‐visited.     

Spatial diversity patterns of Pristimantis frogs in the Tropical Andes

Although biodiversity gradients have been widely documented, the factors governing broad‐scale patterns in species richness are still a source of intense debate and interest in ecology, evolution, and conservation biology. Here, we tested whether spatial hypotheses (species–area effect, topographic heterogeneity, mid‐domain null model, and latitudinal effect) explain the pattern of diversity observed along the altitudinal gradient of Andean rain frogs of the genus Pristimantis. We compiled a gamma‐diversity database of 378 species of Pristimantis from the tropical Andes, specifically from Colombia to Bolivia, using records collected above 500 m.a.s.l. Analyses were performed at three spatial levels: Tropical Andes as a whole, split in its two main domains (Northern and Central Andes), and split in its 11 main mountain ranges. Species richness, area, and topographic heterogeneity were calculated for each 500‐m‐width elevational band. Spatial hypotheses were tested using linear regression models. We examined the fit of the observed diversity to the mid‐domain hypothesis using randomizations. The species richness of Pristimantis showed a hump‐shaped pattern across most of the altitudinal gradients of the Tropical Andes. There was high variability in the relationship between area and species richness along the Tropical Andes. Correcting for area effects had little impact in the shape of the empirical pattern of biodiversity curves. Mid‐domain models produced similar gradients in species richness relative to empirical gradients, but the fit varied among mountain ranges. The effect of topographic heterogeneity on species richness varied among mountain ranges. There was a significant negative relationship between latitude and species richness. Our findings suggest that spatial processes partially explain the richness patterns of Pristimantis frogs along the Tropical Andes. Explaining the current patterns of biodiversity in this hot spot may require further studies on other possible underlying mechanisms (e.g., historical, biotic, or climatic hypotheses) to elucidate the factors that limit the ranges of species along this elevational gradient.     

Disentangling the drivers of continental mammalian endemism

Endemic species and species with small ranges are ecologically and evolutionarily distinct and are vulnerable to extinction. Determining which abiotic and biotic factors structure patterns of endemism on continents can advance our understanding of global biogeographic processes, but spatial patterns of mammalian endemism have not yet been effectively predicted and reconstructed. Using novel null model techniques, we reconstruct trends in mammalian endemism and describe the isolated and combined effects of physiographic, ecological, and evolutionary factors on endemism. We calculated weighted endemism for global continental ecoregions and compared the spatial distribution of endemism to niche-based, geographic null models of endemism. These null models distribute species randomly across continents, simulating their range sizes from their degree of climatic specialization. They isolate the effects of physiography (topography and climate) and species richness on endemism. We then ran linear and structural models to determine how topography and historical climate stability influence endemism. The highest rates of mammalian endemism were found in topographically rough, climatically stable ecoregions with many species. The null model that isolated physiography did not closely approximate the observed distribution of endemism (r² = .09), whereas the null model that incorporated both physiography and species richness did (r² = .59). The linear models demonstrate that topography and climatic stability both influenced endemism values, but that average climatic niche breadth was not highly correlated with endemism. Climate stability and topography both influence weighted endemism in mammals, but the spatial distribution of mammalian endemism is driven by a combination of physiography and species richness. Despite its relationship to individual range size, average climate niche breadth has only a weak influence on endemism. The results highlight the importance of historical biogeographic processes (e.g. centers of speciation) and geography in driving endemism patterns, and disentangle the mechanisms structuring species ranges worldwide.