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.     

Environmental factors, regional body size distributions and spatial variation in body size of local avian assemblages

Aim To determine how well variation in median body size of avian assemblages is predicted by (1) the environmental models usually employed in analyses of Bergmann's rule and (2) random sampling from the regional body size frequency distribution. If body size frequency distributions of local assemblages represent a random sample of a regional frequency distribution, then geographical variation in body sizes of assemblages might be a consequence of the determinants of spatial variation in species richness rather than direct influences on body size per se. Location Southern Africa. Methods Median body masses (as a measure of body size) of avian assemblages were calculated for quarter‐degree grid cells across South Africa and Lesotho. The relationship between median body mass and four environmental variables (minimum and maximum monthly temperatures, precipitation and seasonality in the normalized difference vegetation index, as a measure of seasonality in productivity) was examined using general linear models first without taking spatial autocorrelation into account, and then accounting for it by fitting an exponential spatial covariance structure. Model fit was assessed using the Akaike information criterion and Akaike weights. At each species richness value, random assemblages were sampled by either drawing species randomly from the regional body mass frequency distribution, or drawing species from the regional body mass frequency distribution with a probability proportional to their geographical distribution in the area. The ability of randomizations to predict actual body masses was examined using two‐tailed Fisher exact tests. Results Seasonality in productivity was the only environmental variable that remained a significant predictor of body mass variation in spatially explicit models, though the positive relationship was weak. When species richness was included in the models it remained the only significant predictor of size variation. Randomizations predicted median body mass poorly at low species richness, but well at high richness. Main conclusions Environmental models that have previously been proposed explain little of the variation in body mass across avian assemblages in South Africa. However, much of the variation in the median mass of assemblages could be predicted by randomly drawing species from the regional body mass frequency distribution, particularly using randomizations in which all species were drawn from the regional body mass frequency distribution with equal probability and at high species richness values. This outcome emphasizes the need to consider null expectations in investigations of the geographical variation in body size together with the probable environmental mechanisms underlying spatial variation in average size. Moreover, it suggests that in the South African avifauna, spatial variation in the body sizes of assemblages may be determined indirectly by the factors that influence geographical variation in species richness.     

Reconsidering null models of diversity: Do geometric constraints on species ranges necessarily cause a mid‐domain effect?

Recent null models that place species ranges randomly within a bounded domain have produced controversial results. Many such geometric constraint models predict a peak in species richness in the centre of domains in the absence of underlying environmental gradients or interspecific interactions. We used two-dimensional simulation models to explore different ways that species ranges could interact with the domain boundary. In the rejection model, a randomly generated range that overlaps a domain boundary is removed from the simulation. In the reshaping model, a range that overlaps the domain boundary is reshaped so that the entire range is placed within the domain. The truncation model allows potential ranges to extend across the boundary, but only that portion of the range within the domain is included in the realized range. Both rejection and reshaping models produced a drop in species richness near domain boundaries, though the effect was less pronounced in the reshaping model. Our truncation model did not produce any spatial pattern in species richness. Thus the random placement of species ranges within a bounded domain does not necessarily lead to a mid-domain effect.
  Range truncation is consistent with bioclimate envelope models, which can successfully predict a species range in response to the availability of appropriate climate conditions. We argue that such flexible range sizes are more realistic than the assumption that range size is an unvarying characteristic of a species. Other range characteristics, including size and shape, can change near domain boundaries in the null models, including the truncation model. A broader consideration of range characteristics near domain boundaries could be productive.     

When does diversity fit null model predictions? Scale and range size mediate the mid‐domain effect

Aim  Recently, a flurry of studies have focused on the extent to which geographical patterns of diversity fit mid-domain effect (MDE) null models. While some studies find strong support for MDE null models, others find little. We test two hypotheses that might explain this variation among studies: small-ranged groups of species are less likely than large-ranged species to show mid-domain peaks in species richness, and mid-domain null model predictions are less robust for smaller spatial extents than for larger spatial extents.
Location  We analyse data sets from elevational, riverine, continental and other domains from around the world.
Methods  We use a combination of Spearman rank correlations and binomial tests to examine whether differences within and among studies and domains in the predictive power of MDE null models vary with spatial scale and range size.
Results  Small-ranged groups of species are less likely to fit mid-domain predictions than large-ranged groups of species. At large spatial extents, diversity patterns of taxonomic groups with large mean range sizes fit MDE null model predictions better than did diversity patterns of groups with small mean range sizes. MDE predictions were more explanatory at larger spatial extents than at smaller extents. Diversity patterns at smaller spatial extents fit MDE predictions poorly across all range sizes. Thus, MDE predictions should be expected to explain patterns of species richness when ranges and the scale of analysis are both large.
Main conclusions  Taken together, the support for these hypotheses offers a more sophisticated model of when MDE predictions should be expected to explain patterns of species richness, namely when ranges and the scale of analysis are both large. Thus the circumstances in which the MDE is important are finite and apparently predictable.     

