Abstract

For long-lived species, environmental factors experienced early in life can have lasting effects persisting into adulthood. Large herbivores can be susceptible to cohort-wide declines in fitness as a result of decreases in forage availability, because of extrinsic factors, including extreme climate or high population densities. To examine effects of cohort-specific extrinsic factors on size of adults, we performed a retrospective analysis on harvest data of 450 male black-tailed deer (Odocoileus hemionus columbianus) over 19 years in central California, USA. We determined that population density of females had a more dominant effect than did precipitation on body size of males. Harvest of female deer resulted in increases in the overall size of males, even though a 6-year drought occurred during that treatment period. Body size was most influenced by female population density early in life, while antler size was highly affected by both weather early in life and the year directly before harvest. This study provides insights that improve our understanding of the role of cohort effects in body and antler size by cervids; and, in particular, that reduction in female population density can have a profound effect on the body and antler size of male deer.

1. Introduction

Climate and population density influence organisms throughout their lives and often lead to phenotypic variation among cohorts. Cohort effects occur when individuals born within the same year experience similar environmental conditions, which may alter the life-history characteristics of that group [1]. Large mammalian herbivores, in particular, are strongly influenced by local weather, because of direct effects of weather on forage [24]. Rainfall increases the quality and quantity of vegetation [5], which has a positive influence on body size of some ungulates [68]. Body size is affected more easily by variation in the environment when a population is at moderate to high densities [711], because of increased intraspecific competition for food [12].

Understanding how the interaction between climate and population density influences body mass throughout the lifetime of a large herbivore enables biologists to better predict growth, survival, and reproduction within cohorts [1317]. Extrinsic and intrinsic factors affect body condition, and when resources are limited, important tradeoffs may be made early in life by an individual, or by a mother before parturition, which affects offspring while in utero [18]. Those tradeoffs result in delayed life-history characteristics that can persist into adulthood and even through future generations [4].

In times of severe weather or high population density, tradeoffs are made by the mother when the fetus is in utero [19]. Decreases in maternal investment during periods of low spring precipitation stunt the growth of fetuses, resulting in young that are born small and remain small throughout their lives [2026]. Furthermore, much of the variation in adult size of ungulates has been linked to conditions experienced during the year of birth [27] and environmental conditions and nutrition of their grandmothers 2 generations previously [3, 28].

Antler size has been positively correlated with warmer temperatures and higher precipitation in spring and autumn during the year of harvest [29]. Discerning the influence of extrinsic factors the year preceding harvest with those that occurred during youth is critical for understanding phenotypic variation in size of antlers. Determining the degree to which weather and population density during youth and adulthood affect antler and body size of ungulates is still poorly understood, although body and antler size are related to a number of life-history characteristics among cervids [30].

We conducted a retrospective analysis of harvest data [31] to examine effects of environmental conditions on antler and body size within and across cohorts of a population of black-tailed deer (Odocoileus hemionus columbianus) from 1980 to 1997. We hypothesized that cohort effects would influence body and antler size of male black-tailed deer. Specifically, we predicted that high annual precipitation during gestation, during the first year of growth, and during the year of harvest would be related to body and antler size. Temperatures that were low in summer and high in winter during those periods were predicted to positively affect body size of males. Also, we predicted that a decrease (brought about by harvest) in female population density during gestation, and 1 year following parturition, would result in greater body and antler size of males at the time of harvest. Understanding how a single year of environmental factors affects life-history characteristics of multiple generations can lead to improvements in predicting body growth and antler size, as well as providing insights into conserving and managing those large herbivores.

2. Methods

Research by McCullough [31] was conducted at the Hopland Research and Extension Center (hereafter Hopland) of the University of California (38°58′N, 123°08′W). The 2,149 ha area is located in Mendocino County, about 170 km north of San Francisco, in the central portion of the coastal mountain ranges of California, USA. Hopland is characterized by a Mediterranean climate; mean summer and autumn temperatures range from 20 to 25°C, and winters are mild with occasional frost in valley bottoms and infrequent light snow at high elevations (>800 m). Interannual rainfall was variable among years (1979–1997;  mm; range = 50–191 mm, CV = 37%), with a 6-year drought occurring during the treatment period (1985–1991), where females were harvested (Figure 1).

