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.
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.
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.
BearSilver spoonEnvironmental effectsGPS radiocollarTemporal and spatial heterogeneity
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 . 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 . Likewise, population growth of stoats in New Zealand beech forests is dependent on masting .
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.) . 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 . 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.
Environmental variables used to measure hypothesized environmental drivers of body size patterns in grizzly bears within Alberta, Canada
Hypothesized environmental driver and measurement variable
A. Regional habitat productivity
Temperature (Winter, Spring, Summer)
Precipitation (Winter, Spring, Summer)
B. Inter-annual environments (deviations)
Temperature (Winter, Spring, Summer)
Bt-1, t 0, t+1 & Ct-1, t 0
Precipitation (Winter, Spring, Summer)
Bt-1, t 0, t+1 & Ct-1, t 0
C. Local habitat quality
Shrub habitat (quadratic)
Canopy cover (quadratic)
Variation in canopy cover
Deciduous canopy cover (quadratic)
Forest age (quadratic)
Forest age variation
Regenerating forest habitat (quadratic)
Variation in regen. forest age
Soil wetness (quadratic)
D. Human footprint & activity
Safe harbour habitat
Linear feature density
Distance to human feature
Distance to active energy well
E. Landscape change
Annual rate of habitat change
§ Home ranges estimated by 50% multi-annual kernels; climate variables measured at kernel centroid; † Temporal scales relate to time of measurements; B relates to birth year & C to capture year. For inter-annual variation, 2-yrs prior to and up to 1-yr following birth or 1-yr prior to and the year of capture are considered.
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].
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 , 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 . 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 , 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.
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 . 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  and the Canadian Council on Animal Care  for the safe handling of wildlife.
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 .
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 . 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 . 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 , 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  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 (Ct 0); 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 . 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 ), protected areas and high- or low-risk habitats based on a mortality risk and safe harbour habitat models , 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  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 .
We used the HIREG module  for the software STATA 11 to estimate hierarchical regressions  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 . 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.
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.
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
Block (hypothesized) category and measurement variables
1) Biology and capture effects
Adult Females (AF)
Adult F w/ cubs (AFC)
Male x Age
Number of captures
2) Regional habitat productivity
Spring (May-Jun) temperature
Alpine habitat use (HP)
3) Inter-annual climate variability
Maternal effects (Bt-1):
Summer (Jul-Aug) temperature
Natal effects (Bt0):
Spring (May-Jun) temperature
Summer (May-Oct) temperature
Winter (Dec-Mar) precipitation
Capture effects (Ct):
4) Local habitat quality
Canopy variation (HP)
Regen. forest age variation (HP)
5) Human footprint
6) Landscape change
All measures of habitat use were based on global position system (GPS) telemetry data and relate to a habitat patch (HP) scale of a 30 m pixel (900 m2).
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.
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.
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 . 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 . This is in contrast to brown bears in Sweden that are considered healthy , but where body sizes of adult female bears are inversely related to population density .
Temporal and spatial environmental heterogeneity
Environmental heterogeneity is an important mechanism by which animal populations are regulated . 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 . Since den emergence in grizzly bears in our area typically occurs in April to early May , the amount and timing of spring snowpack is likely a factor affecting the availability of early season food resources such as roots , 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 , Soay sheep , red squirrels  and caribou . 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 , 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 . 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 . For males, average dispersal distances in the province are under 50 kilometers , 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.
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  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].
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 ; 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 , a 50 year examination of regional studies showed that carnivore body sizes have generally increased over the past half century . 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 , 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.
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.
Department of Renewable Resources, University of Alberta
Canadian Cooperative Wildlife Health Centre, University of Saskatchewan
Integrated Ecological Research
Arctos Ecological Services
Department of Geography, University of Calgary
Department of Ecology and Genetics, Uppsala University
Foothills Research Institute
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. 10.1046/j.0021-8790.2001.00532.x.View Article
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. 10.1126/science.269.5227.1112.View ArticlePubMed
Elton C: Animal ecology. 1927, London: Macmillan
Dooley JL, Bowers MA: Demographic responses to habitat fragmentation: experimental tests at the landscape and patch scale. Ecology. 1998, 79: 969-980. 10.1890/0012-9658(1998)079[0969:DRTHFE]2.0.CO;2.View Article
Andren H: Corvid density and nest predation in relation to forest fragmentation: A landscape perspective. Ecology. 1992, 73: 794-10.2307/1940158.View Article
Flaspohler DJ, Temple SA, Rosenfield RN: Effects of forest edges on ovenbird demography in a managed forest landscape. Conserv Biol. 2001, 15: 173-183.View Article
Wiens JA, Chr N, Van Horne B, Ims RA: Ecological mechanisms and landscape ecology. 1993, New York: Oikos
Forman R: Some general principles of landscape and regional ecology. Landscape Ecol. 1995, 10: 133-142. 10.1007/BF00133027.View Article
Fretwell SD, Lucas JRHL: On territorial behavior and other factors influencing habitat distribution in birds I. Theoretical Development. Acta Biotheor. 1970, 19: 16-36.View Article
Pulliam HR, Danielson BJ: Sources, sinks, and habitat selection: a landscape perspective on population dynamics. Am Nat. 1991, 137: S50-S66. 10.1086/285139.View Article
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. 10.1890/07-0234.1.View ArticlePubMed
Schmidt KA, Ostfeld RS: Numerical and behavioral effects within a pulse-driven system: consequences for shared prey. Ecology. 2008, 89: 635-646. 10.1890/07-0199.1.View ArticlePubMed
Yang LH: Pulses of dead periodical cicadas increase herbivory of American bellflowers. Ecology. 2008, 89: 1497-1502. 10.1890/07-1853.1.View ArticlePubMed
Woodward FI, Lomas MR, Quaife T: Global responses of terrestrial productivity to contemporary climatic oscillations. Phil Trans R Soc B. 2008, 363: 2779-2785. 10.1098/rstb.2008.0017.PubMed CentralView ArticlePubMed
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. 10.1038/25328.View Article
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.View Article
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. 10.1098/rstb.2007.2202.View Article
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-View Article
Post E, Stenseth NC: Climatic variability, plant phenology, and northern ungulates. Ecology. 1999, 80: 1322-1339. 10.1890/0012-9658(1999)080[1322:CVPPAN]2.0.CO;2.View Article
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. 10.1046/j.1461-0248.2002.00340.x.View Article
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. 10.1890/08-1709.1.View ArticlePubMed
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. 10.1890/08-1471.1.View ArticlePubMed
Herrera CM: Long-term dynamics of mediterranean frugivorousbirds and fleshy fruits: a 12-year study. Ecol Monogr. 1998, 68: 511-538.
Abrahamson WG, Layne JN: Long-term patterns of acorn production for five oak species in xeric florida uplands. Ecology. 2003, 84: 2476-2492. 10.1890/01-0707.View Article
Curran LM, Leighton M: Vertebrate responses to spatiotemporal variation in seed production of mast-fruiting dipterocarpaceae. Ecol Monogr. 2000, 70: 101-128. 10.1890/0012-9615(2000)070[0101:VRTSVI]2.0.CO;2.View Article
Harrison RD: Drought and the consequences of El Niño in Borneo: a case study of figs. Popul Ecol. 2001, 43: 63-75. 10.1007/PL00012017.View Article
Stenseth NC, Mysterud A, Ottersen G, Hurrell JW, Chan K-S, Lima M: Ecological effects of climate fluctuations. Science. 2002, 297: 1292-1296. 10.1126/science.1071281.View ArticlePubMed
Howe EJ, Obbard ME, Bowman J: Prior reproduction and weather affect berry crops in central Ontario, Canada. Popul Ecol. 2011, 54: 347-356.View Article
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.
Koenig WD, Knops J: Scale of mast-seeding and tree-ring growth. Nature. 1998, 396: 225-226.View Article
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. 10.1890/0012-9658(2002)083[1214:MBENZP]2.0.CO;2.View Article
Peters VS, Macdonald SE, Dale M: The interaction between masting and fire is key to white spruce regeneration. Ecology. 2005, 86: 1744-1750. 10.1890/03-0656.View Article
Kemp GA, Keith LB: Dynamics and regulation of red squirrel (Tamiasciurus hudsonicus) populations. Ecology. 1970, 51: 763-779. 10.2307/1933969.View Article
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. 10.1111/j.1365-2664.2005.01094.x.View Article
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. 10.1046/j.1365-2656.2000.00477.x.View Article
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. 10.1111/j.1365-2656.2006.01079.x.View ArticlePubMed
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. 10.1111/j.1365-2656.2007.01340.x.View ArticlePubMed
Wilkin TA, Sheldon BC: Sex differences in the persistence of natal environmental effects on life histories. Curr Biol. 2009, 19: 1998-2002. 10.1016/j.cub.2009.09.065.View ArticlePubMed
Wittmer HU, Powell RA, King CM: Understanding contributions of cohort effects to growth rates of fluctuating populations. J Anim Ecol. 2007, 76: 946-956. 10.1111/j.1365-2656.2007.01274.x.View ArticlePubMed
Ferguson SH, McLoughlin PD: Effect of energy availability, seasonality, and geographic range on brown bear life history. Ecography. 2008, 23: 193-200.View Article
Nowak RM: Walker’s Mammals of the world. 1999, Baltimore, Maryland: The Johns Hopkins University Press, 6
Gende SM, Quinn TP, Willson MF: Consumption choice by bears feeding on salmon. Oecologia. 2001, 127: 372-382. 10.1007/s004420000590.View Article
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. 10.1139/z04-100.View Article
Mowat G, Heard DC: Major components of grizzly bear diet across North America. Can J Zool. 2006, 84: 473-489. 10.1139/z06-016.View Article
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. 10.1139/z95-082.View Article
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. 10.1644/05-MAMM-A-410R3.1.View Article
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. 10.1644/05-MAMM-A-318R2.1.View Article
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. 10.1644/05-MAMM-A-218R1.1.View Article
Meiri S, Yom-Tov Y, Geffen E: What determines conformity to Bergmann’s rule?. Global Ecol Biogeogr. 2007, 16: 788-794. 10.1111/j.1466-8238.2007.00330.x.View Article
Mattson DJ, Blanchard BM, Knight RR: Yellowstone grizzly bear mortality, human habituation, and whitebark pine seed crops. J Wildl Manage. 1992, 56: 432-444. 10.2307/3808855.View Article
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. 10.2192/1537-6176(2004)015<0010:GBCITG>2.0.CO;2.View Article
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. 10.1139/z98-195.View Article
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. 10.1111/j.1365-2656.2006.01203.x.View ArticlePubMed
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. 10.1644/11-MAMM-A-155.1.View Article
Hamlett G: Delayed implantation and discontinuous development in the mammals. Q Rev Biol. 1935, 10: 432-447. 10.1086/394493.View Article
Spady TJ, Lindburg DG, Durrant BS: Evolution of reproductive seasonality in bears. Mamm Rev. 2007, 37: 21-53. 10.1111/j.1365-2907.2007.00096.x.View Article
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. 10.1016/j.biocon.2008.06.020.View Article
Nielsen SE, Stenhouse GB, Boyce MS: A habitat-based framework for grizzly bear conservation in Alberta. Biol Conserv. 2006, 130: 217-229. 10.1016/j.biocon.2005.12.016.View Article
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. 10.1007/s10980-005-0066-4.View Article
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.View Article
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. 10.1644/08-MAMM-A-095.1.View Article
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.View Article
Cattet M, Caulkett NA, Stenhouse GB: Anesthesia of grizzly bears using xylazine-zolazepam-tiletamine or zolazepam-tiletamine. Ursus. 2003, 14: 88-93.
Stoneberg RP, Jonkel CJ: Age determination of black bears by cementum layers. J Wildl Manage. 1966, 30: 411-414. 10.2307/3797828.View Article
Cattet MRL, Caulkett NA, Obbard ME, Stenhouse GB: A body-condition index for ursids. Can J Zool. 2002, 80: 1156-1161. 10.1139/z02-103.View Article
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. 10.1644/06-MAMM-F-185R1.1.View Article
CANADIAN COUNCIL ON ANIMAL CARE: CCAC guidelines on: the care and use of wildlife. 2003, Ottawa, Ontario, Canada: Canadian Council on Animal Care
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. 10.1644/BRB-134.View Article
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.
Hamann A, Wang TL: Models of climatic normals for genecology and climate change studies in British Columbia. Agr Forest Meteorol. 2005, 128: 211-221. 10.1016/j.agrformet.2004.10.004.View Article
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.
Nielsen SE, Boyce MS, Stenhouse GB: Grizzly bears and forestry. For Ecol Manage. 2004, 199: 51-65. 10.1016/j.foreco.2004.04.014.View Article
Nielsen SE, Munro RHM, Bainbridge EL, Stenhouse GB, Boyce MS: Grizzly bears and forestry. For Ecol Manage. 2004, 199: 67-82. 10.1016/j.foreco.2004.04.015.View Article
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. 10.1016/j.biocon.2010.04.007.View Article
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.
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.View Article
Raudenbush B, Bryk AS: Hierarchical linear models: applications and data analysis methods. 2001, Thousand Oaks, California: Sage Publications, 2
Hosmer DW, Lemeshow S: Applied logistic regression. 2000, New York, New York: John Wiley & Sons, LtdView Article
Wilbur HM: Density-dependent aspects of growth and metamorphosis in Bufo americanus. Ecology. 1977, 58: 196-200. 10.2307/1935122.View Article
Skogland T: The effects of density dependent resource limitation on size of wild reindeer. Oecologia. 1983, 60: 156-168. 10.1007/BF00379517.View Article
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. 10.1098/rspb.2001.1585.PubMed CentralView ArticlePubMed
Benn B, Herrero S: Grizzly bear mortality and human access in Banff and Yoho National Parks, 1971–98. Ursus. 2002, 13: 213-221.
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. 10.1016/j.biocon.2004.02.020.View Article
Zedrosser A, Dahle B, Swenson JE, Gerstl N: Status and management of the brown bear in Europe. Ursus. 2001, 12: 9-20.
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. 10.1038/nature09210.View ArticlePubMed
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. 10.1111/j.1365-2656.2007.01276.x.View ArticlePubMed
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. 10.2192/1537-6176(2005)016[0047:DBADSS]2.0.CO;2.View Article
Coogan S, Nielsen SE, Stenhouse GB: Spatial and temporal heterogeneity creates a “brown tide” in root phenology and nutrition. ISRN Ecology. 2012, 2012: 10-View Article
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. 10.1111/j.1469-7998.1996.tb05449.x.View Article
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. 10.1111/j.1365-2656.2009.01647.x.View ArticlePubMed
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. 10.1111/j.1365-2486.2009.02061.x.View Article
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. 10.1139/B09-013.View Article
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. 10.1139/z04-077.View Article
Meiri S, Guy D, Dayan T, Simberloff D: Global change and carnivore body size: data are stasis. Global Ecol Biogeogr. 2009, 18: 240-247. 10.1111/j.1466-8238.2008.00437.x.View Article
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. 10.1139/z95-197.View Article
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. 10.1890/08-1036.1.View ArticlePubMed
Yom-Tov Y: Body sizes of carnivores commensal with humans have increased over the past 50 years. Func Ecol. 2003, 17: 323-327. 10.1046/j.1365-2435.2003.00735.x.View Article
This article is published under license to 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.