Habitat heterogeneity,temperature, and primary productivity drive elevational gradients in avian species diversity

AimAnticipating and mitigating the impacts of climate change on species diversity in montane ecosystems requires a mechanistic understanding of drivers of current patterns of diversity. We documented the shape of elevational gradients in avian species richness in North America and tested a suite of a priori predictions for each of five mechanistic hypotheses to explain those patterns.LocationUnited StatesMethodsWe used predicted occupancy maps generated from species distribution models for each of 646 breeding birds to document elevational patterns in avian species richness across the six largest U.S. mountain ranges. We used spatially explicit biotic and abiotic data to test five mechanistic hypotheses proposed to explain geographic variation in species richness.ResultsElevational gradients in avian species richness followed a consistent pattern of low elevation plateau‐mid‐elevation peak (as per McCain, 2009). We found support for three of the five hypotheses to explain the underlying cause of this pattern: the habitat heterogeneity, temperature, and primary productivity hypotheses.Main ConclusionsSpecies richness typically decreases with elevation, but the primary cause and precise shape of the relationship remain topics of debate. We used a novel approach to study the richness‐elevation relationship and our results are unique in that they show a consistent relationship between species richness and elevation among 6 mountain ranges, and universal support for three hypotheses proposed to explain the underlying cause of the observed relationship. Taken together, these results suggest that elevational variation in food availability may be the ecological process that best explains elevational gradients in avian species richness in North America. Although much attention has focused on the role of abiotic factors, particularly temperature, in limiting species’ ranges, our results offer compelling evidence that other processes also influence (and may better explain) elevational gradients in species richness.     

Global species–energy relationship in forest plots: role of abundance,temperature and species climatic tolerances

Aim To evaluate the strength of evidence for hypotheses explaining the relationship between climate and species richness in forest plots. We focused on the effect of energy availability which has been hypothesized to influence species richness: (1) via the effect of productivity on the total number of individuals (the more individuals hypothesis, MIH); (2) through the effect of temperature on metabolic rate (metabolic theory of biodiversity, MTB); or (3) by imposing climatic limits on species distributions. Location Global. Methods We utilized a unique ‘Gentry‐style’ 370 forest plots data set comprising tree counts and individual stem measurements, covering tropical and temperate forests across all six forested continents. We analysed variation in plot species richness and species richness controlled for the number of individuals by using rarefaction. Ordinary least squares (OLS) regression and spatial regressions were used to explore the relative performance of different sets of environmental variables. Results Species richness patterns do not differ whether we use raw number of species or number of species controlled for number of individuals, indicating that number of individuals is not the proximate driver of species richness. Productivity‐related variables (actual evapotranspiration, net primary productivity, normalized difference vegetation index) perform relatively poorly as correlates of tree species richness. The best predictors of species richness consistently include the minimum temperature and precipitation values together with the annual means of these variables. Main conclusion Across the world's forests there is no evidence to support the MIH, and a very limited evidence for a prominent role of productivity as a driver of species richness patterns. The role of temperature is much more important, although this effect is more complex than originally assumed by the MTB. Variation in forest plot diversity appears to be mostly affected by variation in the minimum climatic values. This is consistent with the ‘climatic tolerance hypothesis’ that climatic extremes have acted as a strong constraint on species distribution and diversity.     

Microhabitat variations and seed bank-vegetation relationships in a desert wadi ecosystem

