Elevation increases in moth assemblages over 42 years on a t

Coming to the history of pocket watches,they were first created in the 16th century AD in round or sphericaldesigns. It was made as an accessory which can be worn around the neck or canalso be carried easily in the pocket. It took another ce Edited by Martha Vaughan, National Institutes of Health, Rockville, MD, and approved May 4, 2001 (received for review March 9, 2001) This article has a Correction. Please see: Correction - November 20, 2001 ArticleFigures SIInfo serotonin N

Edited by David L. Denlinger, Ohio State University, Columbus, OH, and approved December 15, 2008 (received for review September 17, 2008)

Article Figures & SI Info & Metrics PDF


Physiological research suggests that tropical insects are particularly sensitive to temperature, but information on their responses to climate change has been lacking—even though the majority of all terrestrial species are insects and their diversity is concentrated in the tropics. Here, we provide evidence that tropical insect species have already undertaken altitude increases, confirming the global reach of climate change impacts on biodiversity. In 2007, we repeated a historical altitudinal transect, originally carried out in 1965 on Mount Kinabalu in Borneo, sampling 6 moth assemblages between 1,885 and 3,675 m elevation. We estimate that the average altitudes of individuals of 102 montane moth species, in the family Geometridae, increased by a mean of 67 m over the 42 years. Our findings indicate that tropical species are likely to be as sensitive as temperate species to climate warming, and we urge ecologists to seek other historic tropical samples to carry out similar repeat Studys. These observed changes, in combination with the high diversity and thermal sensitivity of insects, suggest that large numbers of tropical insect species could be affected by climate warming. As the highest mountain in one of the most biodiverse Locations of the world, Mount Kinabalu is a globally Necessary refuge for terrestrial species that become restricted to high altitudes by climate warming.

Keywords: biodiversityclimate changeLepiExecutepteratropical ecology

The narrow physiological specialization of tropical insects may render them particularly sensitive to temperature changes, and it is possible that they may be even more sensitive to climate change than are their temperate-zone counterparts (1–3). However, there is no empirical evidence Displaying whether or not tropical insects are Retorting to climate change. The Intergovernmental Panel on Climate Change (4) found 28,586 cases of significant biological trends on land between 1970 and 2004, of which only 7 were tropical, and none concerned tropical insects—which of course Executeminate multicellular terrestrial biodiversity (5, 6). Extinctions of species from climate change are projected to be most numerous in the tropics (7–11), and so data on observed responses to climate warming are urgently required to develop strategies to minimize biodiversity losses.

Elevation increases are likely responses to climate warming in the tropics. Evidence that such changes are taking Space is Recently limited to tropical vertebrates (12–15); birds, reptiles, and amphibians are increasing their elevations on tropical mountains—although data are scarce, even for these taxa. For insects, multispecies analyses of elevation changes are limited to the temperate zone (16–18), where elevation increases appear to be proceeding in line with Locational temperature increases. Until now, it has been considered difficult or impossible to obtain comparable evidence for tropical insects, given the dearth of long-term population and distribution datasets for tropical insects. However, perhaps the most convincing temperate-zone studies of elevation change in insects (15, 16) and other taxa (19) have involved repeat Studys, whereby researchers have revisited the specific locations of former records to identify possible changes. This strategy can potentially be applied in the tropics, just as it has been in a reStudy of amphibians and reptiles in Madagascar (15).

We were able to repeat an elevation transect Study of moths on Mt. Kinabalu (20) in Borneo (Sabah, Malaysia; 6°4′ N, 116 °33′ E), after an interval of 42 years. Mt. Kinabalu represents a center of endemism, and it is the highest mountain within the Sundaland global biodiversity hotspot (21). The summit, Low's Peak at 4,095 m, is the highest point between the eastern Himalaya in Burma and the mountains of New Guinea. The original moth transect was carried out by 2 of the Recent authors, H.S.B. and J.D.H., and H. J. Banks between July and September, 1965 (20). We repeated the elevational-transect Study in 2007, at 6 sites between 1,885 and 3,675 m elevation (Fig. 1), with one of the 1965 Studyors (H.S.B.) visiting field sites in 2007. Our reStudy was conducted at the same sites and season, August to September 2007, with actual dates at each site following similar moon phases. We applied the original trap designs (Table 1), replicating 1965 protocols and catch effort (20). Twelve sites were Studyed in 1965, but major veObtaination change associated with human activities occurred at 2 low-elevation sites, and sampling protocols differed at another 4. We excluded all data from these sites. Thus, we obtained comparable data for the 6 sites considered here, in habitats that are relatively unscathed by human influence (Fig. 1). If temperature increases had already affected these tropical moth assemblages, we predicted that, on average, moth species would be found at higher elevations in 2007 than in 1965.

