Environmental change and long-term body mass declines in an alpine mammal
Frontiers in Zoologyvolume 11, Article number: 69 (2014)
The Erratum to this article has been published in Frontiers in Zoology 2014 11:88
Climate and environmental change have driven widespread changes in body size, particularly declines, across a range of taxonomic groups in recent decades. Size declines could substantially impact on the functioning of ecosystems. To date, most studies suggest that temporal trends in size have resulted indirectly from climate change modifying resource availability and quality, affecting the ability of individuals to acquire resources and grow.
Here, we investigate striking long-term body mass declines in juvenile Alpine chamois (Rupicapra rupicapra), within three neighbouring populations in the Italian Alps. We find strong evidence that increasing population density and warming temperatures during spring and summer are linked to the mass declines. We find no evidence that the timing or productivity of resources have been altered during this period.
We conclude that it is unlikely that environmental change has driven body size change indirectly via effects on resource productivity or phenology. Instead, we propose that environmental change has limited the ability of individuals to acquire resources. This could be due to increases in the intensity of competition and decreases in time spent foraging, owing to high temperatures. Our findings add weight to a growing body of evidence for long-term body size reductions and provide considerable insight into the potential drivers of such trends. Furthermore, we highlight the potential for appropriate management, for instance increases in harvest size, to counteract the impacts of climate change on body mass.
Climatic and other environmental change has impacted species in a variety of ways, from altering their spatial distributions (e.g. ) to changing the timing of their annual events (e.g. ). Recently, focus has broadened to include the impacts of climate change on life-history traits, body condition and population processes (e.g. -). One emergent generality is that responses to climate change include widespread body size changes, particularly declines, which could have pronounced negative impacts on the functioning and productivity of ecosystems ,. Body size declines driven by climate change have been recorded in the past, for instance in large ungulates during the Pleistocene, and are thought to have led to the extinction of some species .
A variety of climatic drivers of recent body size declines have been proposed. In ectotherms, higher metabolic rates are predicted in warmer environments , so climate change might lead to decreased body size, unless individuals can increase their rate of food intake . In endotherms, in line with Bergmann's rule , it could be beneficial to be smaller (and thus have a larger surface area to volume ratio) in warmer environments, due to a reduced need for heat conservation and a greater need for heat loss . In mammals and birds, intra-specific trends between body size and both latitude , and temperature , have been observed, providing some support for this theory. However, the most frequently cited climatic driver of body mass declines is the indirect link via climatic impacts on resource availability, which has been implicated considerably more frequently than any other mechanism (reviewed in ). Climate and environmental change can alter the timing, quality and quantity of resources, affecting the rate at which individuals acquire resources to invest in growth and energy storage, and ultimately body size (e.g. ,). Non-climatic processes can also drive body size change and could be playing a role alongside climate change. For example, in harvested animal populations, a preference for larger bodied individuals by hunters can drive selection for decreased body size ,.
In ungulates, body mass is an important indicator of fitness , and can respond rapidly to environmental change . High body mass is commonly related to high reproductive success (e.g. ) and survival (e.g. ). As such, changes in body mass can have important effects on population dynamics ,. The indirect link between climate, resources and ungulate body mass is well studied and could be strongly influenced by climate change ,. However, environmental change could also directly affect the ability of ungulates to acquire resources, particularly in seasonal environments, which could lead to temporal body mass change. For instance, changes in population density can alter the intensity of competition for resources, affecting the ability of individuals to forage and grow ,. There is some evidence from other taxa of links between density and temporal body size change . In recent decades, warming climates have led to changing snow cover and depth in some areas , altering the costs of locomotion and foraging for some species . In alpine species, high temperatures in spring and summer can cause heat-stress, reducing the time individuals can spend foraging . As a result, we propose that, in alpine areas, which have been strongly affected by climate warming , reductions in time spent foraging due to higher thermoregulatory costs could lead to reduced body size. Climate change is predicted to drive similar changes in temperature-dependent activity budgets in other taxa ,. Furthermore, temperature induced changes in foraging behaviour have been detected in experimental studies .
Here, we explore long-term variation in the body mass of juveniles in three neighbouring populations of Alpine chamois (Rupicapra rupicapra) in the Italian Alps. The body condition of juvenile ungulates is particularly responsive to environmental variation ,,. This is largely because juveniles invest highly in growth, and not in reproduction, meaning that their body condition is very sensitive to the availability of resources. Experiencing poor environmental conditions in youth can suppress growth rates and result in smaller body size at maturity , which could have important implications for population dynamics. We begin by investigating whether there have been consistent, long-term changes in chamois body mass across sexes and populations. Having identified the existence of temporal trends in body mass, we seek to tease apart the different drivers of these trends. For reasons discussed below (see Methods), we dismiss the role of artificial selection driven by hunting. Consequently, motivated by the findings reviewed above, we seek to evaluate the evidence for three plausible drivers of body size change:
Climate-mediated changes in vegetation productivity or phenology, altering resource quality and availability.
Climate-mediated changes in behaviour, altering access to resources.
Changes in population density, altering per-capita resource availability.
The study area is located in Trento Province in the Central-Eastern Italian Alps (46°02’N, 10°38’E), across three chamois hunting districts: Adamello (area = 373 km2), Presanella (146 km2) and Brenta (263 km2). The area is forested up to the tree-line at about 2,000 m, above which it consists of Alpine meadows, rocky outcrops, scree fields and open rock faces. The average altitude varies among the districts, though with considerable overlap (mean altitude ± SD: Adamello, 1,901 ± 616 m; Presanella, 2,098 ± 540 m; Brenta, 1,594 ± 603 m). Adamello and Presanella are characterised by nutrient-poor siliceous vegetation whilst Brenta is characterised by nutrient-rich calcareous vegetation . Typically, meadows in Adamello and Presanella are dominated by Festuca scabriculmis and Carex curvula, whilst those in Brenta are composed of Sesleria albicans and Carex firma. Throughout the study area, meadows are grazed by small herds of livestock (sheep, goats and cows) during summer, a practice that has been maintained at consistent levels throughout the study period. Several potential predators of chamois were present during the study, including a small, stable population of brown bear (Ursus arctos) in Brenta, a very small number of Eurasian lynx (Lynx lynx) and the golden eagle (Aquila chrysaetos). However, predation on chamois is very rare here (personal communication, Adamello Brenta Nature Park, Trento Province, Italy).
Body mass data
Chamois are hunted every year between mid-September and late-December. Data were collected on the eviscerated body mass and day of shooting of 10,455 yearling (≈1.5 year olds; hereafter juveniles) Alpine chamois (5,762 males and 4,693 females), hunted between 1979 and 2010 (see Additional file 1 for annual breakdowns of sample size). Hunting is heavily regulated and there is little potential for artificial selection by hunters, as chamois can easily detect hunters in the predominantly open habitat and will flee from hunters at particularly large distances ,. Moreover, there is no evidence of hunters preferentially harvesting larger bodied age-classes in these populations . Hunting pressure on yearlings varies among sites (mean proportion of yearlings hunted in census years: Adamello males, 0.40 ± 0.01; Adamello females, 0.32 ± 0.01; Presanella males, 0.32 ± 0.01; Presanella females, 0.24 ± 0.02; Brenta males, 0.37 ± 0.01; Brenta females, 0.31 ± 0.02). In order to account for intra-seasonal variation in body mass, which is not the focus of this study, a published model of seasonal body mass change , which considers inter-annual mass variation, was used to estimate juvenile mass standardised to a specific day of the year. Annual estimates (n = 32) of mean juvenile body mass were produced for each sex, within each site, standardised to day 300 of the year (27th October) (see Figure 1). Body mass was estimated after the vegetation growing season (hereafter `growing season') because body condition at that time will have been influenced by the spring and summer environment, which is thought to have a strong influence on ungulate body mass .
There were clear negative temporal body mass trends in all sexes and sites (Figure 1). In order to examine drivers of deviations from the long-term trends (i.e. years in which mean body mass was particularly high or low, even given the trend), the body mass time series were detrended by fitting linear models and calculating residuals. However, detrending can remove long-term fluctuations related to environmental trends ,, which are of primary interest to us. As a result, we modelled body mass data, to examine drivers of long-term trends, and also modelled body mass residuals, to examine drivers of deviations from the trends.
Environmental and demographic data
A range of climatic and non-climatic factors might be expected to influence chamois body mass. Negative effects of population density on mass are common in ungulates ,. In the absence of natural predation, these effects generally operate through increased intra-specific competition at higher population densities, resulting in lower per-capita food intake, particularly during periods when food is scarce -. To investigate density-dependence in these chamois populations, site-specific population density estimates were used from total population censuses performed in September every year between 1981 and 2009 (with the exception of 1990 and 1991; data from these years were excluded from the analysis). Each year, a set of simultaneous censuses was performed from vantage points across different blocks of each hunting district. It was assumed that density estimates from this time of year would reflect the population density over the previous growing season.
To investigate a possible direct thermoregulatory link between climate and body mass, we calculated yearly site-specific estimates of mean daily growing season temperature between 1982 and 2007 from high-altitude meteorological stations located in each of the three study sites (Data provided by The Forecasts and Organization Office, Civil Protection and Infrastructures Department, Trento Province, Italy). Differences in the elevation of weather stations among sites contributed to inter-site temperature differences (see Figure 2a). However, this did not affect our analysis since the drivers of body mass trends were examined separately in each site and, additionally, temperatures were standardised within sites, along with all other environmental predictors (see `Modelling variation in mass and mass residuals'). The bounds of the growing season were estimated using snow cover data, also from meteorological stations located within each site. The growing season was defined as the period between the snow-melt in spring, when snow cover was reduced to 0% (which generally occurs between late March and early May), and the first significant snowfall in winter that results in new snow settling on the ground (which generally occurs between early November and late December).
To investigate the effect of vegetation productivity and phenology on mass, NDVI (normalised difference vegetation index) data were used as a measure of vegetation productivity, processed by the Global Inventory Modelling and Mapping Studies group (GIMMS; ,). These data are global at a 0.07 degree resolution (approximately 8 km by 8 km) and are available at fortnightly intervals between 1982 and 2006 (thus slightly restricting the yearly data range for analyses). In order to focus on vegetation types utilised by chamois for foraging, such as alpine meadows and sparsely vegetated areas, only NDVI pixels dominated by such vegetation types were considered. To do this the Corine land-cover 2006 data-set at a 100 m resolution  was used to select only NDVI pixels within each site containing less than 25% coniferous woodland. As each of the three sites encompassed a number of these NDVI pixels, mean NDVI from these pixels was calculated for each fortnightly time period, within each site. Previous studies have implicated a number of metrics relating to annual NDVI variation as being important to ungulate body condition (e.g. ,). Here, we seek to derive NDVI metrics in a standardised fashion, despite inherent noise in NDVI estimates caused by factors such as cloud cover, water, snow or shadow . As in previous studies , we used a smoothed function to characterise variation in NDVI with time in a given year. The following function was used (see Additional file 2, for an illustration of the functional form):
Here, is predicted NDVI at time-period t in site s and year y, α s,y and β s,y are minimum and maximum NDVI respectively in site s and year y, σ s,y is a parameter related to the width of the function and z s,y is a parameter describing the shape of the function. Variation in NDVI data, p(s,y,t), about the predicted mean was beta distributed. Thus, the likelihood of the model parameters, θ s , y , given the data, parameterised by and the dispersion coefficient ϕs,y, is
The most parsimonious fitted relationships between mean NDVI and time in each year and site were calculated. Using these relationships, four NDVI metrics, described below, were calculated relating to vegetation productivity and phenology. All four of the metrics selected have been highlighted as important either to juvenile chamois specifically or to other ungulate species ,. Previously, climate-induced changes in spring growing conditions  have been linked to higher juvenile body mass in ungulates, including chamois , due to longer growing-seasons and higher vegetation quality . However, warmer springs have also been linked to negative impacts on body mass as higher temperatures lead to faster rates of vegetation `green-up' and, thus, a shorter period of access to nutrient-rich emergent vegetation associated with early spring ,. Here, the following four metrics were used: maximum rate of spring green-up, growing season duration, maximum NDVI and total growing season NDVI. Maximum rate of spring green-up (NDVIrate) was calculated as the maximum first derivative of (i.e. the maximum rate of NDVI increase). The duration of the growing season (NDVIdur) was calculated as the length of time between the maximum second derivative of (the start date of the growing season; when the rate of NDVI increase is increasing at its maximum rate) and the minimum second derivative of (the end date of the growing season; when NDVI is decreasing most rapidly). Maximum NDVI (NDVImax) was calculated as the maximum value of and total growing-season NDVI (INDVI) as the integral of within the bounds of the growing-season. An illustration of the calculation of these metrics can be seen in Additional file 2.
Modelling variation in mass and mass residuals
Environmental predictors were standardised by z-transformation within each site, as follows: , where z s,y is a z-transformed predictor in year y and site s, x s,y is the untransformed predictor, is the site-specific mean of that predictor and σ s the site-specific standard deviation. General linear models were fitted to examine variation in body mass and body mass residuals using R version 2.12.0 . Juvenile phenotypic quality is thought to be strongly influenced by lagged environmental effects . As such, lagged environmental factors ,, affecting the availability and accessibility of resources to mothers, as well as population density ,, affecting the intensity of competition for resources, can be strong predictors of juvenile body mass. To account for lagged environmental effects on juvenile mass, density, temperature and NDVI data were used not only from year y (the year a cohort was shot), but also from year y-1 (the year of a cohort’s birth). Models were fitted with all possible biologically meaningful combinations of population density, NDVIrate, NDVIdur, NDVImax, INDVI and temperature, from years y and y-1. Models were considered with either the same intercept or different intercepts for males and females. In each site, several of the predictors were highly correlated (Pearson correlation coefficients ≥ 0.6), for instance NDVIdur with INDVI, and densityy with densityy-1. To avoid problems of multicollinearity, highly correlated predictors were not permitted in the same model (but were permitted separately). Given the temporal range of the predictors, and to use the same temporal range of data in each model, data were used for the 19 years between 1983 and 2006, excluding 1990, 1991 and 1992, years for which estimates of densityy and/or densityy-1 are not available (population censuses were not performed in 1990 and 1991). In order to identify the most parsimonious models, we used the two-step model selection process suggested by Richards . First, all models having an AIC within six units of the smallest AIC calculated were selected (i.e. ΔAIC ≤ 6). Second, in order to remove overly complex models, we disregarded those that had a higher AIC than any simpler nested model. To investigate the potential for sex-specific environmental effects, we tested for sex interaction terms with all predictors within models in each site’s top model set. To assess the relative importance of different predictors, Akaike model averaged coefficients were calculated from all models in each site’s top model set . Since we wanted to compare how the relative importance of predictors varied among sites, model averaged coefficients were z-transformed within each site. This standardised coefficients from different sites to the same scale, allowing the most important effects to be readily compared among sites.
Juvenile chamois body masses decreased strongly between 1979 and 2010 in all three populations (Figure 1). The extent of this decrease varied considerably among sexes and sites but decreases in male mass have been more pronounced than decreases in female mass in all sites (slopes of temporal mass trends ± SE: Adamello males, -0.11 ± 0.01; Adamello females, -0.03 ± 0.02; Presanella males, -0.17 ± 0.02; Presanella females, -0.13 ± 0.01; Brenta males, -0.05 ± 0.01; Brenta females, -0.04 ± 0.01). Decreases have been less pronounced in Brenta than in the other two sites. There have been striking increases in growing season temperatures in all sites between 1982 and 2007 (Figure 2a). During the same period, all three populations increased in density substantially, peaking in the mid-1990s before declining slightly in recent years (Figure 2b). This growth coincides with the implementation of stricter controls on hunting in the area (including increases in the number of rangers and a more strictly enforced quota system). In contrast, there have been no pronounced long-term trends in the four NDVI metrics between 1982 and 2006, although growing seasons have tended to be longer (Figure 3b) and more productive (Figure 3d) between 2004 and 2006.
The most parsimonious body mass models fitted the observed data well (Table 1, Figure 4) (R2: Adamello, 0.78; Presanella, 0.83; Brenta, 0.69). In comparison, the fits of body mass residual models were poor (R2: Adamello, 0.44; Presanella, 0.13; Brenta, 0.36; see Additional file 3); thus, our inferences focus on models that describe longer-term changes in body mass, rather than those focused on explaining variation around the long-term trend. Temperature, density and NDVI all appear to play a role in describing long-term variation in juvenile body mass (Table 1, Figure 5). Within each site, the top set of body mass models contains a number of closely competing models (Table 1) but clear and consistent patterns across sites are illustrated by model averaging (Figure 5). Temperature and/or density have the strongest negative effects on mass in all sites (Figures 5 and 6). Furthermore, strong negative effects of temperature and density in the current year are present in all of the top models within Adamello and Presanella (Table 1), providing good evidence for these effects. A slightly weaker negative effect of lagged temperature appears in most of the top models for these sites. In Brenta, there is some evidence for a negative effect of temperature, which appears in the top three models in this site but there is no evidence of a density effect. Despite mass declines in males being consistently stronger than in females (Figure 1), we found no support for any interaction terms between sex and environmental predictors, suggesting that there are not strong differences in the magnitude of environmental effects on body mass between males and females. The lack of evidence for sex-specific effects could stem from the slightly restricted temporal range of data used in the analysis.
The modelled effects of NDVI are much weaker than the effects of density and temperature, and are generally positive (Table 1, Figure 5). Only in Brenta is there strong evidence for an influence of vegetation productivity and phenology; a strong, lagged, positive effect of INDVI is present in all top models and there is some evidence for a lagged, positive effect of maximum NDVI. In Adamello and Presanella, NDVI effects are generally weaker and inconsistent, although there is some support for a positive effect of long growing seasons in Presanella and, surprisingly, for a negative effect of INDVI in Adamello.
Pronounced long-term body mass declines were identified in all three populations. Models of juvenile body mass performed well in explaining these declines (Figure 4). However, the poor performances of residual models show that deviations from the general trends (i.e. years in which juveniles are particularly heavy or particularly light relative to the trend) cannot be explained as accurately. Mass declines appear to be linked most strongly to increases in growing season temperatures and population density; temperature and/or density have the strongest negative effects on body mass and deviations from yearly trend in all sites (Table 1, Figure 5). In comparison, there is no evidence of a strong negative effect of vegetation productivity or phenology on body mass, suggesting that body mass declines were not mediated by vegetation change as indexed by NDVI. Our results provide an interesting contrast to the majority of studies to date, which have implicated changes in food availability and quality as the most likely cause of body size declines .
Growing season temperature appears to be strongly linked to mass declines in all three populations (Table 1, Figures 5 and 6). However, despite significant temperature change (Figure 2a), there has been no consistent change in either vegetation productivity or phenology (Figure 3). Furthermore, there is no relationship between temperature and any NDVI metric. As such, we find no evidence for our first putative driver of mass change: that climate change has affected body mass via effects on resource productivity or phenology. Instead, our results provide support for our second putative driver: that climate change could be directly affecting chamois behaviour or physiology, limiting their ability to acquire resources. In another alpine ungulate, the ibex (Capra ibex), temperatures above 15-20°C result in heat discomfort in males, reducing the time they can spend foraging . In ungulates, behavioural changes, such as allocating less time to foraging, play an important role in thermoregulation because they can be more energetically efficient than autonomic thermoregulation ,. Higher daily temperatures during spring and summer may have led to juvenile chamois spending more time resting and less time foraging than in the past, reducing their ability to store energy reserves and invest in growth. Indeed, chamois, like many ungulates, reduce their feeding activity during the hottest period of the day . Furthermore, chamois spend less time foraging when it is hotter, independent of time of day . Our results suggest that temperature affects the ability of juveniles to acquire resources themselves rather than via lagged effects on their mothers (Figure 5). It is possible that juvenile ungulates are more susceptible to higher temperatures due to having higher relative metabolic demands than adults ,. Whilst we did not find any evidence for sex-specific environmental effects on body mass, that mass declines in juvenile males are consistently more pronounced than in females (Figure 1) suggests that, even in a relatively monomorphic ungulate species, a more energetically demanding growth phase could make males more susceptible to environmental stressors .
Population density also appears to be strongly linked to mass declines in juvenile chamois, providing support for our third putative driver: that changes in density affect per-capita resource availability (Table 1, Figures 5 and 6). The study populations have grown substantially following the implementation of stricter controls on hunting. As a result, increased intra-specific competition for resources may have led to reduced rates of per-capita food intake. Additionally, climate change may have affected chamois survival or fecundity, contributing to this population growth. Density dependence in body mass has been detected in a range of ungulate species (e.g. ,,) and even specifically in juvenile chamois . At high latitudes and in alpine regions, this effect is most likely to manifest during winter, when snow cover reduces forage availability and increases forage patchiness, resulting in higher levels of agonistic interactions among individuals ,. It is also possible that density dependence could be mediated by overgrazing; however, no changes in NDVI consistent with overgrazing, such as decreases in maximum NDVI or INDVI , were detected. Overgrazing can also encourage the colonisation of grazing-tolerant, less palatable species , which might not result in detectable changes in productivity. However, such changes would not alter our finding that increases in density appear to have influenced body mass via changes in per-capita resource availability. As with temperature, the effect of density in the year of shooting appears most important (Figure 5), suggesting that competition for forage following weaning, rather than a lagged influence of competition on mothers in the previous year, has a stronger influence on juvenile mass. Indeed, intra-specific competition for resources is most likely to result in the displacement of subordinate individuals, such as juveniles, from food patches -. Previously, the importance of lagged effects on juvenile condition has been stressed (e.g. ,). However, our findings are consistent with recent work on juvenile chamois showing that environmental conditions during the second year of life, which influence the ability of juveniles to acquire resources directly, have an important effect on investment in growth and energy storage ,. Only in Brenta is there no evidence for density dependence in body mass (Figures 5 and 6). Brenta is calcareous, harbouring a more nutrient-rich plant community , so forage availability could be less limiting here than in the other populations. A hunting regime maintaining this population below its carrying capacity, and limiting the effect of resource competition on mass, could also contribute to the less pronounced mass declines observed in Brenta (Figure 1).
There was no evidence for a strong role of vegetation productivity or phenology in the observed long-term mass declines (Table 1, Figure 5). Given that NDVI is only a proxy for vegetation productivity, it remains possible that other changes in vegetation quality and phenology have had more of an effect than suggested here. As with overgrazing, it is possible that warmer temperatures in spring and summer could lead to increases in the abundance of less palatable species , or reductions in the protein content of the species present , changes which NDVI might not detect. Although they do not appear to play a strong role in mass declines, there is some evidence that vegetation productivity and phenology influence variation in juvenile body mass. This is particularly the case in the calcareous area, Brenta, where there is strong evidence that productive growing seasons (in terms of INDVI and maximum NDVI) have positive effects on juvenile body mass (Table 1, Figure 5). This is consistent with findings from another calcareous Alpine area, where long growing seasons have been linked to higher juvenile chamois body mass . In our study, it is intriguing that the relative importance of environmental factors should vary among neighbouring areas that differ predominantly only in their geological substrate. These differences may arise because calcareous vegetation, whilst being more nutrient-rich, varies more in quality in response to environmental variation than siliceous vegetation ,. The variation in substrate in this study area has previously been shown to mediate differences in reproductive strategy, body mass  and horn length ; this study further highlights the importance of considering geological variation in studies of animal morphology, physiology and life history.
We have detected strong links between recent environmental change and negative temporal body mass trends in juvenile chamois. Increases in both temperature, due to climate change, and population density, due to stricter controls on hunting, appear to be driving the mass declines. Since we found no evidence for an effect of changing resource productivity or phenology on body mass, the observed patterns may not be mediated by changes in vegetation. We propose that heightened thermoregulatory demands and intra-specific competition, to both of which juveniles are particularly susceptible, could be responsible. These findings add an interesting contrast to the large number of studies implicating climate-mediated changes in resource productivity or phenology as drivers of declines in animal body size. Interestingly, our findings highlight that changes in management could ameliorate negative impacts of climate change. For example, future increases in hunting offtake could reduce the intensity of resource competition, counteracting negative impacts of temperature warming on body mass. This study further highlights the importance of considering management when examining the influence of environmental change on species, which is an increasingly important focus of ecological research (e.g. ,,). Future reductions in body size are likely to have far-reaching effects on the functioning and productivity of ecosystems ,; decreases in juvenile body condition could result in reduced survival of juveniles and depressed population growth rates. An increasing focus on the links between climate, management, body condition and demography will shed light on the influence of anthropogenic and climate change on populations.
TM, MA, SW and PS conceived the study. TM analysed the data and, with SW and PS, wrote the manuscript. MA and RC provided the data. All authors contributed to revisions and approved the final manuscript.
Parmesan C, Yohe G: A globally coherent fingerprint of climate change impacts across natural systems. Nature. 2003, 421: 37-42. 10.1038/nature01286.
Menzel A, Sparks TH, Estrella N, Koch E, Aasa A, Ahas R, Alm-Kubler K, Bissolli P, Braslavska O, Briede A, Chmielewski FM, Crepinsek Z, Curnel Y, Dahl A, Defila C, Donnelly A, Filella Y, Jatcza K, Mage F, Mestre A, Nordli O, Penuelas J, Pirinen P, Remisova V, Scheifinger H, Striz M, Susnik A, Van Vliet AJH, Wielgolaski FE, Zach S: European phenological response to climate change matches the warming pattern. Glob Chang Biol. 2006, 12: 1969-1976. 10.1111/j.1365-2486.2006.01193.x.
van de Pol M, Vindenes Y, Saether B-E, Engen S, Ens BJ, Oosterbeek K, Tinbergen JM: Effects of climate change and variability on population dynamics in a long-lived shorebird. Ecology. 2010, 91: 1192-1204. 10.1890/09-0410.1.
Musolin DL: Insects in a warmer world: ecological, physiological and life-history responses of true bugs (Heteroptera) to climate change. Glob Chang Biol. 2007, 13: 1565-1585. 10.1111/j.1365-2486.2007.01395.x.
Ozgul A, Childs DZ, Oli MK, Armitage KB, Blumstein DT, Olson LE, Tuljapurkar S, Coulson T: Coupled dynamics of body mass and population growth in response to environmental change. Nature. 2010, 466: 482-U485. 10.1038/nature09210.
Ozgul A, Tuljapurkar S, Benton TG, Pemberton JM, Clutton-Brock TH, Coulson T: The Dynamics of Phenotypic Change and the Shrinking Sheep of St. Kilda. Science. 2009, 325: 464-467. 10.1126/science.1173668.
Gardner JL, Peters A, Kearney MR, Joseph L, Heinsohn R: Declining body size: a third universal response to warming?. Trends Ecol Evol. 2011, 26: 285-291. 10.1016/j.tree.2011.03.005.
Sheridan JA, Bickford D: Shrinking body size as an ecological response to climate change. Nat Clim Change. 2011, 1: 401-406. 10.1038/nclimate1259.
Guthrie RD: Rapid body size decline in Alaskan Pleistocene horses before extinction. Nature. 2003, 426: 169-171. 10.1038/nature02098.
Bickford D, Howard SD, Ng DJJ, Sheridan JA: Impacts of climate change on the amphibians and reptiles of Southeast Asia. Biodivers Conserv. 2010, 19: 1043-1062. 10.1007/s10531-010-9782-4.
Bergmann C: Ueber die Verhältnisse der Wärmeökonomie der Thiere zu ihrer Grösse (Concerning the relationship of heat conservation of animals to their size). Gottinger Studien. 1847, 3: 595-708.
Mayr E: Geographical character gradients and climate adaptation. Evolution. 1956, 10: 105-108. 10.2307/2406103.
Ashton KG, Tracy MC, de Queiroz A: Is Bergmann’s rule valid for mammals?. Am Nat. 2000, 156: 390-415. 10.1086/303400.
Meiri S, Dayan T: On the validity of Bergmann's rule. J Biogeogr. 2003, 30: 331-351. 10.1046/j.1365-2699.2003.00837.x.
Freckleton RP, Harvey PH, Pagel M: Bergmann's rule and body size in mammals. Am Nat. 2003, 161: 821-825. 10.1086/374346.
Teplitsky C, Mills JA, Alho JS, Yarrall JW, Merila J: Bergmann’s rule and climate change revisited: Disentangling environmental and genetic responses in a wild bird population. Proc Natl Acad Sci U S A. 2008, 105: 13492-13496. 10.1073/pnas.0800999105.
Husby A, Hille SM, Visser ME: Testing Mechanisms of Bergmann’s Rule: Phenotypic Decline but No Genetic Change in Body Size in Three Passerine Bird Populations. Am Nat. 2011, 178: 202-213. 10.1086/660834.
Coltman DW, O'Donoghue P, Jorgenson JT, Hogg JT, Strobeck C, Festa-Bianchet M: Undesirable evolutionary consequences of trophy hunting. Nature. 2003, 426: 655-658. 10.1038/nature02177.
Allendorf FW, England PR, Luikart G, Ritchie PA, Ryman N: Genetic effects of harvest on wild animal populations. Trends Ecol Evol. 2008, 23: 327-337. 10.1016/j.tree.2008.02.008.
Clutton-Brock TH, Guinness FE, Albon SD: Red deer. Behavior and ecology of two sexes. 1982, University of Chicago & Edinburgh University Press, USA & Edinburgh
Gaillard JM, Festa-Bianchet M, Delorme D, Jorgenson J: Body mass and individual fitness in female ungulates: bigger is not always better. Proc R Soc London B Biol Sci. 2000, 267: 471-477. 10.1098/rspb.2000.1024.
Stewart KM, Bowyer RT, Dick BL, Johnson BK, Kie JG: Density-dependent effects on physical condition and reproduction in North American elk: an experimental test. Oecologia. 2005, 143: 85-93. 10.1007/s00442-004-1785-y.
Loison A, Langvatn R, Solberg EJ: Body mass and winter mortality in red deer calves: disentangling sex and climate effects. Ecography. 1999, 22: 20-30. 10.1111/j.1600-0587.1999.tb00451.x.
Pettorelli N, Weladji RB, Holand O, Mysterud A, Breie H, Stenseth NC: The relative role of winter and spring conditions: linking climate and landscape-scale plant phenology to alpine reindeer body mass. Biol Lett. 2005, 1: 24-26. 10.1098/rsbl.2004.0262.
Pettorelli N, Pelletier F, von Hardenberg A, Festa-Bianchet M, Cote SD: Early onset of vegetation growth vs. rapid green-up: Impacts on juvenile mountain ungulates. Ecology. 2007, 88: 381-390. 10.1890/06-0875.
Toïgo C, Gaillard JM, Van Laere G, Hewison M, Morellet N: How does environmental variation influence body mass, body size, and body condition? Roe deer as a case study. Ecography. 2006, 29: 301-308. 10.1111/j.2006.0906-7590.04394.x.
Mysterud A, Yoccoz NG, Stenseth NC, Langvatn R: Effects of age, sex and density on body weight of Norwegian red deer: evidence of density-dependent senescence. Proc R Soc London B Biol Sci. 2001, 268: 911-919. 10.1098/rspb.2001.1585.
Yom-Tov Y, Yom-Tov S, MacDonald D, Yom-Tov E: Population cycles and changes in body size of the lynx in Alaska. Oecologia. 2007, 152: 239-244. 10.1007/s00442-006-0653-3.
Laternser M, Schneebeli M: Long-term snow climate trends of the Swiss Alps (1931-99). Int J Climatol. 2003, 23: 733-750. 10.1002/joc.912.
Parker KL, Robbins CT, Hanley TA: Energy expenditures for locomotion by mule deer and elk. J Wildl Manage. 1984, 48: 474-488. 10.2307/3801180.
Aublet JF, Festa-Bianchet M, Bergero D, Bassano B: Temperature constraints on foraging behaviour of male Alpine ibex (Capra ibex) in summer. Oecologia. 2009, 159: 237-247. 10.1007/s00442-008-1198-4.
Luterbacher J, Dietrich D, Xoplaki E, Grosjean M, Wanner H: European seasonal and annual temperature variability, trends, and extremes since 1500. Science. 2004, 303: 1499-1503. 10.1126/science.1093877.
Sinervo B, Mendez-de-la-Cruz F, Miles DB, Heulin B, Bastiaans E, Cruz MVS, Lara-Resendiz R, Martinez-Mendez N, Calderon-Espinosa ML, Meza-Lazaro RN, Gadsden H, Avila LJ, Morando M, De la Riva IJ, Sepulveda PV, Rocha CFD, Ibarguengoytia N, Puntriano CA, Massot M, Lepetz V, Oksanen TA, Chapple DG, Bauer AM, Branch WR, Clobert J, Sites JW: Erosion of Lizard Diversity by Climate Change and Altered Thermal Niches. Science. 2010, 328: 894-899. 10.1126/science.1184695.
Kearney M, Shine R, Porter WP: The potential for behavioral thermoregulation to buffer “cold-blooded” animals against climate warming. Proc Natl Acad Sci U S A. 2009, 106: 3835-3840. 10.1073/pnas.0808913106.
Stuble KL, Pelini SL, Diamond SE, Fowler DA, Dunn RR, Sanders NJ: Foraging by forest ants under experimental climatic warming: a test at two sites. Ecol Evol. 2013, 3: 482-491. 10.1002/ece3.473.
Kjellander P, Gaillard JM, Hewison AJM: Density-dependent responses of fawn cohort body mass in two contrasting roe deer populations. Oecologia. 2006, 146: 521-530. 10.1007/s00442-005-0188-z.
Chirichella R, Ciuti S, Grignolio S, Rocca M, Apollonio M: The role of geological substrate for horn growth in ungulates: a case study on Alpine chamois. Evol Ecol. 2013, 27: 145-163. 10.1007/s10682-012-9583-1.
Mason THE, Chirichella R, Richards SA, Stephens PA, Willis SG, Apollonio M: Contrasting life histories in neighbouring populations of a large mammal. PLoS One. 2011, 6: e28002-10.1371/journal.pone.0028002.
Mason THE, Stephens PA, Willis SG, Chirichella R, Apollonio M, Richards SA: Intraseasonal variation in reproductive effort: young males finish last. Am Nat. 2012, 180: 823-830. 10.1086/668082.
Oftedal OT: Body size and reproductive strategy as correlates of milk energy output in lactating mammals. Acta Zool Fenn. 1984, 171: 183-186.
Wu ZH, Huang NE, Long SR, Peng CK: On the trend, detrending, and variability of nonlinear and nonstationary time series. Proc Natl Acad Sci U S A. 2007, 104: 14889-14894. 10.1073/pnas.0701020104.
Jonzen N, Hedenstrom A, Hjort C, Lindstrom A, Lundberg P, Andersson A: Climate patterns and the stochastic dynamics of migratory birds. Oikos. 2002, 97: 329-336. 10.1034/j.1600-0706.2002.970303.x.
Skogland T: The effects of density-dependent resource limitations on the demography of wild reindeer. J Anim Ecol. 1985, 54: 359-374. 10.2307/4484.
Fryxell J: Food limitation and demography of a migratory antelope, the white-eared kob. Oecologia. 1987, 72: 83-91. 10.1007/BF00385049.
Sinclair AR, Dublin H, Borner M: Population regulation of Serengeti Wildebeeest: a test of the food hypothesis. Oecologia. 1985, 65: 266-268. 10.1007/BF00379227.
Tucker CJ, Pinzon JE, Brown ME, Slayback DA, Pak EW, Mahoney R, Vermote EF, El Saleous N: An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data. Int J Remote Sens. 2005, 26: 4485-4498. 10.1080/01431160500168686.
Pettorelli N, Vik JO, Mysterud A, Gaillard JM, Tucker CJ, Stenseth NC: Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends Ecol Evol. 2005, 20: 503-510. 10.1016/j.tree.2005.05.011.
Corine land cover 2006 (CLC2006) 100m - version 2/2010. 2010
Beck PS, Atzberger C, Høgda KA, Johansen B, Skidmore AK: Improved monitoring of vegetation dynamics at very high latitudes: A new method using MODIS NDVI. Remote Sens Environ. 2006, 100: 321-334. 10.1016/j.rse.2005.10.021.
Burnham KP, Anderson DR: Model selection and multimodel inference: a practical information-theoretic approach. 2002, Springer, New York
Richards SA: Dealing with overdispersed count data in applied ecology. J Appl Ecol. 2008, 45: 218-227. 10.1111/j.1365-2664.2007.01377.x.
Garel M, Gaillard JM, Jullien JM, Dubray D, Maillard D, Loison A: Population abundance and early spring conditions determine variation in body mass of juvenile chamois. J Mammal. 2011, 92: 1112-1117. 10.1644/10-MAMM-A-056.1.
Linderholm HW: Growing season changes in the last century. Agr Forest Meteorol. 2006, 137: 1-14. 10.1016/j.agrformet.2006.03.006.
Crawley MJ: Herbivory: the dynamics of animal-plant interactions. 1983, Blackwell, Oxford
R: A language and environment for statistical computing. 2011, R Foundation for Statistical Computing, Vienna, Austria
Côté SD, Festa-Bianchet M: Birthdate, mass and survival in mountain goat kids: effects of maternal characteristics and forage quality. Oecologia. 2001, 127: 230-238. 10.1007/s004420000584.
Weladji RB, Holand O: Global climate change and reindeer: effects of winter weather on the autumn weight and growth of calves. Oecologia. 2003, 136: 317-323. 10.1007/s00442-003-1257-9.
Bonenfant C, Gaillard JM, Klein F, Loison A: Sex- and age-dependent effects of population density on life history traits of red deer Cervus elaphus in a temperate forest. Ecography. 2002, 25: 446-458. 10.1034/j.1600-0587.2002.250407.x.
Efron B, Tibshirani R: Statistical Data Analysis in the Computer Age. Science. 1991, 253: 390-395. 10.1126/science.253.5018.390.
Dussault C, Ouellet JP, Courtois R, Huot J, Breton L, Larochelle J: Behavioural responses of moose to thermal conditions in the boreal forest. Ecoscience. 2004, 11: 321-328.
Maloney SK, Moss G, Cartmell T, Mitchell D: Alteration in diel activity patterns as a thermoregulatory strategy in black wildebeest (Connochaetes gnou). J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 2005, 191: 1055-1064. 10.1007/s00359-005-0030-4.
Ruttimann S, Giacometti M, McElligott AG: Effect of domestic sheep on chamois activity, distribution and abundance on sub-alpine pastures. Eur J Wildl Res. 2008, 54: 110-116. 10.1007/s10344-007-0118-y.
Mason THE, Stephens PA, Apollonio M, Willis SG: Predicting potential responses to future climate in an alpine ungulate: inter-specific interactions exceed climate effects.Glob Chang Biol 2014, In press.,
Hudson RJ: Body size, energetics, and adaptive radiation. 1985, CRC Press, Inc., Boca Raton, Florida
Solberg EJ, Jordhoy P, Strand O, Aanes R, Loison A, Saether BE, Linnell JDC: Effects of density-dependence and climate on the dynamics of a Svalbard reindeer population. Ecography. 2001, 24: 441-451. 10.1034/j.1600-0587.2001.d01-200.x.
Forchhammer MC, Clutton-Brock TH, Lindstrom J, Albon SD: Climate and population density induce long-term cohort variation in a northern ungulate. J Anim Ecol. 2001, 70: 721-729. 10.1046/j.0021-8790.2001.00532.x.
Weckerly FW: Social bonding and aggression in female Roosevelt elk. Can J Zool. 1999, 77: 1379-1384. 10.1139/z99-101.
Espmark Y: Social behaviour of roe deer at winter feeding stations. Appl Anim Ethol. 1974, 1: 35-47. 10.1016/0304-3762(74)90006-6.
Hilker T, Natsagdorj E, Waring RH, Lyapustin A, Wang Y: Satellite observed widespread decline in Mongolian grasslands largely due to overgrazing. Glob Chang Biol. 2014, 20: 418-428. 10.1111/gcb.12365.
Fernandez-Gimenez ME, Allen-Diaz B: Testing a non-equilibrium model of rangeland vegetation dynamics in Mongolia. J Appl Ecol. 1999, 36: 871-885. 10.1046/j.1365-2664.1999.00447.x.
Veiberg V, Loe LE, Mysterud A, Langvatn R, Stenseth NC: Social rank, feeding and winter weight loss in red deer: any evidence of interference competition?. Oecologia. 2004, 138: 135-142. 10.1007/s00442-003-1399-9.
Shi J, Dunbar R: Feeding competition within a feral goat population on the Isle of Rum, NW Scotland. J Ethol. 2006, 24: 117-124. 10.1007/s10164-005-0170-6.
Thouless C: Feeding competition between grazing red deer hinds. Anim Behav. 1990, 40: 105-111. 10.1016/S0003-3472(05)80669-4.
Gensac P: Plant and soil groups in the alpine grasslands of the Vanoise Massif, French Alps. Arctic Alpine Res. 1990, 22: 195-201. 10.2307/1551304.
Sage RF, Kubien DS: Quo vadis C4? An ecophysiological perspective on global change and the future of C4 plants. Photosynth Res. 2003, 77: 209-225. 10.1023/A:1025882003661.
Barbehenn RV, Chen Z, Karowe DN, Spickard A: C3 grasses have higher nutritional quality than C4 grasses under ambient and elevated atmospheric CO2. Glob Chang Biol. 2004, 10: 1565-1575. 10.1111/j.1365-2486.2004.00833.x.
Michalet R, Gandoy C, Joud D, Pages JP, Choler P: Plant community composition and biomass on calcareous and siliceous substrates in the northern French Alps: Comparative effects of soil chemistry and water status. Arctic Antarctic Alpine Res. 2002, 34: 102-113. 10.2307/1552514.
Brodie J, Johnson H, Mitchell M, Zager P, Proffitt K, Hebblewhite M, Kauffman M, Johnson B, Bissonette J, Bishop C: Relative influence of human harvest, carnivores, and weather on adult female elk survival across western North America. J Appl Ecol. 2013, 50: 295-305. 10.1111/1365-2664.12044.
Juillet C, Choquet R, Gauthier G, Lefebvre J, Pradel R: Carry-over effects of spring hunt and climate on recruitment to the natal colony in a migratory species. J Appl Ecol. 2012, 49: 1237-1246. 10.1111/j.1365-2664.2012.02199.x.
We are grateful to Trento Province, the Adamello Brenta Nature Park, and the “Associazione Cacciatori Trentini” (The regional hunting association), as well as to the hunting reserve of “Regole Spinale e Manez” for supplying data. We are indebted to M. Rocca, A. Brugnoli and all the game wardens who collected Alpine chamois measurements. TM was funded by a Natural Environment Research Council Doctoral Training Grant (NE/H524573/1) and supervised by SGW, PAS and MA. MA and RC were supported by the Italian Ministry of Education, University and Research (PRIN 2010-2011, 20108 TZKHC). This study complies with all relevant national, regional and provincial Italian laws. The data used in this study are available on request from the corresponding author.
The authors declare that they have no competing interests.
An erratum to this article is available at http://dx.doi.org/10.1186/s12983-014-0088-3.