Due to extreme variability in patterns of rainfall, plant seed banks are an important component of desert habitats. Here I report on effects of standing vegetation and three different microhabitats (channel, bank and terrace) on the soil seed bank of a desert wadi ecosystem in the Eastern Desert of Egypt. A total of 450 soil samples at 45 stands were collected to represent the different wadi microhabitats. The germinable seed bank was estimated by controlled counts of seedling emergence. The floristic composition, functional properties and diversity of the soil seed bank, as well as its similarity with the standing vegetation varied among wadi microhabitats. Such variation could be attributed to differences in disturbance intensity among microhabitats (terrace < bank < channel) and variation of soil factors along the microtopographic gradient. Channel showed the highest species richness and size of soil seed bank, followed by bank and then terrace. Moreover the Shannon index of diversity of the seed bank and its similarity with standing vegetation were significantly greater in both channel and bank microhabitats than in terrace. At the level of plant functional groups, number of seeds of annuals was higher in both channel and bank than in terrace. Shrubs were more abundant in seed banks of channel compared to terrace. The size and species richness of seed bank were increased with the total plant cover, annual/perennial ratio and species richness of the standing vegetation.     

Testing the effectiveness of surrogates for assessing biological diversity of arthropods in cereal agricultural landscapes

Agricultural intensification is altering biodiversity patterns worldwide. Rapid and effective methods are needed to monitor these changes in farmland biodiversity, but it becomes both a cost- and time-prohibitive task, particularly for hyper-diverse groups such as arthropods. We evaluated the effectiveness of surrogates in irrigated and rainfed wheat fields in a Mediterranean farmland in NW Spain in order to get a rapid tool to assess arthropod biodiversity. We studied six groups with different ecological needs (i.e. Aphididae, Aphidiinae, Coccinellidae, Formicidae, Heteroptera and Syrphidae) at species level (147 species), genus (105), family (10, only Heteroptera) and order (19) level. Higher taxa, cross-taxa and subset-taxa or total richness approaches were tested as well as the correlation in composition between levels for the selected groups, and the influence of farming regime. Genus richness was a good surrogate of species richness in all six groups studied (R² = 0.38–0.60), like family and order were for Heteroptera (R² = 0.37 and 0.29, respectively). Cross-taxa analyses showed that Aphididae and Aphidiinae genera (R² = 0.19 and 0.30, respectively) and species (R² = 0.20 and 0.28, respectively) were good surrogates for Aphidiinae and Aphididae species respectively. Coccinellidae genera (R² = 0.26) and species (R² = 0.25) were good surrogates for Heteroptera species. Finally, Aphididae and Coccinellidae both at genera (R² = 0.14 and 0.20, respectively) and at species levels (R² = 0.12–0.22, respectively) were good surrogates for total species richness of all groups. Genera composition was the best surrogate for the species composition within each group. Farming regime had no influence on the relationships between surrogates and species patterns in most cases. Our results suggest that genera level is a useful surrogate for all the studied groups and family is appropriate for Heteroptera. Genus level provided a saving of 15% of identification time in Aphididae and 80% for Coccinellidae. This proves its usefulness to asses and monitor biodiversity in wheat croplands and the possibility to reduce costs.     

Rodent endemism,turnover and biogeographical transitions on elevation gradients in the northwestern Argentinian Andes

The aim of this research is to relate patterns of endemism and turnover along a local elevation gradient in northwest Argentina with continental biogeographical transitions. Specimen based records constituted the principal source of information to infer rodent distribution along the elevation gradient. I assessed elevational variation of richness, endemism and turnover by means of non-linear regression analysis. Then I identified five distributional patterns based on the overlap of species geographic range. Their frequency along elevation was used to validate biogeographical boundaries inferred by turnover rates. Eleven species out of 37 (30%) are endemic to the study area. Species richness and endemism were hump-shaped. The rate of endemism reached its maximum value at the upper limit of the forest (2500 m). By contrast, species turnover was U-shaped, with a small peak at 1500 m and a maximum at 3500 m. The species’ geographic range patterns were not randomly distributed along elevation but agglomerated at specific elevation. Species turnover and chorological analysis suggest two biogeographical boundaries, a weaker at 1500 m and a stronger at 3500 m. The 1500 m boundary marks the transition from assemblages dominated by Lowland-widespread fauna at lower elevation to Montane (Andean eastern slopes) species at middle elevation. This boundary is characterized by moderate species turnover and high species richness. The strong turnover rate at 3500 and the dominance of highland Andean and Andean-Patagonian species above this elevation suggest the occurrence of the transition between the Neotropical and Andean regions; which is characterised by an almost complete species replacement.     

