Skip to main content

Movement behavior in a dominant ungulate underlies successful adjustment to a rapidly changing landscape following megafire

Abstract

Background

Movement plays a key role in allowing animal species to adapt to sudden environmental shifts. Anthropogenic climate and land use change have accelerated the frequency of some of these extreme disturbances, including megafire. These megafires dramatically alter ecosystems and challenge the capacity of several species to adjust to a rapidly changing landscape. Ungulates and their movement behaviors play a central role in the ecosystem functions of fire-prone ecosystems around the world. Previous work has shown behavioral plasticity is an important mechanism underlying whether large ungulates are able to adjust to recent changes in their environments effectively. Ungulates may respond to the immediate effects of megafire by adjusting their movement and behavior, but how these responses persist or change over time following disturbance is poorly understood.

Methods

We examined how an ecologically dominant ungulate with strong site fidelity, Columbian black-tailed deer (Odocoileus hemionus columbianus), adjusted its movement and behavior in response to an altered landscape following a megafire. To do so, we collected GPS data from 21 individual female deer over the course of a year to compare changes in home range size over time and used resource selection functions (RSFs) and hidden Markov movement models (HMMs) to assess changes in behavior and habitat selection.

Results

We found compelling evidence of adaptive capacity across individual deer in response to megafire. Deer avoided exposed and severely burned areas that lack forage and could be riskier for predation immediately following megafire, but they later altered these behaviors to select areas that burned at higher severities, potentially to take advantage of enhanced forage.

Conclusions

These results suggest that despite their high site fidelity, deer can navigate altered landscapes to track rapid shifts in encounter risk with predators and resource availability. This successful adjustment of movement and behavior following extreme disturbance could help facilitate resilience at broader ecological scales.

Background

Movement is a key trait that allows animal species to adjust to changing landscapes and track resources in dynamic ecosystems, such as fire-prone landscapes [1]. This ability has become increasingly critical in an age of constant anthropogenic global change and extreme environmental disturbances, such as increasingly frequent and severe megafires [2]. In fire-prone ecosystems, megafires, defined as wildfires larger than 100 km2 that surpass the size and severity of historical fires, have become increasingly prevalent [3]. Fire has served an important ecological and evolutionary role in many of these ecosystems [4], but historical policy, climate change, and land use change are responsible for the increasing number of unprecedented megafires. Megafires dramatically alter ecosystems at much broader scales than normal fires by rapidly removing resources and triggering rapid conversions in habitat [3]. Though many wild animal species in these fire-prone ecosystems have adaptations to coexist with their historic fire regimes [5, 6], novel megafires may challenge, and even overwhelm, the behaviors and adaptive capacity of individual animals [7]. Behavioral plasticity and movement play an important role in defining an animal’s capacity to adapt and adjust to novel disturbance regimes. Recent work has documented the role movement and behavioral plasticity play in governing the adaptive capacity of species to other forms of global change [8,9,10]. For large-bodied animals, plasticity in movement and behavior allows individuals to adjust to changes in their local environments [11, 12]. For fire specifically, larger-bodied animals may partition their space-use across recently burned landscapes to take advantage of new resources or avoid risky areas [13].

By rapidly altering landscapes, megafire may impact how some animals are able to navigate and use habitat. Burn severity is a characteristic of fire that defines the loss of below and above ground organic matter [14]. High severity burns, even those outside of megafires, can remove important structural resources from landscapes [15] and even cause direct mortality to animals [16]. Changes in structural cover in these systems may alter interspecies interactions, such as predator-prey dynamics, by altering the success of predator hunting strategies and prey predator-avoidance strategies [17]. Specifically, ambush predators, such as mountain lions (Puma concolor), may prefer more unburned areas that maintain cover to successfully ambush prey [18], while more cursorial predators, such as coyote (Canis latrans) may prefer more open areas following fire to find prey [19]. High severity fires may remove important food resources (i.e., forbs, grasses, seeds, etc.), potentially influencing herbivorous species populations [20]. Finally, the short and long-term effects of fire on habitat may be directly related to the dominant vegetation type of that habitat. Oak woodland ecosystems, which are a composite mix of grassland, woodland, and shrubland patches, are typically characterized by a moderate or mixed-severity fire regime in which fires burn patchily at low (grass and tree-understory) and high severities (tree-crowns and shrub-crowns) [21]. Recovery times vary across these different vegetation types as well, with grasslands typically recovering faster than oak trees and shrubs. For example, grassland ecosystems typically recover within the first year following fire [22], whereas shrubs and trees may take 5–10 years to fully regenerate following typical woodland fires [23, 24].

Ungulates serve key ecological roles in many fire-prone ecosystems around the world through their herbivory and by serving as a link between different trophic levels. Changes in their movement and behaviors following fire may have important implications for ecosystem-level processes. Fire may influence patterns of ungulate herbivory across landscapes over space and time [25, 26]. Past work has specifically documented a “magnet effect” across several ungulate species, where individuals select moderately burned areas that have improved forage post-fire [27, 28]. Following more severe fire events, recent work suggests that ungulate behavioral plasticity may buffer the short-term impacts of megafire, as ungulates select covered, woodland habitat and expand their home ranges to compensate for a decrease in foraging resources [29]. Kreling et al. (2021) also ask whether these adjustments could become maladaptive to these populations as megafires become more frequent. The seasonality of fire events may also modulate short- and long-term responses of ungulates, with fires potentially increasing the scarcity of rare vegetation resources during the dry seasons or limiting required resources during energetically costly periods of the year (i.e. spring breeding season) [30]. Ungulate browsers, such as Black-tailed mule deer (Odocoileus hemionus columbianus) across California, depend on forbs, acorns masts, and tree saplings especially during the dry, non-growing seasons [31], and fires that remove these key resources at these points in time may constrain resources until the next growing season.

Behavioral plasticity mediated by movement likely played an important role in shaping the evolution of ungulate populations in dynamic landscapes where plasticity allowed species to adapt to changes in resource availability caused by fire and other disturbances [32]. Variation in behavior across space permits animals to respond accordingly to dynamic landscapes [33, 34]. Behavioral plasticity, therefore, likely continues to play a significant role in influencing the resilience of ungulate species to major wildfire events and changes to local fire regimes. Unlike other large ungulates, behavioral plasticity of mule deer migratory movement specifically has been found to be non-plastic or rigid [35]. Alternatively, previous work has also established that mule deer are capable of efficiently navigating burned landscapes to simultaneously minimize predation risk from a variety of predators whilst identifying and using areas with forage [36]. Therefore, deer may not change the location of their range, even if it is entirely burned in large-scale events such as megafires, but they may alter their space usage within these large burn extents. Understanding the conditions and thresholds under which behavioral plasticity is adopted as an adaptive strategy may be key in tailoring management for this species following major environmental disturbances. Megafires may overwhelm this capacity, and it may take much longer for species to recover to pre-fire conditions due to the dramatic changes imposed on the landscape. It’s important to understand how ungulates adjust their behavior in burned areas over time. This insight is crucial for assessing adaptability to changing fire regimes and guiding management strategies.

In this study, we examined the long-term consequences of megafire on an ecologically and economically important Californian ungulate, the black-tailed mule deer. As a direct follow-up to a study on the initial effects of the 2018 Mendocino Complex Fire on deer behavior [29], we investigate how short-term movement responses of deer to megafire vary over the year following fire. Due to the scale and severity of this megafire, we hypothesized that (1) observed changes in deer behavior (increased home range size and habitat preferences) would persist until the end of the study period. Secondly, we hypothesized that (2) deer would preferentially select habitat that burned at low severity immediately following the fire (“Recently Burned”) to avoid exposure to predators and select areas more likely to have forage remaining. This, however, may not be the case if the activity of mountain lions, the predominant predator of deer in this system, decreases in exposed areas that lack the cover conducive for ambush hunting. In line with the magnet effect, (3) we predicted black-tailed deer would select areas that burned at moderate severities the following growing season (“First Spring”) due to the increased nutritional value of forage in these areas. We hypothesized that the behavioral plasticity of deer would allow them to adjust their movement to minimize risk and maximize access to resources. Specifically, we hypothesized that (4) black-tailed deer would be more likely to travel through severely burned areas to avoid exposure to potential predators, and to rest in low severity burned areas where perceived risk may be lower. Similar to deer habitat selection, changes in deer behavioral modes are likely dependent on whether the activity of their predominant predator, mountain lion, changes in these severely burned areas as well where ambush hunting may be less successful. Large, high severity patches of this fire removed extensive shrub cover and top-killed patches of oak woodland throughout the study site. We predicted these large-scale, structural changes would lead to long-lasting behavioral adjustments in both habitat selection and behavioral modes that would be consistent throughout the study period.