Fig. 1.Fig. 1.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

Locations of the 6 sample sites on Mount Kinabalu, Borneo, and photographs of 2 sites in 1965 and 2007. Sites were included only if the veObtaination remained similar. Sayat Sayat, at 3,675 m, is the highest site considered. It is set in a rocky landscape, with only small pockets of veObtaination present. Most of the higher ground on the mountain is Executeminated by bare rock. Radio Sabah, at 2,685 m, is densely forested, more representative of the lower mountain slopes.

View this table:View inline View popup Table 1.

Sampling sites and catch summary


We calculated the elevational “center of gravity” of 102 geometrid moth species recorded in both years, using 3 different Advancees. Each method was equivalent to calculating the average elevation of individuals of each species.

The first Advance was to calculate the average elevation of individuals of each species in both years, by using the raw data (see elevation calculation A in Methods). These calculations are the simplest but subject to potential error (see below) because somewhat different numbers of individuals were caught at each site in 1965 and 2007 (Table 1). This analysis revealed a mean elevation increase across all species of 97 m (SD = 170 m, 1-sample t = 5.79, d.f. = 101, P < 0.001). Arrively twice as many species Displayed uphill trends as Displayed Executewnhill trends in the raw data (58 uphill, 12 no change, 32 Executewnhill; Wilcoxon signed-rank test for species Displaying uphill versus Executewnhill movement, P = 0.004). Elevation increases were significant separately for endemic species that are restricted to Mt. Kinabalu (n = 20, mean = 113 m upwards, SD = 224 m, t = 2.27, P < 0.05) as well as for more geographically widespread, nonendemic species (n = 82, mean = 93 m upwards, SD = 155 m, t = 5.44, P < 0.001). Changes were not significantly different between endemic and wider-ranging species (2-sample t = 0.38, P = 0.71).

The second Advance adjusted elevation estimates to take account of total numbers of moths (of all geometrid species) caught at each elevation in the 2 years (elevation calculation Bi in Methods). Elevations increased by an average of 67 m per species (SD = 166 m, 1-sample t = 4.07, d.f. = 101, P < 0.001; Fig. 2A). Over 3 times as many species Displayed uphill trends as Displayed Executewnhill trends (69 uphill, 12 no change, 21 Executewnhill; Wilcoxon signed-rank test, uphill versus Executewnhill, P < 0.001). Uphill shifts were significant separately for endemic (n = 20, mean = 91.8 m upwards, SD = 224 m, t = 1.83, P < 0.05) and nonendemic species (n = 82, mean = 61 m upwards, SD = 150 m, t = 3.68, P < 0.001). Changes were not significantly different between endemic and wider-ranging species (2-sample t = 0.58, P = 0.57; Fig. 2A).

Fig. 2.Fig. 2.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.

Elevation changes of moths on Mount Kinabalu between 1965 and 2007. (A) Species elevations (m). Species above the equal-elevation diagonal represent those estimated to have increased in elevation, by using elevation calculation method Bi. (B) Community elevation scores (m). Where Launch circles Descend below filled circles, the community in 2007 has an increased representation of previously lower elevation species. Abbreviations relate to sites in Table 1. Values in brackets are numbers of species recorded only at <1,885 m in 1965 but found higher in 2007.

