Species-range size distributions in Britain

The detailed forms of species-range size distributions in Britain are determined and contrasted for ten taxonomic assemblages (liverworts, vascular plants, molluscs [aquatic and terrestrial], dragonflies, macro-moths. butterflies, birds [breeding and wintering], mammals). All are strongly right-skewed when range sizes are untransformed. A logarithmic transformation fails to normalise the distribution for all but one group, and the distributions for several groups are not readily normalised at all. Taxa with larger median range sizes have species-range size distributions that are less strongly right-skewed. The median observed range sizes of species in each of the taxonomic groups fall, in terms of decreasing range size, in the sequence wintering birds < breeding birds < mammals < butterflies < terrestrial molluscs < dragonflies < aquatic molluscs < vascular plants < moths < liverworts. Despite the difficulties in deriving a simple and sensible mechanistic model for range size distributions, this is likely to be the most important next step towards understanding their forms.     

Species-range size distributions: products of speciation,extinction and transformation

One basic summary of the spatial pattern of biodiversity across the surface of the Earth is provided by a species-range size distribution, the frequency distribution of the numbers of species exhibiting geographic ranges of different sizes. Although widely considered to be approximately lognormal, increasingly it appears that across a variety of groups of organisms this distribution systematically departs from such a form. Whatever its detailed shape, however, the distribution must arise as a product of three processes, speciation, extinction and transformation (the temporal dynamics of the range sizes of species during their life times). Considering the role potentially played by each of these processes necessitates drawing on information from a diverse array of research fields, and highlights the possible role of geographic range size as a common currency uniting them.     

On the heritability of geographic range sizes

Within taxonomic groups, most species are restricted in their geographic range sizes, with only a few being widespread. The possibility that species-level selection on range sizes contributes to the characteristic form of such species-range size distributions has previously been raised. This would require that closely related species have similar range sizes, an indication of "heritability" of range sizes at the species level. Support for this view came from a positive correlation between the range sizes of closely related pairs of fossil mollusc species. We extend this analysis by considering the relationship between the geographic range sizes of 103 pairs of contemporary avian sister species. Range sizes in these sister species show no evidence of being more similar to each other than expected by chance. A reassessment of the mollusc data also suggests that the high correlation was probably overestimated because of the skewed nature of range size data. The fact that sister species tend to have similar life histories and ecologies suggests that any relationship between range sizes and biology is likely to be complicated and will be influenced by historical factors, such as mode of speciation and postspeciation range size transformations.     

Geographical range size heritability: what do neutral models with different modes of speciation predict?

Aim Phylogenetic conservatism or heritability of the geographical range sizes of species (i.e. the tendency for closely related species to share similar range sizes) has been predicted to occur because of the strong phylogenetic conservatism of niche traits. However, the extent of such heritability in range size is disputed and the role of biology in shaping this attribute remains unclear. Here, we investigate the level of heritability of geographical range sizes that is generated from neutral models assuming no biological differences between species. Methods We used three different neutral models, which differ in their speciation mode, to simulate the life‐history of 250,000 individuals in a square lattice of 50 × 50 cells. These individuals can speciate, reproduce, migrate and die in the metacommunity according to stochastic events. We ran each model for 3000 steps and recorded the range size of each species at each step. The heritability of geographical range size was assessed using an asymmetry coefficient between range sizes of sister species and using the coefficient of correlation between the range sizes of ancestors and their descendants. Results Our results demonstrated the ability of neutral models to mimic some important observed patterns in the heritability of geographical range size. Consistently, sister species exhibited higher asymmetry in range sizes than expected by chance, and correlations between the range sizes of ancestor–descendant species pairs, although often weak, were almost invariably positive. Main conclusions Our findings suggest that, even without any biological trait differences, statistically significant heritability in the geographical range sizes of species can be found. This heritability is weaker than that observed in some empirical studies, but suggests that even here a substantial component of heritability may not necessarily be associated with niche conservatism. We also conclude that both present‐day and fossil data sets may provide similar information on the heritability of the geographical range sizes of species, while the omission of rare species will tend to overestimate this heritability.     

