Erosion of Lizard Diversity by Climate Change and Altered Thermal Niches

See allHide authors and affiliations

Science  14 May 2010:
Vol. 328, Issue 5980, pp. 894-899
DOI: 10.1126/science.1184695

Demise of the Lizards

Despite pessimistic forecasts from recent studies examining the effects of global climate change on species, and observed extinctions in local geographic areas, there is little evidence so far of global-scale extinctions. Sinervo et al. (p. 894; see the Perspective by Huey et al.) find that extinctions resulting from climate change are currently reducing global lizard diversity. Climate records during the past century were synthesized with detailed surveys of Mexican species at 200 sites over the past 30 years. Temperature change has been so rapid in this region that rates of adaptation have not kept pace with climate change. The models were then extended to all families of lizards at >1000 sites across the globe, and suggest that climate change-induced extinctions are currently affecting worldwide lizard assemblages.


It is predicted that climate change will cause species extinctions and distributional shifts in coming decades, but data to validate these predictions are relatively scarce. Here, we compare recent and historical surveys for 48 Mexican lizard species at 200 sites. Since 1975, 12% of local populations have gone extinct. We verified physiological models of extinction risk with observed local extinctions and extended projections worldwide. Since 1975, we estimate that 4% of local populations have gone extinct worldwide, but by 2080 local extinctions are projected to reach 39% worldwide, and species extinctions may reach 20%. Global extinction projections were validated with local extinctions observed from 1975 to 2009 for regional biotas on four other continents, suggesting that lizards have already crossed a threshold for extinctions caused by climate change.

Global climate change affects organisms in all biomes and ecosystems. Two natural compensatory responses are possible. Given enough time and dispersal, species may shift to more favorable thermal environments, or they may adjust to new environments by behavioral plasticity, physiological plasticity, or adaptation. Alternatively, failure to adjust or adapt culminates in demographic collapse and extinction. Despite accumulating evidence of contemporary climate change affecting species ranges and phenologies (13), evidence of extinctions at either local or global scales is lacking (46). Moreover, current forecasting models (7, 8) are not calibrated with actual extinctions, but are premised on hypothesized effects of thermal physiology on demography and extinction. Alternatively, models are based on range shifts or species-area relations in mobile species (1), but not extinctions (9). Hence, there is still much uncertainty regarding the expected magnitude of extinctions resulting from climate change (10).

Empirical validation of global extinction forecasts requires three forms of evidence. First, actual extinctions should be linked to macroclimate and validated to biophysical thermal causes arising from microclimate (11). Second, the pace of climate change should compromise thermal adaptation (10), such that evolutionary rates lag behind global warming owing to constraints on thermal physiology (12, 13). Third, extinctions due to climate should be global in extent.

From 2006 to 2008, we resurveyed 48 Sceloporus lizard species at 200 sites in Mexico that were first sampled in 1975 to 1995, and 12% of sites were locally extinct by 2009 (table S1).

Although Sceloporus lizards are heliotherms that bask and require solar radiation to attain physiologically active body temperatures (Tb) (14, 15), activity in hot weather may result in Tb exceeding CTmax, the critical thermal maximum, leading to death. Lizards retreat to cool refuges rather than risk death by overheating. However, hours of restriction (hr) in thermal refuges limit foraging, constraining costly metabolic functions like growth, maintenance, and reproduction, thereby undermining population growth rates and raising extinction risk. Lizards could evolve higher Tb, but this brings them closer to CTmax, which increases risk of overheating. Extinction risk may increase because of other thermal adaptations. For example, viviparity, which is posited to be a thermal adaptation to cold climates (16), may elevate extinction risk because high Tb can compromise embryonic development in utero (17).