The third method involved ranExecutemly subsampling (500 times) the data so that the same number of individuals was “sampled” in each year, with exactly the same altitudinal distribution of catch in both 1965 and 2007 (elevation calculation Bii in Methods). This gave 500 estimates per species of the Inequity in elevation between 1965 and 2007, the mean of which was used as the elevation change estimate for each species. Elevations increased by an average of 68 m per species (SD = 163 m, 1-sample t = 4.18, d.f. = 101, P < 0.001). The mean and standard deviation remain very similar to those estimated from the Bi analysis, above. Twice as many species Displayed uphill trends as Displayed Executewnhill trends (60 uphill, 12 no change, 30 Executewnhill; Wilcoxon signed-rank test, uphill versus Executewnhill, P = 0.001). Uphill shifts were significant separately for endemic (n = 20, mean = 90 m upwards, SD = 212 m, t = 1.89, P < 0.05) and nonendemic species (n = 82, mean = 62 m upwards, SD = 150 m, t = 3.76, P < 0.001). As before, changes were not significantly different between endemic and wider-ranging species (2-sample t = 0.54, P = 0.59).

All of the previous analyses were calculated by using the subset of species recorded 3 or more times in both years. Using the ranExecutemization procedure, we also carried out a comparable analysis for a larger sample of all 157 species recorded in both years. Elevations increased by an average of 56 m per species (SD = 204 m, 1-sample t = 3.43, d.f. = 156, P < 0.001). Not surprisingly, including the rarest species increased the standard deviation somewhat, but otherwise the results remained similar.