Age, area and avian diversification

Using coarse resolution data on the spatial distribution of the entire New World avifauna, we test for phylogenclic patterns in the mean and total geographic range sizes of taxa. The analyses reveal that (i) the species-range size distribution is only approximately normalized, and remains significantly left-skewed, under logarithmic transformation. Most variance in range sizes is explained at the level of species within genera; (ii) there is no effect of the age of taxa on mean clade range size, although older taxa are more likely to have larger total range sizes; (iii) there is some evidence that taxa comprising more species have larger total range sizes; (iv) there is little or no evidence for a relationship between rate of cladogenesis and range size. The results suggest that geographic range size is a labile trait, at least for New World birds, and that the influence of evolutionary history is only weakly detectable in the range size variation of extant taxa, at least at the scale of analysis used here. In addition to these conclusions, two general and important procedural issues emerge.     

Speciation and extinction drive the appearance of directional range size evolution in phylogenies and the fossil record

While the geographic range of a species is a fundamental unit of macroecology and a leading predictor of extinction risk, the evolutionary dynamics of species' ranges remain poorly understood. Based on statistical associations between range size and species age, many studies have claimed support for general models of range evolution in which the area occupied by a species varies predictably over the course of its life. Such claims have been made using both paleontological data and molecular estimates of the age of extant species. However, using a stochastic model, we show that the appearance of trends in range size with species' age can arise even when range sizes have evolved at random through time. This occurs because the samples of species used in existing studies are likely to be biased with respect to range size: for example, only those species that happened to have large or expanding ranges are likely to survive to the present, while extinct species will tend to be those whose ranges, by chance, declined through time. We compared the relationship between the age and range size of species arising under our stochastic model to those observed across 1,269 species of extant birds and mammals and 140 species of extinct Cenozoic marine mollusks. We find that the stochastic model is able to generate the full spectrum of empirical age-area relationships, implying that such trends cannot be simply interpreted as evidence for models of directional range size evolution. Our results therefore challenge the theory that species undergo predictable phases of geographic expansion and contraction through time.     

Geographical range size and host specificity in ectoparasites: a case study with Amphipsylla fleas and rodent hosts

Aim  To identify the factors that determine the geographical range sizes of ectoparasites with different degrees of host specificity.
Location  The study used data on the distributions of fleas of the genus Amphipsylla and their rodent hosts across the Holarctic.
Methods  All known points of occurrence of 32 flea species and 51 species of their rodent hosts were mapped. The shape and size of the geographical range of each species were estimated using a combination of the minimal convex polygon technique and modelling with the garp algorithm. Factors determining the geographical range sizes of the fleas were identified using stepwise multiple regression analysis.
Results  The geographical range size of fleas that are strongly host-specific across their entire ranges correlated positively with the geographical range size of the fleas' principal hosts, and negatively with the geographical range size of the fleas' potential competitors. The geographical range sizes of both (1) fleas that are locally host-specific but that shift their host preferences geographically, and (2) host-opportunistic fleas were positively correlated only with the area of the geographical ranges of their principal hosts. Strongly host-specific fleas occupied 0.2–80.0% of the geographical range of their principal hosts, whereas this figure was 0.9–83.7% in locally host-specific fleas and 16.6–63.7% in host-opportunistic fleas.
Main conclusions  The main determinant of the geographical range size of a flea species is the size of the geographical range of its hosts. The role of potential competitors in determining the geographical range size is stronger in host-specific than in host-opportunistic fleas. Cases in which the geographical range of a parasite is smaller than the geographical range(s) of its host(s) owing to narrower parasite environmental tolerances are much more frequent in host-opportunistic than in host-specific fleas.     

demoniche – an R‐package for simulating spatially‐explicit population dynamics

demoniche is a freely available R‐package which simulates stochastic population dynamics in multiple populations of a species. A demographic model projects population sizes utilizing several transition matrices that can represent impacts on species growth. The demoniche model offers options for setting demographic stochasticity, carrying capacity, and dispersal. The demographic projection in each population is linked to spatially‐explicit niche values, which affect the species growth. With the demoniche package it is possible to compare the influence of scenarios of environmental changes on future population sizes, extinction probabilities, and range shifts of species.     

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.     

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.     