Methods

Study site and fire history

We conducted this study at the Hopland Research and Extension Center (HREC hereafter) in Mendocino County in northern California (39°00′ N, 123°04’ W) (Fig. 1). HREC is composed of a diverse set of vegetation types including chaparral shrublands (Adenostoma fasciculatum), oak woodland savannah (Quercus kelloggii, Quercus douglassi, and Quercus lobata), and a mix of introduced and native, open grassland. HREC is characterized by a Mediterranean climate with cool, wet winters and warm, dry summers. HREC also operates as a working rangeland landscape, containing a sheep farming facility and several agricultural plots throughout the property [37]. Deer hunting occurs seasonally on site generally from August-September. A maximum of 120 hunters are allowed on site annually (20 hunters per day) and an average of 23.4 deer are harvested annually [38]. No hunting was permitted during 2018 because of megafire. Mountain lions are the predominant predator of mule deer on site, but coyote and black bear (Ursus americanus) also occur on site and have been observed to prey on deer and their fawns periodically [39].

Fig. 1
figure 1

Maps of the 2018 Mendocino Complex Fire and the study site, the U.C. Hopland Research and Extension Center (HREC) (39°00′ N, 123°04’ W). Map “a” displays the total burn perimeter of the Mendocino Complex Fire. This fire burned into HREC on July 27, 2018. Map “b” displays the severity of the fire across the HREC property boundary. Sentinel-2 satellite imagery was acquired via Google Earth Engine to calculate fire severity. Fire severity was quantified as the Differenced Normalized Burn Ratio (dNBR). For visualization purposes, dNBR values were binned into categorical values based on those established by US Geological Survey as follows (Unburned = 0–99, Low = 99–269, Moderate-Low = 269–439, Moderate-High = 439–659, High = 659+). Map “c” displays the compositional makeup of dominant vegetation types across HREC. In this map, yellow denotes grassland, green denotes woodland, and brown denotes chaparral shrubland

On July 27, 2018, the River Fire (southern half of the Mendocino Complex Fire), swept through the northern half of HREC, burning approximately 13.76 km2 (65%) of the property. The 2018 Mendocino complex fire burned 1,858 km2 total and is currently the third largest wildfire in California’s recorded history (CALFIRE-FRAP, 2022). Fires in this region typically burn relatively frequently every 5–15 years at relatively low severities in the more open woodland and grassland habitats and more infrequently, but more severely, in the dense shrubland chaparral habitats every 30–60 + years [40, 41]. The River Fire burned a much larger contiguous area and much more severely than recent fires within HREC. Atypical of fires in woodland fire regimes, several oak trees (Q. kelloggii, Q. douglassi, and Q. lobata), whose acorn masting normally provides a key food resource for local deer populations, were top-killed in certain high severity patches of the this fire. Acorn masting remains highly variable across different oak species, as well as between individual trees, but mature trees typically mast every 2–3 years dependent on climatic factors and interspecies interactions [42].

Monitoring black-tailed deer movement and home range estimation

We deployed GPS-collars (Vertex Plus and Lotek Iridum Track M) across 28 female deer between July 2017 and July 2019. These data provided the basis for a natural experiment to observe the effects of megafire on deer movement and behavior. We programmed all collars to record GPS locations once per hour. Deer were captured using Clover traps and were manually restrained to place collars on, without the use of chemical immobilizers (permit #P1680002). We monitored deer remotely post-capture for 1 week to ensure collared deer did not experience capture myopathy and that the distance they traveled within a day was typical of a healthy individual.

To observe how deer movement and behavior changed over time following the megafire, we subset the collected GPS data to only include deer that had GPS points that overlapped the fire perimeter of the Mendocino Complex Fire, excluding seven individual deer. We then subset the collected GPS points temporally into three two-month-long time periods: just after the fire (August 1st – October 1st 2018), the first spring green up following the fire (March 1st – May 1st 2019), and one full year post-fire (August 1st – October 1st 2019) (Additional File 1 - Table S1). We included two additional pre-fire time periods to compare deer home range size before and after the fire and to examine any seasonal differences in home range size that may impact our results. These additional pre-fire time periods included: two separate spring seasons before the fire (March 1st – May 1st 2017 and March 1st – May 1st 2018) and just before the fire (May 25th – July 25th 2018). We collected > 1,000 GPS fixes for most deer within each time period (29 of 38 deer-periods) (Additional File 1 – Table S2). The “Prefire” time period occurs at the intersection of our annual collaring efforts when many collars from the previous year are programmatically dropped off and collars for the coming year are deployed. This resulted in us collecting fewer GPS locations from individual deer during this time period. To facilitate analyses that included this key time period, we therefore included deer-periods that had a minimum of 500 recorded GPS locations, excluding four individual deer-periods. This resulted in 21 unique individuals collared across these five time periods. Thirteen of these animals maintained their collars across two or more time periods, resulting in 38 period-specific deer home ranges (Additional File 1 – Table S2). We removed 10 erroneous, outlier GPS locations that were greater than 2 km from their consecutive points for these deer between hour fixes.

For each deer, and within each time period (deer-period), we used the two months of collected GPS data to estimate individual home range sizes. We used 95% Kernel Utilization Densities (KUD) in the “adehabitatHR” (v.0.4.19) package in “R” (v.4.1.1) to create these home ranges [43,44,45]. To assess whether deer home range sizes continue to change following the megafire, we used paired Welch’s unequal variance t-test to compare deer home range sizes 1) just after fire (“Recently Burned”), 2) the first spring following fire (“First Spring”), 3) one full year post-fire (“1 Year Post-Fire”), 4) the spring seasons before fire (“Pre-spring”), and 5) just before the fire burned (“Pre-fire”). To assess the robustness of this analysis to the smaller sample size of under sampled time periods (i.e. “Prefire”), we randomly sampled 500 GPS locations from each deer-period and repeated this analysis to compare changes in home range size with the rarefied data as well.

Environmental covariates

We compiled fire and other environmental covariates alongside deer movement data to evaluate black-tailed deer movement responses to the megafire over time. We expected that fire severity, predator encounter probability, vegetation type, distance to water, and time since burning would be strong predictors of both deer habitat selection and deer movement during each post-fire time period. We originally planned to include NDVI (Normalized Difference Vegetation Index) as a measure of forage availability across the landscape, but we found measures of NDVI were highly correlated with measures of fire severity, our primary covariate of interest. Therefore, we included fire severity and excluded NDVI. To quantify fire severity on the landscape [46], we calculated the differenced Normalized Burn Ratio (NBR) collected via Sentinel-2 [47] satellite imagery (10 m resolution) and processed in Google Earth Engine [48] from both before (July 25th, 2018) and after (August 25th, 2018) the fire. NBR was calculated using the following equations (14):

∆NBR = NBRprefire- NBRpostfire.

NBR = Near-infrared (NIR) – shortwave infrared (SWIR) / Near-infrared (NIR) + shortwave infrared (SWIR).

We also included a quadratic term for fire severity to examine whether deer may preferentially select for moderately burned areas that, according to the magnet effect, may eventually have more nutritious forage after vegetation regrowth.