Discrete seasons of vegetative growth and senescence were evident between 1979 and 1997, which followed seasonal patterns of precipitation (, ). Most precipitation occurred from October through April; rain seldom fell during May through September, which also was characterized by the warmest average monthly temperatures (Figure 2). Winters are usually mild and wet, encouraging herbaceous growth that extends from November through April. The percentage of vegetative cover changes as wet conditions shift towards the dry season, with an increase in annual grasses and decrease in forbs [2]. Common vegetation types on Hopland include oak woodlands (Quercus spp.), grasslands (e.g., Bromus spp., Erodium spp., and Festuca spp.), and chaparral (Adenostoma spp.) [32].

Domestic sheep (Ovis aries) were pastured on Hopland and maintained at a population size of 900–1,500 females throughout the study (1979–1997) [33]. Black-tailed deer were abundant and hunted annually. The prehunt population before the study began was estimated at 443 deer (25.9 deer/km2) [31]. In 1991, at the end of the treatment period (female and male harvest), the population was approximately 250 deer (14.6 deer/km2). Predators of deer include bobcats (Lynx rufus), coyotes (Canis latrans), mountain lions (Puma concolor), and black bears (Ursus americanus). No attempt was made to control predators during the study.

2.1. Collection and Measurements of Deer

To determine if density-dependent effects influenced antler growth and body size of male black-tailed deer, we reanalyzed data collected by McCullough [31] from 1979 to 1997. No data were used from males sampled in 1979, however, because too few morphometric measurements were collected. During that study, there was a 7-year treatment period (1985–1991) where a mean of 24 females (SD = 4.6 females) were harvested each year in addition to the annual harvest of males. Six years before and after the treatment period of female harvest, only male deer were killed (~30 males per year). The initial research by McCullough [31] was conducted under animal use protocols approved by the University of California, Hopland Research and Extension Center, and a memorandum of understanding and collecting permits from the California Fish and Game Commission.

In our analysis, we used weight, antler measurements, and estimated age from 450 male deer that were harvested during 1980–1997. Males were not included if they lacked data for age or multiple measurements of antlers. When only a single antler measurement (diameter or length of either the right or left antlers) was missing, data were interpolated () from the complete data set. We used the regression function in the R 3.1.2 statistical environment [34, 35] to determine those values. Dressed weight was used in place of whole-body weight because whole weights were not recorded from 1991 to 1997. We used this approach because dressed and whole-body weights [36] of black-tailed deer exhibited a strong positive relationship (, , ). Moreover, carcass weight in mule deer (O. h. hemionus) is strongly correlated with other body measurements, including chest girth, body length, head length, and mandible length [37]. From 1980 to 1997, deer were aged from either cementum annuli or tooth eruption and wear [38]. Data on cementum annuli were used for aging deer during years where both types of data were available. We were able to use both types of data because age determined from cementum annuli and tooth eruption and wear was positively related (, , ). Aging black-tailed deer based on tooth eruption and wear is highly accurate through 2.5 years of age. Approximately 73% of deer 3.5 to ≥5.5 years of age were categorized correctly with a tendency to underestimate ages [39].

To reduce the number of variables, a Principal Component Analysis (PCA) based on the covariance matrix was conducted on antler measurements (right and left diameter and length, and spread) collected by McCullough [31]. Spread and length were measured to the nearest 6.4 mm, and diameters were measured to the nearest 0.5 mm. The number of left and right antler tines (points) was not included in the PCA, because of problems associated with the use of a small range of meristic characters to estimate phenotypic development [40]. The proportion of variance explained by the PCA was nearly identical between the whole raw data set (0.886) and with the addition of interpolated values (0.880); consequently, we included the interpolated values to increase sample size (). PC1 for antler measurements had negative loadings, with larger antlers represented by negative values and smaller antlers by positive values. The proportion of variance explained by PC1 was 0.76. Loadings for antler measurements were similar (Table 2), indicating that PC1 was a reliable measure of antler size [41, 42]. We multiplied the principal components by negative 1 for all analyses so positive values represented larger antlers, making results more logical to follow.

Age is related to antler and body size in male cervids [43, 44]. We fitted the Von Bertalanffy growth curve to examine the relationship between body size and age in male black-tailed deer using the package Fishmethods [45] in R [46, 47]. This procedure also provided the best fit for the relation between age and antler size (PC1), because we had few old individuals that might have resulted in a quadratic curve [42]. Residuals were obtained for both antler and body size, which corrected those size variables for age.

2.2. Local Weather

We determined effects of extrinsic factors on body weight and antler size of males by compiling local weather data. Total precipitation and average temperature were collected at Hopland each month. Monthly average temperatures and total monthly precipitation were compared in a Walter climate graph (Figure 2) to distinguish seasons of vegetative growth (October through April) and times of water deficit (May through September) [48].