We analyzed rate of change in maximum air temperature Embedded Image at 99 Mexican weather stations and constructed climate surfaces (tables S2 and S3, 1973 to 2008; fig. S1). Rate of change in Tmax was greatest for winter-spring (January to May; fig. S1 and table S3A) and increased faster in northern and central México and at high elevation, as evidenced by significant coefficients for fitted climate surfaces. We found a correlation between rate of change in Tmax during winter-spring breeding periods and local extinctions of Sceloporus species (table S3).

Many viviparous species in México are confined to high-elevation “islands,” where climate change has been most rapid. Logistic regression and multiple regression with phylogenetic independent contrasts (18, 19) revealed that extinction risk was significantly related to low latitudinal and altitudinal range limits (Fig. 1, A and B), where thermal physiology and/or ecological interactions limit species (20, 21). Phylogenetic correlation analysis (18) showed that extinction risk of viviparous lizards (18%) was twice that of oviparous lizards (9%, n =10000 bootstrap replications P < 0.001). Moreover, multiple regression based on phylogenetic independent contrasts (PICs; Fig. 1C and table S4) showed that extinction risk of viviparous taxa was significantly related to low Embedded ImageEmbedded Image and cool montane habitats Embedded Image, where climate has changed most rapidly in México.

Fig. 1

(A) Logistic regression of extinction probability (0 = extant, 1 = extinct) of Sceloporus lizards and reproductive mode: χ2 = 7.41, P = 0.025, Δelevation (χ2 = 8.53, P = 0.014), Δlatitude (χ2 = 7.14, P = 0.004), and Δlongitude (not significant), where Δ refers to deviations from species range midpoints. (B) Phylogenetic independent contrasts (PICs) of lineage survival (survival probability of local populations) and Δelevation (t = 2.15, P = 0.03), Δlatitude (t = 3.94, P = 0.0001), and Δlongitude (t = 2.66, P = 0.009). (C) PICs of lineage survival, Tb (t = 2.32, P = 0.02), Tair (t = 2.31, P = 0.02), and reproductive mode (t = −2.92, P = 0.005).

To validate patterns of extinction risk and Tb, we deployed thermal models (22) that record operative temperatures (Te) at two extinct and two persistent Yucatán sites of S. serrifer. Hours of restriction in activity (hr) during reproduction was significantly higher at extinct versus persistent sites (t = 9.26, P < 0.0001). By April 2009, hr at extinct Yucatán sites had become so severe that if S. serrifer were still present, it would have to retreat shortly after emergence (fig. S4A). Daily Tmax was positively correlated with hr assessed by Te (P < 0.001, fig. S4B). The relation between hr as a function of Tmax relative to S. serrifer’s Tb [hr = 6.12 + 0.74 × (Tmax Tb), eq. S2 (23)] is a general formula for predicting extinctions.

We modeled extinct/persistence status based on values for hr at Sceloporus sites derived from eq. S2 (23). The Yucatán ground truth for S. serrifer suggests that extinction occurs when hr exceeds 4. We calibrated this value with extinct/persistent Sceloporus sites. Goodness-of-fit tests of the model indicate that the best fit for observed and predicted extinctions at Sceloporus sites is hr > 3.85. If a species with a given Tb at a given geo-referenced site, subjected to Tmax,i, experienced hr > 3.85 during the 2-month reproductive period (March to April), we assumed that it would go extinct by 2009. Association of predicted and observed extinctions from this physiological model was significant for oviparous (χ2 = 49.0, P < 0.001) and viviparous taxa (χ2 = 4.2, P < 0.04).

As demography of high-elevation taxa becomes compromised due to climate change, species at low elevation that were previously limited by physiology and competition should expand into historically cooler habitat that is now warmer (20, 24), perhaps accelerating extinction of high-elevation forms. For viviparous taxa, six erroneously assigned extinct sites involved six of the eight cases of range expansion by low-elevation taxa, which all invaded from low to high altitudes or latitudes (table S1; significant by sign test, P < 0.001). Adding range shifts of competitors as a factor improved fit significantly between observed and predicted extinctions (Δlog likelihood = 45.37, 1 df, P < 0.0001, logistic regression). Therefore, competitive exclusion by invading low-elevation taxa appears to exacerbate climate-change extinctions of high-elevation taxa.