Beyond Rapoport's rule: evaluating range size patterns of New World birds in a two-dimensional framework

Aim To evaluate Rapoport's rule for New World birds in two‐dimensional geographical space. We specifically test for a topography × climate interaction that predicts little difference in range sizes between lowlands and mountains in cold climates, whereas in the tropics, montane species have narrow ranges and lowland species have broad ranges. Location The western hemisphere. Methods We used digitized range maps of breeding birds to generate mean range sizes in grids of 27.5 × 27.5 km and 110 × 110 km across North and South America. We examined the geographical pattern with respect to range in elevation, mean temperature in the coldest month, their interaction, biome size and continental width, using model II analysis of variance, multiple regression and simple correlation. Results In northern latitudes species have broad ranges in both mountainous and flat areas. However, range sizes in the mountains and lowlands diverge southwards, with the most extreme differences in the tropics. Further, there are minimal differences in range sizes across latitudes in lowlands. The smallest mean ranges occur in the tropical Andes. Mean range sizes in north‐central Canada, Central America and Argentina/Chile are also small, reflecting the narrowing of the continents in these areas. The best regression model explained 51% of the variation in mean range size. Main conclusions The two‐dimensional range size pattern indicates that neither winter temperature nor annual variability in temperature strongly influences the distribution of range sizes directly; rather, climate influences bird range sizes indirectly via effects on habitat size. Also, macroclimate interacts with topographic relief across latitudes, generating sharp mesoscale habitat gradients in tropical mountains but not in high latitude mountains or in lowlands at any latitude. Birds respond to these habitat gradients, resulting in ‘latitudinal’ range size gradients in topographically complex landscapes but not in simple landscapes.     

Assembly of local communities: consequences of an optimal body size for the organization of competitively structured communities

Much empirical evidence suggests that there is an optimal body size for mammals and that this optimum is in the vicinity of l00g. This presumably reflects an underlying fitness function that is greatest at this mass. Here, I combine such a fitness function with an equilibrium model of competitive character displacement to assess the potential influence of a globally optimal body size in structuring local ecological communities. The model accurately predicts the range of body sizes and the average difference in size for species in communities of varying species richness. The model also predicts a uniform spacing of body sizes, rather than the gaps and clumps in the sizes of coexisting species observed in real communities. Alternative explanations for this phenomenon are discussed. The allometric relationships that result in a body size optimum subsume a large number of characteristics associated with the physiological, behavioral, demographic, and evolutionary dynamics of the species. Further integration of the underlying dynamics (e.g. individual energetics) of these relationships into all hierarchical levels of ecology will have to incorporate multiple interactive sites, spatial heterogeneity, and phylogenetic structure, but it has the potential to provide important discoveries into the means by which natural selection operates.     

Wing morphology and migration status,but not body size,habitat or Rapoport's rule predict range size in North‐American dragonflies (Odonata: Libellulidae)

Understanding why species range sizes vary is important for predicting the impact of environmental change on biodiversity. Here we use a multi‐variable approach in a phylogenetic comparative context to understand how four morphological, two ecological, and two eco‐geographical variables are associated with range size, latitudinal range and longitudinal range in 81 species of North‐American libellulid dragonflies. Our results show that: 1) migratory species and species with a more expanded basal hindwing lobe have a larger range size; 2) opposite to Rapoport's rule, latitudinal range is negatively correlated with mid‐range latitude; 3) longitudinal range is predicted by wing morphology and migration; 4) body size and larval habitat are not correlated with range size, latitudinal range or longitudinal range. These results suggest that dispersal‐related traits, such as wing shape and migratory status, are important factors in predicting the range size of libellulid dragonflies. In addition, the reverse Rapoport's rule suggests that more northern‐centred species might be more specialized than more southern‐centred species. We suggest that the variables predicting range size are likely imposed by taxon‐specific morphological, ecological, physiological and behavioural traits. Taxon‐specific knowledge is thus necessary to understand the dynamics of range sizes and is important to implement successful restoration and conservation plans of threatened species.     

