Skip to main content

Naivety dies with the calf: calf loss to human hunters imposes behavioral change in a long-lived but heavily harvested ungulate

Abstract

Background

In prey, patterns of individual habitat selection and movement can be a consequence of an individuals’ anti-predator behavior. Adjustments of anti-predator behavior are important for prey to increase their survival. Hunters may alter the anti-predator behavior of prey. In long-lived animals, experience may cause behavioral changes during individuals’ lifetime, which may result in altered habitat selection and movement. Our knowledge of which specific events related to hunting activity induce behavioral changes in solitary living species is still limited.

Methods

We used offspring loss in a solitary and long-lived ungulate species, moose (Alces alces), as our model system. We investigated whether offspring loss to hunters induces behavioral changes in a species subjected to heavy human harvest but free from natural predation. To test for behavioral change in relation to two proxies for experience (calf fate and age), we combined movement data from 51 adult female moose with data on their offspring survival and female age. We tested for adjustments in females’ habitat selection and movement following calf harvest using Hidden Markov Models and integrated Step Selection Analysis to obtain behavioral state specific habitat selection coefficients.

Results

We found that females with a harvested calf modified habitat selection and movement during the following hunting season. Female moose selected for shorter distance to roads during the night, selected for shorter distance to forests and greater distance to human settlements following calf harvest than females who had not lost a calf. The survival of twins in a given hunting season was related to female age. Older females we more likely to have twins survive the hunting season.

Conclusions

Our findings suggest that losing offspring to human harvest imposes behavioral changes in a long-lived ungulate species, leading to adjustments in females' habitat selection and movement behavior, which may lower the risk of encountering hunters. In our study, female moose that experienced calf loss selected for lower distance to forest and selected for greater distance to human settlements during periods of high hunting pressure compared to females without the experience of calf loss during the previous hunting season. We interpret this as potential learning effects.

Introduction

The behavioral traits of animals are innate or learned through experience [1, 2]. In prey species, anti-predator behavior is vital for increasing both individual and offspring survival, thereby contributing to individuals’ fitness [3]. Efficient anti-predator behavior demands flexibility and development that allows prey to adjust to their major source of mortality. This suggests that an individual’s experiences throughout its lifetime shape anti-predator behavior [4,5,6]. The expression of anti-predator behaviors can differ between individuals across species’ geographical ranges according to predator occurrence, thereby affecting individuals’ likelihood of survival under the risk of predation [7], as well as the species’ evolutionary success [3, 8].

Generally, anti-predator behavior addresses the predation risk of the most abundant predator or the predator causing the most fatalities at a given time and place, including human hunters [4, 9]. Animals across different taxa have been found to adjust their behavior in accordance to varying degrees of predation risk, which has been interpreted as potentially learning from predation attempts [10, 11]. In systems with intensive human harvest, and where humans act as the most dominant predator, the impact of hunters can overshadow the effect of natural predators on prey response [4, 12, 13].

In general, different species have shown indications of learning processes or memory formation (e.g., roe deer (Capreolus capreolus), [14, 15]), as well as in terms of adjusted anti-predator behavior (e.g., in bettongs (Bettongia lesueur) [16] or great tits (Parus major) [17]). Specifically, some behavioral changes suggest learning effects may be triggered by experiencing mortal events of a conspecific, group member, or offspring (e.g., in damselfish (Acanthochromis polyacanthus) [18], crows (Corvus brachyrhynchos) [19] or common goldeneyes (Bucephala clangula) [20]). In species subjected to human harvest or predation, individuals may adjust their behavior in relation to hunting activity [6, 21], as well as towards carnivore species (e.g., Brown bears (Ursus arctos) [22] or wolves (Canis lupus) [23]). Despite the known capability of individuals to adjust their behavior towards the mortality risk during the hunting season [6, 24, 25], our knowledge of the mechanisms underlying this change in anti-predator behavior is still limited, particularly for solitary-living species. Similarly, we are limited in measuring the behavioral responses involved.

Ungulates are among the most harvested species groups worldwide [26, 27]. The risk of human harvest invokes the need for prey to modify their behavioral responses towards human disturbances, especially during the hunting season [9, 28, 29]. For example, hunting of ungulates is usually only allowed during daytime and human hunters often utilize areas near roads [24], making their spatiotemporal distribution predictable. For prey species that are subject to human harvest, behavioral changes have to match the spatiotemporal risk for harvest (i.e. minimizing the risk of encountering hunters by hiding or moving less [12]). For example, human-habituated and gregarious reindeer (Rangifer tarandus) developed new behaviors once they started to be culled [30], and solitary-living roe deer (C. capreolus) modified their behavior in relation to hunting regimes [31, 32]. Furthermore, experiencing mortality events of conspecifics may cause individuals to alter their behavior, which in return may increase their own survival changes [6, 31, 33, 34]. In short, we still lack a comprehensive understanding of which events related to nonlethal experiences that may lead to behavioral adjustments in long-lived solitary species and so far mostly proxies for experience, such as age, have been used [6].

In this study, we investigated whether offspring loss to human harvest results in adjusted anti-predator behavior in terms of habitat selection and movement. As our model system, we used adult Swedish female moose in two populations, which are heavily harvested. In our study areas, hunting accounts for most of the mortality, as large carnivores are functionally extinct [35, 36]. Importantly, female moose live mostly solitary and within the Swedish hunting policy, calves must be harvested before their mothers, potentially generating several opportunities throughout the females’ approximate 20-year lifetime to test whether females adjust their behavior after experiencing calf harvest [37]. Furthermore, this policy framework generates a “life insurance” for females accompanied by calves, and potential double learning events for females with twins.

In ungulates, movement and habitat selection are major factors that can define an individuals’ fate and are generally closely linked to the surrounding environment [38, 39]. Further, the spatial behavior and movement of an individual may depend on the presence of offspring [40,41,42]. Recent research has highlighted that considering behavior-dependent habitat selection (e.g., is the animal foraging or explorative) improves the estimation and prediction of animal habitat preferences and advances our understanding of ecology [43, 44]. To test for individual adjustments in anti-predator behavior, we analyzed movement- and behavior-specific habitat selection in adult female moose following the harvest of their calves. Specifically, we investigated differences in females’ behavioral-state-specific habitat selection and movement in relation to two proxies for experience, as we expected change of anti-predator behavior in response to harvest risk and female experience to be dependent on the behavioral state an animal is in [6, 25].

We tested our hypothesis that female moose that lose a calf to human harvest will adjust their state-specific habitat selection and movement in the following hunting season (in an effort to evade human hunters). Hence, our null hypothesis is that no behavioral adjustment occurs following calf loss to human harvest.

First, we predict (a) that female moose exhibit non-naive anti-predator behavior (i.e., anti-predator behavior is adjusted to evade human hunters) after losing a calf (a specific experience) during the previous hunting season. We expect female to select habitats farther away from human infrastructure (i.e. roads), stay more sheltered and reduce movement following calf loss. Second, we predict (b) that changes in state-specific behavior in relation to age (as a proxy for and measure of experience) will show similar trends as they do for calf fate. Finally, we predict (c) that older females have a greater chance of having their calves survive the hunting season as they are more experienced and have adapted non-naive anti-predator behaviors.

