Email updates

Keep up to date with the latest news and content from BMC Ecology and BioMed Central.

Open Access Research article

Environmental, biological and anthropogenic effects on grizzly bear body size: temporal and spatial considerations

Scott E Nielsen1*, Marc RL Cattet2, John Boulanger3, Jerome Cranston4, Greg J McDermid5, Aaron BA Shafer6 and Gordon B Stenhouse7

Author Affiliations

1 Department of Renewable Resources, University of Alberta, Edmonton, Alberta T6G 2H1, Canada

2 Canadian Cooperative Wildlife Health Centre, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5B4, Canada

3 Integrated Ecological Research, Nelson, BC V1L 5T2, Canada

4 Arctos Ecological Services, Hinton, Alberta T7V 1H9, Canada

5 Department of Geography, University of Calgary, Calgary, Alberta T2N 1N4, Canada

6 Department of Ecology and Genetics, Uppsala University, Uppsala, SE 75240, Sweden

7 Foothills Research Institute, Hinton, Alberta T7V 1X6, Canada

For all author emails, please log on.

BMC Ecology 2013, 13:31  doi:10.1186/1472-6785-13-31

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1472-6785/13/31


Received:3 April 2013
Accepted:6 September 2013
Published:8 September 2013

© 2013 Nielsen et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Individual body growth is controlled in large part by the spatial and temporal heterogeneity of, and competition for, resources. Grizzly bears (Ursus arctos L.) are an excellent species for studying the effects of resource heterogeneity and maternal effects (i.e. silver spoon) on life history traits such as body size because their habitats are highly variable in space and time. Here, we evaluated influences on body size of grizzly bears in Alberta, Canada by testing six factors that accounted for spatial and temporal heterogeneity in environments during maternal, natal and ‘capture’ (recent) environments. After accounting for intrinsic biological factors (age, sex), we examined how body size, measured in mass, length and body condition, was influenced by: (a) population density; (b) regional habitat productivity; (c) inter-annual variability in productivity (including silver spoon effects); (d) local habitat quality; (e) human footprint (disturbances); and (f) landscape change.

Results

We found sex and age explained the most variance in body mass, condition and length (R2 from 0.48–0.64). Inter-annual variability in climate the year before and of birth (silver spoon effects) had detectable effects on the three-body size metrics (R2 from 0.04–0.07); both maternal (year before birth) and natal (year of birth) effects of precipitation and temperature were related with body size. Local heterogeneity in habitat quality also explained variance in body mass and condition (R2 from 0.01–0.08), while annual rate of landscape change explained additional variance in body length (R2 of 0.03). Human footprint and population density had no observed effect on body size.

Conclusions

These results illustrated that body size patterns of grizzly bears, while largely affected by basic biological characteristics (age and sex), were also influenced by regional environmental gradients the year before, and of, the individual’s birth thus illustrating silver spoon effects. The magnitude of the silver spoon effects was on par with the influence of contemporary regional habitat productivity, which showed that both temporal and spatial influences explain in part body size patterns in grizzly bears. Because smaller bears were found in colder and less-productive environments, we hypothesize that warming global temperatures may positively affect body mass of interior bears.

Keywords:
Bear; Silver spoon; Environmental effects; GPS radiocollar; Temporal and spatial heterogeneity

Background

Understanding how spatial and temporal heterogeneity of environments affect life-history traits and the growth of individuals has been a central theme in ecology and population biology [1-3]. Among other measures of phenotype, body size for many species is highly variable across different spatial and temporal scales, which illustrates the importance of environmental heterogeneity on the growth of individuals and populations. Understanding how these spatial and temporal dynamics affect phenotypes is critical to helping identify and prioritize management actions for many species of special concern, especially in today’s rapidly changing world.

There is little argument that spatial heterogeneity of environments shape populations by affecting population density, fitness, dispersal and behaviour [4-6]. Indeed, such relationships are a cornerstone of landscape ecology [7,8] and habitat selection theory [9,10], and form the basis for natural-resource-management. Inter-annual variability in environments creates pulsed-resource dynamics that affect many animal populations [11-13] by affecting primary productivity [14-16] and the frequency and intensity of landscape disturbances [17,18]. For example, climatic oscillations that impact plant productivity will in turn affect primary consumer populations [1,19,20] and thus other trophic levels dependent on primary consumers [21,22]. For consumers that are specialized on fruit (frugivores), which often exhibit supra-annual variation in productivity [23,24], climate conditions can have an important effect on population dynamics and the health of animals. For example, masting events or mast failures are often signalled by climatic conditions [25-28]. On Barro Colorado Island in Panama, warm ENSO events stimulate fruit masting in tropical trees resulting in population increases of frugivore species [14,29]. Likewise, acorn production for many species of oaks in the USA and cones for spruce in Canada are known to mast synchronously across broad spatial scales [30-32] having profound effects on consumer populations [21,33,34].

Increasingly, it appears that such inter-annual variations have long-term effects on individuals, particularly for those experiencing boom or bust conditions during early life. In fact, conditions during in utero or natal periods can be as, or more, important than recent conditions on animal health and fitness [35-37]. This phenomenon is referred to as the “silver-spoon” effect as it emphasizes the importance of being born into “rich” environments [38]. Since resource conditions vary among years for nearly all ecosystems, populations often exhibit cohort effects that structure population dynamics [1,39]. For instance, cone production in white spruce during natal periods and temperature during in utero conditions had long-lasting effects on red squirrel reproductive success in the Yukon of Canada [37]. Likewise, population growth of stoats in New Zealand beech forests is dependent on masting [39].