Effects of plant diversity on microbial biomass and community metabolic profiles in a full-scale constructed wetland

There has been less understanding of relations of microbial community patterns with plant diversity in constructed wetlands. We conducted a single full-scale subsurface vertical flow constructed wetland (SVFCW, 1000 m²) study focusing on domestic wastewater processing. This study measured the size and structure of microbial community using fumigation extraction and BIOLOG Ecoplate? techniques, to examine the effects of macrophyte diversity on microbial communities that are critical in treatment efficiency of constructed wetlands. We also determined the relationship of plant diversity (species richness) with its biomass production under disturbance of the same wastewater supply. Linear regression analysis showed that plant biomass production strongly correlated with plant species richness (R = 0.407, P < 0.001). Increase in plant species richness increased microbial biomass carbon and nitrogen (R = 0.494, P < 0.001; R = 0.465, P < 0.001) and utilization of amino acids on Ecoplates (R = 0.235, P = 0.03), but limited the utilization of amine/amides (R = ?0.338, P = 0.013). Principal components analysis (PCA) showed that the diversity and community-level physiological profiles (CLPP) of microbial community at 168 h of incubation strongly depended on the presence or absence of plant species in the SVFCW system, but not on the species richness. This is the first step toward understanding relations of plant diversity with soil microbial community patterns in constructed wetlands, but the effect of species diversity on microbial community should be further studied.     

Environmental drivers of ant species richness and composition across the Argentine Pampas grassland

Understanding the underlying mechanisms causing diversity patterns is a fundamental objective in ecology and science‐based conservation biology. Energy and environmental‐heterogeneity hypotheses have been suggested to explain spatial changes in ant diversity. However, the relative roles of each one in determining alpha and beta diversity patterns remain elusive. We investigated the main factors driving spatial changes in ant (Hymenoptera, Formicidae) species richness and composition (including turnover and nestedness components) along a 500 km longitudinal gradient in the Pampean region of Argentina. Ants were sampled using pitfall traps in 12 sample sites during the summer. We performed a model selection approach to analyse responses of ant richness and composition dissimilarity to environmental factors. Then, we computed a dissimilarity partitioning of the contributions of spatial turnover and nestedness to total composition dissimilarity. Temporal habitat heterogeneity and temperature were the primary factors explaining spatial patterns of epigean ant species richness across the Pampas. The distance decay in species composition similarity was best accounted by temperature dissimilarity, and turnover had the greatest contribution to the observed beta diversity pattern. Our findings suggest that both energy and environmental‐heterogeneity‐related variables are key factors shaping richness patterns of ants and niche‐based processes instead of neutral processes appear to be regulating species composition of ant assemblages. The major contribution of turnover to the beta diversity pattern indicated that lands for potential reconversion to grassland should represent the complete environmental gradient of the Pampean region, instead of prioritizing a single site with high species richness.     

Modelling potential impacts of climate change on the bioclimatic envelope and conservation of the Maned Wolf (Chrysocyon brachyurus)

Forecasting the influence of climatic changes on the distribution of the Maned Wolf (Chrysocyon brachyurus) is important for the conservation of the species. We explored the environmental characteristics than best explain the current distribution of the species, modelled the past and present distribution, projected the niche model into the future, and identified suitable areas for conservation. Niche modelling was performed using Maxent and 21 environmental variables. For past conditions, we considered the Last Glacial Maximum (LGM) and the mid-Holocene (MH) climates. For future conditions, we used the A2a greenhouse gas emission scenario for 2050. Four General Circulation Models (FGOALS 1.0, HADCM3, IPSL-CM4 and MIROC 3.2) were used. The resulting niche model (AUC = 0.89 ± 0.02) predicts maximum probability of presence at precipitation of 106 mm during the coldest quarter, of 396 mm during the warmest quarter, and in totally flat areas. The suitable area for the Maned Wolf currently covers 4,320,364 km². For the LGM, there were inter-model differences in predicted areas (from 819,324 km² to 6,395,886 km²) and in geographic location. The MH models showed drastic changes with respect to the present and considerable inter-model variation. Predictions for 2050 show significant (at least 33%) reductions in distribution. Only a minor fraction (39%) of the current distribution can be considered stable for the period LGM-2050. The FGOALS model was the best option for projecting species occurrence into the future because it included the three localities known for the Maned Wolf from the late Pleistocene and predicts stable areas that coincide with spatial patterns of genetic diversity. The FGOALS projection for 2050 predicts a 33% reduction in suitable habitats, indicating some stable areas (central South America) that will probably be key sites for the conservation of the species.     

