Tracing the patterns of non-marine turtle richness from the Triassic to the Palaeogene: from origin to global spread

Turtles are key components of modern vertebrate faunas and their diversity and distributions are likely to be affected by anthropogenic climate change. However, there is limited baseline data on turtle taxonomic richness through time or assessment of their past responses to global environmental change. We used the extensive Triassic–Palaeogene (252–223 Ma) fossil record of terrestrial and freshwater turtles to investigate diversity patterns, finding substantial variation in richness through time and between continents. Globally, turtle richness was low from their Triassic origin until the Late Jurassic. There is strong evidence for high richness in the earliest Cretaceous of Europe, becoming especially high following the Cretaceous Thermal Maximum and declining in all continents by the end-Cretaceous. At the K–Pg boundary, South American richness levels changed little while North American richness increased, becoming very high during the earliest Palaeogene (Danian). Informative data are lacking elsewhere for this time period. However, the Selandian–Thanetian interval, approximately 5 myr after the K–Pg mass extinction, shows low turtle richness in Asia, Europe and South America, suggesting that the occurrence of exceptional turtle richness in the post-extinction Paleocene fauna of North America is not globally representative. Richness decreased over the Eocene–Oligocene boundary in North America but increased to its greatest known level for Europe, implying very different responses to dramatic climatic shifts. Time series regressions suggest number of formations sampled and palaeotemperature are the primary influencers of face-value richness counts, but additional factors not tested here may also be involved.     

Influence of light and soil moisture on Sierran mixed-conifer understory communities

Sierra Nevada forests have high understory species richness yet we do not know which site factors influence herb and shrub distribution or abundance. We examined the understory of an old-growth mixed-conifer Sierran forest and its distribution in relation to microsite conditions. The forest has high species richness (98 species sampled), most of which are herbs with sparse cover and relatively equal abundance. Shrub cover is highly concentrated in discrete patches. Using overstory tree cover and microsite environmental conditions, four habitats were identified; tree cluster, partial canopy, gap, and rock/shallow soil. Herb and shrub species were strongly linked with habitats. Soil moisture, litter depth and diffuse light were the most significant environmental gradients influencing understory plant distribution. Herb cover was most strongly influenced by soil moisture. Shrub cover is associated with more diffuse light, less direct light, and sites with lower soil moisture. Herb richness is most affected by conditions which influence soil moisture. Richness is positively correlated with litter depth, and negatively correlated with direct light and shrub cover. Disturbance or management practices which change forest floor conditions, shallow soil moisture and direct light are likely to have the strongest effect on Sierran understory abundance and richness.     

Assessment of the relationships of geographic variation in species richness to climate and landscape variables within and among lineages of North American freshwater fishes

Aim Geographic variation in species richness is a well‐studied phenomenon. However, the unique response of individual lineages to environmental gradients in the context of general patterns of biodiversity across broad spatial scales has received limited attention. The focus of this research is to examine relationships between species richness and climate, topographic heterogeneity and stream channel characteristics within and among families of North American freshwater fishes. Location The United States and Canada. Methods Distribution maps of 828 native species of freshwater fishes were used to generate species richness estimates across the United States and Canada. Variation in species richness was predicted using spatially explicit models incorporating variation in climate, topography and/or stream channel length and stream channel diversity for all 828 species as well as for the seven largest families of freshwater fishes. Results The overall gradient of species richness in North American freshwater fishes is best predicted by a model incorporating variables describing climate and topography. However, the response of species richness to particular climate or landscape variables differed among families, with models possessing the highest predictive ability incorporating data on climate, topography and/or stream channel characteristics within a region. Main conclusions The correlations between species richness and abiotic variables suggest a strong influence of climate and physical habitat on the structuring of regional assemblages of North American freshwater fishes. However, the relationship between these variables and species richness varies among families, suggesting the importance of phylogenetic constraints on the regulation of geographic distributions of species.     

Latitudinal variation in the shape of the species body size distribution: an analysis using freshwater fishes