One species that inhabits highly variable environments with limited resources relative to their dietary needs and large body size are grizzly (brown) bears (Ursus arctos L.) [40]. All the calories necessary to survive and reproduce are acquired in the approximately seven months that they are active prior to about five months of fasting in a den. The importance of limiting resources and phenotypic plasticity is further emphasized by nearly a 10-fold difference in adult body mass across the species’ range [41]. Most often, grizzly bears rely on the seasonal or inter-annual pulsing of high-calorie resources, such as salmon in coastal ecosystems [42-44] or hard and soft mast in interior populations [45-47]. Not surprisingly, body size in bears varies accordingly [48,49], having ramifications to both survival [43,50,51] and reproduction [48,52,53]. Given these resource demands and the existence of environmental uncertainty, grizzly bears have evolved a reproductive mechanism to compensate for these factors – the delayed facultative implantation of the fertilized egg dependent on autumn body condition [54-56]. Understanding body size-environment relations is therefore critical to understanding population processes in grizzly bears, particularly reproductive success and population growth.

Here, we evaluated the importance of six different factors on springtime body size patterns in grizzly bears of Alberta, Canada (see Table 1). The six factors we examined were: (1) regional habitat productivity; (2) inter-annual variability in productivity (e.g. silver-spoon effects); (3) habitat quality; (4) human footprint and activity; (5) rate of landscape change; and (6) density dependence. Our objective was to examine how each of these factors affected body mass, length and condition after accounting for age, sex, offspring dependence and capture effects.

Table 1. Environmental variables used to measure hypothesized environmental drivers of body size patterns in grizzly bears within Alberta, Canada

Methods

Study area

Our study area consisted of a span of 750 km along the eastern slopes of the Canadian Rocky Mountains in Alberta, Canada (Figure 1). Grizzly bears in Alberta are considered ‘interior’ since they lack marine subsidized salmon resources. The area is characterized by cold continental climates without a dry season. Protected areas dominate the mountains, where as the foothills consist largely of multiple resource-use activities of forestry and energy extraction resulting in higher levels of forest fragmentation and human activity [57-59].

thumbnailFigure 1. Grizzly bear capture locations in Alberta, Canada for 112 unique animals across a 750 km distance. Years of capture by population unit indicated along the side of each population unit. Inset map illustrates location within the current range of the species in North America.

Grizzly bear observations

Grizzly bears were captured during the springtime (April to June) from 1999 to 2008 using remote drug delivery (Pneu-Dart Inc. or Dan-Inject) either by helicopter or following restraint by foot snares [60], and since 2004 with culvert traps [61,62]. Capture data, including a breakdown by sex, age and number of individual captures is shown in Figure 2. All bears were anesthetized using a combination of xylazine and zolazepam–tiletamine administered intramuscularly as xylazine at 2 mg/kg and Telazol at 3 mg/kg estimated body weight [63]. We administered atipamezole at 0.15–0.20 mg/kg, half-volume intramuscularly and half-volume intravenously, to reverse the effects of xylazine. Grizzly bears were weighed using a load scale (MSI-7200 Dynalink) and measured for length using a standard tape stretched over the top of the bear from the tip of the nose to the last tail vertebrae. A premolar was collected for aging bears using the number of cementum annuli [64], with adult status considered to be five years of age. For each bear a VHF ear-tag transmitter (Advanced Telemetry Systems, Isanti, MN) was attached and a Televilt Simplex, Tellus or Advanced Telemetry Systems GPS radiocollar fitted. Animal locations were transmitted every 4 hours prior to 2004 and at 1–2 hour intervals since 2004. Here we use data for 112 unique bears (57 female, 55 male) having an average age of 8.0 years (SD = 5.1) and ranging from subadult (2 years of age) to 22 years old.

thumbnailFigure 2. Grizzly bear capture data for 112 animals. a) Percent of animals, by sex, captured at each age class; b) Breakdown of the number of times an individual was captured (by overall percentage).

We used three measures of body size to represent short- to long-term measures of growth: mass; length; and body condition. Body condition was estimated using a body condition index where mass is measured relative to length [65]. Although we had multiple capture events for some animals, we only used the most recent capture because it maximized the range of ages considered. All captures and handling were done on public lands with permits and the capture and handling procedures approved by the University of Saskatchewan’s Committee on Animal Care and Supply (Permit Number: 20010016) following guidelines provided by the American Society of Mammalogists’ Animal Care and Use Committee [66] and the Canadian Council on Animal Care [67] for the safe handling of wildlife.

Statistical analysis

Age, sex and reproductive status (with or without offspring) of each animal was recorded. Number of times captured and density were also considered as response variables for body size measures. The local-population density was indexed as the number of genetically identified individuals surrounding a radiocollared bear [48,53]. Each bear was assigned a single geographic centroid based on their GPS telemetry locations and a buffer around this centroid based on the radius distance of the average daily movement rate of that animal’s sex-age class (4340 m to 10380 m radius). The number of detections of unique bears within each circular buffer was then estimated from DNA hair-snag information collected within 7x7km grids in 2004 to 2008. These counts were divided by the proportion of the buffer overlapping the DNA survey grid, and by the probability of capture (derived from data of the closest observed distance of GPS collared bears to known bait sites – see 67), which varied by the age, sex and reproductive status of the individual being detected, and the DNA survey stratum [68].

Regional environmental productivity was estimated for each bear at their home range centroid location based on monthly temperature and precipitation normal (or average) from 1971 – 2000 estimated with the software ClimateAB [69]. ClimateAB measures of climate normals are downscaled ANUSPLIN-interpolated monthly normal data (2.5 x 2.5 arcmin) using local weather-station data and an elevation lapse-rate adjustment [70]. Monthly climate normals for precipitation and temperature were considered for four seasonal periods (winter, spring, summer and growing season) and for the two individual months of March and July that represented late winter conditions affecting snowpack at high altitudes and peak primary productivity respectively (Table 1). We also considered ecosystem type (i.e., alpine, subalpine and foothills) as a surrogate of regional productivity based on habitat use (exposure at three possible zones of influence) measured from GPS radio-telemetry information. Zones of influence considered around each telemetry location included the local habitat-patch (HP) scale at the 30 m raster resolution, a flight-response (FR) scale of a 300 m radius representing exposure to direct human activity [71], and a landscape-encounter (LE) scale representing the average daily movement rate by sex group (scale or radius buffer).