Relationships between geographical range size, body size, local abundance, and habitat breadth in North American suckers and sunfishes

Aim I examine the relationship between geographical range size and three variables (body size, an index of habitat breadth, and an index of local abundance) within a phylogenetic framework in North American species of suckers and sunfishes. Location North America Methods Regressions after independent contrasts of geographical range size, body size, habitat breadth, and local abundance. Results Species with large range sizes tend to be larger-bodied, be more locally abundant, and have higher habitat breadths. Character reconstructions support the prediction that variables associated with rarity (small geographical range size, low local abundance, low niche breadth, and large body size) evolve in unison, although large body size was associated with the opposite traits in these taxa. Gaston & Blackburn (1996a) suggested using visual identification of the lower boundary of the geographical range-body size relationship to identify extinction-prone species; this resulted in thirteen species that are potentially extinction-prone. Main conclusions Similar evolutionary mechanisms appear to operate on body size and other variables related to rarity, even in distantly related taxa.     

The shape and temporal dynamics of phylogenetic trees arising from geographic speciation

Phylogenetic trees often depart from the expectations of stochastic models, exhibiting imbalance in diversification among lineages and slowdowns in the rate of lineage accumulation through time. Such departures have led to a widespread perception that ecological differences among species or adaptation and subsequent niche filling are required to explain patterns of diversification. However, a key element missing from models of diversification is the geographical context of speciation and extinction. In this study, we develop a spatially explicit model of geographic range evolution and cladogenesis, where speciation arises via vicariance or peripatry, and explore the effects of these processes on patterns of diversification. We compare the results with those observed in 41 reconstructed avian trees. Our model shows that nonconstant rates of speciation and extinction are emergent properties of the apportioning of geographic ranges that accompanies speciation. The dynamics of diversification exhibit wide variation, depending on the mode of speciation, tendency for range expansion, and rate of range evolution. By varying these parameters, the model is able to capture many, but not all, of the features exhibited by birth-death trees and extant bird clades. Under scenarios with relatively stable geographic ranges, strong slowdowns in diversification rates are produced, with faster rates of range dynamics leading to constant or accelerating rates of apparent diversification. A peripatric model of speciation with stable ranges also generates highly unbalanced trees typical of bird phylogenies but fails to produce realistic range size distributions among the extant species. Results most similar to those of a birth-death process are reached under a peripatric speciation scenario with highly volatile range dynamics. Taken together, our results demonstrate that considering the geographical context of speciation and extinction provides a more conservative null model of diversification and offers a very different perspective on the phylogenetic patterns expected in the absence of ecology.     

样本容量和物种特征对BIOCLIM模型模拟物种分布准确度的影响——以12个中国特有落叶栎树种为例

物种分布模型被广泛应用于生态学、生物地理学及保护生物学等领域的研究。由于难于取样或标本记录不完善等原因, 真正能够用于模型预测的物种分布数据非常有限。因此, 有必要搞清楚样本容量和物种特征对模型模拟准确度的影响, 为确定以物种特征为区分条件的最小样本容量奠定基础。为了探讨应用BIOCLIM模型预测中国特有植物种的效果, 以12个落叶栎树种为例, 从不同的样本容量和生态特征两方面研究其对BIOCLIM模型模拟准确度的影响。结果表明: BIOCLIM模型模拟准确度随样本容量的增加在初期几乎呈直线增加趋势至样本容量达到25, 随后渐变平缓至样本容量为75~100时达到最大值。此外, 生态幅窄和环境特化物种比生态幅宽和对环境耐受性强的物种更容易获得较高的准确度。结果说明, BIOCLIM可有效地用于样本数量较小的狭域型物种分布预测。     

The influence of spatial resolution on macroecological patterns of range size variation: a case study using parrots (Aves: Psittaciformes) of the world