Many taxonomic and ecological assemblages of species exhibit a right-skewed body size-frequency distribution when characterized at a regional scale. Although this distribution has been frequently described, factors influencing geographic variation in the distribution are not well understood, nor are mechanisms responsible for distribution shape. In this study, variation in the species body size-frequency distributions of 344 regional communities of North American freshwater fishes is examined in relation to latitude, species richness, and taxonomic composition. Although the distribution of all species of North American fishes is right-skewed, a negative correlation exists between latitude and regional community size distribution skewness, with size distributions becoming left-skewed at high latitudes. This relationship is not an artifact of the confounding relationship between latitude and species richness in North American fishes. The negative correlation between latitude and regional community size distribution skewness is partially due to the geographic distribution of families of fishes and apparently enhanced by a nonrandom geographic distribution of species within families. These results are discussed in the context of previous explanations of factors responsible for the generation of species size-frequency distributions related to the fractal nature of the environment, energetics, and evolutionary patterns of body size in North American fishes.     

One,two and three‐dimensional geometric constraints and climatic correlates of North American tree species richness

The ‘mid‐domain effect’ (MDE) has received much attention recently as a candidate explanation for patterns in species richness over large geographic areas. Mid‐domain models generate a central peak in richness when species ranges are randomly placed within a bounded geographic area (i.e. the domain). The most common terrestrial mid‐domain models published to date have been 1‐D latitude or elevation models and 2‐D latitude‐longitude models. Here, we test 1‐D, 2‐D and 3‐D mid‐domain models incorporating latitude, longitude and elevation, and assess independent and concurrent effects of geometric constraints and climatic variables on species richness of North American trees. We use both the traditional ‘global’ regression models as well as geographically weighted regressions (‘local’ models) to examine local variation in the contribution of MDE and climatic variables to species richness across the domain. Our results show that in some dimensions the contribution of MDE to patterns of species richness can be quite substantial, and we show that in most cases a combination of MDE and climate predicted empirical species richness best in both local and global models. For the North American domain, MDE in the elevation dimension is clearly important in describing patterns of empirical species richness. We also show that the assumption of stationarity in global models is not met in the North American domain and that results of these models mask complex patterns in both the effect of MDE on richness and the response of species richness to climate. In particular we show the increased explanatory role of MDE in predicting species richness as domain edges are approached. Our results support the hypothesis that geometric constraints contribute to species richness patterns and we suggest the mid‐domain effect should be considered alongside more traditional environmental correlates in understanding patterns of species diversity.     

An evolutionary tolerance model explaining spatial patterns in species richness under environmental gradients and geometric constraints

In this study, we developed a simulation model based on the ecological and evolutionary dynamics of geographical ranges, to understand the role of species' environmental tolerances and the strength of the environmental gradient in determining spatial patterns in species richness. Using an one-dimensional space, we present the model and dissect its parameters. Also, we test the ability of the model to simulate richness in complex two-dimensional domains and to fit real patterns in species richness, using South American Tyrannidae as an example. We found that a mid-spatial peak in species richness arises spontaneously under conditions of high environmental tolerances and/or a weak environmental gradient, since this condition causes wide species' geographic ranges, which are constrained by domain's boundary and tend to overlap in the middle. Our model was also a good predictor of real patterns in species richness, especially under conditions of high environmental strength and small species' tolerance. We conclude that this kind of spatial simulation models based on species' physiological tolerance may be an important tool to understand the evolutionary dynamics of species' geographic ranges and in spatial patterns of species richness.     

Can stochastic geographical evolution re‐create macroecological richness–environment correlations?

Aim Richness gradients are frequently correlated with environmental characteristics at broad geographic scales. In particular, richness is often associated with energy and climate, while environmental heterogeneity is rarely its best correlate. These correlations have been interpreted as evidence in favour of environmental determinants of diversity gradients, particularly energy and climate. This interpretation assumes that the expected‐by‐random correlation between richness and environment is zero, and that this is equally true for all environmental characteristics. However, these expectations might be unrealistic. We investigated to what degree basic evolutionary/biogeographical processes occurring independently of environment could lead to richness gradients that correlate with environmental characteristics by chance alone. Location Africa, Australia, Eurasia and the New World. Methods We produced artificial richness gradients based on a stochastic simulation model of geographic diversification of clades. In these simulations, species speciate, go extinct and expand or shift their distributions independently of any environmental characteristic. One thousand two hundred repetitions of this model were run, and the resulting stochastic richness gradients were regressed against real‐world environmental variables. Stochastic species–environment relationships were then compared among continents and among three environmental characteristics: energy, environmental heterogeneity and climate seasonality. Results Simulations suggested that a significant degree of correlation between richness gradients and environment is expected even when clades diversify and species distribute stochastically. These correlations vary considerably in strength; but in the best cases, environment can spuriously account for almost 80% of variation in stochastic richness. Additionally, expected‐by‐chance relationships were different among continents and environmental characteristics, producing stronger spurious relationships with energy and climate than with heterogeneity. Main conclusions We conclude that some features of empirical species–environment relationships can be reproduced just by chance when taking into account evolutionary/biogeographical processes underlying the construction of species richness gradients. Future tests of environmental effects on richness should consider structure in richness–environment correlations that can be produced by simple evolutionary null models. Research should move away from the naive non‐biological null hypotheses that are implicit in traditional statistical tests.     