Species-pool hypothesis: Limits to its testing

The species richness of a community depends both on the pool of available species and on biotic mechanisms that lead to the exclusion of some of the species from a community. The method suggested byPärtel et al.,Oikos 75: 111–117, 1996 to test the effect of species-pool size on the species richness of a community is discussed in this paper. This method is based on the calculation of a correlation between the actual species richness of a community and the actual species-pool size, and compares the correlation found in the data with that simulated by a null model. In the null model, the species richness has a uniform distribution between zero and the size of the species pool. A correlation significantly higher than that in the null model is interpreted as evidence of the greater role of species pool than of local interactions in formation of community species richness (Zobel,Folia Geobot. 36: 3–8, 2001). It is shown that the interpretation of discrepancies between the null model and reality is difficult or impossible, because: (1) a null model with a uniform distribution of species richness is unrealistic, (2) both models based on the random selection of species from a species pool and models that include competitive interactions in the community predict a higher positive correlation of species richness and size of the species pool than the null model, and (3) local species richness might be affected by species-pool size, but a large species-pool size can also be a result of high local species richness. Caution is urged when interpreting the analyses based on the size of the filtered species-pool size.     

Using potential distributions to explore determinants of Western Palaearctic migratory songbird species richness in sub-Saharan Africa

Aim  Although the breeding ranges of most Western Palaearctic migratory passerines are well documented in Europe, their overwintering ranges and patterns of species richness in Africa remain poorly understood. To illustrate potential patterns of species richness despite severely limited data, we extrapolated species ranges from a new and unique data bank of locality records that documents overwintering locations of these birds in Africa.
Location  Sub-Saharan Africa.
Methods  We predicted potential geographical distributions of 60 species of passerine birds based on overwintering records using bioclimatic models. We then combined these predictions to estimate potential species richness and explored response shapes using spatial linear regression. We also evaluated the evidence for a mid-domain effect using a one-dimensional null model.
Results  Spatial linear regression analyses of the species richness pattern revealed non-linear relationships to seasonality in precipitation, minimum net primary productivity, minimum average temperature, habitat heterogeneity, percentage of tree cover, distance from the Sahara Desert and inter-annual variability in net primary productivity. The explanatory power of these variables decreased with geographic range size. The one-dimensional null model of species richness based on distance from the Sahara Desert did not show evidence of a mid-domain effect.
Main conclusions  Distributions of migrants seem generally strongly determined by distance from the Sahara Desert working in concert with climatic effects, but this cannot adequately explain richness patterns of species with small ranges in Africa, many of which are of substantial conservation concern.     

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.     

Seedling and Sapling Dynamics of Treefall Pits in Puerto Rico1

Seedling and sapling dynamics in a Puerto Rican rain forest were compared between forest understory and soil pits created by the uprooting of 27 trees during Hurricane Hugo. Soil N and P, organic matter, and soil moisture were lower and bulk densities were higher in the disturbed mineral soils of the pits than in undisturbed forest soils ten months after the hurricane. No differences in N and P levels were found in pit or forest soils under two trees with N–fixing symbionts (Inga laurina and Ormosia krugii) compared to soils under a tree species without N–fixing sym–bionts (Casearia arborea), but other soil variables (Al, Fe, K) did vary by tree species. Forest plots had greater species richness of seedlings (<10 cm tall) and saplings (10–100 cm tall) than plots in the soil pits (and greater sapling densities), but seedling densities were similar between plot types. Species richness and seedling densities did not vary among plots associated with the three tree species, but some saplings were more abundant under trees of the same species. Pit size did not affect species richness or seedling and sapling densities. Recruitment of young Cecropia schreberiana trees (>5 m tall) 45 months after the hurricane was entirely from the soil pits, with no tree recruitment from forest plots. Larger soil pits had more tree recruitment than smaller pits. Defoliation of the forest by the hurricane created a large but temporary increase in light availability. Recruitment of C. schreberiana to the canopy occurred in gaps created by the treefall pits that had lower soil nutrients but provided a longer–term increase in light availability. Treefall pits also significantly altered the recruitment and mortality of many understory species in the Puerto Rican rain forest but did not alter species richness.     