Periods for which weather data were considered in models included gestation by the mother while the male was in utero (October through April); the first year of growth by the male offspring (May through April); and the weather during the year before harvest (October the previous year through September). For both the first year of growth and the year before harvest, we divided those periods into wet (October through April) and dry (May through September) seasons for a clearer understanding of how temperature and precipitation influenced deer. For instance, weather variables that would be included in data for a male deer harvested at 3.5 years of age in autumn of 1990 would include weather experienced while in utero during gestation by the mother (October 1986 through April 1987), the dry season of the first year of growth by the male (May 1987 through September 1987), the wet season of the first year of growth by the male (October 1987 through April 1988), the wet season of the year preceding harvest (October 1989 through April 1990), and the dry season of the year preceding harvest of the male (May 1990 through September 1990). Although the last third of gestation is the most energetically expensive for large mammalian females [49], weather during full gestation was used in our analyses, because female body condition affects the offspring throughout its time in utero. Furthermore, female deer inhabiting a Mediterranean ecosystem must build sufficient body reserves early in gestation, because parturition occurs at the beginning of the dry season when vegetation starts to senesce [2, 50].

2.3. Statistical Analyses

Intrinsic factors (body and antler size) were modeled as dependent variables in R [34] against extrinsic variables that included treatment period (pretreatment and combined treatment and posttreatment) and total monthly precipitation and average monthly temperature for each time period (i.e., in utero, the dry season during the first year of growth, the wet season during the first year of growth, the wet season before harvest, and the dry season before harvest). Treatment and posttreatment were combined into a single period because of lingering density-dependent factors after 7 years of female harvest. An interaction term between precipitation and population density was added to models. Body weight corrected for age was included as a variable in the antler-size model, because higher body reserves would allow the male to invest more in antler growth [41, 42]. Contrarily, antler size was omitted from the body-size model because antler growth is less likely to affect body size. Colinearity was examined for all variables with variance inflation factors (VIF) in R package Car [51]. All continuous quantitative variables had VIF values ≤5 indicating that multicolinearity did not bias our models [52]. We used the package MuMIn [53] in R to dredge and average the top models based on all models with values ≤2 using Akaike’s information criterion for small samples (AICc) of the top model [54]. Confidence intervals for model-averaged parameters were lowered to 85% to more closely conjoin the information theoretic approach of AICc to most frequentist confidence intervals. In particular, holding likelihood constant, AICc will increase by , where is the number of added parameters in a model. Thus, just as , where is a chi-square distribution with one degree of freedom, it is also true that [55]. We used multiple lines of inquiry to make the most informed decision about which variables were most influential in determining body weight or antler size. We examined 85% confidence intervals, and the relative importance of each variable, calculated as a sum of the Akaike weights over all of the models in which a particular parameter appeared. To further understand how a reduction in female population density affected antler and body size, we performed the median test [56]. We compared the median before and after female harvest began with the median for the entire study.

3. Results

We performed a retrospective analysis on data for 450 black-tailed deer, harvested between 1980 and 1997 at the Hopland Research Extension Center. Male black-tailed deer, on average, weighed 36.9 kg. Asymptotic body mass, determined from Von Bertalanffy growth equation, was 46 kg at 5.5 years of age (), ; Figure 3). To account for the relationship between age and body and antler size, residuals from the Von Bertalanffy growth curve were used in subsequent analyses (Figure 3). Dressed weight and antler size (PC1) also were positively related (Figure 4), but to better understand changes in size of male deer, we examined those variables separately.

Body size of male deer increased during the study (Figure 5(a), Table 1). Body size (i.e., dressed weight corrected for age) increased profoundly after the treatment period but was highly variable among individuals (, , ; Figure 5(a)).

Deer harvested after the treatment period began were significantly larger than those harvested during the pretreatment period. When corrected for age, 180 of 323 male deer harvested after the treatment period had a body weight greater than the median for the entire sample size (; , ). Through the study, growth rate and asymptotic body weight of male deer rose from pretreatment (0.908 per year; 38.7 kg) to posttreatment period (1.22 per year; 446.4 kg) (, ; , , resp.).

Antler size (PC1) increased slightly during the study (, , ; Figure 5(b); Table 1). Overall, antler size increased with age (Figure 6, , ), and when size was corrected for age, the difference in the number of deer with antlers larger than the grand median was significantly greater after female harvest began, compared with pretreatment (, ). When both high precipitation and low female density occurred simultaneously, there were substantial increases in size of antlers, such as those observed in males harvested in 1995 (Figures 1 and 5(b)).

