Patterns of species richness and turnover along the pH gradient in deciduous forests: testing the continuum hypothesis

Question: (i) How do species richness and species turnover change along a pH gradient? (ii) What are possible driving factors behind these patterns? (iii) Can the observed patterns be explained by an individualistic continuum concept that postulates independence of species responses and constant turnover rates? Location: Semi‐natural, deciduous hardwood forests in NW Germany (558 plots). Methods: The instantaneous rate of compositional turnover is measured by the sum of slope angles of modelled response curves (119 understorey species) at any point along the pH gradient. Total turnover rate, positive turnover rate (species increasing in probability of occurrence) and negative turnover rate (species decreasing in probability of occurrence) are calculated separately. Species richness is modelled using Poisson regression and by calculating the sum of predicted probabilities at any gradient point. Turnover rates are compared with those calculated from a null model based on a Gaussian community model. Soil chemical analyses of 49 plots are used to interpret biodiversity patterns. Results: Species richness shows a hump‐shaped relation to pH(CaCl₂) with a minor decline at approximately pH>5.0. The decline is possibly due to the confounding influence of water regime and local species pool effects. Increasing richness from pH 2.5 to 4.7 can be traced back to positive turnover exceeding negative turnover. Peaks in turnover rates, dominated by positive turnover, are located at pH 3.7 and 2.8, where turnover rates considerably exceed rates derived from the null model. The turnover pattern can be related to soil chemical conditions, e.g. decreasing base saturation, Al and H⁺ toxicity and the occurrence of mor. Conclusions: The high turnover rates and the massive imbalance in positive and negative turnover rates found in deciduous forests cannot be explained by the individualistic continuum concept. Physiological constraints at the gradient limits and species pool effects could be responsible for this. Their role should be considered more explicitly in vegetation concepts dealing with the continuum‐discontinuum controversy. The presented approach can be regarded as a comprehensive analytical tool for a better understanding of biodiversity patterns along environmental gradients by linking species richness, turnover and response curve types.     

Grazing effects on the species‐area relationship: Variation along a climatic gradient in NE Spain

Questions: Does grazing have the same effect on plant species richness at different spatial scales? Does the effect of spatial scale vary under different climatic conditions and vegetation types? Does the slope of the species‐area curve change with grazing intensity similarly under different climatic conditions and vegetation types? Location: Pastures along a climatic gradient in northeastern Spain. Methods: In zones under different regimes of sheep grazing (high‐, low‐pressure, abandonment), plant species richness was measured in different plot sizes (from 0.01 to 100 m2) and the slope of the species‐area curves was calculated. The study was replicated in five different locations along a climatic gradient from lowland semi‐arid rangelands to upland moist grasslands. Results: Species richness tended to increase with grazing intensity at all spatial scales in the moist upland locations. On the contrary, in the most arid locations, richness tended to decrease, or remain unchanged, with grazing due to increased bare soil. Grazing differentially affected the slope (z) of the species‐area curve (power function S=c A^z) in different climatic conditions: z tended to increase with grazing in arid areas and decrease in moist‐upland ones. ß‐diversity followed similar pattern as z. Conclusions: Results confirm that the impact of grazing on plant species richness are spatial‐scale dependent. However, the effects on the species‐area relationship vary under different climatic conditions. This offers a novel insight on the patterns behind the different effects of grazing on diversity in moist vs. arid conditions reported in the literature. It is argued that the effect of spatial scale varies because of the different interaction between grazing and the intrinsic spatial structure of the vegetation. Variations in species‐area curves with grazing along moisture gradients suggest also a different balance of spatial components of diversity (i.e. a‐ and ß‐diversity).     

Can environmental constraints determine random patterns of plant species co-occurrence?