To account for mountain lion encounter probability across the landscape for this study, we included a high-resolution mountain lion habitat suitability map produced for the entire State by Dellinger et al. 2020 in our analyses [49]. This habitat suitability modeling effort used a suite of biotic and abiotic variables, including terrain ruggedness, canopy cover, and a rough categorical estimate of deer density. Though predator occurrence is not necessarily associated with predation risk [50], and we do not have data of observed mountain lion kills at the site to create a study-specific risk map, we used this habitat suitability map to serve as a proxy of perceived risk [51] for deer across our study site.

We classified the study site into three broad land cover categories: woodland, shrubland (chaparral), and grassland. To do this, we hand digitized vegetation layers using high-resolution (< 1 m) aerial imagery from the National Agriculture Imagery Program (2014–2015). In 2015, we ground truthed these digitizations by checking 50 randomly generated points across the study site to validate classifications (results were 98% accurate). Our primary interest was to compare the strength of selection and avoidance of broad vegetation types by deer following fire. Therefore, we chose to represent vegetation types as three dominant land cover categories as opposed to continuous covariates. Previous work in this region found that deer prefer woodland habitat following fire [29], so we chose to use woodland as the reference category within our model to compare against the other vegetation types. To calculate the distance between collected GPS points and potential water sources, we obtained stream vector data from the National Hydrography Dataset and used the “sf” (1.0.2) package in R to calculate the distance between collected GPS locations and seasonal streambeds throughout the study site.

We checked the (Variance Inflation Factor) score of covariates to ensure there was no underlying collinearity between modeled covariates (VIF < 3) and qualitatively inspected plotted covariates as well (Additional File 1: Figure S1; Additional File 1: Figure S2).

Resource selection functions

We used Resource Selection Functions (RSFs) to assess black-tailed deer habitat selection across each post-fire time period. Previous work illustrates that point (RSFs) and path (SSFs) selection functions work similarly well in defining habitat selection [52], but given the small home range sizes for most deer in our study and their proportionally large step-lengths, we chose to employ RSFs specifically as most of their defined home range should be “available” to use at any given time step. We collected > 1,000 GPS locations for all deer-periods used in this analysis (26 of 38 deer-periods, 16 unique individuals) (Additional File 1: Table S2). For these RSFs, we used the home ranges defined by the 95% Kernel Utilization Densities (KUDs) [43,44,45]. We modeled habitat selection for all time periods combined to improve interpretability of model results. We also included an interaction term between time period and severity (Severity*BurnLag). For each deer-period home range, we randomly generated four-times as many “available” points from within each deer’s estimated KUD home range (mean available point density across all modeled deer = 8,174.72 available points/km2) [29]. Available points were stratified by time period so that the number of available points had the same ratio across time periods as the true use points. We compared the environmental characteristics of “used” and “available” GPS points using a mixed effects logistic regression via the “lme4” (v.1.1.27.1) package in R [45, 53].

We used an a priori hypothesis-driven approach to select a model to describe deer habitat selection, that included fire severity and its quadratic term (to account for nonlinear effects), encounter probability with mountain lions, vegetation type (chaparral, woodland, or grassland), distance to water, time since burn, and an interaction between severity and time since burn as covariate predictors. We used woodland as the reference vegetation category within these RSFs. We randomly sampled “time since burn” for each available point as a randomly selected date from within its respective time period. Prior to modeling, we standardized each of the included covariates (mean = 0, standard deviation = 1). We included a random intercept of “Deer ID” within our RSF to account for individual differences in behavior and resource availability for each deer (individual deer retained their same “Deer ID” across time periods).

To assess goodness of fit of the RSF model, we used the “performance” (v.0.7.3) package in R [54] to calculate marginal and conditional R2 values for the model and visually inspect overall model fitting.

Hidden markov movement models

While examining habitat selection provides an important opportunity to uncover where animals tend to spend time across landscapes, it is equally important to understand how animals use the time they spend in the areas they are selecting or avoiding as mediated by behavior.

By defining certain movement parameters (i.e. turning angle and step-length), we can use hidden Markov movement models (HMMs hereafter) to predict behavioral states of animals at individual GPS-fixes and compare how the distribution of these states may change in response to environmental covariates across a landscape, such as fire [55]. These behavioral states represent types of responses to an animal’s environment such as “foraging”, “traveling”, or “resting” [56].

To assess how deer behavioral decisions were impacted by megafire, we fit a HMM across the combined, three post-burn time periods within our study (26 deer-periods, 16 individuals) using the “moveHMM” (v.1.8) package within R [57]. We modeled two behavioral states (state 1 = resting, state 2 = traveling) to increase model interpretability and to specifically observe whether deer traveling and resting behavior changes across landscape variables to avoid perceived risks following fire. We calculated step lengths (via von Mises distributions) and turning angles (via gamma distribution) to characterize the two behavioral states. We randomly generated 25 different pairs of starting values from a range of plausible values as defined by the range of each calculated movement parameter (step-length and turning angle) (Additional File 1: Table S3). We ran each of the 25 randomly generated step-length and turning-angle pairs in a model without covariates and compared the negative-log likelihood of each model. We checked that each model had similar maximum log-likelihood values and we selected the best fitting pair of movement parameters based on maximum likelihood [58].

Using these starting values, we fit a single hidden Markov model with a set of a priori selected covariates (Severity + Mountain Lion Encounter Probability + Distance to Water + Time Since Burn + Vegetation Cover + Severity*Time Since Burn) to estimate how the probability of being in a certain behavioral state (i.e. resting vs. traveling) changed as a function of these environmental factors. We then used the “stationary” function of the “moveHMM” package to estimate the probability of each GPS point being in a given behavioral state and used these to create activity budgets by summing the estimated probabilities for being in each state at each recorded GPS point [59]. We used a Chi-squared test to assess whether the proportions of the two behavioral states were significantly different across time periods.

We assessed goodness of fit for the HMM using pseudo-residuals drawn from the fit model. Pseudo-residuals of the step length parameter should be normally distributed given good model fit [60, 61]. Therefore, we visually inspected step length pseudo residuals and used a Shapiro-Wilk normality test using a random subset of pseudo residual values (n = 1000).

Results

Home range comparison across seasons

We found that the average deer home range size across all time periods was 0.75 km2 (sd ± 0.42). Deer home range sizes were largest in the two time periods directly following the fire (“Recently Burned” and “First Spring”) and were smallest in the two pre-fire time periods (“Prespring” and “Prefire”) as well as “1 Year Post Fire”. The average home range size was 0.95 km2 (sd ± 0.46) during the “Recently Burned” period and 1.08 km2 (sd ± 0.35) during the “First Spring” period. During the “1 Year Post Fire” time period the average home range size was 0.38 km2 (sd ± 0.13). Finally, during the pre-fire time periods, the average deer home range size was 0.43 km2 (sd ± 0.11) for the “Prespring” time period and 0.52 km2 (sd ± 0.11) during the “Prefire” time period. (Additional File 1: Table S1; Fig. 2). We found no significant difference between deer home range sizes during the “Recently Burned” and “First Spring Periods (t = -0.72, df = 14.71, p-value = 0.48). We did find a significant difference (p-value < 0.05) in deer home range size between the “Recently Burned” and “1 Year Post Fire” periods (t = 3.52, df = 9.82, p-value < 0.01), as well as between the “First Spring” and “1 Year Post Fire” periods (t = 5.98, df = 13.95, p-value < 0.01). We found no significant differences between the home range sizes of the two pre-fire time periods, “Prespring” and “Prefire” (t = 1.26, df = 6.11, p-value = 0.25). We also found no significant difference in home range size between the “1 Year Post Fire” and “Prespring” periods (t = 1.396, df = 8.787, p-value = 0.20), as well as between the “1 Year Post Fire” and “Prefire” periods (t = -0.83, df = 9.80, p-value = 0.42) (Table S4).

Fig. 2
figure 2

