The distributions of amphibians, birds and mammals have underpinned global and local conservation priorities, and have been fundamental to our understanding of the determinants of global biodiversity. In contrast, the global distributions of reptiles, representing a third of terrestrial vertebrate diversity, have been unavailable. This prevented the incorporation of reptiles into conservation planning and biased our understanding of the underlying processes governing global vertebrate biodiversity. Here, we present and analyse the global distribution of 10,064 reptile species (99% of extant terrestrial species). We show that richness patterns of the other three tetrapod classes are good spatial surrogates for species richness of all reptiles combined and of snakes, but characterize diversity patterns of lizards and turtles poorly. Hotspots of total and endemic lizard richness overlap very little with those of other taxa. Moreover, existing protected areas, sites of biodiversity significance and global conservation schemes represent birds and mammals better than reptiles. We show that additional conservation actions are needed to effectively protect reptiles, particularly lizards and turtles. Adding reptile knowledge to a global complementarity conservation priority scheme identifies many locations that consequently become important. Notably, investing resources in some of the world’s arid, grassland and savannah habitats might be necessary to represent all terrestrial vertebrates efficiently.
Our knowledge of the distributions of a broad variety of organisms has improved greatly in the past decade1,2,3, in turn aiding our efforts to conserve biodiversity4,5,6 and significantly enhance our grasp of broad-scale evolutionary and ecological processes7,8,9,10,11,12. Nevertheless, despite comprising one third of terrestrial vertebrate species, knowledge of reptile distributions remained poor and unsystematic. This represented a major gap in our understanding of the global structure of biodiversity and our ability to conserve nature. Historically, broad-scale efforts towards the protection of land vertebrates (and thus also of reptiles) have been based predominantly on data from plants, birds, mammals and to a lesser degree amphibians13,14,15. Here we present complete species-level global distributions of nearly all reptiles: 10,064 known, extant, terrestrial species for which we could identify precise distribution information. These distributions cover the Sauria (lizards, 6,110 species), Serpentes (snakes, 3,414 species), Testudines (turtles, 322 species), Amphisbaenia (‘worm lizards’, 193 species), Crocodylia (crocodiles, 24 species) and Rhynchocephalia (the tuatara, 1 species).
This dataset completes the global distribution mapping of all described, extant, terrestrial vertebrates (Fig. 1a), providing information that has been missing from much of the global conservation planning and prioritization schemes constructed over the past twenty years4. We use our reptile distribution data to: (a) examine the congruence in general, hotspot and endemism richness patterns across all tetrapod classes and among reptile groups; (b) explore how current conservation networks and priorities represent reptiles; and (c) suggest regions in need of additional conservation attention to target full terrestrial vertebrate representation and highlight current surrogacy gaps, using a formal conservation prioritization technique.
a, Richness of all tetrapods (reptiles, amphibians, birds and mammals). b–e, Species richness of reptile groups: all reptiles (b), lizards (c), snakes (d) and turtles (e). Grey areas denote terrestrial regions devoid of species in a particular group. Blue colours denote regions with few species and red ones denote regions with many species (note that the scale differs between panels). All maps in an equal area, Behrmann projection at a 48.25 × 48.25?km grid-cell resolution), scale bar values represent species richness.
The global pattern of reptile species richness (Fig. 1b) is largely congruent with that of all other terrestrial vertebrates combined (r?=?0.824, e.d.f.?=?31.2, p?«?0.0001; Fig. 2a, Supplementary Fig. 1, Supplementary Table 1). However, the major reptile groups (Figs. 1c–e, 2b,c Supplementary Fig. 1, Supplementary Table 1) show differing degrees of congruence with the other tetrapod taxa. The richness distribution of snakes (Fig. 1d) is very similar to that of other tetrapods (Fig. 2c) in showing pan-tropical dominance (r?=?0.873, e.d.f.?=?30.2, p?«?0.0001). Lizard richness is much less similar to non-reptilian tetrapod richness (r?=?0.501, e.d.f.?=?38.3, p?«?0.001, Fig. 2b). It is high in both tropical and arid regions, and notably in Australia (Fig. 1c, Supplementary Fig. 1). Turtle richness is also less congruent with diversity patterns of the other tetrapods (r?=?0.673, e.d.f.?=?55.2, p?«?0.001), and peaks in the southeastern USA, the Ganges Delta, and Southeast Asia (Fig. 1e).
a–c, Hexagon scatter plots comparing species richness values per grid-cell with binning (black line indicates a LOESS fit, a?=?0.6) of tetrapods without reptiles to all reptiles (a), lizards (b) and snakes (c). d, Map of the ratio of reptile richness to non-reptilian tetrapod richness per grid-cell (note the wide range of values for the top category). Hatched regions designate areas where this proportion is the top 5% (black) and 25% (grey).
Snakes dominate reptile richness patterns due to their much larger range sizes compared with lizards, even though lizards are about twice as speciose (median range size for 3,414 snake species: 62,646?km2; for 6,415 lizard species: 11,502?km2; Supplementary Fig. 2). Therefore snakes disproportionally influence global reptile richness patterns16,17 (Supplementary Table 1, Supplementary Fig. 1).
As with overall richness patterns, hotspots of richness (the richest 2.5%, 5%, 7.5% and 10% of grid-cells) for all reptiles combined, and of snakes, are largely congruent with those of other tetrapod classes. However they are incongruent with hotspots of lizard or turtle richness (Fig. 3; Supplementary Fig. 3).
a–d, In each panel the lightest colour denotes the 10% of 48.25 × 48.25?km grid-cells with the highest numbers of species, and as the colours get darker they represent the top 7.5%, 5% and richest 2.5% cells. Panels show all reptiles (a), lizards (b), snakes (c) and turtles (d).
Congruence in the richness of range-restricted species (those species with the smallest 25% or 10% ranges in each group) between tetrapod groups is lower than the congruence across all species1 (Supplementary Table 1). Endemic lizard and turtle distributions are least congruent with the endemics in other tetrapod classes (Supplementary Table 1). Global hotspots of relative endemism (or range-size-weighted richness, see Methods) for reptiles differ from those of non-reptilian tetrapods (Supplementary Fig. 4). Island faunas in places such as Socotra, New Caledonia and the Antilles are highlighted for reptiles, while hotspots of endemism for non-reptilian tetrapods are more often continental.
Reptiles, like amphibians, are poorly represented in the global network of protected areas (Supplementary Table 2; Supplementary Figs. 5,6). Only 3.5% of reptile and 3.4% of amphibian species ranges are contained in protected areas (median species range overlap per class, with IUCN categories I–IV), compared with 6.5% for birds and 6% for mammals. Within reptile groups, strict protected areas (IUCN category I) overlap less with lizard ranges than with other reptile groups but there are no important differences between taxa for the more permissive protected area types (Supplementary Table 2; Supplementary Fig. 5). Amphibians have the highest proportion of species whose ranges lie completely outside protected areas, when compared to the other tetrapod groups. Lizards, also fare poorly and have the highest proportion of species outside protected areas when compared with the other reptile groups (Supplementary Fig. 6a). Turtles have the lowest proportion of species with at least 10% of their range covered by protected areas (Supplementary Fig. 6b). We suggest that these low overlaps may have been caused by the inability to consider reptile diversity for direct protection, probably arising from ignorance of their distributions.
We explored the coverage of all tetrapods in three global prioritization schemes13,14,18 and a global designation of sites for biodiversity significance15 that have recently used distribution data to highlight regions for targeted conservation. These four global prioritizations/designations cover 6.8–37.4% of the Earth’s land surface with 34–11,815 unique sites. Terrestrial vertebrate groups have 68–98% of their species with at least some range covered by these schemes (Supplementary Fig. 6c). However, reptiles and amphibiansare sampled least well by these global schemes, and within reptiles lizards have the lowest representation (Supplementary Fig. 6c).
Fortunately, reptiles seem better situated in terms of conservation costs compared with other tetrapods. The median conservation opportunity cost19 (using the loss of agricultural revenue as a proxy for land-cost) for reptiles is lower than that for other tetrapods (F3, 31850?=?17.4, p?7). Within reptiles, the opportunity cost is lowest for lizards, and highest for turtles and crocodiles, which could reflect their greater dependence on more valuable freshwater habitats (F3, 10060?=?88.4, p?7b).
Our results suggest that reptiles, and particularly lizards and turtles, need to be better incorporated into conservation schemes. We used relative endemism within a complementarity analysis20 to identify broad areas within which international and local conservation action should reduce reptile extinction risk (Fig. 4, Supplementary Fig. 8), and repeated this analysis to also incorporate conservation opportunity costs19 (Supplementary Fig. 8d,e). Many previously identified priority regions13,14, have been retained with the addition of reptile distributions. These include northern and western Australia; central southern USA and the gulf coast of Mexico; the Brazilian Cerrado; Southeast Asia and many islands.
Cells were ranked in a formal prioritization scheme20, based on complementarity when ranking cells in an iterative manner. Cells were ranked twice: (1) with all tetrapods and (2) with all tetrapods excluding reptiles. a, Patterns per 0.5 degree grid-cell where colours represent the priority ranks for the scheme that included all tetrapods (blue?=?low, red?=?high). The cells that are highlighted with the bold foreground colours are those that pinpoint those regions that gain in conservation importance due to the inclusion of the reptile data. These cells were selected following two rules: (1) they were in the top 10% of increase in rank, when subtracting the ranks of the analysis with reptiles from the ranks of the analysis without them; and (2) were part of statistically significant spatial clusters of rank changes (using local Moran’s I 35). b, The mean change in rank between prioritizations with and without reptiles (using the above method) per ecoregion (red, ecoregions that become more important due to the inclusion of reptile information; blue, ecoregions becoming less important).
Nevertheless, our analyses also reveal many new priority regions, not currently perceived as biodiversity conservation priorities for tetrapods. These priority areas are predominantly arid and semi-arid habitats (see also Supplementary Fig. 8f for mean rank change per biome, for prioritization with and without reptiles). They include parts of northern Africa through the Arabian peninsula and the Levant; around Lake Chad; in inland arid southern Africa; central Asian arid highlands and steppes; central Australia; the Brazilian Caatinga, and the southern Andes. These regions have been previously unrecognized as priorities because their non-reptile vertebrate biotas could be more efficiently represented in other locations. Our analyses show that those locations were poor spatial surrogates for reptile distributions and that conservation efforts in our suggested locations may afford better protection for reptiles while maintaining efficient representation of other vertebrates. We note that many of these locations have low conservation opportunity costs so may be especially attractive for conservation. However, the location of these areas is not primarily driven by conservation opportunity costs. When these costs are incorporated into the analyses, very similar regions are highlighted for special attention due to the inclusion of reptile distributions (Supplementary Fig. 8d,e).
The complete map of tetrapod species richness presented here reveals important and unique properties of reptile diversity, particularly of lizards and turtles (Figs. 1–3). At a regional scale, reptiles have previously been shown to be unusually diverse in arid and semi-arid habitats21,22,23. Here we reveal that this pattern is global, and also show reptile prominence in island faunas (Fig. 2d, Supplementary Fig. 4). Furthermore, we show that the unique diversity patterns of reptiles have important implications for their conservation. Targeted reptile conservation lags behind that of other tetrapod classes, probably through ignorance24,25,26. The distributions provided here could make a vital contribution to bridging this gap. Concentrations of rare species in unexpected locations (Fig. 4) require explicit consideration when planning conservation actions. Highlighting such locations for new taxa could be especially beneficial for resource-constrained planning, especially where land costs are low. The lower global congruence with recognized diversity patterns for reptiles should also serve as a warning sign, contrary to some recent suggestions27, for our ability to use distributions of well-studied groups in order to predict diversity patterns of poorly known taxa. The distinctive distribution of reptiles, and especially of lizards, suggests that it is driven by different ecological and evolutionary processes to those in other vertebrate taxa23,28. The complete distributions of terrestrial tetrapods we now possess could greatly enhance our ability to study, understand and protect nature.
Data collection and assembly was carried out by members of the Global Assessment of Reptile Distributions (GARD) group, which includes all the authors of this paper. Regional specialist group members supervised the integration of geographic data for all species from field guides and books covering the terrestrial reptilian fauna of various regions, as well as revised museum specimen databases, online meta-databases (including the IUCN, GBIF and Vertnet), our own observations and the primary literature. We followed the taxonomy of the March 2015 edition of the Reptile Database29. Source maps were split or joined on that basis. We used the newest sources available to us. Polygonal maps—representing species extent of occurrence—were preferred over other map types. One of the reasons for this preference was that the distributional representations for the other tetrapod classes to which reptiles were compared were of a similar nature. Point locality data were modelled to create polygons representing the extent of occurrence using hull geometries (see Supplementary Information). Gaps in reptile distribution knowledge for particular locations or taxa were filled using de novo polygon and gridded maps created by GARD members specializing in the fauna of particular regions and taxa. These maps and all data obtained from online databases and the primary literature were then internally vetted, in a manner analogous to the IUCN Specialist Group process. Further details on data collection and curation, modelling of point localities and a full list of data sources per species are available in the supplement. Overall we analysed distribution maps for 10,064 extant species, which represent 99% of the species found in the Reptile Database of March 2015. For all analytical purposes we contrasted snakes with the paraphyletic lizards (here defined as lepidosaurs exclusive of snakes).
Polygonal representations of the extent of species’ occurrences, such as we assembled and use in our analyses, are fundamentally important to contemporary conservation planning30. The IUCN’s assessment of the extinction risk of individual species requires (and produces) such data, and both they and many other organizations and researchers have used such data in aggregate and at regional-to-global scales for several decades31. Like any representation of species distributions, polygonal range maps can include errors both of omission and commission. Both kinds of inaccuracy can lead to erroneous conclusions by unwary users and this has led to some controversy over the use of polygonal range maps. Of course, all biogeographic representations—specimen localities, SDM outputs, atlas data, polygonal maps and explorers’ narratives—lie along this omission–commission spectrum, and can equally be misused or found useless32. For global prioritization, we follow a comprehensive recent study33 demonstrating the effectiveness of polygonal range maps in highlighting priority areas, despite errors at the level of individual species. We do, however, recognize that specimen data, if collected, curated and made available (at a suitable scale) remain a gold standard for some uses34.
Our grid-cell analyses were conducted in a Behrmann equal area projection of 48.25?km2 grid-cells (~0.5° at 30° N/S). All analyses were repeated at a grid size of 96.5?km2 (~1° at 30° N/S) and results were qualitatively unchanged. GIS and statistical analyses were carried out in R and PostGIS.
Range-size-weighted richness (rswr) was calculated, for each cell, using the following formula: rswr i = ? j q i j where q ij is the fraction of the distribution of the species j in the cell i.
We used ‘zonation’20 to produce a ranked prioritization amongst cells, assuming equal weight to all species and assuming an equal cost for all cells. Cell value was the maximum proportion of any species range represented in it. Cell priority was calculated by iteratively removing the least valuable cell and updating cell values20. We analysed all tetrapod species combined and tetrapods without reptiles separately, to reveal the change in rank importance induced by adding reptile distributions (See Supplementary Fig. 8). We repeated our prioritization using per-cell agricultural opportunity costs19, and found via rank correlation that our priority regions are fairly insensitive to the use of land costs (Fig. 4, Supplementary Fig. 8).
The reptile distribution data used in this study are available from Dryad (https://doi.org/10.5061/dryad.83s7k).
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Grenyer, R. et al. Global distribution and conservation of rare and threatened vertebrates. Nature 444, 93–96 (2006).
Orme, C. D. L. et al. Global hotspots of species richness are not congruent with endemism or threat. Nature 436, 1016–1019 (2005).
Stuart, S. N. et al. Status and trends of amphibian declines and extinctions worldwide. Science 306, 1783–1786 (2004).
Brooks, T. M. et al. Global biodiversity conservation priorities. Science 313, 58–61 (2006).
Kremen, C. et al. Aligning conservation priorities across taxa in Madagascar with high-resolution planning tools. Science 320, 222–226 (2008).
Wilson, K. A., McBride, M. F., Bode, M. & Possingham, H. P. Prioritizing global conservation efforts. Nature 440, 337–340 (2006).
Holt, B. G. et al. An update of Wallace’s zoogeographic regions of the world. Science 339, 74–78 (2013).
Schipper, J. et al. The status of the world’s land and marine mammals: diversity, threat, and knowledge. Science 322, 225–230 (2008).
Bates, S. T. et al. Examining the global distribution of dominant archaeal populations in soil. ISME J. 5, 908–917 (2011).
Morueta-Holme, N. et al. Habitat area and climate stability determine geographical variation in plant species range sizes. Ecol. Lett. 16, 1446–1454 (2013).
Stuart-Smith, R. D. et al. Integrating abundance and functional traits reveals new global hotspots of fish diversity. Nature 501, 539–542 (2013).
Tittensor, D. P. et al. Global patterns and predictors of marine biodiversity across taxa. Nature 466, 1098–1101 (2010).
Mittermeier, R. A. et al. Hotspots Revisited: Earth’s Biologically Richest and Most Endangered Ecoregions (CEMEX, Mexico City, 2004).
Olson, D. M. & Dinerstein, E. The Global 200: a representation approach to conserving the Earth’s most biologically valuable ecoregions. Conserv. Biol. 12, 502–515 (1998).
Important Bird and Biodiversity Area (IBA) Digital Boundaries. Version 2015 2 (BirdLife Intranational, Cambridge, 2015).
Jetz, W. & Rahbek, C. Geographic range size and determinants of avian species richness. Science 297, 1548–1551 (2002).
Lennon, J. J., Koleff, P., Greenwood, J. J. D. & Gaston, K. J. Contribution of rarity and commonness to patterns of species richness. Ecol. Lett. 7, 81–87 (2004).
Joppa, L. N., Visconti, P., Jenkins, C. N. & Pimm, S. L. Achieving the Convention on Biological Diversity’s goals for plant conservation. Science 341, 1100–1103 (2013).
Naidoo, R. & Iwamura, T. Global-scale mapping of economic benefits from agricultural lands: implications for conservation priorities. Biol. Conserv. 140, 40–49 (2007).
Moilanen, A. et al. Prioritizing multiple-use landscapes for conservation: methods for large multi-species planning problems. Proc. R. Soc. Lond. B 272, 1885–1891 (2005).
Pianka, E. R. in Lacertids of the Mediterranean Region (eds Valakos, E. D., Böhme, W., Pérez-Mellado, V. & Maragou, P.) 121–154 (Hellenic Zoological Society, University of Athens, Athens, 1993).
Lewin, A. et al. Patterns of species richness, endemism and environmental gradients of African reptiles. J. Biogeogr. 43, 2380–2390 (2016).
Powney, G. D., Grenyer, R., Orme, C. D. L., Owens, I. P. F. & Meiri, S. Hot, dry and different: Australian lizard richness is unlike that of mammals, amphibians and birds. Glob. Ecol. Biogeogr. 19, 386–396 (2010).
Böhm, M. et al. The conservation status of the world’s reptiles. Biol. Conserv. 157, 372–385 (2013).
Meiri, S. & Chapple, D. G. Biases in the current knowledge of threat status in lizards, and bridging the ‘assessment gap’. Biol. Conserv. 2014A, 6–15 (2016).
Roll, U. et al. Using Wikipedia page views to explore the cultural importance of global reptiles. Biol. Conserv. 204A, 42–50 (2016).
Bode, M. et al. Cost-effective global conservation spending is robust to taxonomic group. Proc. Natl Acad. Sci. USA 105, 6498–6501 (2008).
Hawkins, B. A. et al. Energy, water, and broad-scale geographic patterns of species richness. Ecology 84, 3105–3117 (2003).
Uetz, P. & Hošek, J. The Reptile Database (2015); http://www.reptile-database.org/
Pouzols, F. M. et al. Global protected area expansion is compromised by projected land-use and parochialism. Nature 516, 383–386 (2014).
Scott, J. M. et al. Gap analysis: a geographic approach to protection of biological diversity. Wildlife Monogr. 123, 3–41 (1993).
Maldonado, C. et al. Estimating species diversity and distribution in the era of Big Data: to what extent can we trust public databases? Glob. Ecol. Biogeogr. 24, 973–984 (2015).
Maréchaux, I., Rodrigues, A. S. L. & Charpentier, A. The value of coarse species range maps to inform local biodiversity conservation in a global context. Ecography https://doi.org/10.1111/ecog.02598 (2016).
Cantú-Salazar, L. & Gaston, K. J. Species richness and representation in protected areas of the Western Hemisphere: discrepancies between checklists and range maps. Divers. Distrib. 19, 782–793 (2013).
Anselin, L. Local indicators of spatial association — LISA. Geogr. Anal. 27, 93–115 (1995).
We thank T. Burbidge, T. Dowe, S. Huang, S. Khela, H.-Y. Lee, K. Tamar, J. Usherwood, M. Hopkins and S. Halle for help in digitizing reptile ranges. We thank librarians and colleagues for help in obtaining relevant literature, G. Bunting and M. Balman for providing IBA polygons and bird species distribution maps from BirdLife International, as well as S. Butchart for insightful comments. A.B. thanks the Gerald M. Lemole endowed Chair funds. G.R.C. thanks CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Conselho Nacional de Desenvolvimento Científico e Tecnológico – CNPq and Fundação de Apoio à Pesquisa do Distrito Federal – FAPDF for financial support. I.D. was supported by a Niche Research Grant Scheme, NRGS/1087/2–13(01). C.N. and M.M. were supported by São Paulo Research Foundation (FAPESP no. 2011/50206-9, no. 2012/19858-2 to C.N.). M.M. acknowledges a research fellowship from CNPq. O.T.C. acknowledges support from SENESCYT. R.G. acknowledges the John Fell Fund of the University of Oxford for support. A.A. and S.M. acknowledge support from a BSF grant no. 2012143.
Present address: Mitrani Department of Desert Ecology, The Jacob Blaustein Institutes for Desert Research, Ben-Gurion University, Midreshet Ben-Gurion, 8499000, Israel
Laurent Chirio is unaffiliated:
U.R., A.F. and M.N. contributed equally to this work.
S.M. and R.G. jointly supervised this work.
A.M.B., R.G., S.M., U.R. conceived the study. R.G., C.D.L.O., U.R. designed the analyses. U.R. conducted the analyses. A.F., S.M., M.N., U.R. complied, designed and curated the dataset. R.G., S.M., U.R. wrote the paper. A.A., A.M.B., M.B., R.B., B.C., F.C.H., L.C., G.R.C., L.D., I.D., T.M.D., A.F., L.L.G., M.H., Y.I., F.K., A.L., M.L., E.M., D.M., M.M., S.M., C.C.N., M.N., Z.T.N., G.P., O.S.G.P., D.P.D., U.R., R.S., O.J.S.T., O.T.C., J.F.T., E.V., P.U., P.W., Y.W. provided, collated and verified underlying data. All authors read and commented on the manuscript.
The authors declare no competing financial interests.
Correspondence to Laurent Chirio or Shai Meiri.
Supplementary Methods, Supplementary Tables 1–4, Supplementary Figures 1–8, Supplementary References and descriptions of Supplementary Tables 3–4
Full list of sources per species and details for each source.
List of species found in the March 2015 version of the Reptile Database (http://www.reptile-database.org/data/reptile_checklist_2015_03.xls.zip) for which we did not present or analyse distributions.