Anuran species richness,complementarity and conservation conflicts in Brazilian Cerrado

Broad-scale correlations between species richness and human population suggest that processes driving species richness, mainly related to high ecological productivity, may also drive human populations. However, it is still under debate if this coincidence implies conflicts between biodiversity conservation and human development. In this paper, we analyzed the relationships among human population size, species richness and irreplaceability in Brazilian Cerrado. We analyzed a dataset with 131 species of anurans distributed in 181 cells with 1° of spatial resolution covering the biome. We found a positive correlation between human population size and anuran species richness (r = 0.46; P = 0.033 with 19.5 geographically effective degrees of freedom, v*), but the irreplaceability of each cell was poorly correlated with human population size (r = 0.075; P = 0.323; v* = 173.9). The 17 cells in the 97 optimal reserve networks contained a total human population ranging from 2942,195 to 4319,845 people, representing on average 11.8% of the human population in the entire Cerrado grid. The comparison of these observed values with 10,000 values from randomly generated networks suggests a relatively high flexibility in optimal complementarity sets for reserve selection. Our results indicated that correlation between richness and human population does not necessarily result in conflicts, given the opportunities for conciliating conservation and development. However, the analyses performed here are initial explorations within the framework of conservation biogeography, so more detailed studies are necessary to establish conservation planning at regional and local scales.     

A test of multiple hypotheses for the species richness gradient of South American owls

Many mechanisms have been proposed to explain broad scale spatial patterns in species richness. In this paper, we evaluate five explanations for geographic gradients in species richness, using South American owls as a model. We compared the explanatory power of contemporary climate, landcover diversity, spatial climatic heterogeneity, evolutionary history, and area. An important aspect of our analyses is that very different hypotheses, such as history and area, can be quantified at the same observation scale and, consequently can be incorporated into a single analytical framework. Both area effects and owl phylogenetic history were poorly associated with richness, whereas contemporary climate, climatic heterogeneity at the mesoscale and landcover diversity explained ca. 53% of the variation in species richness. We conclude that both climate and environmental heterogeneity should be retained as plausible explanations for the diversity gradient. Turnover rates and scaling effects, on the other hand, although perhaps useful for detecting faunal changes and beta diversity at local and regional scales, are not strong explanations for the owl diversity gradient.     

Species richness,equitability, and abundance of ants in disturbed landscapes

Ants are used as indicators of environmental change in disturbed landscapes, often without adequate understanding of their response to disturbance. Ant communities in the southeastern United States displayed a hump-backed species richness curve against an index of landscape disturbance. Forty sites at Fort Benning, in west-central Georgia, covered a spectrum of habitat disturbance (military training and fire) in upland forest. Sites disturbed by military training had fewer trees, less canopy cover, more bare ground, and warmer, more compact soils with shallower A-horizons. We sampled ground-dwelling ants with pitfall traps, and measured 15 habitat variables related to vegetation and soil. Ant species richness was greatest with a relative disturbance of 43%, but equitability was greatest with no disturbance. Ant abundance was greatest with a relative disturbance of 85%. High species richness at intermediate disturbance was associated with greater within-site spatial heterogeneity. Species richness was also associated with intermediate values of the normalized difference vegetation index (NDVI), a correlate of net primary productivity (NPP). Available NPP (the product of NDVI and the fraction of days that soil temperature exceeded 25 °C), however, was positively correlated with species richness, though not with ant abundance. Species richness was unrelated to soil texture, total ground cover, and fire frequency. Ant species richness and equitability are potential state indicators of the soil arthropod community. Moreover, equitability can be used to monitor ecosystem change.     