Plant community ecologists use the null model approach to infer assembly processes from observed patterns of species co‐occurrence. In about a third of published studies, the null hypothesis of random assembly cannot be rejected. When this occurs, plant ecologists interpret that the observed random pattern is not environmentally constrained – but probably generated by stochastic processes. The null model approach (using the C‐score and the discrepancy index) was used to test for random assembly under two simulation algorithms. Logistic regression, distance‐based redundancy analysis, and constrained ordination were used to test for environmental determinism (species segregation along environmental gradients or turnover and species aggregation). This article introduces an environmentally determined community of alpine hydrophytes that presents itself as randomly assembled. The pathway through which the random pattern arises in this community is suggested to be as follows: Two simultaneous environmental processes, one leading to species aggregation and the other leading to species segregation, concurrently generate the observed pattern, which results to be neither aggregated nor segregated – but random. A simulation study supports this suggestion. Although apparently simple, the null model approach seems to assume that a single ecological factor prevails or that if several factors decisively influence the community, then they all exert their influence in the same direction, generating either aggregation or segregation. As these assumptions are unlikely to hold in most cases and assembly processes cannot be inferred from random patterns, we would like to propose plant ecologists to investigate specifically the ecological processes responsible for observed random patterns, instead of trying to infer processes from patterns.     

Diverse responses across soil parent materials during ecological restoration

A major challenge to advancing the science and practice of ecological restoration is working across large landscapes containing diverse sites that may respond differently to restoration. We conducted a 5‐year restoration experiment, replicated across 9 sites spanning 3 soil parent material types within a 9,000‐ha Pinus ponderosa forest landscape. We evaluated plant community response to restoration Pinus thinning, grazing, and aqueous smoke application. We measured vegetation before (2003) and 3 (2006) and 5 (2008) years after treatment. Plant community responses of species richness, cover, and composition were diverse, ranging from increases, decreases, or no change depending on soil parent material, tree thinning, and presence or exclusion of grazing. Restoration outcomes were under hierarchical control: soil parent material constrained response to Pinus thinning, which in turn influenced grazing effects. On limestone‐derived soil, responses included no change in species richness but increased plant cover with Pinus thinning. Both plant richness and cover increased on benmorite soil after thinning, and cover generally increased more without grazing. On rocky, basalt soil, plant richness increased but cover did not after any treatment. Diversity of responses to restoration has implications for: (1) setting goals or monitoring indicators tailored to inherent soil capability; (2) identifying where grazing most affects restoration outcomes; and (3) forecasting responses to restoration across landscapes. Diverse responses to restoration along physiographic gradients such as soil parent material warrant consideration when developing restoration across degraded landscapes.     

Grazing without grasses: Effects of introduced livestock on plant community composition in an arid environment in northern Patagonia

Question: How does grazing intensity affect plant density, cover and species richness in an Patagonian arid ecosystem? Location: Monte steppe ecoregion, SW Argentina. Methods: I analysed the effect of grazing on plant density, cover and species richness using a stocking rate gradient within the same habitat. Six paddocks were used with stocking rates ranging between 0.002 – 0.038 livestock/ha. Plant density, species richness, plant cover and percentage of grazed branches were determined by sampling plots within each paddock. The percentage of grazed branches was used as an independent measurement of grazing intensity. Results: Higher stocking rates were related to lower plant density, species richness and plant cover. The paddock with the lowest grazing intensity had 86% more plants per unit area, 63% more plant cover and 48% higher species richness. The percentage of grazed branches and the quantity of dung increased with stocking rate. Conclusions: Introduced livestock seriously affect native vegetation in the Patagonian Monte. The damage observed in this xerophytic plant community suggests that plant adaptations to aridity do not provide an advantage to tolerate or avoid grazing by vertebrate herbivores in this region. Plant degradation in this arid environment is comparable to the degradation found in more humid ecosystems.     

Community versus single‐species distribution models for British plants