Mammalian species richness and morphological complexity along an elevational gradient in the arid south-west

Aim We examined the relationship between species richness and morphological complexity of terrestrial mammal communities along an elevational gradient. Location The gradient is in the Sonoran Desert in Southern California and extends from a sand dune habitat near sea level to coniferous forest ending at >2600 m. Methods Morphological diversity, characterized by both size and shape of coexisting mammal species, was estimated within and between sites from projections of variables on principal components axes. Similarities among species were calculated as Euclidean distances. To tease apart size and shape, we constructed two principal component analyses: one based on log-transformed original measurements, the other on log-transformed proportional shape variables. To test whether species number accounted for the morphological diversity at each site we designed two null models. The models generated were random communities generated from the forty-two-species pool. Indices of morphological diversity for real communities were compared with the results of 500 simulations of each null model. Results Species richness varied along the gradient, peaking in the mid-elevation agave-ocotillo habitat. Morphological diversity of shapes and sizes correlated strongly with species richness. Locomotor, tooth, and skull traits were all important in distinguishing among species. Main conclusions Two important patterns emerged: (1) diversity of both sizes and of shapes of species within communities correlated positively with species number, and both sets of variables behaved similarly across this gradient; (2) the most species rich sites were not composed of specialists on these best places, but rather, a community of species derived from overlapping faunal groups.     

Geographic range size, life history and rates of diversification in Australian mammals

Abstract What causes species richness to vary among different groups of organisms? Two hypotheses are that large geographical ranges and fast life history either reduce extinction rates or raise speciation rates, elevating a clade's rate of diversification. Here we present a comparative analysis of these hypotheses using data on the phylogenetic relationships, geographical ranges and life history of the terrestrial mammal fauna of Australia. By comparing species richness patterns to null models, we show that species are distributed nonrandomly among genera. Using sister‐clade comparisons to control for clade age, we then find that faster diversification is significantly associated with larger geographical ranges and larger litters, but there is no evidence for an effect of body size or age at first breeding on diversification rates. We believe the most likely explanation for these patterns is that larger litters and geographical ranges increase diversification rates because they buffer species from extinction. We also discuss the possibility that positive effects of litter size and range size on diversification rates result from elevated speciation rates.     

The diversity field of New World leaf‐nosed bats (Phyllostomidae)

Aim To analyse how the patterns of species richness for the whole family Phyllostomidae determine the structure of diversity fields (sets of species‐richness values) within the ranges of individual bat species. Location The range of the family Phyllostomidae in North and South America. Methods We generated a database of the occurrence of 143 phyllostomid bat species in 6794 quadrats, analysing the species‐richness frequency distribution for all sites, and for subsets of sites defined by the geographic ranges of species. Range–diversity plots, depicting simultaneously the size and the mean species richness of ranges, were built to explore the patterns of co‐occurrence in widespread and restricted species. We compared the empirical patterns against two null models: (1) with scattered (non‐cohesive) ranges, and (2) with cohesive ranges modelled with the spreading‐dye algorithm. Diversity fields were analysed with richness maps for individual species and with comparisons of species‐richness frequency distributions. Results Overall richness frequency distribution showed a multimodal pattern, whereas simulated distributions showed lower values of variance, and were unimodal (for model 1) and bimodal (for model 2). Range–diversity plots for the empirical data and for the cohesive‐ranges simulation showed a strong tendency of species to co‐occur in high‐diversity sites. The scattered‐ranges simulation showed no such tendency. Diversity fields varied according to idiosyncratic features of species generating particular geographic patterns and richness frequency distributions. Main conclusions Phyllostomid bats show a higher level of co‐occurrence than expected from null models. That tendency in turn implies a higher variance in species richness among sites, generating a wider species‐richness frequency distribution. The diversity field of individual species results from the size, shape and location of ranges, but also depends on the general pattern of richness for the whole family.     

Challenges in the application of geometric constraint models

Discerning the processes influencing geographical patterns of species richness remains one of the central goals of modern ecology. Traditional approaches to exploring these patterns have focused on environmental and ecological correlates of observed species richness. Recently, some have suggested these approaches suffer from the lack of an appropriate null model that accounts for species ranges being constrained to occur within a bounded domain. Proponents of these null geometric constraint models (GCMs), and the mid-domain effect these models produce, argue their utility in identifying meaningful gradients in species richness. This idea has generated substantial debate. Here we discuss what we believe are the three major challenges in the application of GCMs. First, we argue that there are actually two equally valid null models for the random placement of species ranges within a domain, one of which actually predicts a uniform distribution of species richness. Second, we highlight the numerous decisions that must be made to implement a GCM that lead to marked differences in the predictions of the null model. Finally, we discuss challenges in evaluating the importance of GCMs once they have been implemented.     