Materials and methods

Study area

Our study combined data from two sites in the southern boreal region (56° N, 14° E and 58° N, 17° E) (Fig. 1). At both sites, mixed forests as patches of deciduous forests (e.g., birch (Betula sp.), aspen (Populous sp.), elm (Ulmus glabra), oak (Quercus robur), maple (Acer platanoides)) and coniferous forest such as Scots pine (Pinus sylvestris) and Norway spruce (Picea abies) intermix with agricultural fields in a flat to gently rolling terrain. In our study areas, the annual moose hunting season starts after the moose rut on the second Monday in October and ends in February [45, 46]. Moose hunting is usually carried out in a hunting team and with the use of baying dogs [45, 47]. Hunting pressure is generally highest during the first 3 weeks of the hunting season and decreases thereafter. The majority of the moose harvest occurs from the second Monday of October until the end of December [48].

Fig. 1
figure 1

The two study areas (orange) as the minimum convex polygon around all GPS positions in southern Sweden. Fennoscandia is marked in gray, and Sweden is marked in blue

GPS data

We analyzed 10 years of Global Positioning System (GPS) positions of 51 female moose individuals from 2008 to 2018. Moose were tranquilized from a helicopter with a dart-gun that injected a mixture of etorphine-acepromazine and xylazine in winter (January/February) [49]. Moose were equipped with GPS Plus Collars following our standard protocol [50, 51]. The GPS-collars were programmed to a 3-h sampling rate. We estimated the individual age from tooth wear [37, 50]. To cover major parts of the hunting season and the period before the annual moose hunts start, we analyzed moose GPS data from the 1st of August to the end of December (when hunting pressure decreases significantly), leaving in total 5 months of GPS data to analyze in each year.

Calf survival and reproduction data

We assessed the number of calves born and their summer survival by field observations during two periods, shortly after birth and before the hunting season, following our standard protocol [52]. A calf was defined as lost if the calf was absent from the mother on two consecutive observations before and after the hunting season in a given reproductive year. We verified calf survival during the hunting season with field observations at the end of the hunting season and complemented them by hunters reporting if they shot either the calf of the collared female. Therefore, we considered any calf loss during the hunting season as a calf loss to hunters. To analyze the potential behavioral differences in female moose during the hunting season after calf loss, we matched the GPS data in a given year with information on calf loss in the previous hunting season. We omitted GPS data of moose individuals from the analysis in years when we could not verify during any visual calf check whether she had a calf at heel. We omitted this data to avoid inducing further noise in the analysis, as behavioral patterns of female moose with and without offspring can differ [40,41,42]. We did not consider the amount of calves an individual had lost over its collaring period, as the age at capture strongly varied between individuals and older individuals could have lost many calves before being fitted with a GPS collar, which would introduce bias and uncertainty. Given the lack of information on lost calves a given female may have experienced during her lifetime, we decided to use a more conservative approach. To analyze behavioral changes in relation to experience (calf fate in the previous hunting season), we chose to compare behavior between calf fate groups instead of using the accumulated number of calves lost. This approach enabled us to account for differences in recent experience (as measured by calf fate) in the previous hunting season and compare the difference in behavior among females. In total, we analyzed GPS data from 128 individual moose years of 51 female moose individuals; the average tracking time was 2.51 ± 1.57 years.

Statistical analysis

Segmentation of movement behavior

Habitat selection of ungulates can be behavior-specific (i.e., animals select different habitats based on their behavioral state [6, 43, 53]), making the incorporation of behavioral states into our analysis a vital departure point for understanding animal responses to the risk of predation [43]. To obtain movement-specific habitat selection coefficients and test for individual-based behavioral change following calf harvest, we first fitted a two-state Hidden Markov Model (HMM). We performed a segmentation of our movement data in order to avoid three-way and four-way interactions between movement parameters, habitat, hunting pressure and calf fate or age in later sections of the analysis.

We fitted the two-state HMM with two data streams using the momentuHMM—package [54]. We segmented the movement trajectories of each female moose into two unobserved behavioral states (restricted and exploratory), which we expected to differ in the degree an individual exposes itself and thus may affect its encounter risk with hunters [25, 55, 56] (see Fig. 2a, b). The unobserved states were estimated from step length and turning angles. We considered the restricted behavioral state reflect lower movement rates and thus generally more cautious behavior in terms of habitat selection. Further, we assumed the restricted state to reflect foraging behavior and a low chance of encountering human hunters. In contrast, we expected the exploratory behavioral state to reflect higher movement rates with high directional persistence and transitioning behavior, potentially with higher risk behavior due to higher probability of encountering human hunters [6, 33].

Fig. 2
figure 2

ac Our methodological approach to analyzing the state-specific habitat selection of female moose. a The GPS track of an animal, b the segmentation into behavioral states according to the Hidden Markov Model (HMM) (see “Segmentation of movement behavior” section), where yellow points and lines indicate an exploratory behavioral state and blue points and lines indicate a restricted behavioral state. c The fitted integrated Step Selection Function (iSSF) (see “Analysis of habitat selection” section), where dashed lines and hollow points indicate the alternative steps available to the animal at a given observed step. For visualization purposes, we only show one of 25 available random steps in (c)

Initial step length parameters were chosen according to Neumann et al. [57] to account for bimodal movement patterns of moose over a 3-h time interval. The restricted state was defined with an initial mean step length of 100 m (± 100 m SD and 0.001 zero mass parameter (to account for a low probability of no movement)), a mean turning angle of π, and a turning angle concentration of 0.001. Next, we defined the exploratory behavioral state with longer step lengths and high directional persistence. The exploratory state had an initial mean step length of 450 m (± 450 m SD, 0.001 zero mass parameter), a mean turning angle of 0.1, and a turning angle concentration of 0.99. We modeled step lengths using a gamma distribution and turning angles using a Von Mises distribution. We had to fit a zero mass parameter, as we had few steps in our dataset where animals did not move between GPS—locations [54]. We only defined two states and did not test for evidence of further states (e.g., encamped behavioral state) since we considered the spatial resolution (25 × 25 m) of our environmental covariates too coarse (see “Analysis of habitat selection” section) to model a behavioral state with extremely low movement rates. We closely followed procedures of previous applications of HMMs within ecological research [25] to fit the two-state HMM in our analyses. We then applied the Viterbi algorithm to our dataset to assign each step the most likely behavioral state based on the HMM [54, 56].

Analysis of habitat selection

To analyze state-specific habitat selection and movement of female moose (see “Segmentation of movement behavior” section) we fitted integrated Step Selection Functions (iSSF) [25, 58]. We split the data based on year and behavioral state and fitted separate iSSFs to analyze state-specific habitat selection on an annual basis. For each observed step, we generated 25 alternative steps [58,59,60]. We fitted iSSFs using the amt—package [61], using a gamma distribution to fit step lengths and a Von Mises distribution to fit turning angles to generate alternative steps [58] (Fig. 2c). We added one meter to all step lengths to account for all steps with no movement, as this would have led to errors when fitting the gamma distribution.