Aim Species distribution models are increasingly used to predict the impacts of global change on whole ecological communities by modelling the individualistic niche responses of large numbers of species. However, it is not clear whether this single‐species ensemble approach is preferable to community‐wide strategies that represent interspecific associations or shared responses to environmental gradients. Here, we test the performance of two multi‐species modelling approaches against equivalent single‐species models. Location Great Britain. Methods Single‐ and multi‐species distribution models were fitted for 701 native British plant species at a 10‐km grid scale. Two machine learning methods were used – classification and regression trees (CARTs) and artificial neural networks (ANNs). The single‐species versions are widely used in ecology but their multivariate extensions are less well known and have not previously been evaluated against one another. We compared their abilities to predict species distributions, community compositions and species richness in an independent geographical region reserved from model‐fitting. Results The single‐ and multi‐species models performed similarly, although the community models gave slightly poorer predictive accuracy by all measures. However, from the point of view of the whole community they were much simpler than the array of single‐species models, involving orders of magnitude fewer parameters. Multi‐species approaches also left greater residual spatial autocorrelation than the individualistic models and, contrary to expectation, were relatively less accurate for rarer species. However, the fitted multi‐species response curves had lower tendency for pronounced discontinuities that are unlikely to be a feature of realized niche responses. Main conclusions Although community distribution models were slightly less accurate than single‐species models, they offered a highly simplified way of modelling spatial patterns in British plant diversity. Moreover, an advantage of the multi‐species approach was that the modelling of shared environmental responses resolved more realistic response curves. However, there was a slight tendency for community models to predict rare species less accurately, which is potentially disadvantageous for conservation applications. We conclude that multi‐species distribution models may have potential for understanding and predicting the structure of ecological communities, but were slightly inferior to single‐species ensembles for our data.     

Comparing species interaction networks along environmental gradients

Knowledge of species composition and their interactions, in the form of interaction networks, is required to understand processes shaping their distribution over time and space. As such, comparing ecological networks along environmental gradients represents a promising new research avenue to understand the organization of life. Variation in the position and intensity of links within networks along environmental gradients may be driven by turnover in species composition, by variation in species abundances and by abiotic influences on species interactions. While investigating changes in species composition has a long tradition, so far only a limited number of studies have examined changes in species interactions between networks, often with differing approaches. Here, we review studies investigating variation in network structures along environmental gradients, highlighting how methodological decisions about standardization can influence their conclusions. Due to their complexity, variation among ecological networks is frequently studied using properties that summarize the distribution or topology of interactions such as number of links, connectance, or modularity. These properties can either be compared directly or using a procedure of standardization. While measures of network structure can be directly related to changes along environmental gradients, standardization is frequently used to facilitate interpretation of variation in network properties by controlling for some co‐variables, or via null models. Null models allow comparing the deviation of empirical networks from random expectations and are expected to provide a more mechanistic understanding of the factors shaping ecological networks when they are coupled with functional traits. As an illustration, we compare approaches to quantify the role of trait matching in driving the structure of plant–hummingbird mutualistic networks, i.e. a direct comparison, standardized by null models and hypothesis‐based metaweb. Overall, our analysis warns against a comparison of studies that rely on distinct forms of standardization, as they are likely to highlight different signals. Fostering a better understanding of the analytical tools available and the signal they detect will help produce deeper insights into how and why ecological networks vary along environmental gradients.     

Biodiversity patterns of vascular plant species in mountain vegetation in the Faroe Islands

Biodiversity patterns of vascular plant species were studied along altitudinal gradients in the Faroe Islands. Plants were sampled from five different mountains (150–856 m a.s.l.) at 50 m altitudinal intervals. Included in the study were 107 vascular plant species. In order to compare only altitudes with the same number of plots, three different analyses were carried out. One analysis included five mountains from 250 to 750 m a.s.l., one had three mountains from 150 to 750 m a.s.l., and the last one had two mountains from 750 to 850 m a.s.l. The patterns of biodiversity were evaluated on the basis of species richness as the total number of species at each altitudinal interval, as species turnover between altitudes and in relation to the Shannon‐Wiener index. Similar patterns were found for species richness in the three analyses, although richness was higher along the whole transect when five mountains were included. For the Shannon‐Wiener index, only small differences were found among the three analyses. A maximum was seen at 250 m a.s.l. and again at 500 m a.s.l. both in richness and in the Shannon‐Wiener index. Maximum species turnover was found at mid‐altitudes. Total vegetation cover followed the same pattern as richness. In addition to climate, the altitudinal variation of biodiversity may be affected by grazing.     

Plant species richness responses to grazing protection and degradation history in a low productivity landscape