Lizards cannot evolve rapidly enough to track current climate change because of constraints arising from the genetic architecture of thermal preference (12, 13). A phylogenetic correlation between Tb and CTmax constrains adaptation. PIC regression of CTmax on Tb among Phrynosomatidae suggests that a shift in Tb by 1°C yields only a 0.5°C correlated response in CTmax (table S5 and fig. S7). Thus, CTmax may not evolve fast enough to keep up with evolved change in Tb. Furthermore, adaptive increase in Tb due to climate change is constrained by genetic correlations in which high Tb necessarily requires prolonged activity out of retreat sites (25), further increasing risk of overheating. Genetic trade-offs with energetically costly traits such as growth (25) also constrain adaptation.

The evolutionary response (R = h2s; s is the selection differential) necessary to keep pace with climate change is further constrained by low heritability for Tb, which we previously estimated at h2 = 0.17 for Sceloporus occidentalis in the laboratory (25). We used the physiological model to compute the sustained selection differential at each site j, such that Tb,j + ΔtTb,j evolves to match Tmax,j+ ΔtTmax,j, yielding Δhr,j= 0 and thereby rescuing population j from extinction [Δt computed over 1975 to 2009 (historical), 2009 to 2050, and 2050 to 2080]. We assumed sj = Rj/h2 = ΔtTb,j/h2, and generation times of 1 year versus 2 years (i.e., lowland versus montane Sceloporus, table S1). We expressed these critical levels of adaptive response as surfaces for ssustained, the sustained selection differential (Fig. 2B).

Fig. 2

(A) Sustained selection differentials per year required for Tb to keep pace with global warming. (B) Extinctions of Mexican Sceloporus lizards (1975 to 2009, 2009 to 2050, 2050 to 2080).

We compared the magnitude of selection allowing a species to adapt to climate change with maximum rates sustained under artificial or natural selection (26). Such comparisons are facilitated by dividing each sustained selection differential by the standard deviation (σTb = 1.23 for Tb of Mexican lizards) to obtain i, the standardized intensity of selection (26). Whereas i > 0.4 can be sustained in laboratory artificial selection for nine generations (27), studies in nature (26) indicate that i > 0.4 computed on an annual basis are rare (<5%). We also reference i to other anthropogenic causes of selection. Overfishing of Atlantic cod yielded i = 0.55, among the highest measured, but this selection regime caused demographic collapse of the fishery (28). In México, extinct sites sustained significantly higher i than persistent sites Embedded Image. The relation between intensity of selection and demographic collapse is simple. If sustained for decades, the mortality fraction necessary for selective shifts to new optima compromises population growth rate precipitating local extinction.

If climate change Embedded Image continues unabated in México, 56% of viviparous sites will be extinct by 2050 and 66% by 2080 (Fig. 2B). For oviparous sites, 46% will be extinct by 2050 and 61% by 2080. Based on local extinction of all populations surveyed for species, we project 58% species extinction of Mexican Sceloporus by 2080. Species extinction (58% by 2080) mirrors local population extinction (61 to 66%) because high-elevation endemics will go completely extinct as widespread lowland taxa expand to high elevations.