We measured inter-annual variations in environments using ClimateAB [69] by estimating temperature and precipitation by month at each animal’s home range centroid from the time (year) prior to birth (Bt-1) to the year of capture (Ct0); due to missing data (i.e. locations prior to GPS collaring) and computational considerations, we are making the assumption that home range centroids have not changed over time or if changed that local variation in climates are small (see Discussion). The inter-annual variation (anomalies) was estimated as the absolute deviation in temperature and precipitation from 30 year (1971–2000) climate normals over the range of birth years observed in sampled bears for the same home range centroids again using ClimateAB [69]. By using anomalies rather than actual climate observations, we separated effects associated with regional productivity (climate normal) from inter-annual fluctuations (anomalies). Inter-annual variability was measured for: (1) maternal conditions (one year prior to birth; Bt-1); (2) in-utero and natal conditions (birth year and yearling; B0 and Bt+1); and (3) conditions during or prior to capture (Ct-1 and C0) (Table 1).

Local habitat quality was measured as habitat use (GPS telemetry) at the three scales of exposure (HP, FR and LE) for nine different measures of habitat quality reflecting the association of grizzly bears with disturbed and productive environments [72-74]: canopy cover, variation in canopy cover, deciduous canopy cover, amount of shrub habitat, forest age, forest age variation, amount of regenerating forest, variation in regenerating forest age and terrain soil wetness (Table 1). Non-linear effects were considered for canopy cover, deciduous canopy cover, forest age, amount of regenerating forest used and terrain soil wetness since intermediate amounts of these habitat conditions are normally preferred [72,74,75].

We used regional measures of human footprint and activity including the amount of habitat use associated with private lands (i.e., Alberta’s whitezone; see [75]), protected areas and high- or low-risk habitats based on a mortality risk and safe harbour habitat models [58], density of linear-access features, and distance to nearest human feature or recent energy wells (Table 1). Since we did not expect body size to be affected by human features and recent energy wells beyond local effects (distances), we developed exponential decay functions for each distance variable [75] using parameters of 300 m, 1 km and 3 km. A cost-weighted distance to roads was also considered where cost was defined by terrain ruggedness (a continuous variables accounting for change in elevation) under the assumption that more rugged areas near roads would be less penetrable to humans and thus experience lower human activity. Annual rate of landscape change was measured as the annual change (%) in habitat composition using annual remote sensing of major habitat types and anthropogenic features including roads, clear-cuts and energy well-pads [76].

We used the HIREG module [77] for the software STATA 11 to estimate hierarchical regressions [78] of body size based on the six main hypothesized drivers of growth. This approach was taken in order to partition variances and test for differences among the main hypothesized factors, and account for multiple measurement variables within each hypothesized factor (block) using variable ‘blocking’ approaches. The order of hierarchical regression model considered was: (1) biology effects including density-dependence; (2) regional habitat productivity; (3) inter-annual variation in environments in the form of maternal [(year before birth), in utero (year of birth) and natal (year after birth)] and capture effects (year of and before capture); (4) local habitat quality; (5) human footprint; and (6) landscape change. This order reflects the need to first control for biology before examining residual variance due to environment. We chose more regional measures of environment before inclusion of local measures of environment in the hierarchical order of blocks. No interactions among blocks were considered. For each hierarchical category, we selected predictors (i.e. block of variables) based on a forward step-wise regression procedure of variable blocks using a p < 0.1 significance level [79]. An F-test was used to determine whether changes to the coefficient of variation (R2) among the main hypothesized factors for each block were significant.

Results

Body mass

Biological and environmental factors explained 83.5% of the variation (R2, model F = 50.0, df = 10, 68, p < 0.001) in body mass (Table 2, Figure 3). Age of bears was non-linearly related to mass, and the additive effect of age of male animals explained 63.5% of the total variance in springtime body mass. Regional habitat productivity explained an additional 12.4% of variance (F = 18.8, df = 2, 73, p <0.001), as represented by two regional measures: early spring (March) precipitation and the use (habitat patch scale) of alpine habitats (Table 2). In both cases, body mass of bears was inversely related to regional habitat-productivity measures. In addition to regional productivity, inter-annual climate variability explained an additional 6.6% of the variance (F = 6.5, df = 4, 69, p <0.001) and was associated with silver-spoon (maternal and natal) environments. Specifically, body mass of bears was negatively affected by anomalies in summer (July-August) temperatures in the year prior to birth. During the year of birth, anomalies in summer growing season (May-October) temperatures, winter (December-March) precipitation and August precipitation affected springtime body mass measures: specifically, body mass was higher when summer temperature and winter precipitation during the birth year were above average. Above-average August precipitation resulted in lower observed body masses (Table 2). The maternal-summer temperature was the most important effect among inter-annual climate metrics on body mass. When considering local habitats, canopy variability was inversely related to body mass, although only an additional 1% of final model variance was explained (F = 4.3, df = 1, 68, p = 0.043; Table 2). Body mass was not effected by the presence of cubs, number of captures, or density. Spring capture date did not have an effect (model not shown) and was therefore not included in the mass or subsequent models.

Table 2. Standardized regression coefficients and significance ( p ) of model variables describing body mass (log scale), straight line length (log scale), and body condition measures of springtime grizzly bear captures in Alberta, Canada

thumbnailFigure 3. Model coefficient of determination (R2) for body mass (log[kg]), straight line length (log[SLL]), and body condition index (BCI). Hierarchically blocked variables were partitioned to represent different hypothesized biological or environmental factors. Only significant (p < 0.05) blocked variables are illustrated.

Body length

Biological and environmental factors explained 75.3% of the variation (R2, model F = 39.0, df = 7, 62, p < 0.001) in body length (Table 2, Figfure 3). Similar to body mass, age (as non-linear quadratic function) and sex explained a large amount (61.3%) of the variation in body length. Regional-habitat productivity explained an additional 6.6% of variation in body length (F = 13.3, df = 1, 65, p <0.001) based on average springtime (May-June) temperatures. Bears associated with warmer spring temperatures were more likely to be longer. Inter-annual climate variability – based on maternal and natal effects – explained an additional 4.2% of variance in body length (F = 4.7, df = 2, 63, p <0.001). Body length was positively related to warmer summer (July-August) temperatures during maternal periods and warmer spring temperatures during the year of birth (Table 2). Habitat quality and human footprint were not related to body length, but there was a positive association with landscape change (annual rate of change in habitats associated with human disturbances) adding an additional 3.2% of model variance explained (F = 8.1, df = 1, 62, p <0.001). Density, number of captures and human footprint did not influence body length.