Home range size of black-tailed deer (Odocoileus hemionus columbianus) across five time periods both before and after the 2018 Mendocino Complex Fire in Hopland, California, USA. The Mendocino Complex Fire burned July 27th, 2018. These study periods include: 2017 Spring and 2018 Spring before the fire (“Prespring”), the summer season just before the fire burned (“Prefire”), directly following the fire (“Recently Burned”), the first spring following the fire (“First Spring”), and 1 full year post fire (“1 Year Post Fire”) (from left to right)

We found identical results in our additional analysis using a subset of 500 randomly selected points from each deer-period suggesting that our analysis is robust to differences in sample sizes across time periods (Additional File: Table S5; Additional File: Figure S3).

Resource selection functions

Overall, deer avoided areas that burned at high severity, but this response was nonlinear and probability of use was highest at intermediate severities (Table 1). However, we also found that deer habitat selection of fire burned areas changed over time as an interaction with time since burn. During the “Recently Burned” time period, deer were more likely to avoid high severity areas (Additional File 1: Figure S4). Conversely, deer selected higher severity burned areas during the final “1 Year Post Fire” period (Fig. 3). Deer preferred woodland habitat over grassland and chaparral following the fire (Additional File 1: Figure S4). Deer also avoided areas of high mountain lion encounter probability (mean = -0.09 [-0.08, -0.10]) (Table 1).

Table 1 Listed output estimates for each covariate of the resource selection function for black-tailed deer (O. Hemionus columbianus) following the 2018 Mendocino Complex Fire at the Hopland Research and Extension Center in Mendocino County, CA, USA. Beta-coefficients, standard errors, and p-values are listed for each covariate included in the model. For categorical vegetation types, “woodland” was used as the categorical variable. * indicates statistically significant predicter of habitat selection within the model (p-value < 0.05)
Fig. 3
figure 3

Plotted response curves of deer (Odocoileus hemionus columbianus) habitat selection in response to fire severity and time since fire, as predicted from a resource selection function following the 2018 Mendocino Complex Fire at the Hopland Research and Extension Center, CA, USA. To visualize the interaction, we used the midpoint date of each time period to represent a categorical “Time Since Burn” variable in the plot

Hidden markov model results

We found that the 25 iterations of our null model converged on very similar scores of maximum likelihoods (mean = 266750.40; sd = 956.38). Our best fit hidden Markov model estimated two deer behavioral states: a “resting” state with shorter step-lengths and wider turn angles and a “traveling” behavioral state with longer step-lengths and near 0 turning angles (Additional File 1: Figure S5; Additional File 1: Figure S6). We found a significant difference in the composition of behavioral states between all time periods, with deer spending a greater proportion of time traveling than resting immediately following fire and during the first spring (χ2 = 232.97, df = 2, p-value < 0.001) (Fig. 4; Additional File 1: Table S6) compared to the proportion of time spent in each state during the “1 Year Post Fire” time period.

Fig. 4
figure 4

Behavioral state proportions for black-tailed deer (Odocoileus hemionus columbianus) at the Hopland Research and Extension Center in Mendocino County, California. Behavioral states for deer tracks were estimated for each post-fire time period by the hidden Markov model. State frequencies represented the summed probabilities of each GPS point being in a specific behavioral state

We found that deer behavioral states changed as a function of fire severity. Deer were most likely to be in the “resting” behavioral state in unburned and moderately burned areas across all time periods (Fig. 5a). At high severities, deer were more likely to be in the “traveling” behavioral state. The probability of deer being in the “traveling” behavioral state at high severities was significantly higher during the “Recently Burned” time period than in the “First Spring” and “1 Year Post Fire” time periods (Fig. 5b).

Fig. 5
figure 5

Behavioral state probabilities of black-tailed deer (Odocoileus hemionus columbianus) as a function of fire severity and time periods following the 2018 Mendocino Complex Fire at the Hopland Research and Extension Center in Mendocino County, California. Plot a) displays the probability of deer being in the “resting” behavioral state as a function of fire severity across the three time periods (“Recently Burned”, “First Spring” and “1 Year Post Fire”). Plot b) displays the probability of deer being in the traveling as a function of across the three same time periods. Shaded areas represent 95% confidence intervals. *Note that y-axis scaling of plots a and b are different

Pseudo residuals drawn from the HMM suggested good model fit for the deer track data. Overall, plotted pseudo-residuals of deer step-lengths appeared normally distributed (Additional File 1: Figure S7). We failed to reject the null hypothesis of the Shapiro-Wilks significance test (W = 0.99, p-value = 0.61), suggesting pseudo-residuals were drawn from a normal distribution.

Discussion

In this study, we utilized deer movement data collected opportunistically before, during and after a megafire to examine how deer behavior and space use changes following severe environmental disturbances. We found evidence that refuted our initial hypotheses that deer home range size and habitat selection would remain constant throughout the study period due to the severity and size of the megafire. Instead, estimated deer home range size and habitat selection was shown to change significantly throughout the course of the study. We found evidence supporting our hypotheses that deer preferred low severity burned areas immediately following the megafire and preferred moderate-high severity burned areas later in the study. Finally, we found evidence supporting our hypothesis that deer have some behavioral flexibility in adjusting their movement behaviors (traveling vs. resting) across varying burn severities. The distribution of these behavioral modes across low and high severity burned areas also changed over time. Our results demonstrate the mechanisms in which movement and behavior underpin the capacity of black-tailed deer to effectively adjust to a quickly shifting landscape follow megafire.

Contrary to our original hypotheses, we found that black-tailed deer habitat selection and the composition of movement-inferred behavioral states changed as a function of fire severity and time. As Kreling et al., 2021 found, ungulate home ranges were larger directly following the megafire, but we found that this effect does not persist over time. Deer home range size was significantly higher during the first two time periods following megafire (“Recently Burned” and “First Spring”) compared to the pre-fire time periods (“Prespring” and “Prefire”) and “1 Year Post-Fire”. The scale of the change in home range size observed in this study exceeds what has been previously observed in other studies caused by normal inter-season variation [36, 62], suggesting megafire had a significant effect on deer home range size and space usage. These results corroborate the conclusions drawn from Calhoun et al. 2023, which found that the intensity at which black-tailed deer used recently burned areas decreased during the year of this megafire, but returned to pre-fire conditions one-year following the megafire [63]. Deer space use likely becomes more diffuse as their home ranges expand immediately following fire.

Directly following the megafire, deer strongly avoided areas that burned at high severity, but this effect waned in the initial spring months following fire and inverted by the “1-Year Post-Fire” time period, with deer instead selecting for habitat that burned at higher severities. Similarly, deer were more likely to move than to rest in high severity areas immediately following megafire, but this effect diminished over the course of the year. Our results suggest that the behavioral adjustments made by individuals may be effective coping mechanisms for the potential consequences of megafire. This observed behavioral plasticity may also allow deer to eventually take advantage of the resulting resources that become available over time.

Fire severity was a significant predictor of deer habitat selection, but we found that the direction of selection (against high severity areas vs. towards high severity areas) changed as an interaction with the amount of time that had passed since the fire burned. As observed in previous studies [29], black-tailed deer avoided high severity burned areas in the immediate aftermath following the fire. These observed behaviors may be used to avoid exposure to predators in open areas, but the exact mechanism requires further research. While deer may perceive greater predation risk in these exposed areas, mountain lions, their primary predators, may be paradoxically less likely to use recently burned, uncovered areas that limit ambush success [17, 18]. Alternatively, recent work in this same system has shown that coyote, an alternative predator, increases their occurrence and activity in recently burned areas [63], potentially dissuading deer from using these areas as well. This may be especially pertinent for does with fawns or yearlings that are vulnerable to predation by coyote [39]. Future work that compares the influence of perceived and actual predation risk across these potential predators could better elucidate the mechanism underlying these observed responses.