Habitat occurrence and prey distribution of a multi-species community of shrews in the Siberian taiga

The habitat occurrence and invertebrate prey distribution of nine species of shrew in the mid-taiga of central Siberia were investigated. Species richness ranged from 4–9 shrews per habitat. Sorex araneus and S. caecutiens were numerically dominant in all seven habitats (44 and 36% of the total catch, respectively) while Sorex minutus, S. tundrensis, S. isodon , and S. roboratus each constituted 4–6% and Sorex minutissimus, S. daphaenodon , and Neomys fodiens were rare (< 1% each). There was no overall correlation between abundance of shrews and invertebrate prey, but flood-plain habitats supported the greatest abundance and species richness of shrews, and high density and biomass of prey. Oligochaete-eating shrews were twice as numerous here as in other habitats, coincident with high abundance of oligochaetes. The large, earthworm-feeding Sorex roboratus occurred only here. The more acid, typical taiga habitats had lower adundance and species richness of shrews. They had the lowest density and biomass of prey, particularly oligochaetes, and far fewer oligochaete-eating shrews. The relative paucity of shrews in bush-meadow habitats, despite abundant prey, implied that habitat structure influences shrew distribution. Differential numbers of certain species in the presence or absence of larger congeners also suggested that interspecific competitive effects influence habitat selection by shrews. The high species richness of shrews here in the mid-taiga may be accounted for by the heterogeneous nature of the constituent habitats which provide niches for small and large species of shrew with a range of feeding habits.     

Bird diversity patterns in Neotropical temperate farmlands: The role of environmental factors and trophic groups in the spring and autumn

Productivity, habitat heterogeneity and environmental similarity are of the most widely accepted hypotheses to explain spatial patterns of species richness and species composition similarity. Environmental factors may exhibit seasonal changes affecting species distributions. We explored possible changes in spatial patterns of bird species richness and species composition similarity. Feeding habits are likely to have a major influence in bird–environment associations and, given that food availability shows seasonal changes in temperate climates, we expect those associations to differ by trophic group (insectivores or granivores). We surveyed birds and estimated environmental variables along line‐transects covering an E‐W gradient of annual precipitation in the Pampas of Argentina during the autumn and the spring. We examined responses of bird species richness to spatial changes in habitat productivity and heterogeneity using regression analyses, and explored potential differences between seasons of those responses. Furthermore, we used Mantel tests to examine the relationship between species composition similarity and both the environmental similarity between sites and the geographic distance between sites, also assessing differences between seasons in those relationships. Richness of insectivorous birds was directly related to primary productivity in both seasons, whereas richness of seed‐eaters showed a positive association with habitat heterogeneity during the spring. Species composition similarity between assemblages was correlated with both productivity similarity and geographic proximity during the autumn and the spring, except for insectivore assemblages. Diversity within main trophic groups seemed to reflect differences in their spatial patterns as a response to changes between seasons in the spatial patterns of food resources. Our findings suggest that considering different seasons and functional groups in the analyses of diversity spatial pattern could contribute to better understand the determinants of biological diversity in temperate climates.     

Connectivity and vagility determine spatial richness gradients and diversification of freshwater fish in North America and Europe

