Search tips
Search criteria 


Logo of procbhomepageaboutsubmitalertseditorial board
Proc Biol Sci. 2009 June 7; 276(1664): 1939–1948.
Published online 2009 March 4. doi:  10.1098/rspb.2008.1957
PMCID: PMC2677251

Why tropical forest lizards are vulnerable to climate warming


Biological impacts of climate warming are predicted to increase with latitude, paralleling increases in warming. However, the magnitude of impacts depends not only on the degree of warming but also on the number of species at risk, their physiological sensitivity to warming and their options for behavioural and physiological compensation. Lizards are useful for evaluating risks of warming because their thermal biology is well studied. We conducted macrophysiological analyses of diurnal lizards from diverse latitudes plus focal species analyses of Puerto Rican Anolis and Sphaerodactyus. Although tropical lowland lizards live in environments that are warm all year, macrophysiological analyses indicate that some tropical lineages (thermoconformers that live in forests) are active at low body temperature and are intolerant of warm temperatures. Focal species analyses show that some tropical forest lizards were already experiencing stressful body temperatures in summer when studied several decades ago. Simulations suggest that warming will not only further depress their physiological performance in summer, but will also enable warm-adapted, open-habitat competitors and predators to invade forests. Forest lizards are key components of tropical ecosystems, but appear vulnerable to the cascading physiological and ecological effects of climate warming, even though rates of tropical warming may be relatively low.

Keywords: climate warming, heat stress, body temperature, operative temperature

1. Introduction

Climate warming has profoundly affected terrestrial organisms (Walther et al. 2002; Parmesan 2006), and the magnitude of future effects is expected to increase with latitude (Root et al. 2003), paralleling the latitudinal increase in rate of warming (IPPC 2007). However, biotic impacts depend not only on the rate of warming but also on the number of species at risk, the behavioural and physiological abilities of organisms to buffer climate warming (Huey et al. 2003; Helmuth et al. 2005) and on physiological sensitivity to warming (Huey & Slatkin 1976; Deutsch et al. 2008; Tewksbury et al. 2008). The biodiversity of terrestrial ectotherms (e.g. insects, amphibians, reptiles) is of course vastly greater at lower latitudes. Moreover, many tropical ectotherms are surprisingly intolerant of high temperatures (Ruibal 1961; Parsons 1989) or are relatively sensitive to temperature change (Janzen 1967; van Berkum 1988; Deutsch et al. 2008). Nevertheless, biologists have given limited attention to whether tropical ectotherms are vulnerable to warming (Parsons 1989; Pounds et al. 1999; Root et al. 2003; Parmesan 2007; Colwell et al. 2008; Deutsch et al. 2008; Raxworthy et al. 2008; Tewksbury et al. 2008; Williams et al. 2008; Chen et al. 2009; Kearney et al. 2009).

Reliable projections of the biological impacts of climate change require information on environmental biophysics as well as the body temperature, thermal physiology, behaviour and ecology of species (Huey & Slatkin 1976; Tracy & Christian 1983; Dunham 1993; Buckley 2008; Angilletta 2009; Kearney et al. 2009). Air temperatures recorded by weather stations are often used as a proxy for body temperature (Tb), but actual Tbs are far superior (Chown & Terblanche 2007). Unfortunately, Tb data are lacking for most ectotherms. However, such data (as well as habitats, thermoregulatory behaviours and activity patterns) are available for diverse species of lizards from many latitudes. Moreover, these data are often accompanied by laboratory data on thermal preferences (Tp), which are often co-adapted with optimal performance temperatures (To; Huey & Bennett 1987; Martin & Huey 2008), by critical thermal temperatures (CTmin, CTmax) and sometimes by measures of the temperature sensitivity of physiological performance (Bennett 1980; Hertz et al. 1983; Angilletta et al. 2002). Our knowledge of lizard thermal performance comes largely from studies of sprint speed: accordingly, we compiled data on the thermal dependence of this ecologically relevant trait (Bennett 1980) for more than 70 species of lizards (see table 1 in the electronic supplementary material). The completeness and richness of behavioural and physiological data, as well as information on phylogenetic relationships, are unsurpassed for ectotherms and permit the first physiologically and phylogenetically based assessment of latitudinal patterns in the performance consequences of global warming.

We first conducted a comparative, macrophysiological analysis (Chown et al. 2004) of thermoregulatory behaviour and thermal sensitivity of lizards as a function of latitude. Then we focused on several Puerto Rican lizards. Using field Tb data collected several decades ago as a temporal baseline, as well as data on the thermal sensitivity of sprinting, we simulated how warming should affect both lizard Tb and relative performance during summer (Huey 1983). These complementary analyses use different approaches but lead to the same conclusion: some tropical forest lizards ‘can't take the heat’. Moreover, theoretical models yield this same conclusion (Deutsch et al. 2008; Kearney et al. 2009). Given the huge diversity of ectotherms in the tropics, the implications of these consistent patterns are daunting.

2. Material and methods

(a) Macrophysiological analyses

Data on mean values of CTmin, CTmax, To, Tb (of active lizards), habitat associations, basking behaviour and absolute latitude were compiled from the literature (see tables 1 and 2 in the electronic supplementary material). To minimize methodological noise (Chown et al. 2003; Ives et al. 2007), we analysed only species for which the thermal dependence of sprint speed had been quantified. Most studies (63 out of 70) were done by us, our students or our close colleagues; thus, protocols were uniform by comparative standards.

We include independent field data on habitat association (forest and open), thermoregulatory behaviour (basking and non-basking), Tb and Ta for lizards from 12 neotropical sites (available as table 2b in the electronic supplementary material). Most sites were approximately 1.5×1.5 km. These field data were collected by a single team of investigators (Vitt et al. 2008), and thus methods were consistent.

To evaluate an association of absolute latitude, taxonomic affinity and basking behaviour with the physiological traits, we computed (multiple) regressions involving the three independent variables and their interactions. Phylogenetic relationships were considered, and the tree (and sources) is available in table 3 in the electronic supplementary material.

We computed regressions in three ways (reviews in Garland et al. 2005; Lavin et al. 2008): conventional (non-phylogenetic) ordinary least squares (OLS); phylogenetic generalized least squares (PGLS); and regression in which residuals were modelled as having evolved via an Ornstein–Uhlenbeck process (RegOU), thus mimicking stabilizing selection. These three models form a continuum between assuming a star phylogeny with no hierarchical structure (OLS), a phylogeny as specified by the user (PGLS) and something that can take on values intermediate between the star and the specified hierarchical phylogeny (RegOU) (Garland et al. 2005). Candidate-independent variables included absolute latitude, a 0–1 dummy variable to code for non-basking versus basking and a series of dummy variables to code for ‘family’ membership (Lavin et al. 2008). Details of analysis and model comparisons are described in the electronic supplementary material.