The behavior of big game (i.e., ungulates) hunters makes the spatiotemporal harvest risk for prey predictable in relation to infrastructure (e.g., hunting close to roads [24, 62, 63]) or cover. Based on our experience of spatial behavior of our study species, knowledge of behavior of Swedish hunters and previous literature on anti-predator behavior towards human hunters, we matched the used and available steps of the iSSF with three environmental covariates. We chose distance to roads, human settlements and forests as they play a crucial role in the spatial interactions of human hunters and moose [6, 24, 33, 45]. We calculated distance to forest using the Swedish Land Survey map from a binary raster with a 25 × 25 m spatial resolution [64], updated with data on clear cuts from the Swedish Forest Agency [65]. Distance to human settlements was calculated from the Swedish Land Survey map from a binary raster [64]. Distance to roads was calculated from a binary raster from the Swedish Transport Agency [66]. We expected selection for greater distance to forest to represent riskier behavior. For distance to human settlements and distance to roads it was the opposite, here we considered selecting for shorter distance to roads or human settlements to reflect riskier behavior. Covariates were extracted at the end of each step. Next, we added information on harvest risk as indexed by different hunting pressures (categorical, three factors: no hunting pressure (1st of Aug to the first Sunday of October, the day before the hunting season starts); high hunting pressure (2nd Monday of October (first day of the hunting season) and the following 3 weeks when hunting pressure is highest [48]) and low hunting pressure (after high pressure ends until the end of December)). Further, we added time of day (categorical, two factors: day and night) to account for diurnal patterns in habitat selection and no hunting during the night (Table 1). Time of day was calculated using the time_of_day()—function in the amt—package [61]. All environmental covariates were z-score transformed at the individual level [67].

Table 1 Overview of the covariates used in this study

We analyzed the behavioral state-specific iSSFs by fitting individual conditional logistic regressions for each animal in each year, using the matched sets of used and available steps as the response variable [58, 68]. Conditional logistic regressions were fit using the survival—package [69]. We included one-way interactions between both factors (time of day, hunting pressure) and environmental covariates in all models. We further fitted step length, the natural logarithm of step length and the cosine of the turning angles as covariates to account for the underlying movement process in the iSSF. As we fitted a model for every individual in each year and behavioral state, we obtained 256 models. The β coefficients estimated in models for the restricted behavioral state would therefore reflect habitat selection and movement with a lower probability of human hunter encounters but also more cautious behavior, while coefficients of the exploratory behavioral state would reflect a higher probability of hunter encounters and generally less cautious behavior.

Testing for change in habitat selection and movement

We applied a two-step approach to summarize coefficients of the individual models and draw inference at the population level. The two-step approach allowed us to analyze differences in β coefficients of individual models in relation to age and calf fate in a second modeling step to account for among- and within-individual variation [60, 70]. We chose to account for individual variation using the two-step approach instead of a mixed model approach, as the structure of the mixed model would have been very complex and we were interested in testing for differences in the individuals’ behavior [71]. Further, the two-step approach enabled us to avoid complex models with three-way to four-way interactions between variables that are difficult to interpret (e.g., distance to forest × hunting pressure × calf fate × age). The mixed effect model also likely would have led to convergence issues and would have complicated drawing inference about the same parameter sets and interactions for both behavioral states and both calf fate groups.

To ensure that we could analyze calf fate and female age within the same model and we could draw inference from both variables, we tested whether calf fate was not correlated to moose age using the Mann‒Whitney U test. We summarized β coefficients and tested for effects of calf fate (calf loss vs no calf loss in the previous hunting season) and age on habitat selection at the population level by using inverse-variance weighted (IVW) regression [32, 72]. The IVW regression gives less weight to parameters with higher variance, reducing uncertainty when drawing population level inference. We used the β coefficient sets (βdistance to forest, βdistance to forest × night, etc.) obtained from the conditional logistic regressions in each behavioral state as the response variables and fitted and ranked several competing sets of IVW models using Akaike Information Criterion corrected for small sample sizes (AICc). We averaged models with a Δ—AICc < 2 [73]. We included calf fate and age, as well as a two-way interaction between age and calf fate, as predictors in the IVW regressions to test for behavioral changes (see Table 2 for the ecological interpretations between variables and interactions). This resulted in 24 different sets of competing candidate models to analyze behavioral differences, one for each β coefficient set in each behavioral state.

Table 2 Overview of the fitted models used to test for differences in β coefficients and estimated-observed mean movement from the iSSF

Similarly, we tested for effects of calf fate and age on movement by calculating the estimated step length [74] using the formula \(l_{mean} = \frac{k + \beta \ln (l)}{{\theta^{ - 1} - \beta_{l} }}\), where k and θ determine the scale and shape of the observed gamma distribution, respectively; β and βln(l) are the estimated coefficients for the observed step length and natural logarithm of step length, respectively, of the conditional logistic regression [58]. We then subtracted the mean step length of each individual in each year and behavioral state. Like with habitat selection, we used AICc model selection to select the most parsimonious linear model describing the effects of age and calf loss on the Δ—estimated mean step length—observed mean step length in each behavioral state (Table 2, Appendix. Tables 7, 8).

To investigate whether experience of females increases calf survival in the ongoing hunting season (i.e. prediction (c)), we applied a generalized linear model with a binomial family using moose age and the number of calves before the hunting season as predictors and calf survival in the hunting season as the response (binary). Due to relatively small sample size for a given female (2.89 ± 1.97 SD years), we decided against a mixed modeling approach as we would not have a sufficient sample size for each individual, which would introduce bias [75, 76]. We included a two-way interaction between age and the number of calves. As we had no limitations regarding known calf loss status in the previous season and available data for the current season, the reproduction dataset was larger than the habitat selection dataset and included reproduction and survival data for 76 individual females in 218 moose years.

All the statistical analyses and data handling were conducted in R version 4.1.3 [77].

Results

Out of our 51 females, 13 never experienced calf loss, whereas 14 experienced calf loss every year. The remaining 24 females had both experiences, i.e., years in which they lost a calf and years in which they did not lose a calf to human hunters. The Mann‒Whitney U test showed that calf loss to human hunters occurred at all ages (W = 7815.5, p = 0.833) (see Fig. 1 in Appendix).

Model selection of calf fate and age

Including calf fate as a predictor of behavioral change in females’ habitat selection increased the model performance by 70.8% in all models. Calf fate had a greater impact on habitat selection in a restricted behavioral state than in an exploratory behavioral state (i.e., calf fate as a predictor of behavioral change improved 83.3% of the models for habitat selection in the restricted behavioral state, compared to 58.3% of the models for habitat selection in the exploratory behavioral state) (Table 3). In contrast, adding female age—as an additional proxy for experience—improved the overall model fit for 58.3% of the models. We found a greater impact of age on habitat selection in the exploratory behavioral state (50.0% of the models for habitat selection in restricted behavioral and 66.7% of the models for habitat selection in exploratory behavioral state; however, see Appendix Tables 16).