The latitudinal species richness gradient (LRG) has been the subject of intense interest and many hypotheses but much less consideration has been given to longitudinal richness differences. The effect of postglacial dispersal, determined by connectivity and vagility, on richness was evaluated for the species‐poor European and North American Pacific and species‐rich Atlantic regional freshwater fish faunas. The numbers of species, by habitat, migration and distributional range categories, were determined from regional species lists for these three realms. The current orientation and past connections of drainage channels indicate that connectivity is greatest in the Atlantic and least in the Pacific. With increasing connectivity across realms, endemism decreased and postglacial recolonization increased, as did the LRG slope, with the greatest richness difference occurring between southern Atlantic and Pacific regions. Recolonizing species tended to be migratory, habitat generalists and from families of marine origin. Diversification, as indicated by species/genus ratios, probability of diversification, taxonomic distinctness and endemicity, declined with increasing latitude in all realms and was least in Europe. Richness patterns are consistent with an LRG driven by the time available for postglacial recolonization and by differences in dispersal ability, with richness differences across realms reflecting differences in dispersal and diversification.     

Explaining species richness from continents to communities: the time-for-speciation effect in emydid turtles

Speciation is the process that ultimately generates species richness. However, the time required for speciation to build up diversity in a region is rarely considered as an explanation for patterns of species richness. We explored this "time-for-speciation effect" on patterns of species richness in emydid turtles. Emydids show a striking pattern of high species richness in eastern North America (especially the southeast) and low diversity in other regions. At the continental scale, species richness is positively correlated with the amount of time emydids have been present and speciating in each region, with eastern North America being the ancestral region. Within eastern North America, higher regional species richness in the southeast is associated with smaller geographic range sizes and not greater local species richness in southern communities. We suggest that these patterns of geographic range size variation and local and regional species richness in eastern North America are caused by glaciation, allopatric speciation, and the time-for-speciation effect. We propose that allopatric speciation can simultaneously decrease geographic range size and increase regional diversity without increasing local diversity and that geographic range size can determine the relationship between alpha, beta, and gamma diversity. The time-for-speciation effect may act through a variety of processes at different spatial scales to determine diverse patterns of species richness.     

Species–energy relationships and habitat complexity in bird communities

Species–energy theory is a commonly invoked theory predicting a positive relationship between species richness and available energy. The More Individuals Hypothesis (MIH) attempts to explain this pattern, and assumes that areas with greater food resources support more individuals, and that communities with more individuals include more species. Using a large dataset for North American birds, I tested these predictions of the MIH, and also examined the effect of habitat complexity on community structure. I found qualitative support for the relationships predicted by the MIH, however, the MIH alone was inadequate for fully explaining richness patterns. Communities in more productive sites had more individuals, but they also had more even relative abundance distributions such that a given number of individuals yielded a greater number of species. Richness and evenness were also higher in structurally complex forests compared to structurally more simple grasslands when controlling for available energy.     

Local‐ to continental‐scale variation in the richness and composition of an aquatic food web

Aim We investigated patterns of species richness and composition of the aquatic food web found in the liquid‐filled leaves of the North American purple pitcher plant, Sarracenia purpurea (Sarraceniaceae), from local to continental scales. Location We sampled 20 pitcher‐plant communities at each of 39 sites spanning the geographic range of S. purpurea– from northern Florida to Newfoundland and westward to eastern British Columbia. Methods Environmental predictors of variation in species composition and species richness were measured at two different spatial scales: among pitchers within sites and among sites. Hierarchical Bayesian models were used to examine correlates and similarities of species richness and abundance within and among sites. Results Ninety‐two taxa of arthropods, protozoa and bacteria were identified in the 780 pitcher samples. The variation in the species composition of this multi‐trophic level community across the broad geographic range of the host plant was lower than the variation among pitchers within host‐plant populations. Variation among food webs in richness and composition was related to climate, pore‐water chemistry, pitcher‐plant morphology and leaf age. Variation in the abundance of the five most common invertebrates was also strongly related to pitcher morphology and site‐specific climatic and other environmental variables. Main conclusions The surprising result that these communities are more variable within their host‐plant populations than across North America suggests that the food web in S. purpurea leaves consists of two groups of species: (1) a core group of mostly obligate pitcher‐plant residents that have evolved strong requirements for the host plant and that co‐occur consistently across North America, and (2) a larger set of relatively uncommon, generalist taxa that co‐occur patchily.     

Species richness and soil properties in Pinus ponderosa forests: A structural equation modeling analysis