(b) Focal species analyses of Puerto Rican lizards

Extensive field and laboratory data for several Puerto Rican lizards enable us to describe the thermal biology of these lizards several decades ago, and then to predict how warming will affect their Tb and relative sprint performance. Details of methods are given in the electronic supplementary material. Field data on Tb and Ta, laboratory thermal preference (Tp), CTmax and CTmin for Anolis cristatellus from Punta Salinas were collected in July 1973 (Huey 1974; Huey & Webster 1976). Thermal sensitivities of sprint speed were for lizards from this population collected in 1981 (Huey 1983). The preferred temperature range is arbitrarily delimited by the central 50 per cent of all Tp records (Hertz et al. 1993).

Field data for Tbs and operative temperatures (Te) were obtained for A. cristatellus at San German in August 1983 (Hertz 1992). Tes were obtained using a randomly placed array of 60 hollow electroformed copper models (calibrated against live lizards) that match the size and reflectivity of adults. Thermal preference shows no significant geographical variation (Huey & Webster 1976), and so we used Tps from Punta Salinas.

Field data for Tb and Ta of Anolis gundlachi from El Verde were obtained in July 1972, and Tp were obtained in 1973 (Huey & Webster 1976). Long-term weather records are available for this site (, enabling us to estimate the actual shift in mean July maximum temperature from 1975 to 2008. The predicted shift in Ta was 2.1°C (see figure 1 in the electronic supplementary material).

Field data for lowland populations of diurnal geckos (Sphaerodactylus spp.) were obtained in 1988–1991 (Álvarez 1992; table 2h in the electronic supplementary material). Basic techniques for obtaining Tbs and Tes followed Hertz (1992), except that small polyvinyl chloride cylinders (n=27) were used to estimate Tes.

3. Results and discussion

(a) Natural history, phylogeny and latitude

Basic natural history and phylogenetic issues lay a critical foundation for understanding latitudinal patterns of lizard thermal biology. Over four decades ago, Ruibal (1961) noted that temperate and tropical lizards (diurnal) differ in habitat use, thermoregulatory behaviour and Tb. Almost all temperate-zone lizards are restricted to open habitats where operative temperatures (Tes) are high and sufficiently heterogeneous to enable thermoregulation (Porter et al. 1973; Bakken 1989): most species bask, thermoregulate relatively carefully and maintain high Tb that are often well above associated Ta (Porter et al. 1973; Huey 1982; van Berkum 1988). (The few exceptions are cryptozoic species that infrequently emerge from cover (see below).)

Some tropical and subtropical lizards also live in open habitats and generally share behavioural and physiological characteristics (e.g. basking behaviour, elevated Tb) with high-latitude species living in open habitats (Ruibal 1961; Porter & James 1979; van Berkum 1988; Vitt et al. 1998). However, many other tropical and subtropical lizards live below the canopy inside deeply shaded forests (Inger 1959; Ruibal 1961; Huey 1982), where Tes are relatively low and homogeneous (Hertz 1992), making thermoregulation difficult or expensive (Huey 1974). In fact, these tropical forest lizards rarely bask and are active at relatively low Tbs that are necessarily close to air temperatures (Inger 1959; Ruibal 1961; Huey 1982), because convection dominates heat exchange in such habitats (Bakken 1992).

In the neotropics, non-basking (forest) species often outnumber basking (open-habitat) species at lowland sites with mixed forest and open habitats (see table 2b in the electronic supplementary material). Thus, non-basking species constitute a major component of the diversity of neotropical lizards. In fact, these species constitute an average of 62 per cent of 63 species at 12 sites from Nicaragua to Brazil (see table 2b in the electronic supplementary material).

Close associations among habitat, basking behaviour and Tb are evident for neotropical species (see table 2b in the electronic supplementary material). Of the species found in forests, only 1 out of 40 was basking. Of the species in open or edge habitats, 21 out of 24 were basking. Importantly, activity Tbs of tropical forest lizards averaged 6°C lower than that of tropical open-habitat lizards (t-test, p[double less-than sign]0.001, 28.9±0.31°C (n=38) and 34.9±0.51°C (n=20)).

These linkages appear to have phylogenetic and latitudinal components. Different neotropical lineages tend to comprise either predominately forest-dwelling, non-basking and low Tb species or open-edge dwelling, basking and high Tb species (see table 2a in the electronic supplementary material). (Note: these families show significant heterogeneity in percentages of basking and forest-dwelling species (non-phylogenetic Χ2-tests, both p<0.01).) Moreover, lineages in the tropics that include predominately forest-dwelling and low Tb species (e.g. diurnal Gekkonidae, Gymnophthalmidae, Polychrotidae) are largely restricted to tropical and subtropical regions, whereas lineages with a preponderance of open-habitat species and higher Tb species (e.g. Teiidae) often extend into the temperate zone.

Overall, most temperate-zone lizards are baskers, whereas tropical species can be either baskers on non-baskers. Moreover, these different ‘thermal syndromes’ tend to run in families. Exceptions do exist. A few temperate lineages are cryptozoic and have low Tb (e.g. Anguidae: Kingsbury 1994; Xantusiidae: Mautz et al. 1992; some Scincidae: Bennett & John-Alder 1986; Shine 1999). Also, some tropical families and species (e.g. A. cristatellus) include both types. However, we are confident that these trends are general, at least for New World lizards.

(b) Field body temperatures

Early workers, who apparently expected that tropical lizards would have relatively high Tb simply because lowland tropical sites are always warm, were surprised to find that some tropical species were in fact active at Tbs well below those of most temperate-zone species (Brattstrom 1965). We evaluated field Tbs for species from diverse latitudes (figure 1c). In an ordinary least-squares analysis, which assumes a star phylogeny (Garland et al. 1992), Tb increased with latitude (two-tailed p=0.0145; see table 2c in the electronic supplementary material). However, field Tb showed a strong phylogenetic signal (Blomberg et al. 2003; p=0.007), and the best-fit regression model (RegOU) for Tb excluded latitude but included basking and family membership. Thus, the non-phylogenetic trend probably stems from the many non-basking, low Tb species in the tropics.

Figure 1
Thermal sensitivity of diurnal lizards versus latitude. (a) Heat tolerance (CTmax), (b) optimal sprint temperature (To), (c) mean-field body temperature (Tb) and (d) cold tolerance (CTmin) versus latitude. In phylogenetic analyses, CTmax, To and Tb are ...

(c) Lizard thermal sensitivity

In the OLS regression, CTmin decreased with absolute latitude (p[double less-than sign]0.001; figure 1d; see table 2d in the electronic supplementary material). Phylogenetic signal for CTmin was significant (p<0.001), and absolute latitude was the only variable in the best-fit model (RegOU). As in the non-phylogenetic analysis, CTmin declined with latitude.

Optimal sprint temperature (To; figure 1b) and heat tolerance (CTmax; figure 1d) shared similar patterns. In OLS, both variables increased with latitude (both p[double less-than sign]0.001; see table 2e,f in the electronic supplementary material). Both traits showed strong phylogenetic signal (p[double less-than sign]0.001), and the best-fit models (OLS) showed that To and CTmax are elevated for basking species versus non-baskers and differ among taxa. As with Tb, the non-phylogenetic trend with latitude may be a consequence of the tropical concentration of non-basking taxa.

Note that Tb, To and CTmax appear independent of latitude but are associated with basking behaviour and with taxon, whereas CTmin is independent of basking and taxon and depends only on latitude. These contrasting associations might reflect differential effectiveness of behaviour in buffering climate-associated selection on these thermal traits (van Berkum 1988). When active, a thermoregulating lizard in an open habitat at any latitude can usually achieve some control over its Tb and can simultaneously avoid stressfully high temperatures. Consequently, mean Tb (herein) and also maximal Tb (van Berkum 1988, p. 335) are generally independent of latitude (but see, Clark & Kroll 1974), but do tend to show strong phylogenetic conservatism (Huey 1982; Hertz et al. 1983) and a strong association with habitat and basking behaviour (herein, Ruibal 1961; Clark & Kroll 1974). Because Tb, To and CTmax are likely to be co-adapted traits (Huey & Bennett 1987; Huey & Kingsolver 1993; Angilletta et al. 2006; Martin & Huey 2008; Angilletta 2009), we are not surprised that these thermal traits were closely associated with phylogenetic affinities and basking behaviour.

By contrast, a lizard hibernating at high latitude in a winter refuge has relatively few options for behavioural thermoregulation. Consequently, temperatures of high-latitude lizards in winter (and CTmin) should decline with latitude (van Berkum 1988; Kearney et al. 2009). Thus, thermoregulatory behaviour may buffer selection (Bogert 1949; van Berkum 1988; Huey et al. 2003) on some thermal traits (e.g. Tb, CTmax), but not on all (CTmin).

(d) Some tropical lizards have narrow thermal safety margins

Whether climate warming is detrimental depends not only on the heat tolerance of species (i.e. on To and CTmax) but also on whether current Tes (Roughgarden et al. 1981; Bakken 1992) are already high relative to an ectotherm's To (Tracy & Christian 1983; Huey 1991; Deutsch et al. 2008). If Tes in shaded microenvironments are below an ectotherm's To and if shade is readily accessible, ectotherms confronting climate warming can find thermal refugia and thus remain active (Kearney et al. 2009). However, if Tes in shade are currently equal to or greater than an ectotherm's To, climate warming will force ectotherms to retreat below ground or tolerate activity at a high Tb (Dunham 1993), which will probably induce stress and reduce performance (Huey 1983). If time restrictions and performance declines are substantive, the population might go extinct (Dunham 1993; Kearney et al. 2009).

We quantified potential access to an above-ground thermal refuge by computing a ‘thermal safety margin’, which is the difference between a lizard's To for sprinting and the maximal Te in shade (Deutsch et al. 2008). Ideally, shade Te should be calculated for lizards on a global scale (Kearney et al. 2009), but here we substituted mean maximal daytime air temperature (Ta,max) during the warmest three months of the year (New et al. 2002). Ta,max will closely approximate the maximal average Te for small ectotherms in shade (Rogowitz 1996), because convection dominates heat exchange in shade (Bakken 1992). Of course, Ta,max underestimates the maximal Te in open microhabitats. Nevertheless, the issue here is whether lizards facing climate warming can access thermal refuges (Kearney et al. 2009), and thus Ta,max is a relevant index.

Thermal safety margins showed only a weak phylogenetic signal (p=0.09). OLS regression was the best-fit model (see table 2g in the electronic supplementary material), and safety margin increased with latitude and basking (figure 2). Thermal safety margins are thus small—and sometimes even slightly negative—for most tropical forest lizards. Because these forest species live in already warm environments, but have low Tb, To and CTmax, even a small increase in Te may force them to reduce activity time (Kearney et al. 2009) or to suffer reduced performance in summer. Thus, tropical forest lizards may not be able to ‘take the heat’.

Figure 2
Thermal safety margin (ToTa,max) for diurnal lizards increases with latitude. Open circles, basking species; filled circles, non-basking species. In a phylogenetic analysis, thermal safety margin increases with latitude. The tropical species ...

By contrast, many higher latitude or high-altitude species will not need to retreat if Te in shade increases with climate warming; rather, these species may even benefit (Kearney & Porter 2004; Chamaillé-Jammes et al. 2007; Buckley 2008; Kearney et al. 2009) because warming will improve the thermal quality of their habitat (Tracy & Christian 1983; Huey 1991; Hertz et al. 1993). Of course, some diurnal species (Dunham 1993) living in mid-latitude deserts (and cryptozoic species) also have small safety margins (figure 2), and these may also suffer from warming.

(e) Lizards in lowland Puerto Rico are already at risk

The latitudinal pattern of thermal safety margins (figure 2) suggests that lizards in tropical forests might be relatively vulnerable to climate warming in summer. However, such a macrophysiological approach ignores heterogeneity of microclimates and thermoregulatory behaviour. To add a biophysically and physiologically grounded approach, we analyse data for several lizards from Puerto Rico.

Anolis cristatellus is an abundant and widespread species, and its thermal biology is well known (e.g. Rand 1964; Huey & Webster 1976; Gorman & Hillman 1977; Hertz 1992). In shaded lowland habitats, this species lives on tree trunks, does not bask and is a thermoconformer, such that Tb closely matches Ta, and thus Te (Huey 1974; Hertz 1992; see details of methods in the electronic supplementary material). It will bask in open lowland habitats, but only very early and late in the day (Huey 1974; Hertz 1992). Similar to other anoles (van Berkum 1988), it is relatively intolerant of high Tb (see table 1 in the electronic supplementary material).

Body and air temperatures and behaviour of A. cristatellus were monitored in summer approximately 35 years ago (1972 and 1973) at Punta Salinas on the northern coast (Huey 1974; Huey & Webster 1976), and we use the 1973 data as a temporal baseline. In a dense forest, these lizards were classic thermoconformers (figure 3): they were active from sunrise to sunset; had no opportunity to bask; and had Tbs that averaged only 0.6±0.05°C above Ta (Huey 1974; Huey & Webster 1976). This habitat was thermally suitable for this species: even without thermoregulating, lizards had Tbs that usually fell within the Tp range (figure 3a) and that enabled them to sprint at greater than 90 per cent of their maximal speed from sunrise to sunset (figure 3b).

Figure 3
Impact of climate warming on Tb and performance of a tropical lizard. (a) Body temperature and (b) predicted relative sprint performance of A. cristatellus in a shaded forest at Punta Salinas, Puerto Rico before (grey boxes, 1973) and after (red boxes) ...

To estimate the impact of climate warming on A. cristatellus over the next 100 years, we conservatively assume that Ta will increase by 3°C above levels measured in the early 1970s (Malhi & Wright 2004; IPPC 2007). Because convection dominates heat exchange in forests, Tb will also increase by approximately 3°C. As a result, lizards will have elevated Tb (red boxes in figure 3a) that exceed Tp for most of the day, and their predicted sprint performance should be lower for most of the day (figure 3b).

Additional data (figure 3c) were collected for this species approximately 25 years ago (1983) at San German, a lowland forested site (90 m) on the warmer, southwestern side of Puerto Rico (Hertz 1992). Here, A. cristatellus were again active from sunrise to sunset and were thermoconformers. Their Tbs averaged only 0.6±0.09°C above the mean Te of their habitat. During most of the day, Tb and Te were within the lizard's preferred temperature range, and only 20.8 per cent of Tb (and 7.8% of Te) exceeded that range. If Te increases by 3°C, 68.6 per cent of all predicted Tb records (and 70.6% of all Te) should exceed the lizard's preferred range (figure 3c). In fact, 92.1 per cent of all Te should exceed the preferred range during midday (10.00–14.00 hours). Moreover, relative sprint performance at midday should decline precipitously (figure 3d).

Thermal data are also available for several species of the diurnal ground gecko, Sphaerodactylus, studied in 1988–1991 (Álvarez 1992). These geckos typically inhabit leaf litter. In cool seasons or at high altitude, these geckos sometimes bask; but in lowland forests, they rarely do so (Álvarez 1992). Similar to Anolis, these lizards are relatively intolerant of high temperatures: most have a CTmax below 40°C and a Tp only approximately 30°C (see table 2h in the electronic supplementary material). In 1991, lowland populations of four species were already living in environments that were stressfully warm in summer; in fact, mean Te exceeded mean Tp and mean To by 3.2±1.06 and 5.1±0.70°C, respectively. As a consequence, mean Tb during the day in summer (across species) exceeded mean Tp by 2.1±0.43°C on average and exceeded To for sprinting by 3.9±1.16°C on average (see table 2h in the electronic supplementary material). If Te increases by 3°C because of climate warming, these lizards will probably experience severe heat stress because mean Te over the day in summer will be only 3.5±0.90°C below the CTmax of these lizards.

Although climate warming may well stress Puerto Rican lizards in summer, it should benefit them in winter, as the slightly lower temperatures then depress reproduction (Gorman & Licht 1974; Lister 1981). Therefore, a full demographic model will be necessary to evaluate the overall effects of warming. However, because thermal fitness curves are asymmetric (Gilchrist 1995; Huey & Berrigan 2001), being ‘too hot’ is likely to be much worse physiologically than being ‘too cold’.

(f) Cascading effects of warming on ecological interactions

Climate warming will not only induce thermal stress in tropical forest lizards, but also probably force them to suffer increased competition and predation from warm-adapted, open-habitat lizards. Some open-habitat high-Tb lizards, such as Ameiva festiva (Costa Rica), make periodic forays into shaded forest habitats to search for food. However, they must soon return to the open to bask (van Berkum et al. 1986). As climate warming elevates Te inside forests, A. festiva will cool more slowly, forage longer and thus put increased competitive and predatory pressures on forest lizards (Vitt et al. 1998).

Climate warming may even enable open-habitat species to displace tropical forest species. In the early 1970s at El Verde, Puerto Rico (approx. 350 m), A. gundlachi was a forest-dwelling, non-basking species, whereas A. cristatellus was restricted to open and edge habitats. Relative to A. cristatellus, A. gundlachi had a significantly lower field Tb (figure 4), Tp and CTmax (Huey & Webster 1976; Hertz 1992; Rogowitz 1996).

Figure 4
Predicted impact of recent climate warming on a tropical forest lizard. Observed climate warming between 1972 and 2008 (see figure 1 in the electronic supplementary material) is likely to make forest habitats (El Verde, Puerto Rico) less suitable in summer ...

In July 1972, the thermal environment inside the forest at El Verde was highly suitable for A. gundlachi (Huey & Webster 1976); indeed, 74.7 per cent of all Tbs during the entire day (figure 4) were within 1°C of the Tp range for this species, and only 4.5 per cent of Tbs were 1°C or more above its Tp range. By contrast, forest Tes would have been too cool for A. cristatellus (Gorman & Hillman 1977); 95.5 per cent of all predicted Tbs in summer would have been at least 1°C below its Tp range (figure 4).

The observed 2.1°C climate warming in summer since 1975 (see figure 1 in the electronic supplementary material) should have made the forest less suitable thermally for A. gundlachi in that season. Only 39.6 per cent of its predicted Tbs should now be within 1°C of its Tp range, and 60.3 per cent should be at least 1°C warmer than its Tp range, suggesting likely heat stress. By contrast, 60.4 per cent of predicted Tb of A. cristatellus should now be within 1°C of its Tp range (figure 4). Therefore, if warming continues, A. cristatellus might soon move into the forest, at least during summer, and potentially even displace A. gundlachi.

4. Concluding remarks

Both macrophysiological (figures 1 and and2)2) and focal species (figures 3 and and4)4) approaches suggest that lowland lizards living in neotropical forests are at risk from climate warming; in fact, some species are probably already experiencing body temperatures at or above their physiological optima, at least in summer (figures 3 and and4;4; see table 2h in the electronic supplementary material). The cascading effects of increased competition and predation from open-habitat species will only compound problems for forest species and may have major effects on tropical forest food webs. Moreover, performance of some species may have already been harmed by recent climate warming (figure 4).