In agreement with our initial hypothesis, we found that deer began to select areas that burned at moderate severities during the first spring green-up. However, we did not anticipate that deer would select high severity burned areas during the final time period of the study (“1-Year Post-Fire”). We expected that the high severity burned areas would be depleted of resources for the duration of our study, but our results suggest that once the vegetation in these severely burned areas is able to recover, these areas may attract herbivorous species [64]. This is at least partially supported by a brief, qualitative look at the relationship between forage quality, as represented by EVI (Enhanced Vegetation Index), and fire severity in the latter time periods of the study (Additional File 1: Figure S8) in which forage quality is higher at moderate severities during “First Spring” and increases at high severities during “1 Year Post Fire”. Both ‘First Spring’ and ‘1 Year Post-Fire’ time periods support the ‘magnet effect,’ observed in other studies, where ungulates prefer recently burned areas with enhanced forage [65, 66]. This emphasizes fire’s significant role in ungulate ecosystems, despite broader changes in fire regimes and climate.

Changes in deer behavior over time following fire are likely mediated by how and when vegetation recovers. Though we anticipated that changes in habitat selection would persist longer in severely burned woodland and shrubland areas that take longer to regenerate relative to grasslands [23, 24], our results suggest that these fire-adapted, early-successional plant species can still provide new food resources over much quicker time scales. Alternatively, fire-following forbs and new tree shoots may provide a source of new foraging opportunities for herbivores, particularly for browsing ungulates [67].

The proportion of estimated deer behavioral states changed over the course of the study and the probability of deer being in certain behavioral states (traveling vs. resting) varied significantly with fire severity. However, contrary to our initial hypotheses, we found the probability of deer being in a certain behavioral state as function of severity was not constant across time periods, but instead changed across time. Initially following megafire, we found that deer were more likely to travel through high severity burn areas and spent more time resting in low severity areas. This strategy may allow deer to avoid spending too much time in exposed areas, and to spend more time in the limited areas that contain food and shelter [13]. Spending less time traveling in exposed areas during the late summer may be critical for thermoregulation, specifically to avoid heat stress, and may be especially important for does with fawns [68, 69]. Immediately following megafire, deer spent significantly more time in the “traveling” state than the “resting” state, but this effect waned over time. This, as both a result of and in combination with an overall decrease in resource availability, may potentially result in the initial decreased body condition of ungulates following megafire observed in the study area [29]. Despite this, our findings suggest that black-tailed deer may find benefits in moderate-high severity burned areas over an extended period following megafire that may outweigh these initial costs, but additional research is required to examine how deer body condition changes over time.

Our results suggest that black-tailed deer in this study area have great capacity for short-term behavioral plasticity to allow quick adjustments of their behavior patterns in response to disturbance and vegetation recovery. This mirrors similar findings in nearby deer populations that were found to adjust to seasonal changes in food availability by adjusting their habitat selection in the summer and winter months [70] and may suggest a broader evolutionary adaptation to these types of extreme fire regimes [71]. While the overall extent of this megafire may represent an extreme disturbance event, aspects of this disturbance still resemble key evolutionary pressures to which these animals have adapted to. Behavioral adaptions to dynamic, fire prone landscapes may facilitate deer decisions to maintain high site fidelity in anticipation of the eventual regrowth of increased vegetation resources in burned areas as time passes [72]. While previous work has shown that species traits and fire characteristics play an important part in creating mammalian resilience to fire disturbances [73], this study highlights the key role of behavior underlying these observed effects.

Our results highlight the resilience of black-tailed deer and oak woodland ecosystems overall to megafire via adaptations both have evolved within fire-prone ecosystems. The strategies several plant and animal species have acquired by evolving in fire-prone ecosystems may also be employed when responding to megafire. For large-bodied and vagile animals in particular, movement and behavioral strategies used to adjust within landscapes that burn regularly likely play a significant role in facilitating resilience to changes in resource availability caused by megafire as well [74]. Species in more dynamic landscapes have likely developed some degree of behavioral plasticity to deal with the unpredictability of resources more effectively [75, 76]. This behavioral plasticity may also buffer these species from other types of climatic disturbances [77,78,79], and should be incorporated in future studies that examine the predicted impacts of climate change on species conservation. This study also conveys how movement and behavior are the primary mechanisms by which resilience to megafire may be facilitated as observed in previous work [63]. The success of behavioral plasticity in buffering the effects of megafires likely depends on the degree to which these events depart from the historical fire regimes these species have adapted to [7]. We therefore require more work that can better define the thresholds of extent, severity, and frequency for specific megafires that may overwhelm the behavioral plasticity of ungulates and other fire adapted species.

We found that deer selected for woodland habitat relative to the other vegetation types and strongly avoided shrubland habitat. Deer likely avoided these open shrubland areas to avoid conspicuous encounters with predators [80] or to avoid exposure to temperature extremes in the late summer. Future work that compares these two potential mechanisms using concurrent movement data from ungulates and their predators to asses changes in predation [17], or physiological data related to the thermoregulatory cost of heat exposure could help fill this remaining gap. The limited sample size of our study inhibits our ability to define specific thresholds at which deer habitat selection and movement decisions change across these different environmental gradients (i.e. severity, land cover, and mountain lion encounter probability) or the specific strength of these responses. Despite this, we believe that our findings offer robust insights into the general relationship between deer behavior and these environmental variables within this specific ecosystem type that may guide more targeted studies in the future that are able to define these thresholds across different contexts.

During this study, we observed a preference for burned areas by black-tailed deer in the “First Spring” and “1 Year Post Fire” time periods, potentially highlighting some of the benefits of returning wildfire to fire adapted ecosystems. Whereas megafire is a more extreme example of fire disturbance, more moderate disturbances such as prescribed fire or managed wildfire are known to perform important ecological work in maintaining key ecosystem functioning for local communities [81] and generating improved habitat and resources for wildlife [82], without the more deleterious impacts created initially by megafire. These managed wildfire approaches also serve an important function in reducing the incidence of megafires by promoting landscape heterogeneity and reducing continuous fuel loads [83, 84]. Deer utilization of different habitat patches that have burned at different severities across different portions of the year also further supports the theorized importance of pyrodiversity in these landscapes as well [85]. Thus, utilizing fire management may simultaneously accomplish important wildlife conservation goals (habitat creation and maintenance) and wildfire management goals (megafire prevention) in similar fire-prone ecosystems.

We found evidence to suggest that deer are resilient to the impacts of megafire on a 1-year time scale, but more work is necessary to understand whether these initial responses translate into longer term resilience. The lagged effects of megafire may present novel challenges to species that have adapted to historical fire regimes by altering longer cycles in resource availability [86] as well as interspecies interactions [87]. For example, in oak woodland savannas where acorn masting is a primary food resource for many herbivorous species [88, 89], megafires that top-kill mature oak trees could dramatically alter the availability of these resources until oaks are able to regenerate and begin masting again. These indirect impacts could have powerful effects on future population dynamics through responses like individual fitness and reproduction across the previously burned landscape. Future work that examines how megafire influences the density and demographic trends of mule deer across longer time periods could help assess the resilience of the species across broader temporal scales and help define the potential consequences of megafire on longer-term interspecies interactions such as herbivory and predation.

Conclusion

Black-tailed deer are known generally for their tendency towards non-plastic behaviors and high site fidelity to their home ranges and migration routes [29, 35]. Our study, however, revealed additional nuance in our understanding of deer behavior in response to sudden change. Despite having naturally high site fidelity in the region of our study, we found that this black-tailed deer population had a great deal of adaptive capacity to change their movement and behavior to respond to the impacts and eventual resources following megafire. Climate change and climatic disturbances (such as megafire) may have a more severe impact on species that are unable to adjust their behavior to accommodate sudden changes in their environments [90]. Resilience of dominant herbivores could help facilitate ecological resilience at broader trophic levels following disturbance. We can help facilitate and boost the natural resilience we observed of black-tailed deer and other ungulates through land and fire practices that promote the benefits of fire while simultaneously avoiding the immediate drawbacks of megafire. Identifying the mechanism by which these layers of resilience are produced would not be possible without uncovering the nuances of animal behavior that underly these observed responses.

Data availability