Question: How are the effects of mineral soil properties on understory plant species richness propagated through a network of processes involving the forest overstory, soil organic matter, soil nitrogen, and understory plant abundance? Location: North‐central Arizona, USA. Methods: We sampled 75 0.05‐ha plots across a broad soil gradient in a Pinus ponderosa (ponderosa pine) forest ecosystem. We evaluated multivariate models of plant species richness using structural equation modeling. Results: Richness was highest at intermediate levels of understory plant cover, suggesting that both colonization success and competitive exclusion can limit richness in this system. We did not detect a reciprocal positive effect of richness on plant cover. Richness was strongly related to soil nitrogen in the model, with evidence for both a direct negative effect and an indirect non‐linear relationship mediated through understory plant cover. Soil organic matter appeared to have a positive influence on understory richness that was independent of soil nitrogen. Richness was lowest where the forest overstory was densest, which can be explained through indirect effects on soil organic matter, soil nitrogen and understory cover. Finally, model results suggest a variety of direct and indirect processes whereby mineral soil properties can influence richness. Conclusions: Understory plant species richness and plant cover in P. ponderosa forests appear to be significantly influenced by soil organic matter and nitrogen, which are, in turn, related to overstory density and composition and mineral soil properties. Thus, soil properties can impose direct and indirect constraints on local species diversity in ponderosa pine forests.     

Patterns in the local assembly of Egyptian rodent faunas: Areography and species combinations

We implicate ecological processes in assembly patterns of Egyptian desert rodents using spatial variation of distribution and composition of local assemblages. We also compare our assemblages to prior analyses of North American and Australian small mammal communities in terms of species richness and representation of trophic groups.We studied patterns of occurrence of 29 rodent species among 335 collecting events in 308 sites using museum specimen records resulting from a country-wide survey. The studied taxa vary greatly in their natural histories and ecology. Fully 69% of studied species (20) were localized in 30 or fewer sites (9% of all sites). Site incidence was closely correlated with the geographic range size of species, while local abundance of a species showed little relationship to its geographic range. The species richness of local assemblages ranged from a lone species at 73 sites to 11 species at 3 sites. A total of 214 different combinations were recorded of which 164 (77%) were unique to a single site. G. gerbillus was both the most abundant and most ubiquitous species. Coexistence with other species was positively correlated with incidence and geographic range size. Body mass distribution was remarkably uniform for the fauna as a whole, and influenced the geographic distributions and abundance of individual species.Our sites have low species richness and substantial variability in species composition, which also characterize desert rodent communities elsewhere. Habitat requirements, exclusive distributions of sibling species, low primary productivity, and Egypt's location all influenced species assembly. The Egyptian and Australian deserts supported higher proportions of low richness assemblages compared to North America. As in North America, Egyptian sites were dominated by granivorous species.     

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.     

Broad-scale geographical patterns in local stream insect genera richness

Comprehensive global studies of stream invertebrate assemblages are rare and have produced contradictory results. To address this shortcoming, we compiled data from 495 published estimates of local genera richness for three orders of stream‐dwelling insects (Ephemeroptera, Plecoptera, Trichoptera) from throughout the world and used these data to describe global geographic patterns in stream insect genera richness and to address two questions: 1) does local stream insect richness vary more with regional historical factors or with local ecological factors?, and 2) to what extent have streams converged in the number of taxa they support?
Maximum genera richness varied sharply across the range of latitude examined from the south to north poles for all three orders of aquatic insects. Ephemeroptera richness showed 3 peaks (～30°S, 10°N, and 40°N) with highest richness near 5–10°N and 40°N latitude. Plecoptera richness was distinctly highest at ～40°N latitude with a similar peak at 40°S latitude. Trichoptera richness showed less latitudinal variation than the other taxa but was slightly higher near the equator and at 40°N and S latitude than at other latitudes. Genera richness generally declined with increasing elevation, except for Plecoptera. Maximum genera richness increased steadily with a measure of regional terrestrial net primary production and declined sharply with a measure of hydrologic disturbance for all orders. Richness varied widely among both biogeographical realms and biomes, although ca 2 times as much variation in richness was associated with biome as biogeographic realm. Richness for each order was highest in different biogeographic realms, but all orders had highest richness in broadleaf forest biomes. These latter results imply that spatial variation in local richness of stream insects is more strongly affected by contemporary ecological factors than by historical biogeography and that maintenance of intact forested landscapes may be critical to the conservation of stream invertebrate faunas.     