For models of both antler size and body weight, extrinsic variables incorporated into each model included an index for population density (reflected by treatment period), and temperature and precipitation during critical periods of time for growth of body and antlers in male deer (Tables 2 and 3).

In our model for body size (dressed weight corrected for age), 16 competing models were averaged (Table 2). Two variables were determined to be the best predictors of body weight from 85% confidence intervals and summed Akaike weights (Table 4). The most significant variable was female population density during the gestation by the mother while the male was in utero, but only after female harvest began in 1985 (treatment and posttreatment periods combined). The interaction between weather, specifically precipitation, and treatment period when the male was in utero was important in predicting body size of male deer. Although 85% CI overlapped zero, population density interacted with precipitation during the first year of growth.

We identified and averaged 17 competing models explaining antler size (PC1 corrected for age) (Table 3). Five variables were influential in determining antler size (Table 5). Those variables had a relative importance of 1 when Akaike weights were summed and were different from 0 when we examined 85% confidence intervals. The variable that most positively influenced antler size was body weight corrected for age. Temperature was highly influential during three critical periods of time. All temperature variables exhibited a negative trend with antler size of male deer. The periods in which temperature was most important included, in descending order of significance, the dry season directly before harvest, during the wet season during the first year of growth, and during the wet season the year of harvest. A negative interaction occurred between precipitation and population density (indexed by treatment period) after 1985, which significantly affected antler size. Even though the 85% CI overlapped zero, female population density while the male was in utero had a significant positive trend in determining antler size, but only after female harvest began in 1985 (treatment and posttreatment periods combined).

4. Discussion

We demonstrated that cohort effects could have profound influences on body and antler size of male black-tailed deer, but climatic and density-dependent factors affected those variables differently. We predicted that a reduction in female population density while the male was in utero and during his first year of growth would lead to increases in size of body and antlers of male deer. Following our prediction, treatment period, an index to female population density, was a significant predictor of body size. Specifically, female population density after 1985 (after female harvest had begun—treatment and posttreatment combined) was positively related to body size of male deer at time of harvest. Female density during the first year of growth of the male was an important predictor of the body size of male deer at time of harvest. Notably, younger (~3 years of age) males harvested during posttreatment had larger antlers and obtain heavier body weight than older males (8–12 year olds), which were reared during a time of high female population densities (Figures 3 and 5(a)). Intraspecific competition between females and young likely was high early in the study, reducing the energy available for females to invest in offspring, probably stunting the growth of males early in life. Other studies also have documented that cervids born small remain small as adults, even when conditions are favorable later in life [3].

We predicted that high annual precipitation during critical periods for growth (while in utero, during the first year of growth, and during the year of harvest) would positively affect body and antler size of males. Contrary to our prediction, precipitation by itself had little influence on either body or antler size. Only for the antler size model did precipitation, specifically during the first summer of life, have strong predictive power.

Our results, however, did affirm that there was a significant interaction between precipitation and treatment period, revealed by the 6-year drought during the study. Indeed, that interaction between female population density and precipitation had a significant effect on both antler and body size of males. In the body-size model, there was a negative relation between precipitation during the gestation by the mother and female density after 1985 (after female harvest had begun and female population density was lower). Similarly, that interaction was significant in the model for antler size, but this time involving precipitation during the first year of growth by the male. A reasonable explanation for this result is the increased rainfall during 1994, which also corresponded with some of the lowest female densities following 7 years of female harvest. High precipitation alone did not result in large antler and body sizes documented during 1994 and 1995, because there was even more precipitation in 1982 and 1983, yet deer were significantly smaller during those years. Our models offer evidence that the positive effects of female harvest overrode negative influences of the drought and resulted in greater body weight and antler size, a conclusion also reached by McCullough [31].

As expected, temperature had a negative relation with size of male deer during three critical periods of time—while in utero, during the first year of growth, and the year immediately before harvest. Increases in temperature during the wet season likely reduced snowpack at high elevations and hereby constrained late season forage by decreasing water availability in late spring from snowmelt. Late spring is a critical time for deer in a Mediterranean ecosystem, because females with young must lactate as forage senesces [50]. Unexpectedly, temperature variables only were strong predictors of antler size but not body weight. Male deer grow antlers from March through July [57], which coincides with some of the hottest temperatures, particularly during June and July (Figure 2). The negative relation ostensibly occurred because of effects of hot temperatures on the senescence of forage when antlers were still growing.