Other studies have raised concern for the vulnerability of tropical ectotherms (e.g. Parsons 1989; Pounds et al. 1999; Root et al. 2003; Parmesan 2007; Colwell et al. 2008; Deutsch et al. 2008; Raxworthy et al. 2008; Williams et al. 2008; Kearney et al. 2009). Moreover, population declines of lowland forest lizards in Costa Rica have already been noted (Whitfield et al. 2007). Because lowland tropical forests are centres of biodiversity, these predictions and trends are disturbing.

A reviewer of a draft of this paper felt that we have taken a ‘pessimistic stance’. We disagree. To be sure, physiological acclimation could potentially buffer the impact of climate warming (Chown & Terblanche 2007). However, the few available data suggest that tropical forest ectotherms—including Puerto Rican Anolis (Rogowitz 1996)—show limited acclimation responses (Brattstrom 1968; Feder 1982; Tsuji 1988; Ghalambor et al. 2006). Similarly, genetic adaptation could potentially buffer the impact (Hoffmann & Blows 1993; Travis & Futuyma 1993; Angilletta 2009), and rapid response to temperature in selection experiments is common in invertebrates (Santos et al. 2005). Nevertheless, artificial selection on heat tolerance of a fish was unsuccessful (Baer & Travis 2000). Unfortunately, the heritability of thermal traits in lizards is essentially unstudied (Sinervo 1990; Angilletta 2009). In any case, conservatism seems to be the rule in lizard thermal evolution (Hertz et al. 1983; Huey & Bennett 1987). Thus, adaptive rescue is feasible but unlikely.