Questions: Does species richness and abundance accumulate with grazing protection in low productivity ecosystems with a short evolutionary history of grazing, as predicted by emerging theory? How do responses to grazing protection inform degradation history? Location: Mulga (Acacia aneura) dry forest, eastern Australia, generally considered chronically degraded by livestock grazing. Methods: Three paired exclosures (ungrazed, and macropod‐grazed) were compared with open‐grazed areas after 25 years using quadrats located on either side of the fences. Additionally, the regional flora for mulga dry forest was assessed to identify species that may have declined and could be threatened by grazing. Results: Low herbaceous biomass accumulation (<1.3 t ha^?1) with full grazing protection confirmed a low productivity environment. For most plant life forms the highest species richness was in macropod‐grazed exclosures, an intermediate grazing disturbance that best approximates the evolutionary history of the environment. This was the net outcome of species that both declined and increased in response to grazing. Regeneration and subsequent self‐thinning of mulga was promoted with grazing protection, but did not confound interpretation of species richness and abundance responses. At the regional scale only 11 native species out of 407 comprising the mulga dry forest flora were identified as rare and potentially threatened by grazing. Conclusions: Significant increases in richness or abundance of native plants with grazing protection, persistence of perennial grasses, regeneration of mulga and scant evidence of a major decline in the regional flora are not consistent with established assertions that long‐grazed mulga dry forest has crossed functional thresholds that limit recovery. Further, a peak in species richness under intermediate (macropod) grazing is counter to the shape of the response predicted by emerging theory for recovery of species richness in a low productivity environment. The finding prompts a more thorough understanding of the distinction between environments with inherently low productivity and those degraded by grazing.     

Mapping landscape beta diversity of plants across KwaZulu-Natal,South Africa,for aiding conservation planning

Collective properties of biodiversity, such as beta diversity, are suggested as complementary measures of species richness to guide the prioritisation and selection of important biodiversity areas in regional conservation planning. We assessed variation in the rate of plant species turnover along and between environmental gradients in KwaZulu-Natal, South Africa using generalised dissimilarity modelling, in order to map landscape levels of floristic beta diversity. Our dataset consisted of 434 plots (1000 m²) containing 997 grassland and savanna matrix species. Our model explained 79 % of the null deviance observed in floristic dissimilarities. Variable rates of turnover existed along the major environmental gradients of mean annual temperature, median rainfall in February, and soil cation exchange capacity, as well as along gradients of geographical distance. Beta diversity was highest in relatively warm, drier summer regions and on dystrophic soils. Areas of high beta diversity identify areas that should be included in conservation plans to maximise representation of diversity and highlight areas best suited to protected area expansion. Biome transition areas in high beta diversity areas may be susceptible to climate variability. Including beta diversity turnover rates in regional conservation plans will help to preserve evolutionary and ecological processes that create and maintain diversity.     

Plant species and growth form richness along altitudinal gradients in the southwest Ethiopian highlands

Questions: Do growth forms and vascular plant richness follow similar patterns along an altitudinal gradient? What are the driving mechanisms that structure richness patterns at the landscape scale? Location: Southwest Ethiopian highlands. Methods: Floristic and environmental data were collected from 74 plots, each covering 400 m². The plots were distributed along altitudinal gradients. Boosted regression trees were used to derive the patterns of richness distribution along altitudinal gradients. Results: Total vascular plant richness did not show any strong response to altitude. Contrasting patterns of richness were observed for several growth forms. Woody, graminoid and climber species richness showed a unimodal structure. However, each of these morphological groups had a peak of richness at different altitudes: graminoid species attained maximum importance at a lower elevations, followed by climbers and finally woody species at higher elevations. Fern species richness increased monotonically towards higher altitudes, but herbaceous richness had a dented structure at mid‐altitudes. Soil sand fraction, silt, slope and organic matter were found to contribute a considerable amount of the predicted variance of richness for total vascular plants and growth forms. Main Conclusions: Hump‐shaped species richness patterns were observed for several growth forms. A mid‐altitudinal richness peak was the result of a combination of climate‐related water–energy dynamics, species–area relationships and local environmental factors, which have direct effects on plant physiological performance. However, altitude represents the composite gradient of several environmental variables that were interrelated. Thus, considering multiple gradients would provide a better picture of richness and the potential mechanisms responsible for the distribution of biodiversity in high‐mountain regions of the tropics.     

Comparison of biodiversity and ground cover between a commercial rotationally grazed property and an adjacent nature reserve in semi‐arid rangeland