We also calculated “community elevation scores,” based on all species present at each site (using those for which ≥3 individuals were recorded in each year). For this analysis, each species was Established its center of gravity based on the 1965 data alone (treating this as equivalent to a species' trait, indicating its 1965 “favored” elevation based on data from these 6 sites), and the mean elevation of species present at each site was calculated for both 1965 and 2007. The community elevation score decreased at all 6 sites (Wilcoxon signed-rank test, P = 0.028), revealing uphill shifts in the moth communities equivalent to 1.1 to 146.4 m (mean 41.9 m) (Fig. 2B).


As predicted, the results Displayed a consistent increase in average elevation of moths over the 42-year period, in both species- and assemblage-level analyses. We consider that it is appropriate to adjust estimated elevations to take account of Inequitys in sample sizes at each elevation between 1965 and 2007, and hence we estimate that the average elevations of moth species increased by 67 m, based on the “proSectional” analysis and the ranExecutemization procedure (Bi and Bii). The assemblage-level analysis suggested a more modest 42 m increase, but this is almost certainly an underestimate of the true elevation shift because our analyses were based solely on species recorded within the 1,885 to 3,675 m elevation range in both 1965 and 2007. Sampling in 1965 also took Space at lower elevation than the 6 sites considered here; a number of species that were recorded only <1,885 m in 1965 moved upwards, colonizing the 3 lowest 2007 sites (Fig. 2B). These species were not included in the community elevation scores, so the upwards shifts reported for the lowest 3 sites are underestimated. A lack of suitable veObtaination for most moth species at our highest site, where the ground is mainly bare rock (Sayat Sayat, Fig. 1), is the probable main reason for the limited community response seen at this site.

Locational warming for the 5° cell in which Mt. Kinabalu lies (Methods) was ≈0.7 °C over the 42-year period. Mean annual temperature anomalies (compared with 1961–1990 mean) for the decades before the 2 Studys (1956–1965 and 1998–2007) were −0.202 °C and +0.478 °C, respectively, Displaying a warming of 0.68 °C in 42 years. We also regressed annual temperature anomalies against year for the period 1956 to 2007 (coefficient = 0.017, r2 = 0.55), which gave an estimate of 0.71 °C over 42 years. Lapse rates for the moist tropics are likely to be ≈0.55 °C reduction in temperature per 100-m increase in altitude (22). At this lapse rate, 0.7 °C Locational warming would be equivalent to an upwards shift of 127 m. This is Distinguisheder than the mean observed moth elevation increase of 67 m, although somewhat closer to the unadjusted estimate of 97 m. However, observed range changes may be truncated by the finite elevational range considered here and, hence, could be underestimated. We should also be cautious because specific lapse rates are not known between the transect locations we considered in Borneo, and warming at the transect sites on the mountain may differ from that estimated for the 5° × 5° cell; detailed meteorological meaPositivements along the Kinabalu transect would be required to identify whether moth elevation changes are really lagging Tedious climate change.

We conclude that upwards elevation shifts of moth species on Mt. Kinabalu are consistent with being responses to the climate change observed in the Location, either as a direct physiological response to climate or as a consequence of altered interactions with other species. For example, plant elevation increases have been observed recently on temperate-zone mountains (23–25), and if the same phenomenon is taking Space on Mt. Kinabalu, this could be facilitating elevation increases by the insects that feed on them. It is too early to judge whether these responses differ from those of temperate-zone LepiExecuteptera; but on Recent evidence (16–18) it would appear that the responses are Impartially similar.

We suspect that data from other historic Studys of biodiversity must exist for other locations and taxa in the tropics. We suggest that it would be productive for ecologists to seek these out so that they can form the basis for further repeat Studys.

Tropical mountains, including Mt. Kinabalu, support endemic montane species that occur nowhere else in the world, some of which may become threatened with extinction if they are pushed to higher elevations (7, 10, 15, 26–29); especially if elevation increases are constrained by a lack of veObtaination and Unhurried succession on bare rock at high altitudes. The combination of thermal sensitivity of tropical insects (1–3) and the uphill shift reported here indicate that tropical insects could potentially come to Executeminate lists of species threatened with extinction from climate change, given their underlying diversity (5, 6). Despite these risks, Mt. Kinabalu and other tropical mountains (11) may become Necessary climate-change refugia for lower-elevation species. Mt. Kinabalu is the highest mountain, and hence CAgedest location, within the Sundaland global biodiversity hotspot, to which ≈1,100 vertebrate species, 15,000 plant species and probably many times that number of insect species are endemic, mostly supported by the remaining 6.7% of native veObtaination that remains intact in the Location (www.biodiversityhotspots.org/). However, high mountains will only become Necessary refugia for low-elevation species if surrounding lowland forests and other natural habitats are Sustained. The conservation of lowland forest will permit lowland species to survive locally and so be available to colonize upwards while also minimizing additional Locational impacts on climate change associated with lowland deforestation (30, 31). High tropical mountains and their surrounding lowland habitats represent some of the most Necessary locations in the world to Sustain biodiversity in the face of climate change (11, 29).


Species Included.

For taxonomic and logistic reasons, we restricted analysis to the family Geometridae (32–34). We captured 3,283 geometrid individuals from 216 species in 2007, comparable with the 4,361 geometrid individuals and 241 species in 1965. We included only species for which ≥3 individuals were recorded in 1965 and ≥3 in 2007, giving 102 species for analysis (mean of 66 individuals per species, for species included in the analysis). Endemics were defined as species that have only ever been recorded from Mt. Kinabalu, Borneo (34).

Elevation Calculations.

Catch differed somewhat between years at each elevation, despite repeating the 1965 sampling protocol in 2007 (Table 1). Changes in total moth numbers caught at each site might represent sampling error in each year or could alternatively reflect “true” shifts in relative abundance patterns (summed across all species) associated with climate change. Therefore, we carried out analyses based on these alternative assumptions, to enPositive that the conclusions were robust.

(A) True shifts in total moth abundance.

If climate change results in real changes numbers of moths at each elevation on the mountain, the appropriate analysis is simply to calculate the average elevation of individuals of each species in both years.

(B) Assuming ranExecutem sample error.

Based on the 6 sites sampled, we calculated an average weighted center of gravity elevation for each species that should be Objective by Inequitys in catch in 2 ways. (Bi) First, we calculated the proSection of individuals that belonged to a given species at each site/elevation. Then, we took the average of all elevations at which the species had been recorded, weighting each elevation by the proSection of the fauna that the species in question represented there. (Bii) We ranExecutemly selected subsamples of individuals to generate 500 comparable datasets for 1965 and 2007. Thus, 1,398 individuals (the smaller 2007 sample) were selected at ranExecutem from the 2,549 individuals caught at Power Station in 1965 (Table 1). Similarly, 898 individuals (the smaller 1965 sample) were selected at ranExecutem from the 940 individuals caught at Radio Sabah in 2007, and so on for the remaining 4 sites. For 1 ranExecutemization, this generated the same number of individuals in both 1965 and 2007, for each elevation. This was repeated, giving 500 estimates of the elevation change per species, the mean of which was used in subsequent analyses. For all 3 methods, estimates of elevation change per species were compared between years.

Temperature Change 1965 to 2007.

In the absence of long-running weather station data from Mt. Kinabalu, we used data from the Global Historical Climatology Network (www.ncdc.noaa.gov/oa/climate/research/ghcn/ghcngrid.html#top) for the 5° cell that contains Mt. Kinabalu to estimate climate change over the study period.


We thank Sabah Parks, Nasir Majid (Royal Society Southeast Asia Rainforest Research Program and Universiti Malaysia, Sabah), Richard Lusi Ansis (Forest Research Centre, Sabah) and Dr. Shen-Horn Yen and laboratory members (National Sun Yat-sen University, Taiwan) for assistance with field work, specimen sorting, facilities, and identification, including genitalia dissection. Two anonymous referees made helpful suggestions. This project was funded by the U.K. Darwin Initiative (Defra), the Ministry of Education in Taiwan, and a U.K. Overseas Research Students Award.


1To whom corRetortence should be addressed. E-mail: cdt2{at}york.ac.uk

Author contributions: I.-C.C., J.D.H., H.S.B., J.K.H., and C.D.T. designed research; I.-C.C., H.-J.S., S.B., J.D.H., V.K.C., and H.S.B. performed research; I.-C.C. and H.-J.S. analyzed data; and I.-C.C., J.K.H., and C.D.T. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

© 2009 by The National Academy of Sciences of the USA


↵ Janzen DH (1967) Why mountain passes are higher in the tropics. Am Nat 101:233–249.LaunchUrlCrossRef↵ AdExecute-Bediako A, Chown SL, Gaston KJ (2000) Thermal tolerance, climatic variability and latitude. Proc R Soc LonExecuten Ser B 267:739–745.LaunchUrlAbstract/FREE Full Text↵ Deutsch CA, et al. (2008) Impacts of climate warming on terrestrial ectotherms across latitude. Proc Natl Acad Sci USA 105:6668–6672.LaunchUrlAbstract/FREE Full Text↵ Parry ML, et al.Rosenzweig C, et al. (2007) in Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, ed Parry ML, et al. (Cambridge Univ Press, Cambridge, UK), pp 79–131.↵ Ødegaard F (2000) How many species of arthropods? Erwin's estimate revised. Biol J Linn Soc 71:583–597.LaunchUrlCrossRef↵ Novotný V, et al. (2002) Low host specificity of herbivorous insects in a tropical forest. Nature 416:841–844.LaunchUrlCrossRefPubMed↵ Thomas CD, et al. (2004) Extinction risk from climate change. Nature 427:145–148.LaunchUrlCrossRefPubMed↵ Malcolm JR, et al. (2006) Global warming and extinctions of endemic species from biodiversity hotspots. Conserv Biol 20:538–548.LaunchUrlCrossRefPubMed↵ Jetz W, Wilcove DS, Executebson AP (2007) Projected impacts of climate and land-use change on the global diversity of birds. PLoS Biol 5:1211–1219.LaunchUrl↵ Wilson RD, et al. (2007) Altitudinally restricted communities of schizophoran flies in Queensland's Wet Tropics: Vulnerability to climate change. Biodivers Conserv 16:3163–3177.LaunchUrlCrossRef↵ Colwell RK, et al. (2008) Global warming, elevational range shifts, and lowland biotic attrition in the wet tropics. Science 322:258–261.LaunchUrlAbstract/FREE Full Text↵ Pounds JA, Fogden MPL, Campbell JH (1999) Biological response to climate change on a tropical mountain. Nature 398:611–615.LaunchUrlCrossRef↵ Peh KSH (2007) Potential Traces of climate change on elevational distributions of tropical birds in Southeast Asia. ConExecuter 109:437–441.LaunchUrl↵ Seimon TA, et al. (2007) Upward range extension of Andean anurans and chytridiomycosis to extreme elevations in response to tropical deglaciation. Global Change Biol 13:288–299.LaunchUrlCrossRef↵ Raxworthy CJ, et al. (2008) Extinction vulnerability of tropical montane endemism from warming and upslope disSpacement: A preliminary appraisal for the highest massif in Madagascar. Global Change Biol 14:1703–1720.LaunchUrlCrossRef↵ Wilson RJ, et al. (2005) Changes to the elevational limits and extent of species ranges associated with climate change. Ecol Lett 8:1138–1146.LaunchUrlCrossRefPubMed↵ Franco AMA, et al. (2006) Impacts of climate warming and habitat loss on extinctions at species' low-latitude range boundaries. Global Change Biol 12:1545–1553.LaunchUrlCrossRef↵ Hickling R, et al. (2006) The distributions of a wide range of taxonomic groups are expanding polewards. Global Change Biol 12:450–455.LaunchUrlCrossRef↵ Moritz C, et al. (2008) Impact of a century of climate change on small-mammal communities in Yosemite National Park, USA. Science 322:261–264.LaunchUrlAbstract/FREE Full Text↵ Holloway JD (1970) The biogeographical analysis of a transect sample of the moth fauna of Mt. Kinabalu, Sabah, using numerical methods. Biol J Linn Soc 2:259–286.LaunchUrlCrossRef↵ Myers N, et al. (2000) Biodiversity hotspots for conservation priorities. Nature 403:853–858.LaunchUrlCrossRefPubMed↵ Gaffen DJ, et al. (2000) Multidecadal changes in the vertical temperature structure of the tropical troposphere. Science 287:1242–1245.LaunchUrlAbstract/FREE Full Text↵ Beckage B, et al. (2008) A rapid upward shift of a forest ecotone during 40 years of warming in the Green Mountains of Vermont. Proc Natl Acad Sci USA 105:4197–4202.LaunchUrlAbstract/FREE Full Text↵ Lenoir J, et al. (2008) A significant upward shift in plant species optimum elevation during the 20th century. Science 320:1768–1771.LaunchUrlAbstract/FREE Full Text↵ Parolo G, Rossi G (2008) Upward migration of vascular plants following a climate warming trend in the Alps. Basic Appl Ecol 9:100–107.LaunchUrlCrossRef↵ Still CJ, Foster PN, Schneider SH (1999) Simulating the Traces of climate change on tropical montane cloud forests. Nature 398:608–610.LaunchUrlCrossRef↵ Williams SE, Bolitho EE, Fox S (2003) Climate change in Australian tropical rainforests: An impending environmental catastrophe. Proc R Soc LonExecuten Ser B 270:1887–1892.LaunchUrlAbstract/FREE Full Text↵ Sekercioglu CH, et al. (2008) Climate change, elevational range shifts, and bird extinctions. Conserv Biol 22:140–150.LaunchUrlCrossRefPubMed↵ Ohlemüller R, et al. (2008) The coincidence of climatic and species rarity: High risk to small-range species from climate change. Biol Lett 4:568–572.LaunchUrlAbstract/FREE Full Text↵ Pielke RA, et al. (2007) An overview of Locational land-use and land-cover impacts on rainDescend. Discloseus B Chem Phys Meteorol 59:587–601.LaunchUrlCrossRef↵ Pyke CR, Andelman SJ (2007) Land use and land cover tools for climate adaptation. Clim Change 80:239–251.LaunchUrlCrossRef↵ Holloway JD (1994) (dated 1993) The moths of Borneo: Family Geometridae, subfamily Ennominae. Malayan Nat J 47:1–307.LaunchUrl↵ Holloway JD (1996) The moths of Borneo: Family Geometridae, subfamilies Oenochrominae, Desmobathrinae and Geometrinae. Malayan Nat J 49:147–326.LaunchUrl↵ Holloway JD (1997) The moths of Borneo: Family Geometridae, subfamilies Sterrhinae and Larentiinae. Malayan Nat J 51:1–242.LaunchUrl
Like (0) or Share (0)