Body condition

Biological and environmental factors explained 60.0% of model variation (R2, F = 14.7, df = 7, 68, p < 0.001) in springtime body condition (Table 2, Figure 3). Although body condition represents a standardized mass by length of animal, a non-linear (quadratic) age relationship with body condition was still apparent. Adult females were more likely to have a lower body condition than subadult or adult male bears, and this relationship was more pronounced if a female had cubs. Bears captured multiple times were in lower body condition than bears captured only once. Overall, the biological (including capture effects) base model accounted for 47.7% of the variance in body condition. Unlike mass and length measures, regional productivity did not affect body condition. Effects of inter-annual climate variability were observed with higher-than-normal July precipitation during the year of birth inversely related to body condition (Table 2):, this accounted for an additional 4.5% of the remaining model variance (F = 6.5, df = 1, 69, p = 0.013). Local habitat quality, as measured by use of habitats containing greater variation in regenerating forest age, was positively related to observed body condition (Table 2) and explained an additional 7.8% of model variation (F = 13.3, df = 1, 68, p = 0.001). Density of bears, human footprint, and landscape-change were not related to body condition.

Discussion

Biological factors and body size

Measurements of body mass and length of grizzly bears in Alberta were strongly dependent on intrinsic biological factors: age (positive, non-linear relationship) and sex (males > females). Age, sex and offspring dependence were important factors affecting body condition, which is a short-term measure of growth. Adult females, and especially adult females with cubs of the year, were likely to be in poorer condition than male bears. A negative effect of capture history (number of captures) was also observed for body condition measures which is consistent with previous observations [61]. Although population density (density dependence) is known to inversely affect body-size patterns in animals [80-82], no density dependent effects on body size patterns of grizzly bears were observed in our study. Grizzly bear populations in Alberta are likely to be below carrying capacity given locally high rates of human-caused mortality [83,84], and were recently classified by the province as ‘threatened’ given the low observed population densities [75]. This is in contrast to brown bears in Sweden that are considered healthy [85], but where body sizes of adult female bears are inversely related to population density [48].

Temporal and spatial environmental heterogeneity

Environmental heterogeneity is an important mechanism by which animal populations are regulated [86]. Here, we found that regional heterogeneity in habitat productivity was a moderate predictor of body size patterns of grizzly bears in Alberta. The smallest bears by mass and length occurred in the least-productive and coldest environments as measured by alpine habitat use and home ranges occupying both cool average spring temperatures and high average March precipitation (snowfall). In the Canadian Rocky Mountains, all three of these factors are associated with late timing of spring snowmelt and plant emergence, which are known to affect population dynamics of other alpine mammals [87]. Since den emergence in grizzly bears in our area typically occurs in April to early May [88], the amount and timing of spring snowpack is likely a factor affecting the availability of early season food resources such as roots [89], and generally might restrict access to early spring food resources.

Inter-annual variations in climate during the years’ prior, during and/or just following birth (maternal, in-utero and natal environments, respectively) also affected adult body size. Such silver-spoon effects by which animals that are born into ‘rich’ conditions are favoured throughout life are consistent with observations in other mammals including polar bears [90], Soay sheep [1], red squirrels [37] and caribou [91]. Common among these studies is the importance of winter and spring climate during (natal environments) or just prior (maternal or in utero environments) to the year of birth, which we also observed in this study. Winter and spring climate is related to summer drought conditions in the Canadian Rocky Mountains [92], which suggests that the effect of winter and spring climate may not necessarily be directly associated with the denning period, but rather summer environments when water is limiting. We are unsure, however, how late summer precipitation affects cubs-of-the-year. It may be related to late summer food resources, such as fruit production, or affect food-resource abundance in the following year when bears are yearlings. Further, winter precipitation (December-March) anomalies during the natal birth year were positively related to body mass. We interpreted this as snow cover during winter denning providing energetic benefits (e.g. insulation) in the den for cubs of the year.

During the year prior to birth, late summer (July-August) temperature anomalies were negatively associated with body mass but positively associated with body length in grizzly bears. This late-summer environment may have affected maternal body condition prior to denning and thus subsequent condition of offspring [e.g. 53] or conversely, it may have affected the following years’ food supply during the cub-of-year period, since lag effects in fruit production are caused by weather conditions favourable to flower primordia in the mid-to-late summer period the year prior to fruiting [93]. Although we cannot be certain which factor is more important, the fact that body mass is negatively associated with late-summer temperature anomalies, where as body length is positively associated with late-summer temperature anomalies suggests to us that maternal condition is less likely (as we would expect similar responses in body mass and length if it were solely a maternal effect). Further investigations of mid and late-summer weather on pulsing in food resource abundance the following year are needed, especially in regard to the apparent opposite effects on bear mass and length.

One important consideration to our purported silver spoon effect should be discussed: that is, we have no information on our study animals prior to their first capture. This has two important implications: 1) we cannot account for litter size effects, and 2) the centroid data used to determine natal climatic conditions may not be reflective of the actual natal location. In regards to the former, not accounting for litter size should inflate the variance around our estimates. For the centroid data, this would likely only influence dispersing males, as females are philopatric [94]. For males, average dispersal distances in the province are under 50 kilometers [94], thus still largely reflective of the climate in the centroid of the current home range (differences in climates among bears are mainly regional in effect, not within populations). Further, for this limitation to bias our results, males would consistently have to disperse to poorer environments, again something we deem unlikely. Thus, we argue that the silver spoon pattern is unlikely to be altered by these factors in such away that the statistical pattern would disappear.

Anthropogenic considerations