For some species, thermoregulatory behaviour will effectively buffer the impact of climate warming (Huey et al. 2003; Angilletta 2009; Kearney et al. 2009). For example, many desert lizards readily evade heat stress by retreating underground during warm periods (Porter et al. 1973; Stevenson 1985). By contrast, tropical forest lizards, such as Anolis and Sphaerodactylus, have relatively few thermoregulatory options because forest Tes are relatively homogeneous in space and time (Hertz 1992); the ‘thermal mosaic’ for such forest species is relatively monochromatic. Populations with access to a mountain could potentially migrate uphill (Colwell et al. 2008; Raxworthy et al. 2008; Chen et al. 2009), but the largest remaining tropical rainforest (Amazon basin) has mountains only on one edge.

The biotic pressure of open-habitat species on forest ectotherms will be exacerbated if climate change also opens tropical forest canopies (Clark et al. 2003; Feeley et al. 2007; Whitfield et al. 2007); more solar radiation will penetrate the forest, thus increasing Te. Canopy opening will also reduce humidity and increase wind speed, harming ectotherms sensitive to desiccation (Parsons 1989; Pounds et al. 1999). More importantly, deforestation and fragmentation will directly accelerate extinction rates (Brook et al. 2003; Ferraz et al. 2003) as well as increase hot edge environments (Vitt et al. 1998).