Table 3 Results of the IVW regression analysis of habitat selection of female moose (n = 51) from 2008 to 2018 in an exploratory and restricted behavioral state

Impact of calf fate on state-specific habitat selection

Habitat selection in the restricted behavioral state

Regardless of calf fate, females selected for shorter distance to forests when in a restricted behavioral state (Table 3). In addition, females that had lost a calf selected for shorter distance to forests compared to females that had not lost a calf during low hunting pressure (see also Fig. 2 in the Appendix for visual representation). During high hunting pressure, (the first 3 weeks of the moose hunting season), females that had lost a calf selected for greater distance to human settlements than females that had not lost a calf. Females in both calf fate groups selected for larger distance to roads before the hunting season, with females that had lost a calf selecting to be further away than females with surviving calves. Females who had experienced calf loss selected for shorter distance to roads during the night compared to those that had not lost a calf.

Habitat selection in the exploratory behavioral state

Before the hunting season, in the exploratory behavioral state, females selected for shorter distance to human settlements after a calf loss compared to females that had not lost a calf. This response changed with the onset of the hunting season. During low hunting pressure, females of both calf fate groups selected for greater distance to human settlements. When hunting pressure was low, we found no effect of calf loss on selection for distance to roads, regardless of behavioral state (β coefficient set distance to roads × low pressure in Table 3).

Impact of age on state-specific habitat selection

Habitat selection in the restricted behavioral state

In a restricted behavioral state, we found that age had a general effect on habitat selection in female moose, regardless of the fate of their calf in the previous year. Older females selected for shorter distances to forests. Older females selected for greater distances to human settlements during both day and night. We found only weak adjustments of selection for distance to roads during both high and low hunting pressure.

Habitat selection in the exploratory behavioral state

In an exploratory behavioral state, before the hunting season and during low hunting pressure, females selected for shorter distance to forest as they became older. This pattern reversed under high hunting pressure but was modified by time of day. At night, we found an increase in selection of greater distances to forests for older females in an exploratory behavioral state, compared to daytime. Older females selected to be closer to human settlements when in an exploratory behavioral state and during times of high and low hunting pressure. Age affected the selection of habitats regarding their distance to roads, with older females selecting shorter distances to roads compared to younger females. During times of low hunting pressure, this pattern reversed and older females selected for greater distance to roads. Selection for distance to roads at night was not affected by age, only by calf fate, regardless of behavioral state.

Impact of calf fate and age on state-specific movement

According to the HMM, female moose showed a mean step length of 68.3 m (± 59.7 m SD) in the restricted state and a mean step length of 214 m (± 256 m SD) in the exploratory state (see Fig. 3 in Appendix for state-specific distributions). Age influenced female movement. The most parsimonious model for state-specific movement for both behavioral states included a two-way interaction between calf fate and age (see Appendix Tables 7 and 8). The Δ estimated—observed step length, regardless of calf fate, was negative in a restricted behavioral state, meaning that female moose moved more than expected once movement was accounted for habitat selection. Older females who had lost a calf moved more than expected compared to females who had not lost a calf (Fig. 3). In an exploratory behavioral state, as they aged, females moved less than expected, regardless of calf fate, once movement was accounted for habitat selection.

Fig. 3
figure 3

Effects of calf loss and aging on female (n = 51 females) movement on Δ estimated—observed step length estimated by the most parsimonious linear model for each behavioral state (restricted/exploratory). Red and gray shaded areas indicate estimated standard errors

In-season calf survival

Age alone had no effect on calf survival in the ongoing hunting season (− 0.062 ± 0.049 SE, p = 0.209), indicating that older females were not more successful at avoiding calf harvest during a given season compared to younger females. The number of calves a female had at the onset of the hunting season had a strong negative relationship with offspring survival (− 2.502 ± 0.860 SE, p = 0.004), suggesting that females with twins had a greater risk of losing a calf. Age had a positive effect on the survival of multiple calves (0.232 ± 0.088 SE, p = 0.008), meaning that older females were more successful at keeping twins alive through an ongoing hunting season (Fig. 4).

Fig. 4
figure 4

Marginal effects of age and number of calves on calf survival in female moose (n = 76 females in the ongoing hunting season from 2008 to 2018). The red band shows the interaction between age and twins, and the gray band shows the interaction between age and a single calf. Shaded areas reflect the 95% confidence intervals

Discussion

In this study, we assessed behavioral changes in a solitary ungulate species in relation to two measures that may indicate experience in reproducing females: calf fate and female age. We derive two overarching conclusions from this study.

Change in behavior in relation to experience

First, we found that female moose that had lost their calf during the previous hunting season altered their selection for distance to roads and distance to human settlements during the following hunting season. This suggests that the loss of a calf due to human harvest might be a mechanism that induces behavioral change in female moose, supporting our prediction (a) that females will develop non-naive anti-predator behaviors following calf loss. Most noticeably, behavioral adjustments after calf loss were most prevalent during the hunting season and in the restricted behavioral state. Further, we found that previous calf fate often explained the observed behavioral differences in habitat selection and movement better than age (a common proxy for experience). Here, we found that females who had lost a calf in the previous season selected for greater distance to human settlements and lower distance to forests during the period of low hunting pressure. During the day, we found that female moose selected for shorter distance to forest, suggesting they select for more shelter during the hunting season. Likewise, in years following calf loss, they selected for shorter distance to roads during the night, when hunters in Sweden are not allowed to hunt, suggesting a spatiotemporal avoidance of hunters. Similarly, we found that females that had experienced calf loss strongly selected for greater distance to human settlements during the most intense phase of the hunting season. Lastly, we found that older female moose reduced their movement in the restricted behavioral state following calf loss in the previous hunting season, suggesting an adjustment of anti-predator behavior the following year by moving less and thus possibly reducing exposure to human hunters [6, 33]. We interpret these behavioral changes as a possible indication for learning in this solitary ungulate species. However, we cannot rule out that we missed additional behavioral changes in our analysis that could occur in relation to other habitat parameters, such as forage abundance or visibility [78, 79]. We also acknowledge that we considered a proxy for hunting pressure, not fine-scaled hunting pressure statistics such as number of hunters being present in a given area at a given time in relation to moose GPS—locations. Using this proxy likely added additional variation in our analyses, as females likely may not have been exposed to hunting pressure at any given time and place.

Adjustments of anti-predator behavior in relation to behavioral states