Human footprint did not directly relate to body size patterns of grizzly bears, but human activity indirectly affected body size by influencing habitats. The two most important measures of habitat quality were canopy closure and the age structure of forests. Bears that used habitats associated with higher canopy variability, such as forest/non-forest landscapes in the mountains or expanses of old growth forests with a recent, single-harvest sequence, had lower body masses. Conversely, bears that used forests with higher variability in regenerating forest age had higher body condition. Likewise, body length was positively related to annual landscape change. Taken together, these results suggest that human activities that fragment forests are positively associated with body size measures, although survival of bears in these environments is compromised due to high rates of human-caused mortalities [57,84]. Early successional and highly variable forests are therefore important indicators of improved habitat quality for bears given the relationship to body size patterns reported here, habitat use studies [72] and measures of food resource abundance [73,74]. We hypothesize that positive associations between body size patterns and variability in regenerating forest age are due in part to local landscape patterns in protein availability. For instance, both ungulate and ant resource use in Alberta are associated with disturbed forests [46,74].

Conclusions

While bear body size is largely dictated by age and sex, it only accounted for about 50% of the variation. More consideration of the spatial and temporal patterns of resource availability, including the conditions early in life, is needed to better understand individual performance of animals and population dynamics. For grizzly bears in Alberta, environmental effects on body size are most affected by regional environmental gradients (space) and the environmental conditions animals are born into (time). Local-habitat heterogeneity (particularly young, patchily disturbed forests), and landscape dynamics also had a small influence on body size. It is important to emphasize that while patchily disturbed forests positively affected body size, these areas also have high rates of mortality, which could negate any positive population-level effect.

Worldwide, relationships between carnivore body size and climate warming show ambiguous trends [95]; however, polar bears body sizes have recently declined, which has been attributed primarily to loss in habitat (i.e., sea ice as a platform for hunting; [96,97]). Despite unequivocal global patterns [95], a 50 year examination of regional studies showed that carnivore body sizes have generally increased over the past half century [98]. Given the short season associated with high-alpine environments, such as the Rocky Mountains in Alberta, we hypothesize that individuals with a limited growing season and temperature-limited ecosystems, such as interior grizzly bears, might actually benefit from increases in season length associated with climate change. This prediction is largely consistent with observed body size and seasonality patterns in grizzly bears across North America [40], but may be dependent on sufficient snow cover during the denning period. In conclusion, we have demonstrated a complex interplay of biological, spatial and temporal factors on body size that collectively explained between 60 and 84% of the variation seen in Alberta’s grizzly bears.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: SEN MRLC GJM GBS. Performed the experiments (animal capture and body size measurements): MRLC GBS. Contributed data/analysis tools: SEN MRLC GBS JC GJM. Analyzed the data: SEN JB MRLC JC GBS. Wrote the paper: SEN ABAS. Critically reviewed and revised paper: ABAS MRLC GBS JB JC GJM. All authors read and approved the final manuscript.

Acknowledgements

We thank the Natural Sciences and Engineering Research Council (NSERC) of Canada, Alberta Innovates – Bio Solutions, and partners from the Foothills Research Institute Grizzly Bear Program for financial support. This manuscript was greatly improved from the comments of three anonymous reviewers.