Of course, climate warming would not be physiologically detrimental to all tropical forest lizards. Some species living in cool montane forests should benefit (Huey & Webster 1976; Hertz 1992; Hertz et al. 1993). Also, even lowland species may benefit during cool seasons, when Tes and Tbs are slightly lower (Lister 1981; Hertz 1992) and reproduction is reduced (Gorman & Licht 1974; Lister 1981). Similarly, some high-latitude species may benefit because warmer temperatures will increase potential activity times and growth rates (Kearney & Porter 2004; Chamaillé-Jammes et al. 2007; Buckley 2008; Deutsch et al. 2008; Kearney et al. 2009).

Our prediction that neotropical forest lizards are vulnerable to climate warming can be evaluated in future studies. One predictive approach develops biophysical, physiological and demographic models that transduce predicted climate change onto changes in vital demographic rates (Dunham 1993; Helmuth et al. 2005; Crozier & Dwyer 2006; Buckley 2008). Kearney et al. (2009) have recently made a major step in this direction. Our conclusions are remarkably consistent with theirs: the primary challenge for temperate-zone ectotherms is gaining heat, but that for tropical ectotherms is avoiding it.

A complementary approach involves field studies: if lowland forest populations are indeed currently heat stressed in summer (figures 3 and and4),4), those populations should now show reduced rates of growth and reproduction in summer, relative to cooler seasons (or preferably to earlier decades). Ultimately, time-series monitoring of population ecology (Chamaillé-Jammes et al. 2007; Whitfield et al. 2007), food resources (Buckley 2008), distributions (Colwell et al. 2008; Raxworthy et al. 2008), Tb, habitat usage and microclimates may be necessary to evaluate whether lowland tropical forest lizards are indeed on—and perhaps already even over—the thermal edge. Field data from earlier decades (e.g. figures 3 and and4;4; Chen et al. 2009) provide a critical baseline for monitoring climate-induced changes in temperatures and performance of tropical lizards.


We dedicate this paper to Rodolfo Ruibal, whose 1961 paper pioneered empirical studies of the thermobiology of tropical ectotherms and whose insights into tropical and temperate ectotherms laid the conceptual foundation for this paper. We thank M. Angilletta and M. Kearney for their comments. This study was supported by National Science Foundation grants to P.E.H., R.B.H., J.J.T., L.J.V. and T.G., and by a Program on Climate Change Fellowship to C.A.D. El Verde temperature data were provided courtesy of the Luquillo Experimental Forest (LTER).

Supplementary Material

Data on thermal dependence of lizards:
Summaries of data and comparative analyses:
Phylogenetic information:

Phylogenetic sources and phylogenetic tree