In this study, behavioral change following calf loss was more prevalent in the restricted behavioral state than in the exploratory behavioral state. This finding indicates that when utilizing a limited spatial areal (i.e. the restricted state is defined by shorter steps with low directional persistence, see Fig. 2), female moose select for environmental features that may help them to evade encounters with hunters following harvest of their calf in the previous hunting season. In contrast, in the exploratory behavioral state, we found females to show less behavioral adjustments following calf harvest. Per our definition, this behavioral state includes steps that are longer and more directed. Considering that long scale movements might be necessary in order to avoid human hunters spatially, there might be less room for adjustments of habitat selection or movement in this behavioral state. Additionally, animal responses to mortality risk can be rather short-termed and tortuous, and escapes from human hunters often happen within few minutes [80, 81]. Within this context, it is also important to acknowledge that our analysis was based on 3-h sampling rate of GPS—positions, which might have been too coarse to reflect potential fine-scale adjustments of anti-predator behavior while in the exploratory behavioral state. Alternatively, our chosen environmental covariates and/or their spatial resolution may not have been appropriate to capture the behavioral adjustments in this state. We therefore recommend future studies on behavioral adjustments in relation to experience using more fine-scaled data on hunting pressure, environmental characteristics and at a tighter temporal resolution of GPS—positions.

Our results also suggest that for a female, age and thus likely more experience may result in changed behavioral patterns, supporting our prediction (b) of adjusted behavior over a lifetime. Older female moose selected for habitats providing shelter from human hunters during the hunting season, suggesting that older females counterbalanced the spatiotemporal behavior of human hunters to reduce their harvest risk (i.e., selection for lower distance to forest and selection for greater distance to human settlements), which differed from the selection we observed in younger females. Furthermore, we could link changes in selection for distance to roads, another crucial part of the spatial ecology of human hunters, and the risk of being shot [12, 62] to both calf fate and aging, supporting findings on anti-predator behavior adjustments towards harvest risk in previous research [6, 24]. This suggests an ability to adjust and match behavioral patterns to spatiotemporal variations in harvest risk in a solitary ungulate species like moose. In addition, our results suggest that anti-predator behavior became more refined over time [6, 22, 33]. We want to acknowledge that we could not always link adjusted anti-predator behavior to age or calf fate. This may be due to our study working on a relatively coarse temporal scale (between years and on 3-h time intervals between individual steps where we measured selection). Anti-predator behavior, however, may be expressed at much finer temporal scales (i.e., minutes or seconds) [32, 82]. Lastly, behavioral changes and possible learning processes are complex, emerging from several processes [2, 6, 15], and may occur at points in time not included in this study or may be unrelated to the specific calf loss event that this study addressed. Therefore, performing a similar study in an experimental setting, with detailed information on female behavior and previous experiences, calf fate, and environmental settings could reveal additional information about how anti-predator behavior changes over an individual’s lifetime. Such a setting would help to quantify the impacts of experience on behavior and the adaptive value of the observed behavioral adjustments and would likely remove further unaccounted variation from analyses. Further, it would allow investigating potential accumulation (or the lack thereof) of behavioral change, as these life history parameters are known within such a setting. Within our study setting, we did not know the exact number of calves a female had lost before collaring and our sample size for young females that had yet to reproduce was insufficient to analyze this.

Calf survival in relation to experience

We did not find evidence that general calf survival was greater in older females, suggesting that even though they may adjust their behavior, the alteration does not have a significant effect on the risk of losing an offspring during a given hunting season. However, older females had a greater chance of keeping twins alive throughout the hunting season, supporting our prediction (c) that age (as a proxy for experience) has an effect on calf survival in an ongoing hunting season. We interpret this as a possible learning effect, perhaps due to accumulated experience in older females, particularly as calves must be shot before their mother.

Environmental variation and degree of sociality influence the need for individual learning as well as opportunities for learning [83]. For individuals in non-gregarious prey species, the number of nonlethal predation events likely are more limited than in gregarious species where individuals can learn, to a greater degree, from the fate of conspecifics (e.g., elk [6], or Indian mynahs (Acridotheres tristis) [11]). However, in gregarious species, collection of data on individual offspring loss to human harvest and disentangling relationships between herd members proves to be difficult and require extraordinary amounts of work. In solitary species, this data collection is much simpler since the female–calf relationship is more apparent. Therefore, linking observed behavioral changes in a given female to the survival of her calf is more straightforward.

Previous research has shown that in long-lived species, accumulated experiences allow learning and thus adjustment of behavior over an individuals’ lifetime in relation to changing external conditions, e.g., sea turtles [84] or wandering albatrosses (Diomedea exulans [85]). Moose have a lifespan of up to 20 years and may annually reproduce between three and 14 years of age [35, 50]. Therefore, female moose have up to 10 possible learning opportunities (accounting for singletons only; likely more when considering twins), thus allowing for highly adjusted anti-predator behavior.

Limitations in measuring behavioral responses

Despite largely advancing our understanding of the behavioral ecology of free-ranging animals [86], it is important to note that GPS data do not allow researchers to measure all behavioral adjustments an animal undergoes [87]. For example, some changes in anti-predator behavior are most likely more subtle and expressed on a finer spatiotemporal scale, e.g., short-term tortuous escapes from hunting dogs or decreasing visibility on a local scale by hiding on a minute scale, than the 3-h time interval in this study [88]. Therefore, we may have missed behaviors expressed at finer temporal scales [88]. Moreover, individuals in a non-gregarious species, such as moose, exhibit a considerable degree of variation in movement behavior, likely further increasing variability in anti-predator behavior [89]. We might also have missed behavioral adjustments due to confounding variation between the onset of the hunting season and other ecological factors influencing movement or habitat selection, such as the rutting season. However, this bias applied to all animals and our focus was to investigate the behavioral adjustment at the onset of the hunting season and in relation to diurnal patterns. The use of baying dogs is common for hunting moose in Sweden [47], which may make it difficult even for more experienced individuals to fully undergo the harvest risk. Additionally, the use of hunting dogs and the breed of dogs used may vary among hunting teams and years, which may affect female responses due to different levels of hunting disturbances [80]. Moose may use hiding strategies in the restricted behavioral state that human hunters with baying dogs exploit when releasing the dogs in areas where moose are known to aggregate and track them [80, 88]. Yet, our analyses depart from preys’ perspective, not from the predators. From the preys’ perspective, they still likely aim to expose themselves as little as possible to predators as compared to the explorative behavioral state, which our findings agree on [6, 32, 33].

We look forward to future studies that include data with a higher temporal resolution alongside other data (e.g., acceleration data, heart rate [22, 49, 86, 88]) and perhaps exact habitat quantifiable from cameras on the animal [90] in combination with exact locations and timing of the harvest risk (such as hunters and dogs). Such data may help to detect both stronger and more subtle behavioral changes, and thus increase our understanding of adaptive anti-predator behavior.

Conclusions

Our study highlights that offspring loss to human hunters induces behavioral change in adult females. This behavioral change may reduce exposure to hunters during high hunting pressure. We therefore recommend that future studies include this information when investigating behavioral change in relation to experience, as it might prove useful to compare direct responses to long-term effects that calf loss to human hunters has on female ungulates.

Data availability

Data and code have been archived in the Open Science Framework (OSF) repository and are available via: https://doi.org/10.17605/OSF.IO/NCMZU.

Abbreviations

HMM:

Hidden Markov Model

iSSF:

Integrated Step Selection Function

GPS:

Global Positioning System

IVW Regression:

Inverse variance weighted Regression

References

  1. Sih A, Bell A, Johnson JC. Behavioral syndromes: an ecological and evolutionary overview. Trends Ecol Evol. 2004;19(7):372–8.

    Article  PubMed  Google Scholar 

  2. Dingemanse NJ, Kazem AJ, Reale D, Wright J. Behavioural reaction norms: animal personality meets individual plasticity. Trends Ecol Evol. 2010;25(2):81–9.

    Article  PubMed  Google Scholar 

  3. Sih A. The behavioral response race between predator and prey. Am Nat. 1984;123(1):143–50.

    Article  Google Scholar 

  4. Lima SL. Life in a multipredator environment—some considerations for antipredatory vigilance. Ann Zool Fenn. 1992;29(4):217–26.

    Google Scholar 

  5. Lima SL, Dill LM. Behavioral decisions made under the risk of predation: a review and prospectus. Can J Zool. 1990;68(4):619–40.

    Article  Google Scholar 

  6. Thurfjell H, Ciuti S, Boyce MS. Learning from the mistakes of others: How female elk (Cervus elaphus) adjust behaviour with age to avoid hunters. PLoS ONE. 2017;12(6):e0178082.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Sand H, Wikenros C, Wabakken P, Liberg O. Cross-continental differences in patterns of predation: Will naive moose in Scandinavia ever learn? Proc Biol Sci. 2006;273(1592):1421–7.

    PubMed  PubMed Central  Google Scholar 

  8. Berger J. Carnivore repatriation and holarctic prey: narrowing the deficit in ecological effectiveness. Conserv Biol. 2007;21(4):1105–16.

    Article  PubMed  Google Scholar 

  9. Proffitt KM, Grigg JL, Hamlin KL, Garrott RA. Contrasting effects of wolves and human hunters on elk behavioral responses to predation risk. J Wildl Manag. 2009;73(3):345–56.

    Article  Google Scholar 

  10. Carthey AJR, Banks PB. Naiveté is not forever: responses of a vulnerable native rodent to its long term alien predators. Oikos. 2016;125(7):918–26.

    Article  Google Scholar 

  11. Griffin AS, Boyce HM. Indian mynahs, Acridotheres tristis, learn about dangerous places by observing the fate of others. Anim Behav. 2009;78(1):79–84.

    Article  Google Scholar 

  12. Ausilio G, Wikenros C, Sand H, Wabakken P, Eriksen A, Zimmermann B. Environmental and anthropogenic features mediate risk from human hunters and wolves for moose. Ecosphere. 2022;13(12):e4323.

    Article  Google Scholar 

  13. Wikenros C, Sand H, Wabakken P, Liberg O, Pedersen HC. Wolf predation on moose and roe deer: chase distances and outcome of encounters. Acta Theriol. 2009;54(3):207–18.

    Article  Google Scholar 

  14. Ranc N, Moorcroft PR, Ossi F, Cagnacci F. Experimental evidence of memory-based foraging decisions in a large wild mammal. Proc Natl Acad Sci U S A. 2021;118(15):e2014856118.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Ranc N, Cagnacci F, Moorcroft PR. Memory drives the formation of animal home ranges: Evidence from a reintroduction. Ecol Lett. 2022;25:716–28.

    Article  PubMed  Google Scholar 

  16. West R, Letnic M, Blumstein DT, Moseby KE, Hayward M. Predator exposure improves anti-predator responses in a threatened mammal. J Appl Ecol. 2018;55(1):147–56.

    Article  Google Scholar 

  17. Keen SC, Cole EF, Sheehan MJ, Sheldon BC. Social learning of acoustic anti-predator cues occurs between wild bird species. Proc Biol Sci. 1920;2020(287):20192513.

    Google Scholar 

  18. Manassa RP, McCormick MI. Social learning and acquired recognition of a predator by a marine fish. Anim Cogn. 2012;15(4):559–65.

    Article  PubMed  CAS  Google Scholar 

  19. Swift KN, Marzluff JM. Wild American crows gather around their dead to learn about danger. Anim Behav. 2015;109:187–97.

    Article  Google Scholar 

  20. Dow H, Fredga S. Breeding and natal dispersal of the goldeneye, Bucephala clangula. J Anim Ecol. 1983;52:681–95.

    Article  Google Scholar 

  21. Thiel D, Ménoni E, Brenot J-F, Jenni L. Effects of recreation and hunting on flushing distance of capercaillie. J Wildl Manag. 2007;71(6):1784–92.

    Article  Google Scholar 

  22. Ordiz A, Støen O-G, Sæbø S, Kindberg J, Delibes M, Swenson JE. Do bears know they are being hunted? Biol Cons. 2012;152:21–8.

    Article  Google Scholar 

  23. Ausilio G, Sand H, Månsson J, Mathisen KM, Wikenros C. Ecological effects of wolves in anthropogenic landscapes: the potential for trophic cascades is context-dependent. Front Ecol Evol. 2021;8:577963.

    Article  Google Scholar 

  24. Griesberger P, Obermair L, Zandl J, Stalder G, Arnold W, Hackländer K. Hunting suitability model: a new tool for managing wild ungulates. Wildl Biol. 2022;2022(3):e01021.

    Article  Google Scholar 

  25. Picardi S, Coates P, Kolar J, O’Neil S, Mathews S, Dahlgren D. Behavioural state-dependent habitat selection and implications for animal translocations. J Appl Ecol. 2021;59(2):624–35.

    Article  Google Scholar 

  26. Linnell JDC, Cretois B, Nilsen EB, Rolandsen CM, Solberg EJ, Veiberg V, et al. The challenges and opportunities of coexisting with wild ungulates in the human-dominated landscapes of Europe’s Anthropocene. Biol Cons. 2020;244:108500.

    Article  Google Scholar 

  27. Pascual-Rico R, Morales-Reyes Z, Aguilera-Alcala N, Olszanska A, Sebastian-Gonzalez E, Naidoo R, et al. Usually hated, sometimes loved: a review of wild ungulates’ contributions to people. Sci Total Environ. 2021;801:149652.

    Article  PubMed  CAS  Google Scholar 

  28. Stankowich T. Ungulate flight responses to human disturbance: a review and meta-analysis. Biol Cons. 2008;141(9):2159–73.

    Article  Google Scholar 

  29. Le Grand L, Thorsen NH, Fuchs B, Evans AL, Laske TG, Arnemo JM, et al. Behavioral and physiological responses of Scandinavian brown bears (Ursus arctos) to dog hunts and human encounters. Front Ecol Evol. 2019;7:134.

    Article  Google Scholar 

  30. Mysterud A, Rauset GR, Van Moorter B, Andersen R, Strand O, Rivrud IM, et al. The last moves: the effect of hunting and culling on the risk of disease spread from a population of reindeer. J Appl Ecol. 2020;57(12):2509–18.

    Article  Google Scholar 

  31. Grignolio S, Merli E, Bongi P, Ciuti S, Apollonio M. Effects of hunting with hounds on a non-target species living on the edge of a protected area. Biol Cons. 2011;144(1):641–9.

    Article  Google Scholar 

  32. Picardi S, Basille M, Peters W, Ponciano JM, Boitani L, Cagnacci F. Movement responses of roe deer to hunting risk. J Wildl Manag. 2019;83(1):43–51.

    Article  Google Scholar 

  33. Ciuti S, Muhly TB, Paton DG, McDevitt AD, Musiani M, Boyce MS. Human selection of elk behavioural traits in a landscape of fear. Proc Biol Sci. 2012;279(1746):4407–16.

    PubMed  PubMed Central  Google Scholar 

  34. Dijkgraaf L, Stenbacka F, Cromsigt JP, Ericsson G, Neumann W. Bear in mind! Bear presence and individual experience with calf survival shape the selection of calving sites in a long-lived solitary ungulate. Ecol Evol. 2024;14(3):e11177.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Ericsson G, Wallin K, Ball JP, Broberg M. Age-related reproductive effort and senescence in free-ranging moose, Alces alces. Ecology. 2001;82(6):1613–20.

    Article  Google Scholar 

  36. Swedish Environmental Protection Agency. Population estimates and distribution of large carnivores in Sweden. 2020;2021(06.11).

  37. Ericsson G. Reduced cost of reproduction in moose Alces alces through human harvest. Alces. 2001;37(1):61–9.

    Google Scholar 

  38. Morales JM, Moorcroft PR, Matthiopoulos J, Frair JL, Kie JG, Powell RA, et al. Building the bridge between animal movement and population dynamics. Philos Trans R Soc B Biol Sci. 2010;365(1550):2289–301.

    Article  Google Scholar 

  39. Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, Saltz D, et al. A movement ecology paradigm for unifying organismal movement research. Proc Natl Acad Sci. 2008;105(49):19052–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Dussault C, Ouellet JP, Courtois R, Huot J, Breton L, Jolicoeur H. Linking moose habitat selection to limiting factors. Ecography. 2005;28(5):619–28.

    Article  Google Scholar 

  41. Bjørneraas K, Solberg EJ, Herfindal I, Moorter BV, Rolandsen CM, Tremblay J-P, et al. Moose Alces alces habitat use at multiple temporal scales in a human-altered landscape. Wildl Biol. 2011;17(1):44–54.

    Article  Google Scholar 

  42. Singh NJ, Ericsson G. Changing motivations during migration: linking movement speed to reproductive status in a migratory large mammal. Biol Lett. 2014;10(6):20140379.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Beumer LT, Schmidt NM, Pohle J, Signer J, Chimienti M, Desforges JP, et al. Accounting for behaviour in fine-scale habitat selection: a case study highlighting methodological intricacies. J Anim Ecol. 2023;92(10):1937–53.

    Article  PubMed  Google Scholar 

  44. Paterson JT, Johnston AN, Ortega AC, Wallace C, Kauffman M. Hidden Markov movement models reveal diverse seasonal movement patterns in two North American ungulates. Ecol Evol. 2023;13(7):e10282.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Neumann W, Ericsson G. Influence of hunting on movements of moose near roads. J Wildl Manag. 2018;82(5):918–28.

    Article  Google Scholar 

  46. Swedish County Board. 2021. Available from: https://www.lansstyrelsen.se/.

  47. Heberlein TA. The gun, the dog and the thermos: culture and hunting in Sweden and the United States. Sweden America. 2000;13:24.

    Google Scholar 

  48. Swedish Association For Hunting and Wildlife Management. 2020. Available from: https://www.viltdata.se/.

  49. Kreeger T, Arnemo J. Handbook of wildlife chemical immobilization Laramie. Wyoming: International Wildlife Veterinary Services; 2007.

    Google Scholar 

  50. Ericsson G, Wallin K. Senescence in a northern ungulate: age and sex-specific patterns of mortality in moose, Alces alces. Ecoscience. 2001;8:157–63.

    Article  Google Scholar 

  51. Neumann W, Ericsson G, Dettki H. The non-impact of hunting on moose Alces alces movement, diurnal activity, and activity range. Eur J Wildl Res. 2008;55(3):255–65.

    Article  Google Scholar 

  52. Neumann W, Singh NJ, Stenbacka F, Malmsten J, Wallin K, Ball JP, et al. Divergence in parturition timing and vegetation onset in a large herbivore-differences along a latitudinal gradient. Biol Lett. 2020;16(6):20200044.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Passoni G, Coulson T, Ranc N, Corradini A, Hewison AM, Ciuti S, et al. Roads constrain movement across behavioural processes in a partially migratory ungulate. Mov Ecol. 2021;9:1–12.

    Article  Google Scholar 

  54. McClintock BT, Michelot T. momentuHMM: R package for generalized hidden Markov models of animal movement. Methods Ecol Evol. 2018;9(6):1518–30.

    Article  Google Scholar 

  55. Edelhoff H, Signer J, Balkenhol N. Path segmentation for beginners: an overview of current methods for detecting changes in animal movement patterns. Mov Ecol. 2016;4(1):21.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Zucchini W, MacDonald IL, Langrock R. Hidden Markov models for time series. New York: Chapman and Hall/CRC; 2017.

    Book  Google Scholar 

  57. Neumann W, Ericsson G, Dettki H, Bunnefeld N, Keuler NS, Helmers DP, et al. Difference in spatiotemporal patterns of wildlife road-crossings and wildlife-vehicle collisions. Biol Cons. 2012;145(1):70–8.

    Article  Google Scholar 

  58. Avgar T, Potts JR, Lewis MA, Boyce MS, Börger L. Integrated step selection analysis: bridging the gap between resource selection and animal movement. Methods Ecol Evol. 2016;7(5):619–30.

    Article  Google Scholar 

  59. Avgar T, Lele SR, Keim JL, Boyce MS. Relative selection strength: quantifying effect size in habitat- and step-selection inference. Ecol Evol. 2017;7(14):5322–30.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Thurfjell H, Ciuti S, Boyce MS. Applications of step-selection functions in ecology and conservation. Mov Ecol. 2014;2(1):4.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Signer J, Fieberg J, Avgar T. Animal movement tools (amt): R package for managing tracking data and conducting habitat selection analyses. Ecol Evol. 2019;9(2):880–90.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Stedman R, Diefenbach DR, Swope CB, Finley JC, Luloff AE, Zinn HC, et al. Integrating wildlife and human-dimensions research methods to study hunters. J Wildl Manag. 2004;68(4):762–73.

    Article  Google Scholar 

  63. Bonnot N, Morellet N, Verheyden H, Cargnelutti B, Lourtet B, Klein F, et al. Habitat use under predation risk: hunting, roads and human dwellings influence the spatial behaviour of roe deer. Eur J Wildl Res. 2013;59(2):185–93.

    Article  Google Scholar 

  64. Swedish Land Survey. Digital data maps of land cover, elevation and infrastructure 2002. Available from: http://www.lantmateriet.se.

  65. Swedish Forest Agency. skogsstyrelsen 2019. Available from: www.skogsstyrelsen.se.

  66. Trafikverket. Digital data maps on roads and settlements, vector 2014. Available from: www.trafikverket.se.

  67. Zuur AF, Ieno EN, Elphick CS. A protocol for data exploration to avoid common statistical problems. Methods Ecol Evol. 2010;1(1):3–14.

    Article  Google Scholar 

  68. Fortin D, Beyer HL, Boyce MS, Smith DW, Duchesne T, Mao JS. Wolves influence elk movements: behavior shapes a trophic cascade in Yellowstone National Park. Ecology. 2005;86(5):1320–30.

    Article  Google Scholar 

  69. Therneau T, Lumley T. R survival package. R Core Team. 2013;523.

  70. Fieberg J, Matthiopoulos J, Hebblewhite M, Boyce MS, Frair JL. Correlation and studies of habitat selection: Problem, red herring or opportunity? Philos Trans R Soc B Biol Sci. 2010;365(1550):2233–44.

    Article  Google Scholar 

  71. Muff S, Signer J, Fieberg J. Accounting for individual-specific variation in habitat-selection studies: efficient estimation of mixed-effects models using Bayesian or frequentist computation. J Anim Ecol. 2020;89(1):80–92.

    Article  PubMed  Google Scholar 

  72. Beardsworth CE, Whiteside MA, Laker PR, Nathan R, Orchan Y, Toledo S, et al. Is habitat selection in the wild shaped by individual-level cognitive biases in orientation strategy? Ecol Lett. 2021;24(4):751–60.

    Article  PubMed  Google Scholar 

  73. Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer; 2002.

    Google Scholar 

  74. Viana DS, Granados JE, Fandos P, Perez JM, Cano-Manuel FJ, Buron D, et al. Linking seasonal home range size with habitat selection and movement in a mountain ungulate. Mov Ecol. 2018;6:1.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Hox J, Moerbeek M, Van de Schoot R. Multilevel analysis: techniques and applications. Routledge; 2017.

    Book  Google Scholar 

  76. Snijders TA, Bosker R. Multilevel analysis: an introduction to basic and advanced multilevel modeling; 2011.

  77. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, 2021. 2022.

  78. Felton AM, Hedwall P-O, Felton A, Widemo F, Wallgren M, Holmström E, et al. Forage availability, supplementary feed and ungulate density: associations with ungulate damage in pine production forests. For Ecol Manag. 2022;513:120187.

    Article  Google Scholar 

  79. Zong X, Wang T, Skidmore AK, Heurich M. Estimating fine-scale visibility in a temperate forest landscape using airborne laser scanning. Int J Appl Earth Obs Geoinf. 2021;103:102478.

    Google Scholar 

  80. Graesli AR, Le Grand L, Thiel A, Fuchs B, Devineau O, Stenbacka F, et al. Physiological and behavioural responses of moose to hunting with dogs. Conserv Physiol. 2020;8(1):coaa122.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Neumann W, Martinuzzi S, Estes AB, Pidgeon AM, Dettki H, Ericsson G, et al. Opportunities for the application of advanced remotely-sensed data in ecological studies of terrestrial animal movement. Mov Ecol. 2015;3(1):8.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Meisingset EL, Gusevik J, Skjørestad A, Brekkum Ø, Mysterud A, Rosell F. Impacts of human disturbance on flight response and habitat use of red deer. Ecosphere. 2022;13(11):e4281.

    Article  Google Scholar 

  83. Whitehead H. Learning, climate and the evolution of cultural capacity. J Theor Biol. 2007;245(2):341–50.

    Article  PubMed  Google Scholar 

  84. Heithaus MR, Wirsing AJ, Thomson JA, Burkholder DA. A review of lethal and non-lethal effects of predators on adult marine turtles. J Exp Mar Biol Ecol. 2008;356(1–2):43–51.

    Article  Google Scholar 

  85. Weimerskirch H, Cherel Y, Delord K, Jaeger A, Patrick SC, Riotte-Lambert L. Lifetime foraging patterns of the wandering albatross: life on the move! J Exp Mar Biol Ecol. 2014;450:68–78.

    Article  Google Scholar 

  86. Wilmers CC, Nickel B, Bryce CM, Smith JA, Wheat RE, Yovovich V. The golden age of bio-logging: how animal-borne sensors are advancing the frontiers of ecology. Ecology. 2015;96(7):1741–53.

    Article  PubMed  Google Scholar 

  87. Fraser KC, Davies KT, Davy CM, Ford AT, Flockhart DT, Martins EG. Tracking the conservation promise of movement ecology. Front Ecol Evol. 2018;6:150.

    Article  Google Scholar 

  88. Ericsson G, Neumann W, Dettki H. Moose anti-predator behaviour towards baying dogs in a wolf-free area. Eur J Wildl Res. 2015;61(4):575–82.

    Article  Google Scholar 

  89. Singh NJ, Börger L, Dettki H, Bunnefeld N, Ericsson G. From migration to nomadism: movement variability in a northern ungulate across its latitudinal range. Ecol Appl. 2012;22(7):2007–20.

    Article  PubMed  Google Scholar 

  90. Kaczensky P, Khaliun S, Payne J, Boldgiv B, Buuveibaatar B, Walzer C. Through the eye of a Gobi khulan–application of camera collars for ecological research of far-ranging species in remote and highly variable ecosystems. PLoS ONE. 2019;14(6):e0217772.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We acknowledge the field staff and all the contributing members of the public, managers, and organizations of the hunters. We would especially like to thank Mark Boyce and Jörg Fabian Knufinke for their valuable and constructive contributions and assessments of earlier versions of this work. We would further like to express our gratitude to the two reviewers for their work and comments, which have helped to improve this work significantly.