Soundscape analysis and acoustic monitoring document impacts of natural gas exploration on biodiversity in a tropical forest

Natural resource extraction is increasing rapidly in tropical forests, but we lag behind in understanding the impacts of these disturbances on biodiversity. In high diversity tropical habitats, acoustic monitoring is an efficient tool for sampling a large proportion of the fauna across varied spatial and temporal scales. We used passive acoustic monitoring in a pre-montane forest in Peru to investigate how soundscape composition and richness of acoustic frequencies varied with distance from a natural gas exploratory well and with operational phase (construction and drilling). We also evaluated how anuran and avian species richness and vocal activity varied with distance and between phases. Soundscape analyses showed that acoustic frequency similarity was greatest among sites closer to (≤250 m) and farther from (≥500 m) the platform. Soundscapes revealed more frequencies were used during construction and showed a weak trend of increasing frequency richness with increasing distance from the disturbance. Avian species richness and detections increased with distance from the platform, but anuran richness and detections declined with distance. Operational phase did not play a significant role in overall richness or activity patterns of either group. Among birds, insectivore detections increased with distance from the platform, and nectarivores were detected more frequently during the drilling phase. Results demonstrate that acoustic monitoring and soundscape analyses are useful tools for evaluating the impact of development activity on the vocalizing community, and should be implemented as a best practice in monitoring biodiversity and for guiding specific mitigation strategies.     

An exploration of factors influencing lotic insect species richness

An understanding of factors that influence species richness of lotic insects is generally lacking. We present comparative data on aquatic insect species richness from several North American and other streams. Factors such as large sample numbers and drainage area (species area relationships) are not significant predictors of species richness across the streams we examined. We explore several hypotheses regarding the origins and maintenance of species richness using Upper Three Runs Creek (UTR), South Carolina, USA, as a reference stream. UTR has the highest species richness of any stream in the Western Hemisphere. Hypotheses examined included historical, regional and local processes such as: (1) Evolutionary time, (2) disturbance regime/environmental variability, (3) temperature/evolutionary-speed, (4) productivity, and (5) habitat heterogeneity. Of these hypotheses, we suggest that productivity and habitat heterogeneity appear to contribute most to the high species richness found in UTR. We believe that multidisciplinary analysis of other streams is necessary because without this crucial information our knowledge of, and desire to protect biodiversity in streams will be wanting.     

A habitat quality indicator for common birds in Europe based on species distribution models

The EU 2020 Biodiversity Strategy requires the gathering of information on biodiversity to aid in monitoring progress towards its main targets. Common species are good proxies for the diversity and integrity of ecosystems, since they are key elements of the biomass, structure, functioning of ecosystems, and therefore of the supply of ecosystem services. In this sense, we aimed to develop a spatially-explicit indicator of habitat quality (HQI) at European level based on the species included in the European Common Bird Index, also grouped into their major habitat types (farmland and forest). Using species occurrences from the European Breeding Birds Atlas (at 50 km × 50 km) and the maximum entropy algorithm, we derived species distribution maps using refined occurrence data based on species ecology. This allowed us to cope with the limitations arising from modelling common and widespread species, obtaining habitat suitability maps for each species at finer spatial resolution (10 km × 10 km grid), which provided higher model accuracy. Analysis of the spatial patterns of local and relative species richness (defined as the ratio between species richness in a given location and the average richness in the regional context) for the common birds analysed demonstrated that the development of a HQI based on species richness needs to account for the regional species pool in order to make objective comparisons between regions. In this way, we proved that relative species richness compensated for the bias caused by the inherent heterogeneous patterns of the species distributions that was yielding larger local species richness in areas where most of the target species have the core of their distribution range. The method presented in this study provides a robust and innovative indicator of habitat quality which can be used to make comparisons between regions at the European scale, and therefore potentially applied to measure progress towards the EU Biodiversity Strategy targets. Finally, since species distribution models are based on breeding birds, the HQI can be also interpreted as a measure of the capacity of ecosystems to provide and maintain nursery/reproductive habitats for terrestrial species, a key maintenance and regulation ecosystem service.     