Continuous livestock grazing can have negative effects on biodiversity and landscape function in arid and semi‐arid rangelands. Alternative grazing management practices, such as rotational grazing, may be a viable option for broad‐scale biodiversity conservation and sustainable pastoral management. This study compared ground cover, plant species composition and floristic and functional diversity along gradients of grazing intensity between a pastoral property rotationally grazed by goats and an adjacent nature reserve (ungrazed by commercial livestock) in semi‐arid south‐eastern Australia. Understorey plant species composition differed significantly between the rotationally grazed property and the nature reserve, with a greater proportion and frequency of palatable species recorded in the nature reserve. Understorey plant species richness, diversity, functional biodiversity measures and ground cover declined with increasing grazing pressure close to water points under commercial rotational grazing management. However, at a whole‐paddock scale, there were few differences in plant biodiversity and ground cover between the rotationally grazed property and the nature reserve, despite differences in overall plant species composition. Flexible, adaptive, rotational grazing should be investigated further for its potential to achieve both socio‐economic and biodiversity conservation outcomes in semi‐arid rangelands to complement existing conservation reserves.     

Effect of sampling grain on patterns of species richness and turnover in Amazonian forests

Grain (size of sampling units) affects the spatial resolution at which ecological patterns can be observed and analysed, and potentially has an important effect on the results of broad‐scale studies on diversity gradients. Here we examine the effect of grain on patterns of species richness and turnover in lowland rainforests of western Amazonia (Peru and Colombia). We inventoried pteridophytes (ferns and lycophytes), melastomes (Melastomataceae) and palms (Arecaceae) in four line transects of 22–29 km length. Different grains were obtained by aggregating original 100‐m‐long sampling units into larger segments up to 19.2 km long. With any given grain and plant group, local species richness varied considerably both within and among transects, and a transect segment that was species‐rich with one grain could be relatively species‐poor with another. Which transect had the highest vs lowest mean species richness per sampling unit (α richness) differed among plant groups. It also varied to some degree with grain, as transects differed in how rapidly local species richness increased with increasing grain. Patterns of species turnover were more consistently correlated among plant groups than patterns of species richness were, and NMDS ordinations were rather similar with all grains and plant groups. Floristic heterogeneity within the Amazonian terra firme rainforest seems to contain a general compositional pattern that is sufficiently robust to be detectable with various sampling schemes, but patterns of species richness appear more case‐specific. Therefore, using one plant group as an indicator for patterns in other plant groups can be expected to work better for species turnover than for species richness.     

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.     

Effects of grazing intensity in grasslands of the Espinal of central Chile

Question: How is grazing intensity associated with species and morpho‐functional traits (MFTs) composition, productivity and richness of annual dominated grasslands? Have native and exotic species similar associations to this gradient? Location: Anthropogenic grassland in the Espinal vegetation in the sub‐humid area of the mediterranean type climate region of Chile (35°58’ S, 72°17’ W). Methods: Data were obtained from a long‐term (eight years) experiment with six stocking rates (1 to 3.5 sheep/ha). Detrended Correspondence Analysis (DCA) and regression analysis were used to determinate the relationship between grazing intensity and biomass, richness, abundance and traits of the species. Results: The first DCA axis was related to grazing intensity and explained most of the floristic variation (69.3%); the abundance of some non‐native species, e.g. Vulpia megalura were highly correlated with this axis. In the DCA for MFTs the first axis explained 87% of the variance and was also related to grazing intensity; the abundance of small size plants and shallow roots increased with grazing intensity. The relative abundance of grasses and composites, but not of legumes, changed with stocking rate: as grazing intensity increased composites became the predominant species to the detriment of grasses. The above‐ground biomass measured in exclusion cages declined with increasing grazing pressure. The richness of exotic species was greater compared to native ones at low stocking rates, but they converge to similar values at higher stocking rates. However, the relative abundance of exotic species was greater than 75% in all stocking rates. Conclusions: Grazing intensification has large effects in the structure of grassland in central Chile. With grazing intensities greater than 1 sheep/ha species characteristics change; evolving in a few years (6–8) towards a similar community regardless of the stocking rate. The overgrazed community has more native than exotic species richness, possibly due to greater defence traits against herbivory of this group of species.     

Impact of shrub canopies on understorey vegetation in western Eurasian tundra