We used the model to derive global extinction projections (Fig. 3) for 34 lizard families (Table 1) with 1216 geo-referenced Tb records (table S6). Our data include heliotherms that bask and thermoconformers that do not bask, but track ambient air and surface temperature. Tmax was obtained from the WorldClim database (29) at 10–arc min resolution (1975, 2020, 2050, and 2080). We used distributional limits of heliothermic lizards of the world in 1975 to calibrate hr by family, which if exceeded at a given site would precipitate extinction. The extinction model is easily adapted to thermoconformers that maintain Tb close to Tair or retreat when Tair > Tpreferred. Assuming a sine wave for Tair between Tmin and Tmax (24-hour period), if the cumulative hours that Tair > Tb for a thermoconformer at a given geo-referenced site (table S6) exceeded the hr of a given lizard family, we assumed it would go extinct. Given Tmax Tb at each geo-referenced site, we computed the hr each species sustained in 1975, and for each family we used the upper 95% confidence level of hr (Table 1) as the extinction threshold (iteratively estimated, given global climate surfaces). Calibration with these 1975 distributional limits for Sceloporus yields hr = 3.9, which was cross-validated by hr = 3.85 computed from observed extinctions in México (1975 to 2009), and hr = 4, which was estimated directly from Te at persistent S. serrifer sites on the verge of extinction.

Fig. 3

Contour plots of global levels of local extinction for heliothermic lizards (1975 to 2009, 1975 to 2050, 1975 to 2080), assuming Embedded Image = 4.55 (23) and various Tb values.

Table 1

Sample size, Tb range, Embedded Image, Embedded Image, hr, and nspecies for 34 lizard families. Local extinction rates are based on geo-referenced Tb data and a physiological model of extinction. We also validated model predictions of local extinction risk in 2080 for six families: 57% (±3, n = 200) for Mexican Phrynosomatidae, 13% (±2, n = 3155) for South American Liolaemidae, 56% (±5, n = 117) for European Lacertidae (L. vivipara), 13% (±2, n = 1438) for African Cordylidae + Gerrhosauridae, 57% (±4, n = 125) on Madagascar, and 10% (±1, n = 2841) for Australian Egernia Group lizards species. Estimates of species extinctions in each family are derived from the relationships for extinction of all local populations for these six families (table S8).

View this table:

As in the validation of Mexican Sceloporus extinction, we computed hr for temperate lizards over 2 critical reproductive months, but were conservative in modeling critical months required for hr to be exceeded in the equatorial zone (±12° latitude) where lizards potentially breed year-round (hr exceeded over 12 months), and in the wet-dry tropical zone (±12° to 24° latitude: hr exceeded for 5 to 6 months).

Geo-referenced Tb samples indicate that current (2009) local extinctions average 4% worldwide (Table 1). Global averages will increase fourfold to 16% by 2050 and nearly eightfold to 30% by 2080, while equatorial extinctions will reach 23% by 2050 and 40% by 2080. Assuming reproduction shifts 1 month earlier in temperate zones [h2 = 1.0 lay date (30)] and proportionately less to the trade zones (i.e., no shift), 2080 global extinctions jump to 38% because spring seasons are warming faster across the globe. Our model is robust to plasticity in Tb (table S7) and initial assumptions made for reproductive periods in the tropics. If hr for equatorial taxa is computed over the 9 hottest months of reproduction, rather than the conservative assumption of 12 months, global extinctions increase to 39% by 2080.

The global generality of our model is verified by concordant distributions of current observed and predicted local extinctions of lizard biotas from four other continents (table S7). Our model pinpoints exact locations of two Liolaemid species going extinct in South America (Liolaemus lutzae, Phymaturus tenebrosus: χ2 = 32.1, P < 0.0001). In addition, the model predicts recent (2009) extinctions among 24 resurveyed populations of L. lutzae2 = 8.8, P = 0.003). In Europe, our resurvey of Lacerta vivipara revealed 14 extinct sites out of 46 (30%), which are predicted quite precisely by the model (χ2 = 24.4, P < 0.001). In Australia, the model pinpoints 2009 extinctions of Liopholis slateri2 = 17.8, P < 0.00001) and 2009 extinctions of Liopholis kintorei2 = 3.93, P = 0.047). In Africa, analysis of Gerrhosauridae and Cordylidae at 165 sites predicts <1% extinctions, and yet the model pinpoints the single extinction reported by 2009 (exact P-value = 0.006). We temper this value with extinction projections of 23% for 2009 at Malagasy Gerrhosauridae sites, which is validated by the observed 21% levels of local extinction across several lizard families in Madagascar nature reserves (23).