Aim To assess the extent to which the resolution at which geographical range sizes are measured influences macroecological patterns in this variable. Location Global. Methods Data on the geographical ranges of parrot species were digitized, and a Geographic Information System used to produce nine range size estimates for each species using different degrees of spatial resolution. The inter‐correlation of these estimates was then compared, together with their patterns of covariation with population size, body mass and migratory behaviour (across species and controlling for phylogeny), their pattern of phylogenetic correlation, and the frequency distributions of the different measures. Results Strong correlations exist among all nine range size measures across species, albeit that measures of similar spatial resolution are more strongly correlated. All measures show similar patterns of covariation with population size, body mass and migratory behaviour, and similar patterns of phylogenetic correlation. The skewness of frequency distributions increases towards zero as the resolution of the range size measure declines. Main conclusions The results of macroecological analyses are little affected by the resolution with which geographical range sizes are calculated, at least for the parrots of the world. Previously published studies based on crude measures of range size would be unlikely to have produced markedly different conclusions had they used more refined range size metrics.     

Incorporating the temporal autocorrelation of demographic rates into structured population models

Population dynamics are typically temporally autocorrelated: population sizes are positively or negatively correlated with past population sizes. Previous studies have found that positive temporal autocorrelation increases the risk of extinction due to ‘inertia’ that prolongs downward fluctuations in population size. However, temporal autocorrelation has not yet been analyzed at the level of life cycle transitions. We developed an R package, colorednoise, which creates stochastic matrix population projections with distinct temporal autocorrelation values for each matrix element. We used it to analyze long-term demographic data on 25 populations from the COMADRE and COMPADRE databases and simulate their stochastic dynamics. We found a broad range of temporal autocorrelation across species, populations and life cycle stages. The number of stage-classes in the matrix strongly affected the temporal autocorrelation of the growth rate. In the plant populations, reproduction transitions had more negative temporal autocorrelation than survival transitions, and matrices dominated by positive temporal autocorrelation had higher extinction risk, while in animal populations transition type was not associated with noise color. Our results indicate that temporal autocorrelation varies across life cycle transitions, even among populations of the same species. We present the colorednoise package for researchers to analyze the temporal autocorrelation of structured demographic rates.     

Halfway up the trophic chain: development of parasite communities in the sparid fish Boops boops

We examined the patterns of composition and structure of parasite communities in the Mediterranean sparid fish Boops boops along a gradient of fish sizes, using a large sample from a single population. We tested the hypothesis that species forming the core of the bogue parasite fauna (i.e. species which have a wide geographical range and are responsible for recognizable community structure) appear early in the fish ontogeny. The sequential community development observed supported the prediction that core species appear in the fish population earlier than rare and stochastic species. There was also a strong correlation between the order of 'arrival' of the species and their overall prevalence. Six key species were responsible for recognizable community structure across size/age cohorts; the addition to this baseline community of key parasite species resulted in a nested structure that is linked to differential species abundance rather than fish size. Information on the life-cycles, distribution and host range of the parasites is used to explain the observed patterns of parasite community structure. We conclude that the small mouth size of B. boops coupled with suction feeding may provide a setting for passive sampling as a mechanism leading to non-random parasite community structure.     

Effects of deterministic and random refuge in a prey-predator model with parasite infection

Most natural ecosystem populations suffer from various infectious diseases and the resulting host-pathogen dynamics is dependent on host's characteristics. On the other hand, empirical evidences show that for most host pathogen systems, a part of the host population always forms a refuge. To study the role of refuge on the host-pathogen interaction, we study a predator-prey-pathogen model where the susceptible and the infected prey can undergo refugia of constant size to evade predator attack. The stability aspects of the model system is investigated from a local and global perspective. The study reveals that the refuge sizes for the susceptible and the infected prey are the key parameters that control possible predator extinction as well as species co-existence. Next we perform a global study of the model system using Lyapunov functions and show the existence of a global attractor. Finally we perform a stochastic extension of the basic model to study the phenomenon of random refuge arising from various intrinsic, habitat-related and environmental factors. The stochastic model is analyzed for exponential mean square stability. Numerical study of the stochastic model shows that increasing the refuge rates has a stabilizing effect on the stochastic dynamics.