Question: How does the composition and species richness of understorey vegetation associate with changing abundance of deciduous shrub canopies? What are the species‐specific associations between shrubs and understorey plants? Location: Tundra habitats along an over 1000‐km long range, spanning from NW Fennoscandia to the Yamal Peninsula in northwest Russia. Methods: The data from 758 vegetation sample plots from 12 sites comprised cover estimates of all plant species, including bryophytes and lichens, and canopy height of deciduous shrubs. The relationships between shrub volume and cover of plant groups and species richness of vegetation were investigated. In addition, species‐specific associations between understorey species and shrub volume were analysed. Results: Shrub abundance was shown to be associated with the composition of understorey vegetation, and the association patterns were consistent across the study sites. Increased forb cover was positively associated with shrub volume, whereas bryophyte, lichen, dwarf shrub and graminoid cover decreased in association with increasing volume of deciduous shrubs. The total species richness of vegetation declined with increasing shrub volume. Conclusions: The results suggest that an increase of shrubs – due to climatic warming or a decrease in grazing pressure – is likely to have strong effects on plant–plant interactions and lead to a decrease in the diversity of understorey vegetation.     

A field test of the stress‐gradient hypothesis along an aridity gradient

Aims: The stress‐gradient hypothesis (SGH) predicts how plant interactions change along environmental stress gradients. We tested the SGH in an aridity gradient, where support for the hypothesis and the specific shape of its response curve is controversial. Location: Almería, Cáceres and Coimbra, three sites in the Iberian Peninsula that encompass the most arid and wet habitats in the distribution range of a nurse shrub species –Retama sphaerocarpa L. (Boiss) – in Europe. Methods: We analysed the effect of Retama on its understorey plant community and measured the biomass and species richness beneath Retama and in gaps. We estimated the frequency (changes in species richness), importance and intensity of the Retama effects, and derived the severity–interaction relationship pattern, analysing how these interaction indices changed along this aridity gradient. Results and conclusions: The intensity and frequency of facilitation by Retama increased monotonically with increasing environmental severity, and the importance tended to have a similar pattern, overall supporting the SGH. Our data did not support predictions from recent revisions of the SGH, which may not apply to whole plant communities like those studied here or when interactions are highly asymmetrical. Facilitation by Retama influenced community composition and species richness to the point that a significant fraction of species found at the most arid end of the gradient were only able to survive beneath the nurse shrub, whereas some of these species were able to thrive in gaps at more mesic sites, highlighting the ecological relevance of facilitation by nurse species in mediterranean environments, especially in the driest sites.     

Associations among spring-dependent aquatic assemblages and environmental and land use gradients in a Mojave Desert mountain range

In many arid landscapes, springs provide the only reliable source of water. Accordingly, both native species and human land uses, including diversion of water, livestock grazing, and recreation, tend to concentrate around springs and spring‐fed riparian areas. We examined whether species richness and composition of aquatic macroinvertebrates at 45 springs in the Spring Mountains, an isolated mountain range in the eastern Mojave Desert (Nevada, USA), could be predicted using readily measured environmental gradients and estimates of disturbance intensity. The Spring Mountains is a focus of regional conservation planning, and managers are charged with prioritizing its springs for conservation and rehabilitation. Our results suggested that species richness of aquatic macroinvertebrates in the Spring Mountains system may be greatest at intermediate levels of natural and human disturbance. Discharge and springbrook length appeared to be only weakly correlated with species richness, whereas neither elevation, nor water temperature, nor electrical conductance was significantly associated with species richness. Nestedness analyses demonstrated that species present in relatively depauperate locations tended to be subsets of the species present in locations that were richer in species, but that pattern did not appear to be driven by either disturbance intensity or by the environmental variables we measured. Disturbance intensity was not associated with the extent to which species presences and absences were predictable. Although our results should not be interpreted to mean that major environmental gradients and disturbance intensity have no effect on distributional patterns of aquatic invertebrates in the Spring Mountains, the ability of these variables to serve as predictors of species richness and composition may be relatively low. Springs and other wetlands in arid landscapes are characterized by isolation and unpredictable disturbances, and faunal responses to environmental gradients may tend to be individualistic and taxon‐specific.