References

  1. Forchhammer MC, Clutton-Brock TH, Lindström J, Albon SD: Climate and population density induce long‒term cohort variation in a northern ungulate.

    J Anim Ecol 2001, 70:721-729. Publisher Full Text OpenURL

  2. Krebs C, Boutin S, Boonstra R, Sinclair ARE, Smith JNM, Dale M, Turkington R: Impact of food and predation on the snowshoe hare cycle.

    Science 1995, 269:1112-1115. PubMed Abstract | Publisher Full Text OpenURL

  3. Elton C: Animal ecology. London: Macmillan; 1927. OpenURL

  4. Dooley JL Jr, Bowers MA: Demographic responses to habitat fragmentation: experimental tests at the landscape and patch scale.

    Ecology 1998, 79:969-980. Publisher Full Text OpenURL

  5. Andren H: Corvid density and nest predation in relation to forest fragmentation: A landscape perspective.

    Ecology 1992, 73:794. Publisher Full Text OpenURL

  6. Flaspohler DJ, Temple SA, Rosenfield RN: Effects of forest edges on ovenbird demography in a managed forest landscape.

    Conserv Biol 2001, 15:173-183. OpenURL

  7. Wiens JA, Chr N, Van Horne B, Ims RA: Ecological mechanisms and landscape ecology. New York: Oikos; 1993. OpenURL

  8. Forman R: Some general principles of landscape and regional ecology.

    Landscape Ecol 1995, 10:133-142. Publisher Full Text OpenURL

  9. Fretwell SD, Lucas JRHL: On territorial behavior and other factors influencing habitat distribution in birds I. Theoretical Development.

    Acta Biotheor 1970, 19:16-36. OpenURL

  10. Pulliam HR, Danielson BJ: Sources, sinks, and habitat selection: a landscape perspective on population dynamics.

    Am Nat 1991, 137:S50-S66. Publisher Full Text OpenURL

  11. Anderson WB, Wait DA, Stapp P: Resources from another place and time: responses to pulses in a spatially subsidized system.

    Ecology 2008, 89:660-670. PubMed Abstract | Publisher Full Text OpenURL

  12. Schmidt KA, Ostfeld RS: Numerical and behavioral effects within a pulse-driven system: consequences for shared prey.

    Ecology 2008, 89:635-646. PubMed Abstract | Publisher Full Text OpenURL

  13. Yang LH: Pulses of dead periodical cicadas increase herbivory of American bellflowers.

    Ecology 2008, 89:1497-1502. PubMed Abstract | Publisher Full Text OpenURL

  14. Woodward FI, Lomas MR, Quaife T: Global responses of terrestrial productivity to contemporary climatic oscillations.

    Phil Trans R Soc B 2008, 363:2779-2785. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Tian H, Melillo JM, Kicklighter DW, McGuire AD, Helfrich JVK, Moore B, Vörösmarty CJ: Effect of interannual climate variability on carbon storage in Amazonian ecosystems.

    Nature 1998, 396:664-667. Publisher Full Text OpenURL

  16. Black TA, Chen WJ, Barr AG, Arain MA, Chen Z, Nesic Z, Hogg EH, Neumann HH, Yang PC: Increased carbon sequestration by a boreal deciduous forest in years with a warm spring.

    Geophys Res Lett 2012, 27:1271-1274. OpenURL

  17. Macias Fauria M, Johnson EA: Climate and wildfires in the North American boreal forest.

    Philos Trans of the R Soc B: Biol Sci 2008, 363:2315-2327. Publisher Full Text OpenURL

  18. Macias Fauria M, Johnson EA: Large‒scale climatic patterns and area affected by mountain pine beetle in British Columbia, Canada.

    J of Geo Res: Biogeosci 2009, 114:G01012. OpenURL

  19. Post E, Stenseth NC: Climatic variability, plant phenology, and northern ungulates.

    Ecology 1999, 80:1322-1339. Publisher Full Text OpenURL

  20. Sæther BE, Smith FM, Cooper EJ, Wookey PA: The Arctic Oscillation predicts effects of climate change in two trophic levels in a high‒arctic ecosystem.

    Ecol Lett 2002, 5:445-453. Publisher Full Text OpenURL

  21. Haynes KJ, Liebhold AM, Fearer TM, Wang G, Norman GW, Johnson DM: Spatial synchrony propagates through a forest food web via consumer-resource interactions.

    Ecology 2009, 90:2974-2983. PubMed Abstract | Publisher Full Text OpenURL

  22. Warne RW, Pershall AD, Wolf BO: Linking precipitation and C3-C4 plant production to resource dynamics in higher-trophic-level consumers.

    Ecology 2010, 91:1628-1638. PubMed Abstract | Publisher Full Text OpenURL

  23. Herrera CM: Long-term dynamics of mediterranean frugivorousbirds and fleshy fruits: a 12-year study.

    Ecol Monogr 1998, 68:511-538. OpenURL

  24. Abrahamson WG, Layne JN: Long-term patterns of acorn production for five oak species in xeric florida uplands.

    Ecology 2003, 84:2476-2492. Publisher Full Text OpenURL

  25. Curran LM, Leighton M: Vertebrate responses to spatiotemporal variation in seed production of mast-fruiting dipterocarpaceae.

    Ecol Monogr 2000, 70:101-128. Publisher Full Text OpenURL

  26. Harrison RD: Drought and the consequences of El Niño in Borneo: a case study of figs.

    Popul Ecol 2001, 43:63-75. Publisher Full Text OpenURL

  27. Stenseth NC, Mysterud A, Ottersen G, Hurrell JW, Chan K-S, Lima M: Ecological effects of climate fluctuations.

    Science 2002, 297:1292-1296. PubMed Abstract | Publisher Full Text OpenURL

  28. Howe EJ, Obbard ME, Bowman J: Prior reproduction and weather affect berry crops in central Ontario, Canada.

    Popul Ecol 2011, 54:347-356. OpenURL

  29. Wright SJ, Carrasco C, Calderón O, Paton S: The El Niño Southern Oscillation, variable fruit production, and famine in a tropical forest.

    Ecology 1999, 80:1632-1647. OpenURL

  30. Koenig WD, Knops J: Scale of mast-seeding and tree-ring growth.

    Nature 1998, 396:225-226. OpenURL

  31. Schauber EM, Kelly D, Turchin P, Simon C, Lee WG: Masting by eighteen New Zealand plant species: the role of temperature as a synchronizing cue.

    Ecology 2002, 83:1214-1225. Publisher Full Text OpenURL

  32. Peters VS, Macdonald SE, Dale M: The interaction between masting and fire is key to white spruce regeneration.

    Ecology 2005, 86:1744-1750. Publisher Full Text OpenURL

  33. Kemp GA, Keith LB: Dynamics and regulation of red squirrel (Tamiasciurus hudsonicus) populations.

    Ecology 1970, 51:763-779. Publisher Full Text OpenURL

  34. Bieber C, Ruf T: Population dynamics in wild boar Sus scrofa: ecology, elasticity of growth rate and implications for the management of pulsed resource consumers.

    J Appl Ecol 2005, 42:1203-1213. Publisher Full Text OpenURL

  35. Madsen T, Shine R: Silver spoons and snake body sizes: prey availability early in life influences long‒term growth rates of free‒ranging pythons.

    J Anim Ecol 2000, 69:952-958. Publisher Full Text OpenURL

  36. Van de Pol M, Bruinzeel LW, Heg D, Van der Jeugd HP, Verhulst S: A silver spoon for a golden future: long-term effects of natal origin on fitness prospects of oystercatchers (Haematopus ostralegus).

    J Anim Ecol 2006, 75:616-626. PubMed Abstract | Publisher Full Text OpenURL

  37. Descamps S, Boutin S, Berteaux D, McAdam AG, Gaillard J-M: Cohort effects in red squirrels: the influence of density, food abundance and temperature on future survival and reproductive success.

    J Anim Ecol 2008, 77:305-314. PubMed Abstract | Publisher Full Text OpenURL

  38. Wilkin TA, Sheldon BC: Sex differences in the persistence of natal environmental effects on life histories.

    Curr Biol 2009, 19:1998-2002. PubMed Abstract | Publisher Full Text OpenURL

  39. Wittmer HU, Powell RA, King CM: Understanding contributions of cohort effects to growth rates of fluctuating populations.

    J Anim Ecol 2007, 76:946-956. PubMed Abstract | Publisher Full Text OpenURL

  40. Ferguson SH, McLoughlin PD: Effect of energy availability, seasonality, and geographic range on brown bear life history.

    Ecography 2008, 23:193-200. OpenURL

  41. Nowak RM: Walker’s Mammals of the world. 6th edition. Baltimore, Maryland: The Johns Hopkins University Press; 1999. PubMed Abstract | Publisher Full Text OpenURL

  42. Gende SM, Quinn TP, Willson MF: Consumption choice by bears feeding on salmon.

    Oecologia 2001, 127:372-382. Publisher Full Text OpenURL

  43. Boulanger J, Himmer S, Swan C: Monitoring of grizzly bear population trends and demography using DNA mark–recapture methods in the Owikeno Lake area of British Columbia.

    Can J Zool 2004, 82:1267-1277. Publisher Full Text OpenURL

  44. Mowat G, Heard DC: Major components of grizzly bear diet across North America.

    Can J Zool 2006, 84:473-489. Publisher Full Text OpenURL

  45. Mclellan BN, Hovey FW: The diet of grizzly bears in the Flathead River drainage of southeastern British Columbia.

    Can J Zool 1995, 73:704-712. Publisher Full Text OpenURL

  46. Munro RHM, Nielsen SE, Price MH, Stenhouse GB, Boyce MS: Seasonal and diel patterns of grizzly bear diet and activity in west-central Alberta.

    J Mamm 2006, 87:1112-1121. Publisher Full Text OpenURL

  47. Naves J, Fernández-Gil A, Rodríguez C, Delibes M: Brown bear food habits at the border of its range: a long-term study.

    J Mamm 2006, 87:899-908. Publisher Full Text OpenURL

  48. Zedrosser A, Dahle B, Swenson JE: Population density and food conditions determine adult female body size in brown bears.

    J Mamm 2006, 87:510-518. Publisher Full Text OpenURL

  49. Meiri S, Yom-Tov Y, Geffen E: What determines conformity to Bergmann’s rule?

    Global Ecol Biogeogr 2007, 16:788-794. Publisher Full Text OpenURL

  50. Mattson DJ, Blanchard BM, Knight RR: Yellowstone grizzly bear mortality, human habituation, and whitebark pine seed crops.

    J Wildl Manage 1992, 56:432-444. Publisher Full Text OpenURL

  51. Gunther KA, Haroldson MA, Frey K, Cain SL, Copeland J, Schwartz CC: Grizzly bear–human conflicts in the Greater Yellowstone ecosystem, 1992–2000.

    Ursus 2004, 15:10-22. Publisher Full Text OpenURL

  52. Hilderbrand GV, Schwartz CC, Robbins CT, Jacoby ME, Hanley TA, Arthur SM, Servheen C: The importance of meat, particularly salmon, to body size, population productivity, and conservation of North American brown bears.

    Can J Zool 1999, 77:132-138. Publisher Full Text OpenURL

  53. Zedrosser A, Bellemain E, Taberlet P, Swenson JE: Genetic estimates of annual reproductive success in male brown bears: the effects of body size, age, internal relatedness and population density.

    J Anim Ecol 2007, 76:368-375. PubMed Abstract | Publisher Full Text OpenURL

  54. Robbins CT, Ben-David M, Fortin JK, Nelson OL: Maternal condition determines birth date and growth of newborn bear cubs.

    J Mamm 2012, 93:540-546. Publisher Full Text OpenURL

  55. Hamlett G: Delayed implantation and discontinuous development in the mammals.

    Q Rev Biol 1935, 10:432-447. Publisher Full Text OpenURL

  56. Spady TJ, Lindburg DG, Durrant BS: Evolution of reproductive seasonality in bears.

    Mamm Rev 2007, 37:21-53. Publisher Full Text OpenURL

  57. Nielsen SE, Stenhouse GB, Beyer HL, Huettmann F, Boyce MS: Can natural disturbance-based forestry rescue a declining population of grizzly bears?

    Biol Conserv 2008, 141:2193-2207. Publisher Full Text OpenURL

  58. Nielsen SE, Stenhouse GB, Boyce MS: A habitat-based framework for grizzly bear conservation in Alberta.

    Biol Conserv 2006, 130:217-229. Publisher Full Text OpenURL

  59. Linke J, Franklin SE, Huettmann F, Stenhouse GB: Seismic cutlines, changing landscape metrics and grizzly bear landscape use in alberta.

    Landscape Ecol 2005, 20:811-826. Publisher Full Text OpenURL

  60. Cattet MR, Christison K, Caulkett NA, Stenhouse GB: Physiologic responses of grizzly bears to different methods of capture.

    J Wildlife Dis 2003, 39:649-654. OpenURL

  61. Cattet M, Boulanger J, Stenhouse G, Powell RA, Reynolds-Hogland MJ: An evaluation of long-term capture effects in ursids: implications for wildlife welfare and research.

    J Mamm 2008, 89:973-990. Publisher Full Text OpenURL

  62. Cattet M, Stenhouse G, Bollinger T: Exertional myopathy in a grizzly bear (Ursus arctos) captured by leghold snare.

    J Wildlife Dis 2008, 44:973-978. OpenURL

  63. Cattet M, Caulkett NA, Stenhouse GB: Anesthesia of grizzly bears using xylazine-zolazepam-tiletamine or zolazepam-tiletamine.

    Ursus 2003, 14:88-93. OpenURL

  64. Stoneberg RP, Jonkel CJ: Age determination of black bears by cementum layers.

    J Wildl Manage 1966, 30:411-414. Publisher Full Text OpenURL

  65. Cattet MRL, Caulkett NA, Obbard ME, Stenhouse GB: A body-condition index for ursids.

    Can J Zool 2002, 80:1156-1161. Publisher Full Text OpenURL

  66. Gannon WL, Sikes RS: Guidelines of the American Society of Mammalogists for the use of wild mammals in research.

    J Mamm 2007, 88:809-823. Publisher Full Text OpenURL

  67. CANADIAN COUNCIL ON ANIMAL CARE: CCAC guidelines on: the care and use of wildlife. Ottawa, Ontario, Canada: Canadian Council on Animal Care; 2003. PubMed Abstract | Publisher Full Text OpenURL

  68. Boulanger J, Stenhouse G, Munro R: Sources of heterogeneity bias when dna mark-recapture sampling methods are applied to grizzly bear (Ursus arctos) populations.

    J Mamm 2004, 85:618-624. Publisher Full Text OpenURL

  69. Wang T, Hamann A, Mbogga M:

    ClimateAB v3.21: a program to generate projected, normal, decadel, annual, seasonal and monthly interpolated climate data for Alberta. OpenURL

  70. Hamann A, Wang TL: Models of climatic normals for genecology and climate change studies in British Columbia.

    Agr Forest Meteorol 2005, 128:211-221. Publisher Full Text OpenURL

  71. Archibald WR, Ellis R, Hamilton AN: Responses of grizzly bears to logging truck traffic in the Kimsquit River Valley, British Columbia.

    Bears: Their Biol and Manag 1987, 7:251-257. OpenURL

  72. Nielsen SE, Boyce MS, Stenhouse GB: Grizzly bears and forestry.

    For Ecol Manage 2004, 199:51-65. Publisher Full Text OpenURL

  73. Nielsen SE, Munro RHM, Bainbridge EL, Stenhouse GB, Boyce MS: Grizzly bears and forestry.

    For Ecol Manage 2004, 199:67-82. Publisher Full Text OpenURL

  74. Nielsen SE, McDermid G, Stenhouse GB, Boyce MS: Dynamic wildlife habitat models: Seasonal foods and mortality risk predict occupancy-abundance and habitat selection in grizzly bears.

    Biol Conserv 2010, 143:1623-1634. Publisher Full Text OpenURL

  75. Nielsen SE, Cranston J, Stenhouse GB: Identification of priority areas for grizzly bear conservation and recovery in Alberta, Canada.

    J of Conserv Plann 2009, 5:38-60. OpenURL

  76. Linke J, McDermid GJ, Laskin DN, McLane AJ, Pape A, Cranston J, Hall-Beyer M, Franklin SE: A disturbance-inventory framework for flexible and reliable landscape monitoring.

    Photogramm Eng Remote Sens 2009, 75:981-996. OpenURL

  77. Bern PH:

    HIREG: Stata module for hierarchial regression. 2005. [Statistical software components]

    http://EconPapers.repec.org/RePEc:boc:bocode:s432904 webcite

    OpenURL

  78. Raudenbush B, Bryk AS: Hierarchical linear models: applications and data analysis methods. 2nd edition. Thousand Oaks, California: Sage Publications; 2001. OpenURL

  79. Hosmer DW, Lemeshow S: Applied logistic regression. New York, New York: John Wiley & Sons, Ltd; 2000. OpenURL

  80. Wilbur HM: Density-dependent aspects of growth and metamorphosis in Bufo americanus.

    Ecology 1977, 58:196-200. Publisher Full Text OpenURL

  81. Skogland T: The effects of density dependent resource limitation on size of wild reindeer.

    Oecologia 1983, 60:156-168. Publisher Full Text OpenURL

  82. 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 Biol Sci 2001, 268:911-919. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  83. Benn B, Herrero S: Grizzly bear mortality and human access in Banff and Yoho National Parks, 1971–98.

    Ursus 2002, 13:213-221. OpenURL

  84. Nielsen SE, Herrero S, Boyce MS, Mace RD, Benn B, Gibeau ML, Jevons S: Modelling the spatial distribution of human-caused grizzly bear mortalities in the Central Rockies ecosystem of Canada.

    Biol Conserv 2004, 120:101-113. Publisher Full Text OpenURL

  85. Zedrosser A, Dahle B, Swenson JE, Gerstl N: Status and management of the brown bear in Europe.

    Ursus 2001, 12:9-20. OpenURL

  86. 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-485. PubMed Abstract | Publisher Full Text OpenURL

  87. Morrison SF, Hik DS: Demographic analysis of a declining pika Ochotona collaris population: linking survival to broad-scale climate patterns via spring snowmelt patterns.

    J Anim Ecol 2007, 76:899-907. PubMed Abstract | Publisher Full Text OpenURL

  88. Ciarniello LM, Boyce MS, Heard DC, Seip DR: Denning behavior and den site selection of grizzly bears along the Parsnip River, British Columbia, Canada.

    Ursus 2005, 16:47-58. Publisher Full Text OpenURL

  89. Coogan S, Nielsen SE, Stenhouse GB: Spatial and temporal heterogeneity creates a “brown tide” in root phenology and nutrition.

    ISRN Ecology 2012, 2012:10. OpenURL

  90. Atkinson SN, Stirling I, Ramsay MA: Growth in early life and relative body size among adult polar bears ( Ursus maritimus).

    J Zool 1996, 239:225-234. Publisher Full Text OpenURL

  91. Hegel TM, Mysterud A, Ergon T, Loe LE, Huettmann F, Stenseth NC: Seasonal effects of Pacific-based climate on recruitment in a predator-limited large herbivore.

    J Anim Ecol 2010, 79:471-482. PubMed Abstract | Publisher Full Text OpenURL

  92. Meyn A, Taylor SW, Flannigan MD, Thonicke K, Cramer W: Relationship between fire, climate oscillations, and drought in British Columbia, Canada, 1920–2000.

    Glob Chang Biol 2010, 16:977-989. Publisher Full Text OpenURL

  93. Krebs CJ, Boonstra R, Cowcill K, Kenney AJ: Climatic determinants of berry crops in the boreal forest of the southwestern Yukon.

    Botany 2009, 87:401-408. Publisher Full Text OpenURL

  94. Proctor MF, Mclellan BN, Strobeck C, Barclay RMR: Gender-specific dispersal distances of grizzly bears estimated by genetic analysis.

    Can J Zool 2004, 82:1108-1118. Publisher Full Text OpenURL

  95. Meiri S, Guy D, Dayan T, Simberloff D: Global change and carnivore body size: data are stasis.

    Global Ecol Biogeogr 2009, 18:240-247. Publisher Full Text OpenURL

  96. Derocher AE, Stirling I: Temporal variation in reproduction and body mass of polar bears in western Hudson Bay.

    Can J Zool 1995, 73:1657-1665. Publisher Full Text OpenURL

  97. Rode KD, Amstrup SC, Regehr EV: Reduced body size and cub recruitment in polar bears associated with sea ice decline.

    Ecol Appl 2010, 20:768-782. PubMed Abstract | Publisher Full Text OpenURL

  98. Yom-Tov Y: Body sizes of carnivores commensal with humans have increased over the past 50 years.

    Func Ecol 2003, 17:323-327. Publisher Full Text OpenURL