Deer GPS data and novel code generated and analyzed for this study are archived in the Dryad repository: Calhoun, Kendall et al. (2023), Movement behavior in a dominant ungulate underlies successful adjustment to a rapidly changing landscape following megafire, Dryad, Dataset, doi:https://doi.org/10.6078/D12H83. Deer collar data are stored as a CSV spreadsheet of all GPS locations for each deer throughout the study period. We utilize novel code to perform the analyses of this project. Sentiel-2 satellite imagery (https://sentinel.esa.int/web/sentinel/missions/sentinel-2) was obtained and processed using Google Earth Engine (https://earthengine.google.com/) to create fire severity rasters. The Hopland Research and Extension Center Shapefile Boundary and water streams shapefile used to extract satellite imagery and calculate distance to water, respectively, was originally downloaded from Hopland REC (https://hrec.ucanr.edu).

Abbreviations

HREC:

Hopland Research and Extension Center

RSF:

Resource Selection Function

HMM:

Hidden Markov Movement Model

NDVI:

Normalized Difference Vegetation Index

EVI:

Enhanced Vegetation Index

References

  1. Abrahms B, Aikens EO, Armstrong JB, Deacy WW, Kauffman MJ, Merkle JA. Emerging perspectives on Resource Tracking and Animal Movement Ecology. Trends Ecol Evol. 2021;36(4):308–20.

    Article  PubMed  Google Scholar 

  2. Stott P. How climate change affects extreme weather events. Science. 2016;352(6293):1517–8.

    Article  CAS  PubMed  Google Scholar 

  3. Linley GD, Jolly CJ, Doherty TS, Geary WL, Armenteras D, Belcher CM, et al. What do you mean, ‘megafire’? Global Ecol Biogeogr. 2022;31(10):1906–22.

    Article  Google Scholar 

  4. McLauchlan KK, Higuera PE, Miesel J, Rogers BM, Schweitzer J, Shuman JK et al. Fire as a fundamental ecological process: Research advances and frontiers. Durigan G, editor. J Ecol. 2020;108(5):2047–69.

  5. Jones GM, Kramer HA, Whitmore SA, Berigan WJ, Tempel DJ, Wood CM, et al. Habitat selection by spotted owls after a megafire reflects their adaptation to historical frequent-fire regimes. Landsc Ecol. 2020;35(5):1199–213.

    Article  Google Scholar 

  6. Pausas JG, Parr CL. Towards an understanding of the evolutionary role of fire in animals. Evol Ecol. 2018;32(2–3):113–25.

    Article  Google Scholar 

  7. Nimmo DG, Jolly CJ, Carthey AJR. Megafire: the darwinian guillotine? Australian Zoologist. 2022;42(2):217–22.

    Article  Google Scholar 

  8. Hammond TT, Palme R, Lacey EA. Ecological specialization, variability in activity patterns and response to environmental change. Biol Lett. 2018;14(6):20180115.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Riddell EA, Odom JP, Damm JD, Sears MW. Plasticity reveals hidden resistance to extinction under climate change in the global hotspot of salamander diversity. Sci Adv. 2018;4(7):eaar5471.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Schell CJ, Young JK, Lonsdorf EV, Santymire RM, Mateo JM. Parental habituation to human disturbance over time reduces fear of humans in coyote offspring. Ecol Evol. 2018;8(24):12965–80.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Suraci JP, Gaynor KM, Allen ML, Alexander P, Brashares JS, Cendejas-Zarelli S, et al. Disturbance type and species life history predict mammal responses to humans. Glob Change Biol. 2021;27(16):3718–31.

    Article  CAS  Google Scholar 

  12. Gaynor KM, Hojnowski CE, Carter NH, Brashares JS. The influence of human disturbance on wildlife nocturnality. Science. 2018;360(6394):1232–5.

    Article  CAS  PubMed  Google Scholar 

  13. Nimmo DG, Avitabile S, Banks SC, Bliege Bird R, Callister K, Clarke MF, et al. Animal movements in fire-prone landscapes. Biol Rev. 2019;94(3):981–98.

    Article  PubMed  Google Scholar 

  14. Keeley JE. Fire intensity, fire severity and burn severity: a brief review and suggested usage. Int J Wildland Fire. 2009;18(1):116.

    Article  Google Scholar 

  15. Steel ZL, Fogg AM, Burnett R, Roberts LJ, Safford HD. When bigger isn’t better—Implications of large high-severity wildfire patches for avian diversity and community composition. Archibald S, editor. Divers Distrib. 2021;ddi.13281.

  16. Jolly CJ, Dickman CR, Doherty TS, Eeden LM, Geary WL, Legge SM et al. Animal mortality during fire. Global Change Biology. 2022;gcb.16044.

  17. Doherty TS, Geary WL, Jolly CJ, Macdonald KJ, Miritis V, Watchorn DJ et al. Fire as a driver and mediator of predator–prey interactions. Biol Rev. 2022;20.

  18. Blakey RV, Sikich JA, Blumstein DT, Riley SPD. Mountain lions avoid burned areas and increase risky behavior after wildfire in a fragmented urban landscape. Curr Biol. 2022;32(21):4762–e47685.

    Article  CAS  PubMed  Google Scholar 

  19. Geary WL, Doherty TS, Nimmo DG, Tulloch AIT, Ritchie EG. Predator responses to fire: a global systematic review and meta-analysis. Barrios‐Garcia MN, editor. J Anim Ecol. 2020;89(4):955–71.

    Article  PubMed  Google Scholar 

  20. Rickbeil GJM, Hermosilla T, Coops NC, White JC, Wulder MA. Barren-ground caribou (Rangifer tarandus groenlandicus) behaviour after recent fire events; integrating caribou telemetry data with Landsat fire detection techniques. Glob Change Biol. 2017;23(3):1036–47.

    Article  Google Scholar 

  21. Perry DA, Hessburg PF, Skinner CN, Spies TA, Stephens SL, Taylor AH, et al. The ecology of mixed severity fire regimes in Washington, Oregon, and Northern California. For Ecol Manag. 2011;262(5):703–17.

    Article  Google Scholar 

  22. Halofsky JE, Donato DC, Hibbs DE, Campbell JL, Cannon MD, Fontaine JB, et al. Mixed-severity fire regimes: lessons and hypotheses from the Klamath-Siskiyou Ecoregion. Ecosphere. 2011;2(4):art40.

    Article  Google Scholar 

  23. Schriver M, Sherriff RL, Varner JM, Quinn-Davidson L, Valachovic Y, Age. and stand structure of oak woodlands along a gradient of conifer encroachment in northwestern California. Ecosphere [Internet]. 2018 Oct [cited 2023 Jun 22];9(10). https://onlinelibrary.wiley.com/doi/https://doi.org/10.1002/ecs2.2446.

  24. Keeley JE, Fotheringham CJ, Baer-Keeley M, DEMOGRAPHIC PATTERNS OF, POSTFIRE REGENERATION IN MEDITERRANEAN-CLIMATE SHRUBLANDS OF CALIFORNIA. Ecol Monogr. 2006;76(2):235–55.

    Article  Google Scholar 

  25. Roerick TM, Cain JW, Gedir JV. Forest restoration, wildfire, and habitat selection by female mule deer. For Ecol Manag. 2019;447:169–79.

    Article  Google Scholar 

  26. Cherry MJ, Chandler RB, Garrison EP, Crawford DA, Kelly BD, Shindle DB, et al. Wildfire affects space use and movement of white-tailed deer in a tropical pyric landscape. For Ecol Manag. 2018;409:161–9.

    Article  Google Scholar 

  27. Allred BW, Fuhlendorf SD, Engle DM, Elmore RD. Ungulate preference for burned patches reveals strength of fire-grazing interaction: strength of Fire-Grazing Interaction. Ecol Evol. 2011;1(2):132–44.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Archibald S, Bond WJ, Stock WD, Fairbanks DHK. SHAPING THE LANDSCAPE: FIRE–GRAZER INTERACTIONS IN AN AFRICAN SAVANNA. Ecol Appl. 2005;15(1):96–109.

    Article  Google Scholar 

  29. Kreling SES, Gaynor KM, McInturff A, Calhoun KL, Brashares JS. Site fidelity and behavioral plasticity regulate an ungulate’s response to extreme disturbance. Ecol Evol. 2021;11(22):15683–94.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Proffitt KM, DeVoe J, Barker K, Durham R, Hayes T, Hebblewhite M, et al. A century of changing fire management alters ungulate forage in a wildfire-dominated landscape. Forestry: Int J for Res. 2019;92(5):523–37.

    Article  Google Scholar 

  31. Monteith KL, Bleich VC, Stephenson TR, Pierce BM, Conner MM, Kie JG, et al. Life-history characteristics of mule deer: effects of nutrition in a variable environment: Nutrition and Life history of mule deer. Wild Mon. 2014;186(1):1–62.

    Article  Google Scholar 

  32. Buskirk JV. Behavioural plasticity and environmental change. In: Candolin U, Wong BBM, editors. Behavioural Responses to a Changing World [Internet]. Oxford University Press; 2012 [cited 2023 Dec 7]. pp. 145–58. https://academic.oup.com/book/9434/chapter/156303289.

  33. Peters NM, Beale CM, Bracebridge C, Mgumba MP, Kendall CJ. Combining models for animal tracking: defining behavioural states to understand space use for conservation. J Biogeogr. 2022;49(11):2016–27.

    Article  Google Scholar 

  34. Chimienti M, Beest FM, Beumer LT, Desforges J, Hansen LH, Stelvig M et al. Quantifying behavior and life-history events of an Arctic ungulate from year‐long continuous accelerometer data. Ecosphere [Internet]. 2021 Jun [cited 2023 Feb 8];12(6). https://onlinelibrary.wiley.com/doi/https://doi.org/10.1002/ecs2.3565.

  35. Sawyer H, Merkle JA, Middleton AD, Dwinnell SPH, Monteith KL. Migratory plasticity is not ubiquitous among large herbivores. Chapman J, editor. J Anim Ecol. 2018;1365-2656.12926.

  36. Ganz TR, DeVivo MT, Kertson BN, Roussin T, Satterfield L, Wirsing AJ, et al. Interactive effects of wildfires, season and predator activity shape mule deer movements. J Anim Ecol. 2022;91(11):2273–88.

    Article  PubMed  Google Scholar 

  37. Gomes WR, UC Research and Extension Centers. Statewide system provides local answers to local needs. Calif Agr. 2001;55(6):3–6.

    Article  Google Scholar 

  38. Gaynor KM, McInturff A, Brashares JS. Contrasting patterns of risk from human and non-human predators shape temporal activity of prey. J Anim Ecol. 2022;91(1):46–60.

    Article  PubMed  Google Scholar 

  39. Wittmer HU, Forrester TD, Allen ML, Marescot L, Casady DS. Black-tailed Deer Population Assessment In The Mendocino National Forest, California, 2009–2013. 2014 [cited 2023 Feb 16]; http://rgdoi.net/https://doi.org/10.13140/2.1.3440.9764.

  40. Syphard AD, Keeley JE. Mapping fire regime ecoregions in California. Int J Wildland Fire. 2020;29(7):595.

    Article  Google Scholar 

  41. van Wagtendonk JW, Sugihara NG, Fites-Kaufman J, Shaffer KE, Thode AE. Fire in California’s ecosystems. Berkeley, CA: University of California Press; 2018.

    Book  Google Scholar 

  42. Koenig WD, Mumme RL, Carmen WJ, Stanback MT. Acorn Production by oaks in Central Coastal California: variation within and among years. Ecology. 1994;75(1):99–109.

    Article  Google Scholar 

  43. Barker KJ, Mitchell MS, Proffitt KM. Native forage mediates influence of irrigated agriculture on migratory behaviour of elk. Hoye B, editor. Journal of Animal Ecology. 2019;88(7):1100–10.

  44. Calenge C. The package adehabitat for the R software: a tool for the analysis of space and habitat use by animals. Ecol Model. 2006;197(3–4):516–9.

    Article  Google Scholar 

  45. R Core Team. R: A Language and Environment for Statistical Computing [Internet]. Vienna, Austria: R Foundation for Statistical Computing. 2021. https://www.R-project.org/.

  46. Eidenshink J, Schwind B, Brewer K, Zhu ZL, Quayle B, Howard S. A Project for Monitoring trends in burn severity. fire ecol. 2007;3(1):3–21.

    Article  Google Scholar 

  47. Sentinel Hub. Sinergise Ltd. 2022. https://www.sentinel-hub.com.

  48. Gorelick N, Hancher M, Dixon M, Ilyushchenko S, Thau D, Moore R. Google Earth Engine: planetary-scale geospatial analysis for everyone. Remote Sens Environ. 2017;202:18–27.

    Article  Google Scholar 

  49. Dellinger JA, Cristescu B, Ewanyk J, Gammons DJ, Garcelon D, Johnston P, et al. Using Mountain Lion Habitat Selection in Management. Jour Wild Mgmt. 2020;84(2):359–71.

    Article  Google Scholar 

  50. Moll RJ, Redilla KM, Mudumba T, Muneza AB, Gray SM, Abade L et al. L Prugh editor 2017 The many faces of fear: a synthesis of the methodological variation in characterizing predation risk. J Anim Ecol 86 4 749–65.

    Article  PubMed  Google Scholar 

  51. Gaynor KM, Brown JS, Middleton AD, Power ME, Brashares JS. Landscapes of fear: spatial patterns of risk perception and response. Trends Ecol Evol. 2019;34(4):355–68.

    Article  PubMed  Google Scholar 

  52. Zeller KA, Jennings MK, Vickers TW, Ernest HB, Cushman SA, Boyce WM. Are all data types and connectivity models created equal? Validating common connectivity approaches with dispersal data. Bolliger J, editor. Diversity and Distributions. 2018;24(7):868–79.

  53. Bates D, Mächler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models using lme4 [Internet]. arXiv; 2014 [cited 2023 Feb 16]. http://arxiv.org/abs/1406.5823.

  54. Lüdecke D, Ben-Shachar M, Patil I, Waggoner P, Makowski D. Performance: an R Package for Assessment, comparison and testing of statistical models. JOSS. 2021;6(60):3139.

    Article  Google Scholar 

  55. McClintock BT, Langrock R, Gimenez O, Cam E, Borchers DL, Glennie R et al. T Coulson editor 2020 Uncovering ecological state dynamics with hidden Markov models. Ecol Lett 23 12 1878–903.

  56. Clontz LM, Pepin KM, VerCauteren KC, Beasley JC. Behavioral state resource selection in invasive wild pigs in the Southeastern United States. Sci Rep. 2021;11(1):6924.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Michelot T, Langrock R, Patterson TA. moveHMM: an R package for the statistical modelling of animal movement data using hidden Markov models. McInerny G, editor. Methods Ecol Evol. 2016;7(11):1308–15.

  58. Michelot T, Langrock R, Bestley S, Jonsen ID, Photopoulou T, Patterson TA. Estimation and simulation of foraging trips in land-based marine predators. Ecology. 2017;98(7):1932–44.

    Article  PubMed  Google Scholar 

  59. Langrock R, King R, Matthiopoulos J, Thomas L, Fortin D, Morales JM. Flexible and practical modeling of animal telemetry data: hidden Markov models and extensions. Ecology. 2012;93(11):2336–42.

    Article  PubMed  Google Scholar 

  60. Farhadinia MS, Michelot T, Johnson PJ, Hunter LTB, Macdonald DW. Understanding decision making in a food-caching predator using hidden Markov models. Mov Ecol. 2020;8(1):9.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Patterson TA, Basson M, Bravington MV, Gunn JS. Classifying movement behaviour in relation to environmental conditions using hidden Markov models. J Anim Ecol. 2009;78(6):1113–23.

    Article  PubMed  Google Scholar 

  62. Forrester TD, Casady DS, Wittmer HU. Home sweet home: fitness consequences of site familiarity in female black-tailed deer. Behav Ecol Sociobiol. 2015;69(4):603–12.

    Article  Google Scholar 

  63. Calhoun KL, Goldstein BR, Gaynor KM, McInturff A, Solorio L, Brashares JS. Mammalian resistance to megafire in western U.S. woodland savannas. Ecosphere. 2023;14(7):e4613.

    Article  Google Scholar 

  64. Funk KA, Koenig WD, Knops JMH. Fire effects on acorn production are consistent with the stored resource hypothesis for masting behavior. Can J Res. 2016;46(1):20–4.

    Article  Google Scholar 

  65. Raynor EJ, Joern A, Briggs JM. Bison foraging responds to fire frequency in nutritionally heterogeneous grassland. Ecology. 2015;96(6):1586–97.

    Article  Google Scholar 

  66. Gureja N, Owen-Smith N. Comparative use of burnt grassland by rare antelope species in a lowveld game ranch. South Afr. 2002;32(1):8.

    Google Scholar 

  67. Lecomte X, Caldeira MC, Catry FX, Fernandes PM, Jackson RB, Bugalho MN. Ungulates mediate trade-offs between carbon storage and wildfire hazard in Mediterranean oak woodlands. Mukul S, editor. Journal of Applied Ecology. 2019;56(3):699–710.

  68. Van Moorter B, Gaillard JM, McLoughlin PD, Delorme D, Klein F, Boyce MS. Maternal and individual effects in selection of bed sites and their consequences for fawn survival at different spatial scales. Oecologia. 2009;159(3):669–78.

    Article  PubMed  Google Scholar 

  69. Parker KL. Effects of heat, cold, and rain on coastal black-tailed deer. Can J Zool. 1988;66(11):2475–83.

    Article  Google Scholar 

  70. Bose S, Forrester TD, Casady DS, Wittmer HU. Effect of activity states on habitat selection by black-tailed deer: Activity States and Habitat Selection. Jour Wild Mgmt. 2018;82(8):1711–24.

    Article  Google Scholar 

  71. Romme WH, Boyce MS, Gresswell R, Merrill EH, Minshall GW, Whitlock C, et al. Twenty years after the 1988 Yellowstone fires: lessons about disturbance and ecosystems. Ecosystems. 2011;14(7):1196–215.

    Article  Google Scholar 

  72. Morrison TA, Merkle JA, Hopcraft JGC, Aikens EO, Beck JL, Boone RB et al. Drivers of site fidelity in ungulates. Loison A, editor. J Anim Ecol. 2021;90(4):955–66.

  73. Pocknee CA, Legge SM, McDonald J, Fisher DO. Modeling mammal response to fire based on species’ traits. Conserv Biol. 2023;cobi.14062.

  74. Batista EKL, Figueira JEC, Solar RRC, De Azevedo CS, Beirão MV, Berlinck CN, et al. In case of fire, escape or die: a trait-based Approach for identifying Animal Species threatened by fire. Fire. 2023;6(6):242.

    Article  Google Scholar 

  75. Morse D. Behavioral mechanisms in ecology. Cambridge, MA: Harvard University Press; 1982.

    Google Scholar 

  76. Komers PE. Behavioural plasticity in variable environments. Can J Zool. 1997;75(2):161–9.

    Article  Google Scholar 

  77. Buchholz R, Banusiewicz JD, Burgess S, Crocker-Buta S, Eveland L, Fuller L. Behavioural research priorities for the study of animal response to climate change. Anim Behav. 2019;150:127–37.

    Article  Google Scholar 

  78. Harmon JP, Barton BT. On their best behavior: how animal behavior can help determine the combined effects of species interactions and climate change. Ann N Y Acad Sci. 2013;1297(1):139–47.

    Article  PubMed  Google Scholar 

  79. Abernathy HN, Chandler RB, Crawford DA, Garrison EP, Conner LM, Miller KV, et al. Behavioral responses to ecological disturbances influence predation risk for a capital breeder. Landsc Ecol. 2022;37(1):233–48.

    Article  Google Scholar 

  80. Pierce BM, Bowyer RT, Bleich VC, HABITAT SELECTION BY MULE. DEER: FORAGE BENEFITS OR RISK OF PREDATION? Krausman, editor. J Wildl Manage. 2004;68(3):533–41.

    Article  Google Scholar 

  81. Sangha KK, Evans J, Edwards A, Russell-Smith J, Fisher R, Yates C, et al. Assessing the value of ecosystem services delivered by prescribed fire management in Australian tropical savannas. Ecosyst Serv. 2021;51:101343.

    Article  Google Scholar 

  82. Connor T, Tripp E, Tripp B, Saxon BJ, Camarena J, Donahue A, et al. Karuk ecological fire management practices promote elk habitat in northern California. J Appl Ecol. 2022;59(7):1874–83.

    Article  Google Scholar 

  83. Coppoletta M, Merriam KE, Collins BM. Post-fire vegetation and fuel development influences fire severity patterns in reburns. Ecol Appl. 2016;26(3):686–99.

    Article  PubMed  Google Scholar 

  84. Stephens SL, Burrows N, Buyantuyev A, Gray RW, Keane RE, Kubian R, et al. Temperate and boreal forest mega-fires: characteristics and challenges. Front Ecol Environ. 2014;12(2):115–22.

    Article  Google Scholar 

  85. Jones GM, Tingley MW. Pyrodiversity and biodiversity: A history, synthesis, and outlook. Divers Distrib. 2021;ddi.13280.

  86. Abella SR, Fornwalt PJ. Ten years of vegetation assembly after a north American mega fire. Glob Change Biol. 2015;21(2):789–802.

    Article  Google Scholar 

  87. Nimmo DG, Carthey AJR, Jolly CJ, Blumstein DT. Welcome to the pyrocene: animal survival in the age of megafire. Glob Change Biol. 2021;27(22):5684–93.

    Article  CAS  Google Scholar 

  88. Schnurr JL, Ostfeld RS, Canham CD. Direct and indirect effects of masting on rodent populations and tree seed survival. Oikos. 2002;96(3):402–10.

    Article  Google Scholar 

  89. McShea WJ, Schwede G. Variable Acorn crops: responses of White-tailed deer and other mast consumers. J Mammal. 1993;74(4):999–1006.

    Article  Google Scholar 

  90. Ensbey M, Legge S, Jolly CJ, Garnett ST, Gallagher RV, Lintermans M, et al. Animal population decline and recovery after severe fire: relating ecological and life history traits with expert estimates of population impacts from the Australian 2019-20 megafires. Biol Conserv. 2023;283:110021.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to thank Janelle Dorcy, Alison Smith, Troy McWilliams, and the rest of the Hopland Research and Extension Center staff for their support in making the fieldwork component of this study possible. The California Department of Fish and Wildlife and the University of California—Berkeley Animal Care and Use Committee approved all capture procedures. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Funding

Deployment of GPS collars to deer was funded and supported by the California Department of Fish and Wildlife (CDFW) through permit #P1680002. Kendall L. Calhoun, Kaitlyn M. Gaynor, Amy Van Scoyoc, and Alex McInturff were supported by the NSF Graduate Research Fellowship Program (GRFP) during various stages of the study.

Author information

Authors and Affiliations

Authors

Contributions

KC contributed to collecting field data, conceptualization of the study, analyses, and writing of the manuscript. TC contributed to the conceptualization of the study, performing analyses, and writing the manuscript. KG contributed to collecting field data, analyses, and writing of the manuscript. AVS contributed to collecting field data and writing of the manuscript. AM contributed to collecting field data and writing of the manuscript. SK contributed to collecting field data and writing of the manuscript. JB contributed to colleting field data and conceptualization of the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kendall L. Calhoun.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

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/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Calhoun, K.L., Connor, T., Gaynor, K.M. et al. Movement behavior in a dominant ungulate underlies successful adjustment to a rapidly changing landscape following megafire. Mov Ecol 12, 53 (2024). https://doi.org/10.1186/s40462-024-00488-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40462-024-00488-4

Keywords