We faced several challenges in conducting this study. We were constrained in some analyses because of issues concerning missing data. First, we attempted to examine effects of the normalized difference vegetation index (NDVI) on body and antler growth of deer. Unfortunately, too few images were available because of cloud cover or missing scenes to make that assessment possible across years. Nevertheless, we observed positive correlations (data not shown) between weather variables and available NDVI scenes. During some years, there also was incomplete data on antler measurements of deer [31]. We dealt with this shortcoming by interpolating antler metrics for no more than one measurement for 49 of 450 deer, with few observable effects on results. Data from 1979 were too sparse and therefore not included in analyses. Moreover, data were not collected annually on population size, so we used an index for female population density by dividing the study into two categories—pretreatment compared with the combined treatment (female harvest) and posttreatment periods. Indeed, population density with respect to carrying capacity () is reflected by female condition and reproduction [58]. We have no data on how selectivity of hunters for large males may have changed over the course of the study, but there is no obvious reason to believe it did.

Although female population size was not recorded during the study, three density-dependent patterns allowed us to infer that female density decreased sufficiently to allow for improvements in physical condition during the treatment period. First, females <2 years old were mistaken for older females because of their large body size [31], indicating that population densities were below and females were in good nutritional condition [17, 59, 60]. Pregnancy of females <2 years old also increased markedly during the treatment period [31]. Cementum annuli, induced by periods of annual stress, were more obscure after female harvest began [61]. McCullough [31] noted that those results indicated a higher-quality diet, likely leading to improvements in body condition of females, which increased recruitment of young as well as the size of male offspring.

Common conditions experienced by a cohort can result in distinguishable differences in phenotype when compared with other cohorts of the same age [31, 62, 63]. Therefore, identifying the long-term effects of environmental factors during early stages of life is critical to understanding patterns of growth and reproduction in populations of large herbivores [10, 60]. Unfortunately, the interaction between weather and population density is not always well understood [58]. Even on an optimal diet, offspring originally from an area with high population density and low resource availability exhibited suppressed growth for 1 generation in domestic sheep (Ovis aries) [64] and did not recover until 2 generations in white-tailed deer (Odocoileus virginianus) [3]. In white-tailed deer, survival was directly related to the body condition of the mother during the previous two winters [65]. In contrast, other studies have reported that deer were able to recover from periods of poor nutrition and overcome cohort effects [6668]. Consequently, understanding how a cohort can be affected by weather or population density has immense implications for the management and conservation of ungulates.

We were able to construct models for predicting body and antler size in male black-tailed deer living in a Mediterranean ecosystem. Our research offers new insights into the interaction between extrinsic and intrinsic factors, and effects of maternal and postnatal conditions on overall size of males at time of harvest. Biologists can evaluate effects of harvest on trends of antler size in relation to extrinsic factors during gestation by the mother and the first year of growth of male deer. Even minor changes in quality and availability of forage have resulted in relatively large fluctuations in rate of parturition and recruitment of young, because of the high energetic requirements necessary for gestation and lactation [60, 6972]. In cervids, a cohort-wide decrease in body mass could have lasting consequences, affecting the reproductive performance of a population, because of the close association between size of males and reproductive success [8, 7375]. We have created a valuable model for understanding and examining effects of population density in a variable climate. Moreover, our study provides additional information about interactions between extrinsic and intrinsic factors in large mammalian herbivores and its effects on antler size [25, 26, 29, 44, 67, 76, 77]. Our research is crucial for understanding the role of environmental factors in studies on body size and size of horn-like structures in relation to harvest of males [78, 79]. In addition to considering effects of harvest on the size of horn-like structures, our research demonstrates that density-dependent influences and cohort effects need to be taken into account. Indeed, there is increasing evidence that density-dependent mechanisms may override beneficial and detrimental effects of climate [11, 31, 58].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank many individuals that helped collect data during the initial study and were already acknowledged in McCullough (2001). For this study, the authors thank Z. Simpson for data processing and D. M. Walcott and K. M. Stewart for helpful reviews of their paper.

Supplementary Materials

Supplementary Table 1: Summary statistics on variables measured on harvested black-tailed deer during 1980–1997 from Hopland, California, USA.

Supplementary Table 2: Principal Component Analysis for antler size of male black-tailed deer harvested between 1980 and 1997 in Hopland, California, USA.

Supplementary Table 3: Mean dressed weight of male deer per age group for three treatment periods from 1980 to 1997.

  1. Supplementary Tables