Mapping multiple plant species abundance patterns - A multiobjective optimization procedure for combining reflectance spectroscopy and species ordination

Nature conservation and ecological restoration crucially depends on the knowledge about spatial patterns of plant species that control habitat conversion and disturbance regimes. Especially, species abundances are capable of indicating early development tendencies for setting habitat management strategies. This study demonstrates the transfer of field spectroscopy to hyperspectral imagery to map multiple plant species abundances in an open dryland area using two imaging spectrometers in two different phenological phases. We show that species abundances can partially be described by multiple gradients forming different coordinates in a contour map. For this purpose, species abundances were projected into an ordination space using non-metric multidimensional scaling and subsequent spatial interpolation. It was demonstrated that different gradients can be modeled in a Partial Least Squares regression framework resulting in distinct spectral features for certain gradient directions. We combine both objectives in a multiobjective NSGA-II procedure to maximize the quantitative determination of species abundance in ordination and spectral predictability in related field spectra, simultaneously. NSGA-II was finally used to select optimal spectral models for n = 35 single species that were transferred to hyperspectral imagery for mapping purpose. We can show that abundance predictabilities can be evaluated on the basis of individual model performances that hold different spectral features for each species in a designated phenological phase. Finally, we present spatially explicit multi-species maps for the best n = 18 and abundance maps for n = 8 models that could be linked to patterns of species richness, coexistence, succession stages and habitat type conditions.     

Environmental determinants of amphibian and reptile species richness in China

Understanding the factors that regulate geographical variation in species richness has been one of the fundamental questions in ecology for decades, but our knowledge of the cause of geographical variation in species richness remains poor. This is particularly true for herpetofaunas (including amphibians and reptiles). Here, using correlation and regression analyses, we examine the relationship of herpetofaunal species richness in 245 localities across China with 30 environmental factors, which include nearly all major environmental factors that are considered to explain broad-scale species richness gradients in such theories as ambient energy, water–energy dynamics, productivity, habitat heterogeneity, and climatic stability. We found that the species richness of amphibians and reptiles is moderately to strongly correlated with most of the environmental variables examined, and that the best fit models, which include explanatory variables of temperature, precipitation, net primary productivity, minimum elevation, and range in elevation, explain ca 70% the variance in species richness for both amphibians and reptiles after accounting for sample area. Although water and temperature are important explanatory variables to both amphibians and reptiles, water variables explain more variance in amphibian species richness than in reptile species richness whereas temperature variables explain more variance in reptile species richness than in amphibian species richness, which is consistent with different physiological requirements of the two groups of organisms.     

Phreatophytic vegetation response to climatic and abstraction-induced groundwater drawdown: Examples of long-term spatial and temporal variability in community response

The influence of climatic drought and groundwater abstraction on phreatophytic vegetation dynamics was investigated in the southwest of Western Australia. Two contrasting examples of long-term phreatophytic plant community response to reduced water availability are presented. Multivariate analysis of vegetation and hydrological parameters determined depth to watertable as the dominant biophysical driver of floristic spatial and temporal patterns. Under lower rates of watertable drawdown (9 cm year^?1), a progressive change in floristic composition was observed over a 33-year period. The abundance of species with a preference for wetter sites was significantly reduced, whereas that of more drought-tolerant species increased. Higher rates of drawdown (50 cm year^?1) where groundwater abstraction exacerbated climatic drought resulted in a threshold response in vegetation and 33% dissimilarity to pre-abstraction floristics in 12 years. In the context of an ecohydrological state and transition conceptual model, it is suggested higher rates of groundwater drawdown result in a threshold breach and subsequent transition to an alternative ecohydrological state, whilst lower rates result in a progressive floristic transition.     