Macroevolution of hoverflies (Diptera: Syrphidae): the effect of using higher‐level taxa in studies of biodiversity,and correlates of species richness

We test a near‐complete genus level phylogeny of hoverflies (Diptera: Syrphidae) for consistency with a null model of clade growth having uniform probabilities of speciation and extinction among contemporaneous species. The phylogeny is too unbalanced for this null model. Importantly, the degree of imbalance in the phylogeny depends on whether the phylogeny is analysed at the genus level or species level, suggesting that genera ought not to be used uncritically as surrogates for species in large‐scale evolutionary analyses. Tests for a range of morphological, life‐history and ecological correlates of diversity give equivocal results, but suggest that high species‐richness may be associated with sexual selection and diet breadth. We find no correlation between species‐richness and either body size or reproductive rate.     

Anthropogenic modification disrupts species co‐occurrence in stream invertebrates

The question of whether species co‐occurrence is random or deterministic has received considerable attention, but little is known about how anthropogenic disturbance mediates the outcomes. By combining experiments, field surveys and analysis against null models, we tested the hypothesis that anthropogenic habitat modification disrupts species co‐occurrence in stream invertebrates across spatial scales. Whereas communities in unmodified conditions were structured deterministically with significant species segregation, catchment‐scale conversion to agriculture and sediment deposition at the patch‐ or micro‐habitat scale apparently randomized species co‐occurrences. This shift from non‐random to random was mostly independent of species richness, abundance and spatial scale. Data on community‐wide life‐history traits (body size, dispersal ability and predatory habits) and beta‐diversity indicated that anthropogenic modification disrupted community assembly by affecting biotic interactions and, to a lesser extent, altering habitat heterogeneity. These data illustrate that the balance between predictable and stochastic patterns in communities can reflect anthropogenic modifications that not only transcend scales but also change the relative forces that determine species coexistence. Research into the effects of habitat modification as a key to understanding global change should extend beyond species richness and composition to include species co‐occurrence, species interactions and any functional consequences.     

Latitudinal patterns of range size and species richness of New World woody plants

Aim  Relationships between range size and species richness are contentious, yet they are key to testing the various hypotheses that attempt to explain latitudinal diversity gradients. Our goal is to utilize the largest data set yet compiled for New World woody plant biogeography to describe and assess these relationships between species richness and range size.
Location  North and South America.
Methods  We estimated the latitudinal extent of 12,980 species of woody plants (trees, shrubs, lianas). From these estimates we quantified latitudinal patterns of species richness and range size. We compared our observations with expectations derived from two null models.
Results   Peak richness and the smallest- and largest-ranged species are generally found close to the equator. In contrast to prominent diversity hypotheses: (1) mean latitudinal extent of tropical species is greater than expected; (2) latitudinal extent appears to be decoupled from species richness across New World latitudes, with abrupt transitions across subtropical latitudes; and (3) mean latitudinal extents show equatorial and north temperate peaks and subtropical minima. Our results suggest that patterns of range size and richness appear to be influenced by three broadly overlapping biotic domains (biotic provinces) for New World woody plants.
Main conclusions  Hypotheses that assume a direct relationship between range size and species richness may explain richness patterns within these domains, but cannot explain gradients in richness across the New World.     

Effects of repeated burning on species richness in a Florida pine savanna: A test of the intermediate disturbance hypothesis

Abstract. We studied the effect of burning frequency on the density and species richness of understory flowering stems in a Florida sandhill. Flowering stems were censused weekly for 54 weeks in six sites that had been burned one to six times in the previous 16 years. We concurrently measured overstory characteristics such as species composition, density and basal area. We used maximum likelihood and Akaike's Information Criterion to compare linear, quadratic, saturating, and null models of community response to repeating burning. We did not find a relationship between species richness, diversity or flowering stem density and fire frequency. Tree density was related to fire frequency and may represent an indirect pathway for fire effects on understory characteristics. While we found no support for the Intermediate Disturbance Hypothesis, an analysis of our experimental design indicated that we had low statistical power. We develop the hypothesis that a saturating model of response to fire best describes understory species richness in our system. We test this hypothesis using the most extensive published fire data set we are aware of and find support for a saturating model.     

Niche differentiation depends on body size in a cichlid fish: a model system of a community structured according to size regularities