Funding

Open access funding provided by Swedish University of Agricultural Sciences. The thematic program Wildlife and Forestry at the Swedish University of Agricultural Sciences, the Södra Skogsägarna Foundation for Research, the Swedish Environmental Protection Agency, the Swedish Association for Hunting and Wildlife Management, and SLU funded the data collection.

Author information

Authors and Affiliations

Authors

Contributions

LG, WN and HT analyzed the data. LG wrote the manuscript as a lead with the input of WN, HT and GE. GE and WN collected the data and provided funding. All the authors read the most recent draft and provided consent for publication.

Corresponding author

Correspondence to Lukas Graf.

Ethics declarations

Ethics approval and consent to participate

Marking and handling of study animals were approved by the Animal Care Committee for Northern Sweden in Umeå (Dnr A12-12, A14-15, A116-09, A124-05, A44-16, A50-12, A77-06) and was carried out in accordance with the Swedish laws concerning animal research ethics. All personnel were certified according to the standards of the Swedish Animal Welfare Agency and the Swedish Board of Agriculture.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Graf, L., Thurfjell, H., Ericsson, G. et al. Naivety dies with the calf: calf loss to human hunters imposes behavioral change in a long-lived but heavily harvested ungulate. Mov Ecol 12, 66 (2024). https://doi.org/10.1186/s40462-024-00506-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40462-024-00506-5

Keywords