Diatom biomonitoring of streams: Reliability of reference sites and the response of metrics to environmental variations across temporal scales

Benthic diatoms are widely used indicators of human impacts on stream ecosystems because they are very responsive to changing environmental conditions. However, little research has explicitly focused on their reliability with regards to temporal variation in assemblage structure and environmental conditions. We examined variability in diatom-environment relationships at bi-weekly, monthly, and yearly time scales from 7 reference, 7 agricultural, and 2 acid mine drainage (AMD)-impacted streams, and how nutrient and pH fluctuations may affect the interpretation of diatom metrics and the Diatom Model Affinity (DMA) index. Reference streams had less bi-weekly variability in NO₃-N concentrations than non-reference streams. The % eutraphentic diatoms and DMA scores were more strongly correlated with seasonal means of NO₃-N and PO₄-P concentrations than with same day concentrations. Most nutrient indicator metrics had strong correlations with watershed land use. All 14 non-AMD streams experienced substantial increases in NO₃-N and decreases in temperature from November to May, which were associated with high species turnover, substantial changes in community structure, reduced diversity and richness, increased relative abundances of high nutrient diatoms, and decreases in low nutrient diatoms and DMA scores. The % acidophilic diatoms and DMA scores were significantly correlated with increased pH associated with greater precipitation at AMD sites from December to April (r = ?0.77, r = 0.62, respectively; P < 0.01). Yearly, DMA scores for all reference streams were consistently in the minimally impaired category, whereas scores for non-reference streams varied among impairment categories. Reference sites serve as reliable benchmarks for diatom ecological integrity during the summer. In this region, June to October is a recommended time period for diatom sampling in monitoring programs because subsequent shifts in hydrologic regimes, nutrients, and diatom assemblages occurred, affecting all sites and masking among stream differences attributable to agricultural land uses.     

Squamate richness in the Brazilian Cerrado and its environmental–climatic associations

We investigate patterns of species richness of squamates (lizards, snakes, and amphisbaenians) in the Brazilian Cerrado, identifying areas of particularly high richness, and testing predictions of large‐scale richness hypotheses by analysing the relationship between species richness and environmental climatic variables. We used point localities from museum collections to produce maps of the predicted distributions for 237 Cerrado squamate species, using niche‐modelling techniques. We superimposed distributions of all species on a composite map, depicting richness across the ecosystem. Then, we performed a multiple regression analysis using eigenvector‐based spatial filtering (Principal Coordinate of Neighbour Matrices) to assess environmental–climatic variables that are best predictors of species richness. We found that the environmental–climatic and spatial filters multiple regression model explained 78% of the variation in Cerrado squamate richness (r² = 0.78; F = 32.66; P < 0.01). Best predictors of species richness were: annual precipitation, precipitation seasonality, altitude, net primary productivity, and precipitation during the driest quarter. A model selection approach revealed that several mechanisms related to the different diversity hypothesis might work together to explain richness variation in the Cerrado. Areas of higher species richness in Cerrado were located mainly in the south‐west, north, extreme east, and scattered areas in the north‐west portions of the biome. Partitioning of energy among species, habitat differentiation, and tolerance to variable environments may be the primary ecological factors determining variation in squamate richness across the Cerrado. High richness areas in northern Cerrado, predicted by our models, are still poorly sampled, and biological surveys are warranted in that region. The south‐western region of the Cerrado exhibits high species richness and is also undergoing high levels of deforestation. Therefore, maintenance of existing reserves, establishment of ecological corridors among reserves, and creation of new reserves are urgently needed to ensure conservation of species in these areas.