1. Communities of different species are often structured according to niche differentiation associated with competitive interactions. We show that similar principles may apply on an ecological time-scale when individuals of a species having a wide size variation compete for resources, using the Lake Tanganyika cichlid Lobochilotes labiatus (5-30 cm). This species has a mouth especially adapted to suck up invertebrates from rock crevices. 2. Individuals defended feeding territories against similar-sized conspecifics, but not against different-sized ones. Thus, territories of similar-sized fish rarely overlapped, but up to a total of seven individuals (of seven size-ranks) had broadly overlapping territories with dissimilar-sized individuals. Comparison with expectation from the null model demonstrated clearly that observed size ratios between adjacent size rank were determined non-randomly regardless of sexual combinations. 3. Larger individuals took larger prey types of larger average size, but more importantly used wider rock crevices from which to suck food than smaller individuals. We calculated pairwise values of Schoener's index of diet overlap C(d) and the values of Levin's index of diet breadth B(d) (prey type and prey size) and the same for the width of the rock crevices used for foraging (C(r) and B(r)). C(d) remained high among all combinations of the seven ranks. In contrast, C(r) declined strongly in combinations of adjacent ranks (to 0.27), and was low or zero among further different size ranks. This shows that fish with overlapping territories divided the food resources largely through foraging site partitioning. Accordingly, B(d) did not depend on the size difference to the nearest two coinhabiting fish, whereas B(r) did. 4. We conclude that this L. labiatus community is structured non-randomly: body size-dependent effects on foraging site usage result in competition with, and territorial exclusion of, similar-sized individuals, but not of dissimilar-sized individuals that were accepted as coinhabitants. Accordingly, mean body size ratios (large/small) between two adjacent ranks were consistently approximately 1.28 [standard deviation (SD) = 0.07, n = 104], while approximately 1.34 from the null model (SD = 0.34, n = 10 400 simulations). We discuss our results as an example of Hutchinson's rule, applied originally to size ratios of different species.     

Worldwide patterns in species richness of Falconiformes: analytical null models, geometric constraints, and the mid-domain effect.

Recently, the hypothesis that the geographic distribution of species could be influenced by the shape of the domain edges, the so-called Mid-Domain Effect (MDE), has been included as one of the five credible hypotheses for explaining spatial gradients in species richness, despite all the unsuccessful current attempts to prove empirically the validity of MDE. We used data on spatial worldwide distributions of Falconiformes to evaluate the validity of MDE assumptions, incorporated into two different sorts of null models at a global level and separately across five domains/landmasses. Species richness values predicted by the null models of the MDE and those values predicted by Net Primary Productivity, a surrogate variable expressing the effect of available energy, were compared in order to evaluate which hypothesis better predicts the observed values. Our tests showed that MDE continues to lack empirical support, regardless of its current acceptability, and so, does not deserve to be classified as one possible explanation of species richness gradients.     

Testing multiple assembly rule models in avian communities on islands of an inundated lake,Zhejiang Province,China

Aim A fundamental question in community ecology is whether general assembly rules determine the structure of natural communities. Although many types of assembly rules have been described, including Diamond’s assembly rules, constant body‐size ratios, favoured states, and nestedness, few studies have tested multiple assembly rule models simultaneously. Therefore, little is known about the relative importance of potential underlying factors such as interspecific competition, inter‐guild competition, selective extinction and habitat nestedness in structuring community composition. Here, we test the above four assembly rule models and examine the causal basis for the observed patterns using bird data collected on islands of an inundated lake. Location Thousand Island Lake, China. Methods  We collected data on presence–absence matrices, body size and functional groups for bird assemblages on 42 islands from 2007 to 2009. To test the above four assembly rule models, we used null model analyses to compare observed species co‐occurrence patterns, body‐size distributions and functional group distributions with randomly generated assemblages. To ensure that the results were not biased by the inclusion of species with extremely different ecologies, we conducted separate analyses for the entire assemblage and for various subset matrices classified according to foraging guilds. Results The bird assemblages did not support predictions by several competitively structured assembly rule models, including Diamond’s assembly rules, constant body‐size ratios, and favoured states. In contrast, bird assemblages were highly significantly nested and were apparently shaped by extinction processes mediated through area effects and habitat nestedness. The nestedness of bird assemblages was not a result of passive sampling or selective colonization. These results were very consistent, regardless of whether the entire assemblage or the subset matrices were analysed. Main conclusions Our results suggest that bird assemblages were shaped by extinction processes mediated through area effects and habitat nestedness, rather than by interspecific or inter‐guild competition. From a conservation point of view, our results indicate that we should protect both the largest islands with the most species‐rich communities and habitat‐rich islands in order to maximize the number of species preserved.     

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.