Evidence for gene flow in parasitic nematodes between two host species of shrews

We describe the genetic structure of populations of the intestinal nematode Longistriata caudabullata (Trichostrongyloidea: Heligmosomidae), a common parasite of short-tailed shrews (genus Blarina, Insectivora: Soricidae). Parasites and hosts were collected from a transect across a contact zone between two species of hosts, Blarina brevicauda and B. hylophaga, in central North America. An 800-base pairs (bp) fragment of the ND4 mitochondrial DNA (mtDNA) gene was sequenced for 28 worms and a 783-bp fragment of the mtDNA control region was analysed for 16 shrews. Phylogenetic analyses of mtDNA sequences revealed reciprocal monophyly for the shrew species, concordant with morphological diagnosis, and supported the idea that the transect cuts through a secondary contact zone between well-differentiated B. brevicauda and B. hylophaga. In contrast to this pattern, the parasitic nematode mtDNA phylogeny was not subdivided according to host affiliation. Genealogical discordance between parasite and host phylogenies suggests extensive gene flow among parasites across the host species boundary.     

The environmental basis of North American species richness patterns among Epicauta (Coleoptera: Meloidae)

Understanding regional variability in species richness is necessary for conservation efforts to succeed in the face of large-scale environmental deterioration. Several analyses of North American vertebrates have shown that climatic energy provides the best explanation of contemporary species richness patterns. The paucity of analyses of insect diversity patterns, however, remains a serious obstacle to a general hypothesis of spatial variation in diversity. We collected species distribution data on a North American beetle genus, Epicauta (Coleoptera: Meloidae) and tested several major diversity hypotheses. These beetles are generally grasshopper egg predators as larvae, and angiosperm herbivores as adults. Epicauta richness is highest in the hot, dry American southwest, and decreases north and east, consistent with the species richness-energy hypothesis. Potential evapotranspiration, which is also the best predictor of richness patterns among North American vertebrates, explains 80.2% of the variability in Epicauta species richness. Net primary productivity and variables measuring climatic heat energy only (such as PET) are not generally comparable, though they are sometimes treated as if they were equivalent. We conclude that the species richness-energy hypothesis currently provides a better overall explanation for Epicauta species richness patterns in North America than other major diversity hypotheses. The observed relationship between climatic energy and regional species richness may provide significant insight into the response of ecological communities to climate change.     

Spatial Dynamics of Body Size Frequency Distributions for North American Squamates

Scale dependent patterns of body size frequency distributions (BSFDs) have been explained by competition and an evolutionarily optimal body size in mammals. We test these ideas in a vertebrate group that is a model for competition and evolutionary studies by assessing the scale-dependence of BSFDs. BSFDs (body size defined as maximum total length) of North American squamates were assembled for the entire continent, biomes within the continent and local habitat patches within biomes. We described these distributions using skewness, kurtosis, interquartile range (IQR), and an index of evenness. We compared these parameters among spatial scales using Kolmogorov–Smirnov tests and bootstrap simulations. We assessed the relationship between body size and species richness using correlations (Pearsons and Spearmans R). The North American BSFD is bimodal, with a primary mode (240 mm) corresponding to lizards and small snakes and a secondary mode (912 mm) to snakes. Squamate BSFDs varied in a scale dependent fashion for some biomes and local habitat patches for kurtosis (12% of local patches and 10% of biomes more platykurtic), skewness (30% of biomes skewed to the right) and IQR (12% of patches increased). The index of evenness of BSFDs did not vary with spatial scale. Body size of biomes and local habitat patches closely resembles the North American BSFD as species richness increases. We found limited statistical support for the scale-dependency of North American squamate BSFDs (only 12–30% of patches or biomes conformed to the predicted pattern). These results suggest that the mechanisms implicated in scale-dependent patterns of BSFDs for mammals, geographic turnover of modal-sized species and competition within local assemblages may be of diminished importance in squamates. As geographic turnover of modal-sized species is theoretically linked to an evolutionarily optimal body size, this may suggest that optimal size theory is not adequate to predict spatial scaling of BSFDs in squamates.