Details of methods:
Climate warming in Puerto Rico:


  • Álvarez H.J. University of Puerto Rico; San Juan, PR: 1992. Thermal characteristics of Sphaerodactylus species in Puerto Rico and their implications for the distribution of species in Puerto Rico.
  • Angilletta M.J. Oxford University Press; Oxford, UK: 2009. Thermal adaptation: a theoretical and empirical synthesis.
  • Angilletta M.J., Jr, Niewiarowski P.H., Navas C.A. The evolution of thermal physiology in ectotherms. J. Therm. Biol. 2002;27:249–268. doi:10.1016/S0306-4565(01)00094-8
  • Angilletta M.J., Jr, Bennett A.F., Guderley H., Navas C.A., Seebacher F., Wilson R.S. Coadapation: a unifying principle in evolutionary thermal biology. Physiol. Biochem. Zool. 2006;79:282–294. doi:10.1086/499990 [PubMed]
  • Baer C.F., Travis J. Direct and correlated responses to artificial selection on acute thermal stress tolerance in a livebearing fish. Evolution. 2000;54:238–244. doi:10.1111/j.0014-3820.2000.tb00024.x [PubMed]
  • Bakken G.S. Arboreal perch properties and the operative temperature experienced by small animals. Ecology. 1989;70:922–930. doi:10.2307/1941359
  • Bakken G.S. Measurement and application of operative and standard operative temperatures in ecology. Am. Zool. 1992;32:194–216. doi:10.1093/icb/32.2.194
  • Bennett A.F. The thermal dependence of lizard behaviour. Anim. Behav. 1980;28:752–762. doi:10.1016/S0003-3472(80)80135-7
  • Bennett A.F., John-Alder H.B. Thermal relations of some Australian skinks (Sauria: Scincidae) Copeia. 1986;1986:57–64. doi:10.2307/1444888
  • Blomberg S.P., Garland T., Jr, Ives A.R. Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution. 2003;57:717–745. doi:10.1111/j.0014-3820.2003.tb00285.x [PubMed]
  • Bogert C.M. Thermoregulation in reptiles, a factor in evolution. Evolution. 1949;3:195–211. doi:10.2307/2405558 [PubMed]
  • Brattstrom B.H. Body temperatures of reptiles. Am. Midl. Nat. 1965;73:376–422. doi:10.2307/2423461
  • Brattstrom B.H. Thermal acclimation in anuran amphibians as a function of latitude and altitude. Comp. Biochem. Physiol. 1968;24:93–111. doi:10.1016/0010-406X(68)90961-4 [PubMed]
  • Brook B.W., Sodhi N.S., Ng P.K.L. Catastrophic extinctions follow deforestation. Nature. 2003;424:420–423. doi:10.1038/nature01795 [PubMed]
  • Buckley L.B. Linking traits to energetics and population dynamics to predict lizard ranges in changing environments. Am. Nat. 2008;171:E1–E8. doi:10.1086/523949 [PubMed]
  • Chamaillé-Jammes S., Massot M., Aragón P., Clobert J. Global warming and positive fitness response in mountain populations of common lizards Lacerta vivipara. Glob. Change Biol. 2007;12:392–402. doi:10.1111/j.1365-2486.2005.01088.x
  • Chen I.-C., Shiu H.-J., Benedick S., Holloway J.D., Chey V.K., Barlow H.S., Hill J.K., Thomas C.D. Elevation increases in moth assemblages over 42 years on a tropical mountain. Proc. Natl Acad. Sci. USA. 2009;106:1479–1483. doi:10.1073/pnas.0809320106 [PubMed]
  • Chown S.L., Terblanche J.S. Physiological diversity in insects: ecological and evolutionary contexts. Adv. Insect Physiol. 2007;33:50–152. doi:10.1016/S0065-2806(06)33002-0 [PMC free article] [PubMed]
  • Chown S.L., Addo-Bediako A., Gaston K.J. Physiological diversity: listening to the large-scale signal. Funct. Ecol. 2003;17:562–572. doi:10.1046/j.1365-2435.2003.07622.x
  • Chown S.L., Gaston K.J., Robinson D. Macrophysiology: large-scale patterns in physiological traits and their ecological implications. Funct. Ecol. 2004;18:159–167. doi:10.1111/j.0269-8463.2004.00825.x
  • Clark D.R., Jr, Kroll J.C. Thermal ecology of anoline lizards: temperate versus tropical strategies. Southwest. Nat. 1974;19:9–19. doi:10.2307/3669786
  • Clark D.A., Piper S.C., Keeling C.D., Clark D.B. Tropical rain forest tree growth and atmospheric carbon dynamics linked to interannual variation during 1984–2000. Proc. Natl Acad. Sci. USA. 2003;100:5852–5857. doi:10.1073/pnas.0935903100 [PubMed]
  • Colwell R.K., Brehm G., Cardelús C.L., Gilman A.C., Longino J.T. Global warming, elevational range shifts, and lowland biotic attrition in the wet tropics. Science. 2008;322:258–261. doi:10.1126/science.1162547 [PubMed]
  • Crozier L., Dwyer G. Combining population-dynamic and ecophysiological models to predict climate-induced insect range shifts. Am. Nat. 2006;167:853–866. doi:10.1086/504848 [PubMed]
  • Deutsch C.A., Tewksbury J.J., Huey R.B., Sheldon K.S., Ghalambor C.K., Haak D.C., Martin P.R. Impacts of climate warming on terrestrial ectotherms across latitude. Proc. Natl Acad. Sci. USA. 2008;105:6668–6672. doi:10.1073/pnas.0709472105 [PubMed]
  • Dunham A.E. Population responses to environmental change: physiologically structured models, operative environments, and population dynamics. In: Kareiva P.M., Kingsolver J.G., Huey R.B., editors. Biotic interactions and global change. Sinauer Associates; Sunderland, MA: 1993. pp. 95–119.
  • Feder M.E. Environmental variability and thermal acclimation of metabolism in tropical anurans. J. Therm. Biol. 1982;7:23–28. doi:10.1016/0306-4565(82)90015-8
  • Feeley K.J., Wright S.J., Supardi M.N.N., Kassim A.R., Davies S.J. Decelerating growth in tropical forest trees. Ecol. Lett. 2007;10:461–469. doi:10.1111/j.1461-0248.2007.01033.x [PubMed]
  • Ferraz G., Russell G.J., Stouffer P.C., Bierregaard R.O., Jr, Pimm S.L. Rapid species loss from Amazonian forest fragments. Proc. Natl Acad. Sci. USA. 2003;100:14069–14073. doi:10.1073/pnas.2336195100 [PubMed]
  • Garland T., Jr, Harvey P.H., Ives A.R. Procedures for the analysis of comparative data using phylogenetically independent contrasts. Syst. Biol. 1992;41:18–32. doi:10.2307/2992503
  • Garland T., Jr, Bennett A.F., Rezende E.L. Phylogenetic approaches in comparative physiology. J. Exp. Biol. 2005;208:3015–3035. doi:10.1242/jeb.01745 [PubMed]
  • Ghalambor C.K., Huey R.B., Martin P.R., Tewksbury J.J., Wang G. Are mountain passes higher in the tropics? Janzen's hypothesis revisited. Integr. Comp. Biol. 2006;46:5–17. doi:10.1093/icb/icj003 [PubMed]
  • Gilchrist G.W. Specialists and generalists in changing environments. I. Fitness landscapes of thermal sensitivity. Am. Nat. 1995;146:252–270. doi:10.1086/285797
  • Gorman G.C., Hillman S. Physiological basis for climatic niche partitioning in two species of Puerto Rican Anolis (Reptilia, Lacertilia, Iguanidae) J. Herpetol. 1977;11:337–340. doi:10.2307/1563246
  • Gorman G.C., Licht P. Seasonality in ovarian cycles among tropical Anolis lizards. Ecology. 1974;55:360–369. doi:10.2307/1935223
  • Helmuth B., Kingsolver J.G., Carrington E. Biophysics, physiological ecology, and climate change: does mechanism matter? Annu. Rev. Physiol. 2005;67:177–201. doi:10.1146/annurev.physiol.67.040403.105027 [PubMed]
  • Hertz P.E. Temperature regulation in Puerto Rican Anolis lizards: a field test using null hypotheses. Ecology. 1992;73:1405–1417. doi:10.2307/1940686
  • Hertz P.E., Huey R.B., Nevo E. Homage to Santa Anita: thermal sensitivity of sprint speed in agamid lizards. Evolution. 1983;37:1075–1084. doi:10.2307/2408420
  • Hertz P.E., Huey R.B., Stevenson R.D. Evaluating temperature regulation by field-active ectotherms: the fallacy of the inappropriate question. Am. Nat. 1993;142:796–818. doi:10.1086/285573 [PubMed]
  • Hoffmann A.A., Blows M.W. Evolutionary genetics and climate change: will animals adapt to global warming? In: Kareiva P.M., Kingsolver J.G., Huey R.B., editors. Biotic interactions and global change. Sinauer Associates; Sunderland, MA: 1993. pp. 165–178.
  • Huey R.B. Behavioral thermoregulation in lizards: importance of associated costs. Science. 1974;184:1001–1003. doi:10.1126/science.184.4140.1001 [PubMed]
  • Huey R.B. Temperature, physiology, and the ecology of reptiles. In: Gans C., Pough F.H., editors. Biology of the Reptilia. Physiology (C) vol. 12. Academic Press; London, UK: 1982. pp. 25–91.
  • Huey R.B. Natural variation in body temperature and physiological performance in a lizard. In: Rhodin A.G.J., Miyata K., editors. Advances in herpetology and evolutionary biology: essays in honor of Ernest E. Williams. Museum of Comparative Zoology; Cambridge, MA: 1983. pp. 484–490.
  • Huey R.B. Physiological consequences of habitat selection. Am. Nat. 1991;137:S91–S115. doi:10.1086/285141
  • Huey R.B., Bennett A.F. Phylogenetic studies of coadaptation: preferred temperatures versus optimal performance temperatures of lizards. Evolution. 1987;41:1098–1115. doi:10.2307/2409194
  • Huey R.B., Berrigan D. Temperature, demography, and ectotherm fitness. Am. Nat. 2001;158:204–210. doi:10.1086/321314 [PubMed]
  • Huey R.B., Kingsolver J.G. Evolution of resistance to high temperature in ectotherms. Am. Nat. 1993;142:S21–S46. doi:10.1086/285521
  • Huey R.B., Slatkin M. Costs and benefits of lizard thermoregulation. Q. Rev. Biol. 1976;51:363–384. doi:10.1086/409470 [PubMed]
  • Huey R.B., Webster T.P. Thermal biology of Anolis lizards in a complex fauna: the cristatellus group on Puerto Rico. Ecology. 1976;57:985–994. doi:10.2307/1941063
  • Huey R.B., Hertz P.E., Sinervo B. Behavioral drive versus behavioral inertia: a null model approach. Am. Nat. 2003;161:357–366. doi:10.1086/346135 [PubMed]
  • IPPC. Cambridge University Press; Cambridge, UK: 2007. Climate change 2007: the physical science basis.
  • Inger R.F. Temperature responses and ecological relations of two Bornean lizards. Ecology. 1959;40:127–136. doi:10.2307/1929930
  • Ives A.R., Midford P.E., Garland T., Jr Within-species variation and measurement error in phylogenetic comparative methods. Syst. Biol. 2007;56:252–270. doi:10.1080/10635150701313830 [PubMed]
  • Janzen D.H. Why mountain passes are higher in the tropics. Am. Nat. 1967;101:233–249. doi:10.1086/282487
  • Kearney M., Porter W.P. Mapping the fundamental niche: physiology, climate, and the distribution of a nocturnal lizard. Ecology. 2004;85:3119–3131. doi:10.1890/03-0820
  • Kearney M., Shine R., Porter W.P. The potential for behavioral thermoregulation to buffer ‘cold-blooded’ animals against climate warming. Proc. Natl Acad. Sci. USA. 2009 doi:10.1073/pnas.0808913106 [PubMed]
  • Kingsbury B.A. Thermal constraints and eurythermy in the lizard Elgaria multicarinata. Herpetologica. 1994;50:266–273.
  • Lavin S.R., Karasov W.H., Ives A.R., Middleton K.M., Garland T., Jr Morphometrics of the avian small intestine, compared with that nonflying mammals: a phylogenetic approach. Physiol. Biochem. Zool. 2008;81:536–550. doi:10.1086/590395 [PubMed]
  • Lister B.C. Seasonal niche relationships of rain forest anoles. Ecology. 1981;62:1548–1560. doi:10.2307/1941511
  • Malhi Y., Wright J. Spatial patterns and recent trends in the climate of tropical rainforest regions. Phil. Trans. R. Soc. B. 2004;359:311–329. doi:10.1098/rstb.2003.1433 [PMC free article] [PubMed]
  • Martin T.L., Huey R.B. Why suboptimal is optimal: Jensen's inequality and ectotherm thermal preferences. Am. Nat. 2008;171:E102–E118. doi:10.1086/527502 [PubMed]
  • Mautz W.J., Daniels C.B., Bennett A.F. Thermal dependence of locomotion and aggression in a xantusiid lizard. Herpetologica. 1992;48:271–279.
  • New M., Lister D., Hulme M., Makin I. A high-resolution data set of surface climate over global land areas. Clim. Res. 2002;21:1–25. doi:10.3354/cr021001
  • Parmesan C. Ecological and evolutionary responses to recent climate change. Annu. Rev. Ecol. Syst. 2006;37:637–669. doi:10.1146/annurev.ecolsys.37.091305.110100
  • Parmesan C. Influence of species, latitudes and methodologies on estimate of phenological response to global warming. Glob. Change Biol. 2007;13:1860–1872. doi:10.1111/j.1365-2486.2007.01404.x
  • Parsons P.A. Conservation and global warming: a problem in biological adaptation to stress. Ambio. 1989;18:322–325.
  • Porter W.P., James F.C. Behavioral implications of mechanistic ecology II: the African rainbow lizard, Agama agama. Copeia. 1979;1979:594–619. doi:10.2307/1443867
  • Porter W.P., Mitchell J.W., Beckman W.A., DeWitt C.B. Behavioral implications of mechanistic ecology. Thermal and behavioral modeling of desert ectotherms and their microenvironment. Oecologia. 1973;13:1–54. doi:10.1007/BF00379617
  • Pounds J.A., Fogden M.P.L., Campbell J.H. Biological response to climate change on a tropical mountain. Nature. 1999;398:611–615. doi:10.1038/19297
  • Rand A.S. Ecological distribution in anoline lizards of Puerto Rico. Ecology. 1964;45:745–752. doi:10.2307/1934922
  • Raxworthy C.J., Pearson R.G., Rabibisoa N., Rakotondrazafy A.M., Ramanamanjato J.-B., Raselimanana A.P., Wu S., Nussbaum R.A., Stone D.A. Extinction vulnerability of tropical montane endemism from warming and upslope displacement: a preliminary appraisal for the highest massif in Madagascar. Glob. Change Biol. 2008;14:1703–1720. doi:10.1111/j.1365-2486.2008.01596.x
  • Rogowitz G.L. Evaluation of thermal acclimation and altitudinal variation of metabolism in the Neotropical lizard, Anolis gundlachi. Copeia. 1996;1996:535–542. doi:10.2307/1447517
  • Root T.L., Price J.T., Hall K.R., Schneider S.H., Rosenzweig C., Pounds J.L. Fingerprints of global warming on wild animals and plants. Nature. 2003;421:37–42. doi:10.1038/nature01333 [PubMed]
  • Roughgarden J., Porter W.P., Heckel D. Resource partitioning of space and its relationship to body temperatures in Anolis lizard populations. Oecologia. 1981;50:256–264. doi:10.1007/BF00348048
  • Ruibal R. Thermal relations of five species of tropical lizards. Evolution. 1961;15:98–111. doi:10.2307/2405846
  • Santos M., Céspedes W., Balanyà J., Trotta V., Calboli F.C.F., Fontdevila A., Serra L. Temperature-related genetic changes in laboratory populations of Drosophila subobscura: evidence against simple climatic-based explanations for latitudinal clines. Am. Nat. 2005;165:258–273. doi:10.1086/427093 [PubMed]
  • Shine R. Egg-laying reptiles in cold climates: determinants and consequences of nest temperatures in montane lizards. J. Evol. Biol. 1999;12:918–926. doi:10.1046/j.1420-9101.1999.00093.x
  • Sinervo B. Evolution of thermal physiology and growth rate between populations of the western fence lizard (Sceloporus occidentalis) Oecologia. 1990;83:228–237. doi:10.1007/BF00317757 [PubMed]
  • Stevenson R.D. The relative importance of behavioral and physiological adjustments controlling body temperature in terrestrial ectotherms. Am. Nat. 1985;126:362–386. doi:10.1086/284423
  • Tewksbury J.J., Huey R.B., Deutsch C.A. Putting the heat on tropical animals. Science. 2008;320:1296–1297. doi:10.1126/science.1159328 [PubMed]
  • Tracy C.R., Christian K.A. Ecological relations among space, time and thermal niche axes. Ecology. 1983;67:609–615. doi:10.2307/1937684
  • Travis J., Futuyma D.J. Global change: lessons from and for evolutionary biology. In: Kareiva P.M., Kingsolver J.G., Huey R.B., editors. Biotic interactions and global change. Sinauer Associates; Sunderland, MA: 1993. pp. 251–263.
  • Tsuji J.S. Thermal acclimation of metabolism in Sceloporus lizards from different latitudes. Physiol. Zool. 1988;61:241–253.
  • van Berkum F.H. Latitudinal patterns of the thermal sensitivity of sprint speed in lizards. Am. Nat. 1988;132:327–343. doi:10.1086/284856
  • van Berkum F.H., Huey R.B., Adams B. Physiological consequences of thermoregulation in a tropical lizard (Ameiva festiva) Physiol. Zool. 1986;59:464–472.
  • Vitt L.J., Avila-Pires T.C.S., Caldwell J.P., Oliveira V.R.L. The impact of individual tree harvesting on thermal environments of lizards in Amazonian rain forest. Conserv. Biol. 1998;12:654–664. doi:10.1046/j.1523-1739.1998.96407.x
  • Vitt L.J., Shepard D.B., Vieira G.H.C., Caldwell J.P., Colli G.R., Mesquita D.O. Ecology of Anolis nitens brasiliensis in cerrado woodlands of Cantão. Copeia. 2008;2008:142–151. doi:10.1643/CP-06-251
  • Walther G.-R., Post E., Convey P., Menzel A., Parmesan C., Beebee T.J.C., Fromentin J.-M., Hoegh-Guldberg O., Bairlein F. Ecological responses to recent climate change. Nature. 2002;416:389–395. doi:10.1038/416389a [PubMed]
  • Whitfield S.M., Bell K.E., Philippi T., Sasa M., Bolaños F., Chaves G., Savage J.M., II, Donnelly M.A. Amphibian and reptile declines over 35 years at La Selva, Costa Rica. Proc. Natl Acad. Sci. USA. 2007;104:8352–8356. doi:10.1073/pnas.0611256104 [PubMed]
  • Williams S.E., Isaac J.L., Shoo L.P. The impact of climate change on the biodiversity and ecosystem functions of the wet tropics. In: Stork N.E., Turton S.M., editors. Living in a dynamic tropical forest landscape. Blackwell Publishing; Malden, MA: 2008. pp. 282–294.

Articles from Proceedings of the Royal Society B: Biological Sciences are provided here courtesy of The Royal Society