There is strong correlative evidence that human-induced climate warming is contributing to changes in the timing of natural events. Firm attribution, however, requires cause-and-effect links between observed climate change and altered phenology, together with statistical confidence that observed regional climate change is anthropogenic. We provide evidence for phenological shifts in the butterfly Heteronympha merope in response to regional warming in the southeast Australian city of Melbourne. The mean emergence date for H. merope has shifted −1.5 days per decade over a 65-year period with a concurrent increase in local air temperatures of approximately 0.16°C per decade. We used a physiologically based model of climatic influences on development, together with statistical analyses of climate data and global climate model projections, to attribute the response of H. merope to anthropogenic warming. Such mechanistic analyses of phenological responses to climate improve our ability to forecast future climate change impacts on biodiversity.
There is a strong signal that biological events have happened progressively earlier over the past few decades (Walther et al. 2002; Rosenzweig et al. 2008), including emergence in butterflies, migration in birds and fruiting and flowering in plants (Hughes 2000; Parmesan & Yohe 2003; Root et al. 2003; Chuine et al. 2004; Parmesan 2006; Cleland et al. 2007). These shifts have occurred concurrently with increases in air temperature at regional and global scales and in the direction expected if temperature was playing a causal role (Rosenzweig et al. 2008).
Despite the geographical and taxonomic generality of these signals, the magnitudes and rates of the observed phenological shifts have rarely been tied directly to known physiological responses. This leaves some uncertainty as to the driving forces behind these shifts and also limits the capacity for predictions into the future. Here, we examine historical phenological change in the common brown butterfly Heteronympha merope (Nymphalidae) and test whether (i) the phenological shift could be quantitatively explained by the influence of air temperature change on known physiological processes and (ii) the associated climate change could be attributed to human influence.
2. Material and methods
Heteronympha merope is an annual species that is abundant and widely distributed across southeastern Australia. Males and females emerge and mate in late spring/early summer (November/December), the males then die and the females become dormant until late summer/early autumn (February/March) after which they begin to oviposit on a wide range of grasses (Edwards 1973). Observed emergence times for H. merope in the vicinity of Melbourne, Victoria (37.60–38.54° S, 144.17–145.48° E) were based on opportunistically collected museum records and privately collected data courtesy of Kelvyn L. Dunn. We used the 10-year average of the earliest observed record per year from 1941–2005 as the emergence date. While the use of opportunistically collected data probably adds considerable noise to any signal of phenological shift, we accounted for any chronological bias by including the number of museum records as a covariate in the analysis of the trend through time.
The thermal dependence of development rate for eggs, larvae and pupae was determined from the offspring of 10 field-collected females, raised individually on the grass Ehrharta erecta (Poaceae) from egg to emergence in glass vials (height 130 mm, diameter 25 mm) sealed with a moistened foam stopper. Temperature was controlled at 8°C, 12°C, 15°C, 25°C or 30°C with a 12 : 12 light : dark photoperiod. We observed the animals daily and recorded time to completion of each life cycle stage, the inverse of this being the developmental rate. We then fitted polynomial regressions to the relationship between development rate and temperature for each life stage, using the adjusted r2 method (Quinn & Keough 2002).
We used historical weather data for 1945–2007 (Bureau of Meteorology, Australia) from Laverton (37.86° S, 144.76° E), a rural site close to Melbourne, to model the physiological response of H. merope to temperature. This station is a ‘high-quality’ site, unaffected by changes in exposure, urbanization, instrumentation, etc., during the study period. Weather records (mean monthly maximum and minimum air temperature, wind speed and cloud cover) were translated into microclimates experienced by immature H. merope using biophysical modelling software (Niche Mapper, http://www.zoology.wisc.edu/faculty/Por/Por.html#niche; Porter & Mitchell 2006). This software includes a microclimate model that translates weather-station records into near-surface air temperature and wind speed profiles and a soil surface energy balance (Porter et al. 1973). The model includes a first-principles model of solar radiation (McCullough & Porter 1971) and uses cloud cover records to predict long- and short-wave radiation loads. Individuals were assumed to be 3 cm above the ground in grass tussocks (75% shade), based on our observations of larval behaviour. Microclimate model predictions compared well with field microclimate measurements (figure S1, electronic supplementary material).
Animals were assumed to be at the 3 cm air temperature. Predicted daily cycles in animal body temperature were converted to constant temperature equivalents (CTEs) (as implemented in Mitchell et al. 2008), which were then used as independent variables in the fitted equations for thermal dependence of development rate. There was a good agreement between modelled CTE-based development times and the development times measured under natural, fluctuating thermal conditions in grass tussocks in Melbourne over the winter of 2007 (table S1, electronic supplementary material). We encoded the equations for development rate and CTE calculations into the Ectotherm model of Niche Mapper. The Ectotherm program tracked developmental stage on an hourly basis given a specified laying date (11 April—females lay from around March to April in Melbourne; Pearse 1978) and outputted the time of completion of development based on observed meteorological data from 1944 to 2005. We calculated separate emergence dates for each year and then took the average value for each 10-year window as the final emergence date. Note that we did not solve an energy balance specifically for the caterpillars with the Ectotherm model but rather assumed that they were at the shaded 3 cm air temperature.
We compared observed temperature trends from the high-quality weather station in Laverton with output from extended climate model simulations, both including and excluding anthropogenic climate forcing for the single-model grid box overlying Melbourne and Laverton. Anthropogenic climate forcing included observed increases in greenhouse gases and estimated variations of anthropogenic aerosols, whereas natural external climate forcing included estimated changes in solar irradiance and volcanic aerosols. Multi-member ensembles of simulations from four different climate models with prescribed changes in both anthropogenic and natural external climate forcing were used to provide regional temperature data for 1944–2007 (table S2, electronic supplementary material). The range of possible trends owing to natural internal climate variability was estimated using the variability of regional temperature from extended control model simulations (including only natural climate variation with no changes in external forcing).
First spring observations for H. merope in the vicinity of Melbourne show a trend in earliest emergence date of −1.4 days per decade over the last 65 years (figure 1a). This trend is statistically significant (linear regression: R2 = 0.766, year p = 0.022, museum record number p = 0.046) and occurs coincident with a regional air temperature increase of 0.16°C per decade during the developmentally important period (R2 = 0.893, p = 0.001; figure 1b).
The thermal response of development rate varies considerably across the life stages (no fifth instars survived at constant temperatures at or above 25°C) (figure 1c). Combining the physiological functions of thermal sensitivity with the microclimate model and historical monthly climate data, we predicted a shift in emergence date of −1.5 days per decade (R2 = 0.866, p = 0.024), close and statistically indistinguishable from the observed rate of −1.4 days per decade (heterogeneity of slopes test, p > 0.944; figure 1a). Importantly, the physiological model provided a better prediction of emergence date (R2 = 0.652, p = 0.028) than did air temperature (R2 = 0.426, p = 0.112). While the modelled H. merope temperatures are mainly a function of radiation, wind speed and air temperature, sensitivity analyses showed that most (80.6%) of the predicted shift in emergence date was driven by changes in air temperature. Very little (less than 1%) was independently attributable to wind speed and cloud cover, but the predicted effects of the latter two variables were inversely correlated (r = −0.65, p = 0.01).
Simulations from multi-member ensembles of four different general circulation models for the Melbourne area show that the air temperature trend coincident with the phenological shift in H. merope was very unlikely to be a result of natural internal climate fluctuations (figure 1d). Changes in natural external forcing, owing to changes in solar irradiance and volcanic aerosols, are likely to have led to a small cooling over this period (Karoly & Braganza 2005) and are very unlikely to explain the observed warming. The observed regional warming trend is consistent with the modelled climate response in this region to increasing greenhouse gases and other anthropogenic climate forcing (figure 1d).
The immature life stages of H. merope differ strongly in their thermal sensitivities, corresponding to seasonal changes in temperature during the immature period. Combining these physiological responses with historical weather records via the microclimate model, we predicted a shift towards earlier emergence of the same magnitude as that observed over the same time period. This very strongly suggests that the observed phenological shift in H. merope in the Melbourne area is a direct result of climate change impacts on the development rate.
In our model, egg, caterpillar and pupal temperatures were a function of air temperature and also of radiation (via cloud cover) and wind speed. Cloud cover and wind speed variation did explain part of the predicted shift in emergence date when considered in isolation but had no net effect because they acted in opposite directions. As a consequence, air temperature shifts were seen to dominate the predicted phenological response in this particular case, but it should be clear that this need not necessarily be so. In other cases, concurrent changes in other environmental variables may suppress or exacerbate the effects of air temperature and such complexities can only be considered through a mechanistic analysis of microclimatic impacts on phenology.
The observed shift in air temperature of 0.16°C per decade in the vicinity of Melbourne can very likely be explained through the effects of greenhouse gases emitted by humans (figure 1d). Our analysis thus provides direct causal linkages between the emission of greenhouse gases by humans, a shift in local air temperature, and the physiological response of a butterfly resulting in earlier spring emergence.
Phenological shifts may represent significant environmental challenges, particularly for species of low vagility, or provide opportunities for a species' persistence (Harrington et al. 1999; Visser & Holleman 2001). Our study illustrates how it is possible to attribute phenological shifts to anthropogenic climate change based on physiological data, a model of how the organism experiences climatic conditions from a microclimatic perspective and the likelihood that locally observed climate change is human-caused. A greater understanding of the physiological links between changes in climate and phenology will strengthen our interpretation of past phenological change. It will also allow better forecasts of future shifts in phenology and their consequences (Ibáñez et al. 2006; Bradshaw & Holzapfel 2008; Kearney & Porter 2009).
We thank Chris Thomas, Lesley Hughes, K. McClellan, Cheryl O'Dwyer, Nigel Stork, Anna Lister and Neil Murray for discussion/comments on the manuscript. We especially thank Mark Hendrickx whose critique of a previous version of the article significantly improved the final version. This study was supported by the Australian Research Council grants DP0772837 and FF0668679, Parks Victoria and Department of Sustainability and Environment, and by an Australian Department of Climate Change grant to N.J.B.
- Received January 19, 2010.
- Accepted February 26, 2010.
- © 2010 The Royal Society