Thermoconforming lizards have been posited (31) to be more vulnerable to climate change relative to heliotherms. Even though Embedded Image of thermoconformers (27.5°C ± 1.8°) is significantly less thanEmbedded Image of heliotherms (33.5ºC ± 1.3, t = 2.66, P < 0.02, n = 34 families; Table 1), PICs show that extinction risk was unrelated to thermoregulatory mode (fig. S8), but was significantly increased by low Embedded Image, low hr, and high Embedded Image. The similar level of local extinctions in 2009 for Malagasy thermoconformers (21%, n = 63) and heliotherms [21%, n = 34; (23)] supports this view. Evolved changes in thermoregulatory mode, Tb, hr, lay date, and habitat preference set risk as Tmax rises, but owing to trade-offs, Tb and hr cannot be simultaneously maximized, hence extinction risk is independent of mode (fig. S8). Moreover, extinction risk is not higher for conformers because heliotherms inhabit equatorial regions (i.e., sub-Saharan Africa) that are unavailable to thermoconformers [a factor not considered by (31) or other models (10)], and these areas are warming rapidly (Fig. 3).

Our model, based on Tb, hr in activity during reproduction, and timing of breeding, assesses salient adaptations that affect thermal extinctions. Concordant verification of 2009 levels of local lizard extinction in North and South America, Europe, Africa, and Australia confirm that extinctions span tropical, temperate, rainforest, and desert habitats. Estimates of evolutionary rates required to keep pace with global change indicate that sustained and intense selection compromises population growth rates, precipitating extinctions. Probability of local extinction is projected to result in species extinction probabilities of 6% by 2050 and 20% by 2080 (table S8). Range shifts only trivially offset losses, because widespread species with high Tb shift to ranges of endemics, thereby accelerating their demise. Although global efforts to reduce CO2 may avert 2080 scenarios, 2050 projections are unlikely to be avoided; deceleration in Embedded Image lags atmospheric CO2 storage by decades (4). Therefore, our findings indicate that lizards have already crossed a threshold for extinctions.

Supporting Online Material

Materials and Methods

Figs. S1 to S9

Tables S1 to S8


References and Notes

  1. Materials and Methods are available as supporting material on Science Online.
  2. Research of B.S. was funded by the National Geographic Society, UC Mexus, UCSC Committee-On-Research, NSF awards (DEB 0108577, IBN 0213179, LTREB DEB 051597), CNRS fellowships, and visiting professorships (Museum Nationale d’Histoire Naturelle, Université Paris 6, Université Paul Sabatier Toulouse III), PAPIIT-UNAM IN213405 and 224208 to F.M.-C., a Université Paul Sabatier Toulouse III Visiting Professorship to D.B.M., CONACYT grants (4171N and 52852Q) to M.V.-S.C., grant CONACYT-SEP (43142-Q) to H.G., a CONACYT fellowship to R.N.M.-L., CNRS funding to B.H., and M.M., Biodivera: Tenlamas and from ANR Blanche: DIAME to J.C., CONICET grants to L.J.A. and M.M., FONDYCET 1090664 grants to P.V.S., CGL2005-03156 and CLG2008-04164 grant from SMSI to I.J.R., APCT-PICT1086 grant to N.I., scholarships and grants from Universidad Nacional Autónoma de México and American Museum of Natural History to M.V.-S.C., Academy of Finland grant (108955) to T.A.O., Australian Research Council grants to D.G.C., NSF awards DEB 0515909 and 0844523 to A.M.B., NSF award OISE 0530267, PIRE-Patagonia grant to J.W.S., L.J.A., M.M., and P.V.S. and Brigham Young University funding (Biology Department, Kennedy Center for International Studies, Bean Life Science Museum) to J.W.S.
View Abstract

Navigate This Article