- Open Access
Quantifying effects of snow depth on caribou winter range selection and movement in Arctic Alaska
Movement Ecology volume 9, Article number: 48 (2021)
Caribou and reindeer across the Arctic spend more than two thirds of their lives moving in snow. Yet snow-specific mechanisms driving their winter ecology and potentially influencing herd health and movement patterns are not well known. Integrative research coupling snow and wildlife sciences using observations, models, and wildlife tracking technologies can help fill this knowledge void.
Here, we quantified the effects of snow depth on caribou winter range selection and movement. We used location data of Central Arctic Herd (CAH) caribou in Arctic Alaska collected from 2014 to 2020 and spatially distributed and temporally evolving snow depth data produced by SnowModel. These landscape-scale (90 m), daily snow depth data reproduced the observed spatial snow-depth variability across typical areal extents occupied by a wintering caribou during a 24-h period.
We found that fall snow depths encountered by the herd north of the Brooks Range exerted a strong influence on selection of two distinct winter range locations. In winters with relatively shallow fall snow depth (2016/17, 2018/19, and 2019/20), the majority of the CAH wintered on the tundra north of the Brooks Range mountains. In contrast, during the winters with relatively deep fall snow depth (2014/15, 2015/16, and 2017/18), the majority of the CAH caribou wintered in the mountainous boreal forest south of the Brooks Range. Long-term (19 winters; 2001–2020) monitoring of CAH caribou winter distributions confirmed this relationship. Additionally, snow depth affected movement and selection differently within these two habitats: in the mountainous boreal forest, caribou avoided areas with deeper snow, but when on the tundra, snow depth did not trigger significant deep-snow avoidance. In both wintering habitats, CAH caribou selected areas with higher lichen abundance, and they moved significantly slower when encountering deeper snow.
In general, our findings indicate that regional-scale selection of winter range is influenced by snow depth at or prior to fall migration. During winter, daily decision-making within the winter range is driven largely by snow depth. This integrative approach of coupling snow and wildlife observations with snow-evolution and caribou-movement modeling to quantify the multi-facetted effects of snow on wildlife ecology is applicable to caribou and reindeer herds throughout the Arctic.
Caribou as a snow-adapted species
Snow dominates Arctic landscapes for more than two thirds of the year , and therefore, snow, in its multiple forms, affects Rangifer species residing in northern regions globally. Caribou (Rangifer tarandus ssp.) are chionophiles; their morphology, physiology, and behavior evolved to enable survival in snowy winter environments [2,3,4,5]. Relatively long legs and low “foot loadings” are unique adaptations that help Rangifer cope with snow by providing efficient locomotion and establishment of feeding craters in snowy conditions . Skin that is well-insulated by woolly fur overlaid with hollow guard hairs  provide extreme insulation and protection from wind [8, 9]. In addition to their sharp edges, caribou hooves also have microstructures of high skid and abrasion resistance, which allow caribou to travel on icy surfaces [10,11,12,13]. These traits enable caribou to remain mobile year-round, even during winters with extremely deep snow . Consequently, their existence in the Arctic has been sustained for more than 50,000 years [5, 15, 16].
Despite the specialization of caribou to snowy conditions, observational research indicates that snow significantly affects caribou winter ecology . For example, snow characteristics determine forage accessibility (e.g., [3, 4, 18]), energy expenditure in locomotion , daily activity budget (i.e., the amount of time spent cratering, resting, and moving ), and migration timing (e.g., ). The energetic costs and benefits related to these processes and activities likely influence where caribou spend the winter and their movement paths through snow-covered landscapes. However, the degree of control that snow exerts on caribou movement patterns and rates across different spatial scales throughout winter, and how that may vary among herds and among individuals, has received limited attention. Pruitt  recognized the large spatial and temporal variations in snow characteristics across caribou winter ranges. He argued that for the purpose of explaining caribou winter-range conditions, “...conventional snow data as reported by the existing net of meteorological observatories are not only insufficient but may actually be misleading”. Past studies have typically relied on snow field observations or snow remote sensing datasets that can be incompatible with the dynamic nature of wildlife movement data in terms of temporal frequency, spatial resolution and coverage, or the snow variables being represented; often, these datasets do not characterize the most ecologically meaningful snow properties at the most appropriate scales for wildlife applications . While Pruitt advocated for the importance of incorporating spatiotemporal snow variability in caribou studies over 60 years ago, more recently Boelman et al.  argued it remains a persistent issue and data-gap in wildlife ecology. Here, we present snow datasets that can be integrated with spatiotemporally dynamic caribou movement data. The winter environmental modeling tools we used, SnowModel (presented in section 2.2.2; ), are capable of producing snow information that captures caribou- and other wildlife-relevant snow properties and their evolution across space and time; here, we utilized snow depth. With such information, it is now possible to more deeply investigate and quantify the role of snow in caribou winter ecology .
Fall migration cues
Most barren-ground caribou conduct fall and spring migrations that bracket the Arctic winter; these can represent the longest-distance terrestrial migrations on the planet . For migratory species, the timing and destination of these long, directed movements can be correlated to environmental cues such as vegetation phenology, photoperiod, previously-used migration trails, snow accumulation and melt, or temperature (e.g., [19, 20, 24,25,26]). An animal perceives cues from its surrounding environment that are indicative of conditions at different locations and/or in the future, which may be suitable for individuals to strategically carry out life-cycle events, such as migration, with minimal risk and maximum benefit [27, 28]. The perception of these cues affects an animal’s decision to initiate preparatory steps for migration, and hence, affects the seasonal timing of this life-cycle event .
Snow conditions can be a significant trigger for migratory species inhabiting seasonally variable environments (e.g., [29,30,31]). For instance, the first heavy snowfall after rut can trigger fall migration of caribou [32, 33]. In addition to initiating migration, fall environmental cues may influence selection of winter range location. Based on mostly anecdotal evidence, this appears to be the case among Arctic Alaska caribou herds, for which the first snowfall likely acts as a stimulus for fall migration [34, 35], and the cues for selection of winter range location likely include thresholds of snow depth, timing of first snowfall event, and olfactory sensing of lichen below the snow . In the past, snow data at adequate spatial and temporal resolutions and extents, to capture changes in snow conditions that may trigger such behaviors, have been non-existent or inadequate  and, therefore, quantification of the cueing effect of snow has been limited. Such spatially and temporally explicit snow data, describing, e.g., snow accumulation, need to: 1) cover the entire migration route, including initiation location, migration path, stopover sites, and finally, winter range; 2) have adequate temporal frequency and detail to match the temporal scale at which a definitive change in caribou movement behavior is detected (e.g., daily); and 3) contain spatial and temporal information at sufficient detail to realistically capture the variability in snow accumulation within the initiation location of each individual caribou’s migration. Until now, such dynamic, spatially and temporally comprehensive snow information has been largely unavailable for wildlife applications .
Winter range location and winter movement
At a regional scale, caribou wintering locations can vary considerably from year to year and include tundra and boreal forest habitats [19, 37,38,39]. For many herds across North America, their precise annual winter range location may depend largely on snow conditions and snow’s control on forage accessibility [3, 11, 40, 41]; winter forage distribution, particularly of lichen species [42,43,44,45,46]; fire history [42, 47]; and coincident predator abundance [33, 48].
At a landscape scale within their winter range, field observations suggest that snow depths tolerated by caribou differ between tundra and boreal forest habitats [4, 49]. However, in general, caribou select for forage areas with shallow and/or soft snow cover at both the landscape and patch scales (i.e., within tens of meters; [3, 4, 18, 36, 50, 51]). In tundra habitats, this selection typically results in winter movement patterns where caribou avoid depressions in the landscape, e.g., cut banks of rivers that are more likely full of deep, wind-compacted snow, and instead stay on snow-free ground, such as rounded hilltops with wind-swept slopes where ablation (sublimation and wind erosion) has significantly reduced the snowpack . Caribou wintering in a boreal forest habitat typically move single-file through the characteristically deep, but soft, snow cover ; crater for forage below the snowpack or feed on arboreal lichen, which is more easily accessed than terrestrial lichen when the snowpack is deep ; and they may rest on lakes or seek out tree-less ridges that enable better detection of incoming predators [3, 48]. Such selection patterns translate into an ideal caribou winter range being characterized by snow conditions that provide easy access to forage, minimize energy costs related to mobility , and possibly enhance predator avoidance [55, 56].
The Central Arctic Caribou Herd
During the past four decades, the barren-ground caribou (Rangifer tarandus granti) of the Central Arctic Herd (CAH) in Arctic Alaska have been studied extensively, particularly in summer because their calving grounds are located in oil and gas extraction areas near the Beaufort Sea coast; (Fig. 1 ; [57,58,59,60,61,62,63,64,65]). The CAH population size was estimated at approximately 6000 animals in 1978 and grew steadily to its peak of 68,000 caribou in 2010 (with a slight decline during the early to mid-1990s), and was estimated to be roughly 30,000 animals in 2019 [60, 61, 66]. From their annual gathering on the coastal calving grounds in summer, CAH caribou generally disperse during early fall (~end of August – mid September) across the region extending from the Beaufort Sea coast, south to the Brooks Range (BR), and by October, they typically migrate further south to spend the winter either north or south of the Continental Divide (CD; running west to east along the crest of the BR; Fig. 1). For nearly two decades, from 2002 through 2020, Alaska Department of Fish and Game (ADFG) has annually monitored the general CAH winter distribution by estimating the proportion (%) of collared caribou located south of the CD in mid-March ([66,67,68,69,70]; Table 1). These historical surveys suggest that during the last 19 years, the majority of collared CAH caribou (≥50%) were distributed south of the CD in mid-March. In four of the 19 years, less than 10% of collared CAH caribou were located south of the CD, thus, in those years, the majority of the herd wintered north of the CD, typically in the northern BR foothills both east and west of the Dalton Highway (Table 1; Fig. 1; [57, 59, 67, 71,72,73]).
The winter ranges south of the CD include areas dominated by coniferous forest and areas of non-forested mountain slopes; this habitat type is hereafter called ‘mountainous boreal forest’. Winter ranges north of the CD generally include the broad, north-facing valleys and the northern foothills of the BR. The tundra habitat there includes wind-blown, rocky ridges, and rolling hills intersected by braided river deltas with shrub tundra and tussock tundra. Particularly during the six most recent winters (2014–2020), the CAH winter range location differed markedly from year-to-year between the tundra in the north and the mountainous boreal forest in the south (Table 1). We expected that an Arctic winter environmental attribute with the capability of changing dramatically from one year to the next influenced this winter-range variability, prompting our focus on evaluating the role of snow in CAH winter range selection. While we acknowledge that other climate variables, forage abundance, predation, disturbance, and neighboring herds’ winter range selection have been known to influence caribou winter range variability and winter movement, the effect of snow depth remains unquantified.
Specifically, we quantified the influence of snow depth on caribou winter range selection over a large region (~ 100s of kilometers) and evaluated the effect of snow depth and winter forage on daily movement at a landscape-scale (~ 100s of meters). Because this mountainous, Arctic Alaska study region is remote, nearly road-less, and highly inaccessible in winter, these investigations required a combination of snow field observations and snow modeling tools capable of representing the physics associated with snowpack evolution and distribution at regional and landscape scales to produce the necessary temporally and spatially explicit snow depth data. At the regional scale, these snow depth data were coupled with long-term CAH winter distribution records to study winter range location selection. At the landscape scale, these snow depth data and winter forage data were coupled with Global Positioning System (GPS) locations in animal movement models to quantify the effect of snow depth on daily caribou selection and movement rate (see animation of daily snow depth data and caribou location data in Additional file 1). We addressed the following questions:
Does snow depth encountered by CAH caribou in the fall correspond to observed difference in winter range location?
Once on winter range, how do snow depth and winter forage (lichen) abundance affect CAH caribou daily resource selection and movement rate?
Do the effects of snow depth on daily resource selection and movement rate differ between the two CAH wintering habitats: tundra north of the CD and mountainous boreal forest south of the CD?
Further, given the increased energy expenditure associated with movement through deeper snow , and the importance of lichen as a winter forage resource for barren-ground caribou [19, 44], we hypothesized that on a daily temporal scale, caribou (1) select for areas with lower snow depth than otherwise available in the surrounding landscape, (2) select for areas with higher lichen proportions than otherwise available in the surrounding landscape, and (3) caribou movement rates decrease with increasing snow depth. Finally, we expected this selection and movement behavior to be more pronounced when caribou winter in mountainous boreal forest than on tundra, because the boreal forest snowpack is typically deeper than the tundra snowpack.
Our study area encompassed all CAH GPS locations used in this study. It included parts of the traditional lands of Iñupiat, Nunamiut, Gwich’in, and Koyukon Athabascans, which is a section of the Brooks Range, the highest mountain range in the Arctic, as well as its foothills and the northern-most coastal plain of Alaska (Fig. 1). Specifically, the study area extended from the Beaufort Sea coast in the north, 460 km south to the southern foothills of the BR, and between the northward-flowing rivers, Killik River in the west, and the Hulahula River, 430 km to the east. Two main habitat types compose this area and were central to this research: tundra and mountainous boreal forest. The tundra vegetation is dominated by wet sedges, tussocks, and dwarf shrubs near the coast [74, 75], and sedges, shrubs, and tussocks in the northern foothills of the BR [76, 77]. The mountain slopes of the BR are dominated by low alpine tundra and barren or sparsely vegetated areas, while the valleys south of the CD are covered by more erect shrubs and coniferous boreal forest dominated by black spruce communities . These two habitats are associated with two distinctly different types of seasonal snow cover: Arctic tundra snow is generally shallow (mean snow depth = 45 ± 21 cm; March – April 2018 and 2019 (this study) ;) with windblown top layers that are likely capable of supporting a caribou without collapsing. In contrast, boreal forest snow is shallow to moderately deep (mean snow depth = 76 ± 17 cm; March – April 2018 and 2019 (this study); [78,79,80]) with a more fragile structure of weakly bonded, large snow grains (depth hoar), making it more likely that caribou sink deeper into the snowpack. The main reasons for these snowpack differences are: (1) a windier climate on the tundra resulting in increased redistribution of snow compared to the boreal forest, where the trees inhibit redistribution by wind; and (2) air temperature gradients across the region, defined by topographic elevation and distance to the coast and the ocean with variable sea-ice coverage [81,82,83,84,85].
To characterize the spatial snow depth variability within these two different snow types and between winters, we conducted six-week snow field campaigns March–April in 2018, 2019, and 2020 at 82 snow sites (Fig. 1; measurement protocols and snow site selection detailed by Pedersen et al. ). These 2018–2020 snow observations provided three winters of field-validation for the 2014–2020 snow depth datasets that were paired with caribou location data (see sections on SnowModel and Caribou data). The fieldwork on the tundra snow was completed by snowmobile traverses. The mountainous boreal forest snow sites were accessed by first landing on snow-covered lakes with a fixed-wing aircraft on skis and then traveling to sites on snowshoes. At each snow site, we measured snow depth at approximately every three meters along 200–800 m transects across the landscape using an automated snow-depth probe (MagnaProbe; ). We dug snow pits to identify snowpack stratigraphy and measured bulk and layer-specific snow density according to protocols detailed by Pedersen et al. . To estimate a bulk snow water equivalent (SWE) at each snow site, we paired the bulk density (weighted by layer thicknesses) from each snow pit profile with the average snow depth measured in the nearby transects. Both site-specific SWE and density were assimilated in SnowModel (described in the section below). The individual snow depth measurements made along the transects were used for comparison with the modeled snow depths.
To investigate the effects of snow depth on caribou winter range selection and winter movement patterns across space and time, we used snow depth data produced by SnowModel (Fig. 2; [22, 87]). This system allows users to generate fit-for-purpose snow information that is relevant for unique investigations of animal movement behaviors over appropriate, application-specific, spatial and temporal resolutions. For this application, we produced daily snow depth data at 90 m by 90 m spatial resolution, over the 430 km by 460 km study area (~ 200,000 km2; Fig. 1), for 1 September 2001–31 May 2020. The most relevant spatial resolution for each unique snow dataset should (1) resolve the physical processes that define the variability and distribution of the snow variable of interest across the landscape (e.g., snow redistribution by wind); (2) resolve the scale of the observed, dynamic movements of the study animal; and (3) balance (1) and (2) with the available computational resources while producing manageable output file sizes. We selected a 90-m spatial resolution because it resolved the variability in snow depth encountered over the approximate observed distance (~ 2 km) traveled by our wintering GPS-collared caribou in one day (i.e., our temporal resolution). The 90-m data resolved dominant, local, topography-specific snow-depth differences between valley bottoms and hilltops, and the low snow depth levels observed in the Brooks Range mountain valleys . Lastly, these data dimensions were reasonable given our computational processing capabilities. SnowModel is the core of a suite of physically based snowpack-evolution, snow (re)distribution, and process modeling tools (Fig. 2; [22, 87]). Included in SnowModel is SnowAssim (Fig. 2; ), a submodel designed to assimilate field-based snow observations, such as snow density and SWE. SnowAssim works with SnowModel to simulate physics-based snow-landscape interactions and generate spatially and temporally distributed snow information between snow measurements in space and time. Through this assimilation of observed snow properties, the resulting snow information is physically realistic at measured and unmeasured locations and times across a landscape at user-defined spatial and temporal scales (e.g., [88, 90]). Hence, SnowModel provided spatially and temporally continuous distributions of snow properties that were consistent with our observational datasets.
A high-resolution meteorological model called MicroMet  provides meteorological data inputs to the SnowModel process models (Fig. 2). MicroMet spatially distributes meteorological information from local weather stations (in this application we used 15 weather stations; Fig. 1; [92, 93]) and gridded atmospheric reanalysis data (MERRA-2; ) over the study area. MicroMet requires data for five variables: air temperature, relative humidity, precipitation, wind speed, and wind direction.
SnowModel inputs also included two temporally invariant environmental data layers (Fig. 2): the National Elevation Dataset  was used to describe topography, and North American Land Change Monitoring System (NALCMS; ) provided land-cover data to define the spatial distribution of land-cover types across the region. The topography dataset is used by MicroMet to define temperature and precipitation lapse rates and wind speed acceleration and deceleration across the landscape . This information, in turn, is used to define snow evolution features such as snowfall and snowmelt distributions, and the redistribution (i.e., transport) of snow by wind. The land-cover information defines the variability in vegetation heights used in SnowModel; the heights are used to control whether snow is available for redistribution by wind or is, instead, immobilized by vegetation [97, 98]. We generated our required snow depth datasets by running SnowModel and SnowAssim for a 19-year simulation, while assimilating snow density and SWE observations (together these define snow depth; ) when available across the study area and period. For 2001–2017, we used SWE data from the Imnaviat Creek snow monitoring site in the northern BR foothills (Fig. 1; [100, 101]), and for 2017–2020, we used SWE and density observations from our field expedition/traverse datasets (Fig. 1; ).
Our field snow-property observations  and other historical snow-free date observations  provided evidence for a consistent, relatively thin or non-existent snow cover in the major northward-facing valleys of the BR, north of the CD. During our field campaigns, we observed the greatest number of caribou foraging in these valley areas, and Nicholson et al.  also observed a frequent use of these valleys in earlier years (2003–2007). High wind speeds in these valleys are responsible for increased sublimation and wind-erosion of the snowpack, leaving them relatively snow-free during winter compared to the surrounding landscape . To ensure our models accurately reproduced the observed shallow snow depth and high wind speeds in these valleys, we developed a wind-speed correction mask using annual snow-free date datasets for the region . This adjustment was required because the weather stations used to provide meteorological inputs to SnowModel (Fig. 2) were located too far from these valleys to detect and record these locally increased wind speeds. In the tundra areas covered by the correction mask, the wind speed was increased in SnowModel by a factor of 1 to 7. The largest correction (factor 7) was applied to wind speed in these major north-facing valleys, while wind speed in areas between the valleys received little to no adjustment. This correction served to increase snowpack wind erosion and sublimation in these northward-facing valleys and was required for SnowModel and SnowAssim to reproduce our observed snow depths.
Winter forage data
Lichen constitutes a dominant proportion (60–70%) of barren-ground caribou winter diet (e.g., [19, 44]). Lichen species occur throughout the CAH winter range [46, 104], and we detected lichen types (Cladonia spp. and Cetraria spp.) below the snow in 65% of our snow pits during fieldwork in 2017–2020 (S. H. Pedersen 2020, unpublished data). The patchy distribution of lichen makes mapping of its spatial distribution and biomass challenging. While recent methodologies have been used to successfully map lichen distributions at the landscape scale (~ 30 m; [105,106,107,108]), large portions of the CAH winter range are not yet covered by these lichen abundance datasets. Therefore, we mapped caribou winter forage across our study area using the proportion of lichen-rich vegetation based on the 30-m NALCMS land-cover map . The lichen proportion was calculated as a continuous variable, and defined as the fraction of the nine 30-m by 30-m grid cells (i.e., three grid cells by three grid cells or a 90 m by 90 m area total) classified as erect shrub tundra or alpine tundra (non-tussock) vegetation, which are both defined as containing a minimum of 20% lichen cover . This lichen-proportion dataset produced values that ranged from 0.0 to 1.0; higher lichen proportions (i.e., 0.78–1.0) were more common north of the CD than south. In this study, the proportion of lichen was assumed to be constant through time.
Caribou capture and monitoring
During the months of March, April, or June in years 2001–2019, adult (≥2 years old) female and male caribou were captured by ADFG for collaring using a hand-held net-gun fired from a low-flying helicopter (model R-44, Robinson Helicopter Co., Torrance, California). Caribou were fitted with either Very High Frequency (VHF) or GPS equipped satellite collars (Telonics, Mesa, AZ, USA).
We used caribou GPS location data representing male (11) and female (65) caribou collected over nearly six years (September 2014–May 2020, excluding the months of June, July, and August; Table 2). To obtain equal time intervals (fix rates) throughout the winter season for each individual, we resampled all caribou locations at various fix rates (2, 8, and 12 h) to produce a complete time series from fall to spring using the Animal Movement Tools (amt) R package . The resampling yielded three location datasets of 8-hourly, 12-hourly, and daily location frequency, where 8- and 12-hourly fix rates existed for only a subset of the monitored animals and varied in availability throughout our six study years (Table 2). While daily location data represents the coarsest temporal resolution applicable in many types of movement analyses , we chose to use daily locations (24-h fix rates) for the analyses, because it was the fix rate resulting in the highest number of individuals with adequate data throughout the winter season (Table 2), i.e., providing the maximum number of individuals for our analyses.
Fall snow depth
The CAH caribou typically experience the first snowfall in early fall (August–September) and the bulk of the tundra snowpack is established during the ‘early cold’ season of sustained, below-freezing air temperature, usually late September through early December . The typical early-fall distribution of the CAH is from the Beaufort Sea coast to the northern foothills of the BR (e.g., [67, 71,72,73]). To evaluate whether the CAH use snow depth encountered across this area during early fall as a winter-range selection cue, we used all available GPS locations for the month of September (8-h, 12-h, and 24-h fix rates in years 2014–2019; n = 43,680) to represent this scattered, early fall distribution. Located north of the BR, this distribution extended from 153.9°–146.0° W and from 68.4°–70.2° N (Fig. 1). To generate an annual average fall snow depth, we extracted the daily, modeled snow depth in these GPS locations from 1 September through 31 December and averaged all values per year 2001/02–2019/20. To evaluate the relationship between long-term, observed CAH winter range location (i.e., ADFG mid-March distributions that we assumed approximate their winter range location north or south of the CD; Table 1) and average fall snow depth from 2001/02–2019/20, we used two methods. We first used a simple Pearson’s correlation test to estimate the correlation between the two variables. Secondly, we used a mixed-effects logistic regression with binomial errors  to conduct a deeper examination of the relationship between these two variables and how they may have changed over time. The mixed-effects logistic regression used the total number of individuals detected during annual ADFG surveys that wintered south of CD and the total number of individuals detected during annual surveys that wintered north of CD (Table 1) as the response variable and the annual value of fall snow depth and a linear year trend as continuous explanatory variables. The linear trend in wintering location with year was included to test whether the wintering location changed consistently over the study period. We included random intercepts for each year to account for the interannual variability of the 19 years and to account for overdispersion caused by unexplained variability, likely due to other factors that arguably contribute to the herd’s winter-range selection (see the Discussion section) but are not included in our model. Hence, our mixed-effects logistic regression included fall snow depth and a linear year trend as fixed effects and random intercepts for individual years. We calculated odds ratio (i.e., e, the natural logarithm base, raised to the exponent of the model estimate per explanatory variable, eestimate) to evaluate the effect of the variables included in the model . We used the ‘lme4’ R package for these statistical analyses .
In preparation for our movement analyses, we identified missing data in our CAH caribou location time-series dataset, and placed missing values (NA) at times when a location should have been recorded but was missing. GPS collars are programmed to record a location at a fixed hour, but the location dataset can contain imprecisions, because the actual time of location recording may vary seconds to minutes from that specified hour. Because a constant time lag between successive location measurements was required for our movement analyses, we rounded the times based on a reference date (e.g., with the hour 00:00), so that our locations were exactly 24 h apart. As a consequence of rounding the timestamps, an individual’s position may also need to be corrected; for example, if the time of a location is corrected to occur 30 min later than the actual recording occurred, then the position of the individual also needs a correction, which is made under the assumption that an animal moves forward with a constant speed . We used the functions ‘setNA’ and ‘sett0’ in the R package adehabitatLT to regularize our trajectories .
To objectively identify the wintering period (i.e., determine when animals were on winter range) for each GPS-collared caribou in each year (defined as the date at which each caribou arrived on winter range through the date that it departed winter range), we first calculated the movement rate for each relocation of each animal. We then used the method of Lavielle [116, 117] to partition the animal trajectory into segments characterized by homogenous behavior, in this case, sustained, relatively low movement rates. We identified a maximum of ten segments for each winter, with each segment containing a minimum of ten relocations. Typically, for each individual, this method identified four general segments per year, during the time 1 September–31 May: (1) A period during early fall (prior to fall migration) of intermediate movement rates (~ 5 km per day), where the CAH caribou typically inhabited the coastal plain and the region north of the BR; (2) a short-term period of fall migration with abruptly elevated and maximized movement rates (up to ~ 50 km per day); (3) following the fall migration, a period of minimum, sustained, relatively low winter movement rates (~ 1.8–3.1 km per day; Table 2; this low-movement rate period was typically the longest; ~ 4–7 months, and we selected this longest segment as the winter range period for each individual for the analyses of this study); (4) finally, for individuals with GPS locations available throughout April and/or May, the final segment of elevated movement rates (up to ~ 50 km per day) was identified and defined as spring migration. Six individuals had inadequate data and were excluded from the analyses. The total number of individuals used in the analyses is in Table 2. All preparation of the GPS location data and movement analyses were done in R version 3.6.0 .
Integrated step selection analysis
We used integrated step selection analysis (iSSA)  to test the hypothesis of whether movement patterns of wintering caribou were influenced by snow depth and the proportion of lichen-rich vegetation in tundra or mountainous boreal forest winter habitats, i.e., north and south of the CD. Additionally, elevation was tested as a potential covariate because of the topographic differences between these two habitat types and its likely importance for caribou movement. Step selection functions [110, 120] are widely applied in animal movement ecology to investigate resource selection at the spatial and temporal scale of a movement step, which is defined to be a relocation of an animal between two points: a starting point and endpoint at fixed time intervals. These analyses operate under the assumption that locations closer to an animal are more likely to be used than locations farther away [121, 122]. By fitting a conditional logistic regression model, step selection functions account for the fact that for each observed step/relocation (i.e., used point), there exists alternate locations (i.e., available points), unique to that observed time, that the animal could have selected instead . Therefore, these analyses allow relatively fine-scale assessment of selection as an animal moves through the landscape. We randomly generated ten available points for each used point; the locations of available points were based on each animal’s individual distributions of observed turning angle and step length (the Euclidian distance traveled per day) from their own used steps recorded during each individual winter. Snow-depth and lichen-proportion covariate values were extracted at the endpoint of each used and available step to describe the snow depths and lichen proportions that caribou moved towards, in order to test these covariates' effect on selection. In addition to the characteristics of a step selection function, iSSA allows the user to include movement parameters, such as turning angle and step length, in the regression model, thus producing an estimate of movement and selection parameters in the same model . We included the interaction between snow depth and step length as a covariate in the iSSA to quantify whether snow depth had a significant effect on caribou movement rate (i.e., step length). For this covariate, we used snow depth values that were extracted from our spatially and temporally continuous snow depth dataset at the starting point of each step, rather than the endpoint to evaluate how the initially encountered snow depth influenced subsequent step length, i.e., daily movement rate.
Each of the six winters in this study were categorized according to the primary winter distribution defined by where the majority of the collared caribou (≥50%; Table 1) were wintering. We grouped the three winters 2016/17, 2018/19, and 2019/20, when the CAH primarily wintered on tundra north of the CD and grouped the three winters 2014/15, 2015/16, and 2017/18, when the CAH primarily wintered in mountainous boreal forest south of the CD. We conducted the movement analyses with the caribou location data grouped by this primary winter distribution, and identified each step by individual and winter. Consequently of this habitat-specific (tundra vs. mountainous boreal forest) grouping of the winters, all GPS-collared caribou (Table 2) were included in the movement modeling, i.e., contributed to the population estimates of selection coefficients, despite some of those caribou wintered elsewhere than the habitat designation of a given winter (n = 2 in 2015/16, n = 14 in 2017/18, n = 4 in 2018/19, Table 1). An iSSA was fitted to each individual caribou per winter with the resulting estimates averaged (non-weighted) to produce a population-level estimate per winter range location north and south of the CD (using the R package amt; ). To facilitate easier comparison of effect sizes (i.e., the relative strength of selection) between covariates, the snow depth and lichen proportion values used in the iSSAs were centered by subtracting their mean and were then scaled by dividing the centered values by their standard deviation.
The iSSA results were reported as relative selection strength and probability of use, following methods and recommendations by Avgar et al. . We conducted separate iSSAs using 8-hourly, 12-hourly, and daily location data for the two winter ranges. To avoid collinearity [123, 125] between covariates included in iSSAs, we generated pair-wise Pearson’s correlation coefficients (r) for all variables (snow depth, lichen proportion, and elevation), and considered variables with r ≥ 0.6 to be highly correlated . If two variables were highly correlated, we only retained the most informative variable in our model. We included only snow depth and lichen proportion in the iSSAs; we excluded elevation because it was significantly correlated with snow depth (r ≥ 0.6) in two of the six winters (r = 0.33, 0.65, 0.72, 0.30, 0.31, 0.52 for 2014/15–2019/20, respectively).
Our snow observations collected in 2018, 2019, and 2020 included bulk SWE, snow density, and more than 1500 snow depth measurements at the 82 snow sites (2018 n = 17,743; 2019 n = 10,091; and 2020 n = 1543; Fig. 1). The SnowModel and SnowAssim regional snow depth distributions across our study area described the observed, spatial snow-depth variability among snow site locations (Fig. 3a, b, c; the map shading is modeled snow depth and the filled circles are average snow depth per snow site). The average snow depths ranged from 5 cm to 111 cm (Fig. 4a; the axes scales show the snow-depth range among snow sites and the grey ±1 standard deviation error bars show the variability of observed snow depth within each snow site). The modeled snow depths extracted at transect locations at each snow site were significantly correlated with snow depths observed along those transects for each snow site during the three years, based on Pearson’s correlation coefficients (r = 0.76 (2018); r = 0.80 (2019); r = 0.87 (2020); Fig. 4a). Generally, the range in modeled snow depth captured the observed range with a mean residual of 0.3 cm between mean modeled and observed snow depth (Fig. 4b).
Fall snow depth and winter range selection
In early fall, the CAH is typically distributed north of the CD. The timing of the first snow accumulation in this area coincided with the timing of CAH caribou selection of winter range location (Fig. 5; animation in Additional file 1). The first day on the winter range (i.e., when exhibiting wintering behavior of relatively low, daily movement rate) for GPS-collared individuals of this study was between 7 October and 7 December. Because a substantial portion of the snow cover is typically formed during the months September–December (Fig. 5), the snow depths that caribou encounter at that time are generally a consistent predictor of the snow depths they will experience during the subsequent winter months (r = 0.88, p < 0.001). Though the fall snowpack evolution does vary across our large study region and between years [85, 127], a potential explanation for this predictive quality of the fall snow depth is that often a wind event during the fall months creates a wind slab on top of the snowpack. This compact top layer/crust can act as a lid and prevent the underlying snowpack layers from being eroded by wind transport during the remainder of the winter, i.e., the snowpack will remain at an approximately similar depth after this fall wind event .
The CAH distribution and movement varied throughout the winters 2014/15–2019/20 (Fig. 5). In September, GPS-collared CAH caribou were distributed from the Beaufort Sea coast south to the BR (Fig. 5); these early-fall distributions were similar across the six years. The fall migration, identified as more directed movements, was initiated in October during most years. The GPS data indicated that by November and December, caribou reduced their movement, i.e., exhibited wintering behavior, and remained in a final winter range location (Fig. 5). We evaluated whether fall snow depth (September–December average) in the area north of the CD (marked by red circles in Fig. 5), acted as a cue for caribou to either move south through the BR to winter in mountainous boreal forest, or stay on the north side of the CD to winter primarily on tundra. Two winters that clearly exemplified both cases were 2014/15 and 2018/19. The deeper average snow depths (~ 20–35 cm) north of the CD in October–November 2014 were followed by movement to winter ranges primarily south of the CD, while shallower average snow depths (~ 0–15 cm) in October–November 2018 were followed by use of winter ranges north of the CD. After arrival on the winter range, most GPS-collared caribou remained in that general location until March, or when they initiated spring migration (Fig. 5).
The monthly and annual distribution of caribou provides a visualization of the ways that snow depth may play a role in defining CAH caribou selection of winter range (Fig. 5). During years with relatively shallow fall snow depths north of the CD, which typically equated to lower-snow winters overall, the majority of the CAH stayed north of the CD (e.g., winters 2016/17, 2018/19, and 2019/20; Fig. 6). Conversely, deeper fall snow depths north of the CD, likely signaling a snow-rich winter, corresponded with the majority of CAH caribou wintering south of the CD (as was observed in 2014/15, 2015/16, and 2017/18; Fig. 6). During fall seasons with intermediate snow depths, a large proportion, but not all, of the collared CAH animals wintered south of the CD, as they did during years with deeper fall snow depths (Fig. 6). The modeled snow depths were used to quantify this relationship. Our results indicated that fall snow depths correlate strongly with CAH winter range locations during 2001/02–2019/20 (Fig. 6; Table 1), yielding a significant positive correlation between winter range location and fall snow depth (r = 0.59, p = 0.008, n = 19). The modeled estimates of fall snow depth and linear trend in wintering location with year in the mixed-effects logistic regression model (including random intercepts per year) were both statistically significant. These results confirmed the Pearson’s correlation test result; the modeled estimate of change in wintering location with fall snow depth (estimate 0.39; 0.22 and 0.58, lower and upper 97.5% confidence interval [CI], respectively) was significantly different from zero. This result indicates that proportionately more CAH caribou winter south of the CD in years with higher mean fall snow depth. In other words, a single unit (= 1 cm) increase in fall snow depth increases the odds of wintering south of the CD by a factor of 1.5 (1.2 and 1.8, lower and upper odds ratio 97.5% CI, respectively), when holding all other variables constant. Furthermore, the proportion of caribou wintering south of the CD declined over time. However, within this annual trend in winter distribution, fall snow depth exerted a strong influence on the winter range selection. Based on examination of the residuals, the annual pattern in winter distribution was largely driven by an abrupt change in later years; caribou were more likely to winter south of the CD in years before 2015/16, than in the years after 2015/16.
Snow depth effects on movement rate and selection
In addition to evaluating the relationship between fall snow depths and CAH caribou winter range location, we tested the hypotheses of whether snow depth affects movement rates and selection once the caribou were on winter range. We extracted daily, in situ snow depths and lichen proportion associated with caribou GPS locations (Fig. 5) for the 2014/15–2019/20 winters. The annual average snow depths at the endpoints of both available and used steps ranged from 22 ± 16 cm to 43 ± 29 cm (mean ± SD) in winters primarily spent north of CD (2016/17, 2018/19, and 2019/20) and from 35 ± 15 cm to 64 ± 28 cm in winters primarily spent south of CD (2014/15, 2015/16, and 2017/18).
Our iSSA results revealed that snow depth affected CAH caribou movement and selection on the winter range, but the quantified relationship differed between tundra and mountainous boreal forest habitats. During winters when the CAH primarily remained north of the CD on tundra, at the scale of our analyses (90-m and daily data), snow depth was not a statistically significant driver of caribou selection (Fig. 7a). Therefore, caribou did not display positive or negative selection for snow depth, indicating neither selection for nor avoidance of deeper snow. In contrast, in years when the CAH primarily wintered in the mountainous boreal forest south of the CD, caribou showed an increased probability of use with decreased snow depth (Fig. 7c) indicative of negative selection for snow depth, i.e., avoidance of areas with deeper snow (Fig. 7a).
To account for the role of forage in caribou winter selection patterns, we also included lichen proportion in the iSSA models. The presence of lichen-rich habitat (i.e., lichen proportion = 1.0) was higher north than south of the CD; in endpoints of available and used steps, high lichen proportions occurred more frequently for winters primarily spent north of CD (2016/17, 2018/19, and 2019/20) than south (2014/15, 2015/16, and 2017/18). During all six winters, regardless of winter range habitat type, GPS-collared caribou selected for areas with higher lichen proportions, as indicated by a statistically significant positive selection coefficient (Fig. 7a). Probability of use increased in areas with a higher lichen proportion in both the tundra and mountainous boreal forest environments (Fig. 7b). We found no statistically significant correlation between snow depth and lichen proportion for the six winters (r = − 0.047, − 0.24, − 0.37, − 0.18, − 0.37, − 0.095 for 2014/15–2019/20, respectively). Hence, at this scale of analysis and using our lichen proportion dataset, the lichen abundance was not significantly higher in areas with shallow snow. These results suggest that caribou movement rate and selection within the winter range is driven independently by snow depth and lichen proportion.
Finally, iSSA results suggested that on winter ranges both on tundra and in mountainous boreal forest, caribou movement rates were significantly and inversely proportional to the snow depths encountered at the starting point of each daily step (Fig. 7a). Hence, caribou movement rates were significantly lower when animals experienced greater initial snow depths. These results (based on daily location data from 2014 to 2020) were consistent with the iSSA results from separate analyses using 8-hourly and 12-hourly location data. Furthermore, the average movement rates were significantly different and higher for caribou wintering primarily north (2.6 km per day) than south (2.1 km per day; t (74) = 2.5, p = 0.01) of the CD.
The importance of snow is widely recognized in Arctic wildlife research, and snow information that is relevant and appropriate to wildlife applications is of great interest . While the role of snow in caribou winter ecology has been observed and discussed as being important in movement and definition of winter range characteristics (e.g., [57, 128]), research that integrates snow and animal location data to quantify the direct impacts of snow depth on this snow-adapted species has been rare. This is mainly due to the lack of snow data at wildlife-relevant spatiotemporal scales (e.g., hourly to daily and meters instead of kilometers resolutions) and the lack of wildlife-relevant snow property information (e.g., snow depth vs SWE, which is often inaccurately used as proxy for snow depth). To overcome these issues and incorporate germane snow-property data at appropriate scales, we used snow depth datasets that spanned the entire CAH fall and winter range; these datasets were produced at sufficient spatial and temporal resolution to account for local-scale snow-depth variability and evolution to match the spatiotemporally dynamic character of caribou winter movement. Combined with multiple years of CAH locations, these snow depth datasets enabled us to evaluate how caribou adjusted their regional winter range location, and day-to-day movement patterns once on winter range, in response to regional and local snow depth distributions, respectively.
Fall snow depth as a cue to winter range location
Mountains present little obstacle to migratory caribou [19, 129] and CAH caribou are no exception. In 15 of the recent 19 winters (2001/02–2019/20), the majority of collared CAH caribou migrated through and over the BR mountains to access winter ranges south of the CD (Fig. 5; Fig. 6; Table 1), though the specific winter ranges varied longitudinally among winters (this study; [57, 58, 67]). In that context, the four winters 2007/08, 2016/17, 2018/19, and 2019/20 that the CAH spent primarily north of the CD represent an abrupt change to a long-term pattern of using southern winter ranges. The ADFG survey (Table 1) showed that prior to and including winter 2015/16, the CAH winter range selection was more consistent among consecutive years, but after 2015/16, the selection became more variable among years, and the herd tended to winter north of the CD more frequently. This selection pattern may explain the significant linear trend in winter range selection reported in the mixed-effects logistic regression results. However, in the years after 2015/16 the fall snow depth still seemed to exert a strong influence on winter range selection. Such abrupt deviations from one year to the next are also documented in the neighboring Teshekpuk Caribou Herd (TCH) and other Arctic caribou herds [130, 131]. We found a strong significant correlation (r = 0.59, p = 0.001) between the fall snow depth and the proportion of collared caribou wintering south of the CD, and while this relationship was not completely driven by extreme snow years, snow depth had a significant effect on CAH winter range selection. While both our correlation test and logistic regression results suggest that caribou respond to fall snow depth, it is important to emphasize that fall snow depth levels reported herein (< 25 cm) are considered a cue and may have limited biological significance, e.g., in inhibiting caribou travel given their relatively long legs and efficient locomotion in relatively deep snow [6, 19]. However, from a snow-science perspective, fall snow depth is likely to be a reasonable predictor of snow levels that may inhibit movement or forage access in subsequent winter months (January–March) (r = 0.88). North of the CD, the majority of the snowpack is typically established in September–December [85, 111]. Therefore, fall snow depth is often a reasonable predictor of snow-depth levels throughout the remaining winter.
Snow acting as a trigger or cue for migration and winter range selection has been reported previously for the CAH. Roby  noted that, following the first snow in the fall, groups of CAH caribou made a southward movement into the BR, while others stopped at the northern edge of the foothills. Similar examples exist for other Alaska caribou herds, e.g., an extraordinarily heavy August–September snowfall near the Arctic coast and in Anaktuvuk Pass in 1960 and 1961 triggered migration of Western Arctic Herd (WAH) caribou south of the BR . In interior Alaska, the Denali Herd caribou made an unusual movement 200 km north of their regular winter range following a record September snowfall in 1992 . In Canada, Le Corre et al.  found that caribou arrived earlier on winter range in years with early fall snow. They suggested that caribou use the amount of October precipitation to predict snow abundance during migration, and that animals adjust their migration timing to limit the energetic costs of moving through deep snow . Our temporally and spatially explicit snow data can be used in individual-based quantification of such potential cueing effects of snow depth (or other derived snow variables, e.g., rate of accumulation, timing of first snowfall, depth variability, and snow-coverage of forage species) and other relevant environmental variables produced by SnowModel (e.g., air temperature, wind speed, wind direction, and day length) on caribou winter range selection and migration timing.
In reality, migratory caribou winter-range selection is influenced by a suite of environmental and biological factors and herd-specific demographic parameters (e.g., population size, herd density, and sex and age structure; [47, 134, 135]). The importance of these other factors may vary between years and between herds, and they may act together in complex ways. While our data are coarse (i.e., annual surveys with varying sample sizes and areal average snow depth) and our analyses focus solely on estimating the role of fall snow depth in winter range selection over 19 years, these analyses highlight the importance of snow in influencing winter range selection. Furthermore, our results demonstrate that snow should be included in future analyses, along with a range of other factors documented as being influential for caribou wintering location, e.g., overgrazing of slow-growing lichen species [42, 136]. Macander et al.  reported that caribou in interior Alaska prefer lichen habitats in winter until they become overgrazed, and consequently, the herd shifts wintering area within the boreal forest region. The CAH winter range variation may be explained by consistent use of winter ranges south of CD from 2001/02 to 2015/16, interrupted only when the very low fall snow depths allowed caribou to winter north of CD in 2007/08 (Fig. 6). It is possible that by 2016/17, the southern winter ranges could have been overgrazed and the caribou began wintering north, except that they were ‘pushed’ south by record snow in 2017/18 (Fig. 6). The lichen proportion dataset used in our study is invariable in time and provides no information on lichen biomass, and it is therefore inadequate to test this overgrazing hypothesis. Because the TCH, the Porcupine Caribou Herd (PCH), and occasionally the WAH may utilize parts of the same winter ranges as CAH [39, 61] these neighboring and/or larger caribou herds could also exacerbate overgrazing south of the CD.
Although not included in our analyses, predator avoidance may be an additional factor influencing CAH winter range selection, since predation pressure is likely different between the two general wintering areas. During snow-free months, wolves (Canis lupus), grizzly bears (Ursus arctos), and golden eagles (Aquila chrysaetos) prey on CAH caribou, but in winter, wolves are expected to be their primary predator [137, 138]. Wolf surveys report higher wolf densities in the BR mountains (6 wolves per 1000 km2; ) than further north on the tundra of the coastal plain (2–4 wolves per 1000 km2; ). Further research is needed to understand how forage availability and predation pressure, and their interaction with snow, may affect the CAH selection of winter range location and local movement once on winter range.
Finally, herd size is suggested to impact winter range variability; the larger the herd, the larger the required winter range [57, 129, 140]. During the period 2008–2013, when the CAH was near a recent maximum size (estimated population size > 50,000 animals; ), the majority of the herd (≥80% of collared animals; Fig. 6, Table 1) wintered south of the CD, which supports this supposition. However, according to both the strong correlation between winter distribution and fall snow depth and our mixed-effects modeling results suggesting that fall snow depth acts as a cue for winter range selection, and the fact that the winters 2008/09–2012/13 had fall snow-depth levels north of the CD above average for our 19-year time series (> 12.7 cm; Fig. 6), snow depth is also a likely explanation. Moreover, in winter 2017/18, the CAH size was estimated to be ~ 23,000 animals, a recent low, yet most of the collared caribou wintered south of the CD that winter, while fall snow depths were at a record high level on the north side of CD. We propose that the herd’s selection of southern winter ranges may be triggered by relatively deep fall snow depths in addition to herd size. While Nicholson et al.  mapped CAH winter range location and discussed potential drivers of its year-to-year variability, they recommended further research of environmental factors such as inter-annual differences in winter weather and snow cover. Here we investigated one possible mechanism, suggesting that fall snow depth is an important cue, and may likewise be important for caribou and reindeer winter range selection elsewhere in the Arctic outside of our Alaska study domain.
Snow depth and lichen affect caribou movement within winter range
Early scientific accounts describe caribou in Alaska as year-round ‘wanderers’, with directed movements occurring only during fall and spring migration . The mid-winter months are characterized by the least amount of nomadism, and caribou reduce their movement rate to the lowest annual levels during this time (this study; [34, 38, 57, 141, 142]). Our results show that once CAH caribou have selected their winter range, they remain in the same general area until March or the initiation of spring migration (Fig. 5). This pattern is consistent with observations made of other neighboring caribou herds [19, 39, 142].
Decreased movement rate with increasing snow depth
This minimal winter movement may be caused by snow levels. Murie  observed that heavy snowfall can act as a movement barrier for wintering caribou in Alaska. Our results suggest that both on tundra and in mountainous boreal forest, CAH caribou move more slowly when in areas of deeper snow (Fig. 7a). The 24-h fix rate of our GPS locations does not allow for identification of different daily activities (e.g., moving, feeding, resting). However, the daily winter activities of PCH caribou have been found to vary by season, day length, and snow conditions, and these daily activities were dominated by feeding and resting . Russell et al.  found that the least time was spent moving (running, trotting, and walking) and most time was spent lying (ruminating) and feeding during periods of adverse snow conditions, including increased snow depth. These observations of decreased movement and increased lying and feeding intensity can translate to an overall lower movement rate (i.e., shorter distance traveled per day) in deep snow, which our results support. Further, based on our knowledge of this snow-covered landscape, its differential snowpack properties (tundra vs. mountainous boreal forest snow), and informal observations of the caribou behavior and movement in these landscapes during our field campaigns, we speculate that there are several plausible explanations for the decreased movement rate with increasing snow depth (Fig. 7a), and that these are different for the mountainous boreal forest and tundra winter habitats.
In the mountainous boreal forest, a snowpack of more than 70 cm deep can reduce caribou mobility  and increase energetic costs, likely because this depth threshold exceeds their approximate chest height (~ 50 cm; adult domesticated reindeer; ) or leg length (~ 50–60 cm; adult CAH caribou; hind foot length; ). In our mountainous boreal forest snow sites, average snow depth ranged from 50 ± 8 cm to 111 ± 13 cm (March/April 2018 and 2019 average snow depth = 76 ± 17 cm for all mountainous boreal forest snow sites, n = 18), and the average modeled snow depth of 49 cm suggests that many of the starting locations of used steps exceed or are near the snow-depth threshold for impeding caribou movement. These observations support the iSSA result of decreased movement rates when caribou encountered increased snow depths. Caribou may choose to remain in a smaller area due to the higher energetic costs of movement in areas of deep snow. Such reduction in movement rate with increased snow depth is also observed for PCH caribou during October in the central BR mountains , where the CAH and PCH winter ranges overlap . Conversely, the iSSA result can also be interpreted as increased movement rate with decreasing snow depth. On the winter ranges south of the CD, the snow depth is relatively shallow in tree-less alpine tundra areas or on lakes large enough for the wind to erode the snowpack. This relatively shallow snow may promote fast movement, for example, to avoid predators. Finally, it takes more time to crater through deep snow than shallow snow. Hence, feeding caribou in the relatively deep snow of the mountainous boreal forest often stay longer, dig one large crater, and move less than in areas with shallow snow .
On the tundra, the snow depths are generally more shallow (average 45 ± 21 cm for all tundra snow sites; n = 53; ) and present a limited impeding effect on movement, particularly in the frequently-used, wind-blown north-facing valleys with shallow snow depth (5 ± 5 cm – 22 ± 10 cm; n = 5; ). Furthermore, the GPS-collared caribou moved significantly faster in winters spent primarily north of the CD on tundra, than when south of the CD. Hence, the increased movement rate with decreasing snow depth may be the most plausible explanation for the iSSA result for tundra (Fig. 7a). Though, in more sheltered areas of the otherwise windy tundra environment, e.g., areas between the broad, wind-blown valleys, the snow was relatively deep (ranging 42 ± 13 cm – 90 ± 10 cm; n = 5; ). Our snow observations show that these areas experienced limited wind erosion, thus the snow there is both softer and deeper, and possibly capable of hindering caribou movement because of these properties. Finally, the generally shallow tundra snowpack also allows faster feeding, faster movement, and the ability to travel longer distances, because less time is spent digging and cratering .
Selection of snow depth
Our results further suggest that CAH caribou wintering south of the CD in the mountainous boreal forest selected for areas with less snow than was available in the surrounding area (Fig. 7a, c). Arctic Alaska caribou are often found on wind-blown ridges and mountain slopes, where forage is more readily available because of snow-removal by wind . While they do not hesitate to descend into the spruce forest and to crater through considerable snow depositions for lichen [19, 143], caribou wintering in boreal forest systems tend to select feeding sites in more shallow snow . Significant selection of more shallow snow depths may also be explained by the inherent characteristic of the boreal forest snowpack; it is typically fully comprised of relatively large, weakly bonded snow grains (depth hoar). When cratered by caribou, the structure of the snowpack fails and the disturbed snow subsequently metamorphoses into a more dense and well-bonded snow cover [145, 146]. Hence, reworking relatively deep snow that has previously been cratered and solidified is more energy consuming, and for this reason the likelihood of reusing previously cratered snow is low . Selecting areas with shallow snow for feeding may also be advantageous for energy conservation . The hardening process of cratered and disturbed snow also occurs in tundra snow, but the areal extent of undisturbed and shallow snow makes this a negligible factor on selection within the tundra landscape. This may explain our result showing that GPS-collared caribou did not show any significant selection for snow depth during winters spent north of the CD on the tundra.
Selection of lichen proportion
The wandering nature of caribou is undoubtedly influenced by the snow depths they encounter during winter, but the scattered spatial distribution of lichen  may also contribute to their movement behavior. We found that, irrespective of winter range location, caribou selected for areas with higher lichen proportion than otherwise available in the surrounding landscape, indicating the importance of this winter forage resource (Fig. 7a, b). This result agrees with observations by Roby , who found that lichen was consistently selected by the CAH caribou in feeding areas, despite its sparse distribution. The highest relative selection strength for lichen proportion was seen in CAH caribou wintering north of the CD (Fig. 7a), where fruticose lichen is abundant on surfaces where the wind has eroded the snowpack . The effects of both forage distribution, and the snow depth on the accessibility of that forage, are considered to be related and inseparable factors influencing caribou winter movement and foraging behavior and should be assessed together [32, 35, 36]. While the dataset of distributed lichen proportion applied in our movement models does not provide information on lichen biomass, and admittedly represents a coarse measure of winter caribou forage abundance, this study is one of few studies quantifying the combined effects of snow depth and winter forage on caribou movement. New methodologies for mapping and modeling distributions of lichen and other caribou forage species are emerging [105,106,107,108, 147], however, these improved lichen and forage distribution estimates are not available across our entire study area. Future research using such data, e.g., in combination with SnowModel snow depth and caribou location data, could help quantify the effects of snow depth on caribou winter forage accessibility by cratering (e.g., [50, 148]). Key to such investigations are location data of sufficient temporal detail to determine activity states (e.g., foraging/cratering, running, walking, resting/ruminating) of each individual (e.g., [149, 150]). Hence, our results are likely the first of many investigations to combine snow and forage data to advance our understanding of caribou winter foraging behavior across space and time.
We acknowledge that, in addition to snow depth, other aspects of the snowpack can be critical to caribou winter ecology. The strength of the snow surface and individual snow or ice layers buried within or below the snowpack is likely important for caribou locomotion and foraging accessibility in both tundra and mountainous boreal forest wintering habitats. Quantitative models of strength properties of snow and ice across space and time do not yet exist, and essential to their development is the ability to represent the direct effects of ice layers formed by, e.g., rain-on-snow (ROS) and mid-winter melt-freeze events . ROS events are still rare in Arctic Alaska, but predicted to increase in frequency , and have previously caused a major change in the TCH winter distribution [131, 152]. Furthermore, ROS events are hypothesized to account for the majority of PCH population fluctuations , and the icing resulting from ROS events influence ungulate forage access [153, 154]. Across the Arctic, investigations of the relationships between climate variables, including some winter and snow-related parameters, and long-term time series of caribou and reindeer life-history event observations (e.g., calving date and timing of migration), indicate that winter processes are highly influential factors (e.g., [155, 156]). We encourage further investigation of the underlying mechanisms linking snow to caribou life-history events to understand the effects that snowpack depth and strength can exert on caribou nutritional condition through, e.g., regulating forage accessibility and the energetic expense of locomotion. Both are likely to have consequences for caribou fitness at the individual and population level through parameters including survival, parturition or birth rate, recruitment of calves to the herd, or reproduction timing [157,158,159]. Spatially and temporally explicit snow information such as the data acquired and used in the research presented herein are highly applicable to such investigations, because they enable a quantification of and accounting for: (1) inter-seasonal carry-over effects (e.g., [158, 160]); (2) short-lived weather events that have significant long-term population effects (e.g., mid-winter icing or ROS events [64, 161, 162]); (3) snow phenology throughout winter (i.e., temporal evolution of the snowpack and timing of snowfall events ); and (4) differential effects of snow on caribou and their predators (e.g., ). Additionally, such studies would benefit from interdisciplinary research teams of snow and wildlife professionals collaboratively working to gain a more complete understanding of the aspects of caribou ecology dependent on, and influenced by, snow.
Across Arctic Alaska, a long history of use, dependence, and research, through observations (e.g., [7, 129]) and Traditional Ecological Knowledge (e.g., [165, 166]), has laid the foundation for our current understanding of caribou ecology. In our research, we integrated methodologies from snow and wildlife sciences. We used SnowModel and field observations to provide snow depth data that resolved the snow-depth variability across both the entire CAH home range and at a more detailed scale representing the snow-depth distributions caribou encountered on a day-to-day basis. We found that fall snow depth on the north side of the CD likely is an important cue that either, if shallow, motivates CAH caribou to stay on the tundra north of the CD, or if deep, promotes migration to winter ranges in the mountainous boreal forest south of the CD. When wintering south of the CD, caribou avoided areas with deeper snow and selected areas with higher lichen abundance. Both on the tundra and in the mountainous boreal forest, caribou movement rates decreased with increased snow depth. For snow properties, this research focused solely on snow depth, but we anticipate that numerous other snow properties are key in governing wildlife movement and health. Using SnowModel in future investigations will enable answering increasingly complex questions involving other snow properties (such as strength), icing events, forage accessibility, and winter mobility that will greatly advance our knowledge of caribou winter ecology.
Availability of data and materials
The datasets generated and analyzed during the current study are available and permanently archived at NSF Arctic Data Center: Stine Højlund Pedersen. 2021. Observed and modeled snow depth in Arctic Alaska 2001–2020. Arctic Data Center. https://doi.org/10.18739/A2WM13V3Q. R scripts generated during this study are included in Additional file 2. The caribou GPS location datasets analyzed during the current study are not publicly available due to Alaska law [AS 16.05.815(d)]. These laws allow the release of such information to parties that have been authorized to perform specific activities as long as the parties agree to use the information only for purposes as provided under an agreement with the Alaska Department Fish and Game, contact Elizabeth A. Lenart.
Alaska Department of Fish and Game
Central Arctic Herd
Global Positioning System
integrated step selection analysis
kernel density estimate
Kilometer per day
North American Land Change Monitoring System
Porcupine Caribou Herd
Rain on snow
Snow water equivalent
Teshekpuk Caribou Herd
Very high frequency
Western Arctic Herd
Liston GE, Hiemstra CA. The changing cryosphere: pan-Arctic snow trends (1979-2009). J Clim. 2011;24(21):5691–712. https://doi.org/10.1175/JCLI-D-11-00081.1.
Formozov AN. Snow cover as an integral factor of the environment and its importance in the ecology of mammals and birds. English ed. Edmonton: Boreal Institute for Northern Studies, University of Alberta; 1946:153.
Pruitt WO. Snow as a factor in the winter ecology of the barren ground caribou (Rangifer arcticus). Arctic. 1959;12(3):158–79. https://doi.org/10.14430/arctic3723.
LaPerriere AJ, Lent PC. Caribou feeding sites in relation to snow characteristics in northeastern Alaska. Arctic. 1977;30(2):101–8. https://doi.org/10.14430/arctic2690.
Telfer ES, Kelsall JP. Adaptation of some large north American mammals for survival in snow. Ecology. 1984;65(6):1828–34. https://doi.org/10.2307/1937779.
Fancy SG, White RG. Energy expenditures for locomotion by barren-ground caribou. Can J Zool. 1987;65(1):122–8. https://doi.org/10.1139/z87-018.
Murie OJ. Alaska-Yukon Caribou. North American Fauna No. 54. Washington, DC: United States Department of Agriculture, Bureau of Biological Survey; 1935. p. 1–92.
Cuyler C, Øritsland NA. Effect of wind on Svalbard reindeer fur insulation. Rangifer. 2002;22(1):93–9. https://doi.org/10.7557/188.8.131.524.
Cuyler C, Øritsland NA. Do seasonal changes in Svalbard reindeer fur have relevance for heat transfer? Rangifer. 2002;22(2):133–42. https://doi.org/10.7557/184.108.40.2062.
Nieminen M. Hoof and foot loads for reindeer (Rangifer tarandus). Rangifer. 1990;10(3):249–54. https://doi.org/10.7557/220.127.116.115.
Kelsall JP. The migratory barren-ground caribou of Canada. 1. Ottawa: Queen’s Printer; 1968:340.
Fancy SG, White RG. The incremental cost of activity. In: Hudson RJ, White RG, editors. Bioenergetics of wild herbivores. Boca Raton: CRC Press; 1985:143–59.
Zhang R, Qiao Y, Ji Q, Ma S, Li J. Macro-microscopic research in reideer (Rangifer tarandus) hoof suitable for efficient locomotion on complex grounds. J Vet Res. 2017;61(2):223–9. https://doi.org/10.1515/jvetres-2017-0029.
Brown WK, Theberge JB. The effect of extreme snowcover on feeding-site selection by woodland caribou. J Wildl Manag. 1990;54(1):161–8. https://doi.org/10.2307/3808916.
Lorenzen ED, Nogués-Bravo D, Orlando L, Weinstock J, Binladen J, Marske KA, et al. Species-specific responses of Late Quaternary megafauna to climate and humans. Nature. 2011;479(7373):359–64. https://doi.org/10.1038/nature10574.
Flagstad Ø, Røed KH. Refugial origin of reindeer (Rangifer tarandis L.) inferred from metochondrial DNA sequences. Evolution. 2003;57(3):658–70. https://doi.org/10.1111/j.0014-3820.2003.tb01557.x.
Tyler NJC. Climate, snow, ice, crashes, and declines in populations of reindeer and caribou (Rangifer tarandus L.). Ecol Monogr. 2010;80(2):197–219. https://doi.org/10.1890/09-1070.1.
Ihl C, Klein DR. Habitat and diet selection by muskoxen and reindeer in Western Alaska. J Wildl Manag. 2001;65(4):964–72. https://doi.org/10.2307/3803045.
Russell DE, Martell AM, Nixon WAC. Range ecology of the porcupine caribou herd in Canada. Rangifer. 1993;13(5):1–168. https://doi.org/10.7557/18.104.22.1687.
Duquette LS. Snow characteristics along caribou trails and within feeding areas during spring migration. Arctic. 1988;41(2):143–4. https://doi.org/10.14430/arctic1706.
Boelman NT, Liston GE, Gurarie E, Meddens AJH, Mahoney PJ, Kirchner PB, et al. Integrating snow science and wildlife ecology in Arctic-boreal North America. Environ Res Lett. 2019;14(1):e010401. https://doi.org/10.1088/1748-9326/aaeec1.
Liston GE, Elder K. A distributed snow-evolution modeling system (SnowModel). J Hydrometeorol. 2006;7(6):1259–76. https://doi.org/10.1175/JHM548.1.
Joly K, Gurarie E, Sorum MS, Kaczensky P, Cameron MD, Jakes AF, et al. Longest terrestrial migrations and movements around the world. Sci Rep. 2019;9(1):15333. https://doi.org/10.1038/s41598-019-51884-5.
Merkle JA, Sawyer H, Monteith KL, Dwinnell SPH, Fralick GL, Kauffman MJ. Spatial memory shapes migration and its benefits: evidence from a large herbivore. Ecol Lett. 2019;22(11):1797–805. https://doi.org/10.1111/ele.13362.
Middleton AD, Merkle JA, McWhirter DE, Cook JG, Cook RC, White PJ, et al. Green-wave surfing increases fat gain in a migratory ungulate. Oikos. 2018;127(7):1060–8. https://doi.org/10.1111/oik.05227.
Bauer S, Nolet B, Giske J, Chapman J, Akesson S, Hedenstrom A, et al. Cues and Decision rules in animal migration. In: Milner-Gulland EJ, Fryxell JM, ARE S, editors. Animal migration - A synthesis. New York: Oxford University Press; 2011:269.
Visser ME, Caro SP, van Oers K, Schaper SV, Helm B. Phenology, seasonal timing and circannual rhythms: towards a unified framework. Philos Trans R Soc B Biol Sci. 2010;365(1555):3113–27. https://doi.org/10.1098/rstb.2010.0111.
Dingle H, Drake VA. What is migration? Bioscience. 2007;57(2):113–21. https://doi.org/10.1641/B570206.
Oliver RY, Mahoney PJ, Gurarie E, Krikun N, Weeks BC, Hebblewhite M, et al. Behavioral responses to spring snow conditions contribute to long-term shift in migration phenology in American robins. Environ Res Lett. 2020;15(4):45003. https://doi.org/10.1088/1748-9326/ab71a0.
Boyce MS. Migratory behavior and management of elk (Cervus elaphus). Appl Anim Behav Sci. 1991;29(1):239–50. https://doi.org/10.1016/0168-1591(91)90251-R.
Igota H, Sakuragi M, Uno H, Kaji K, Kaneko M, Akamatsu R, et al. Seasonal migration patterns of female sika deer in eastern Hokkaido. Japan Ecol Res. 2004;19(2):169–78. https://doi.org/10.1111/j.1440-1703.2003.00621.x.
Bergerud AT. Relative abundance of food in winter for Newfoundland caribou. Oikos. 1974;25(3):379–87. https://doi.org/10.2307/3543960.
Bergerud AT, Luttich SN, Camps L. The Return of Caribou to Ungava. Montreal: McGill-Queen’s University Press; 2008:656.
Lent PC. The caribou of northwestern Alaska. In: Wilimovsky NJ, Wolfe JN, editors. Environment of the Cape Thompson region, Alaska. Washington: United States Atomic Energy Commission; 1966:481–517.
Eastland WG. Ph.D. Thesis: Influence of weather on movements and migrations of caribou. Fairbanks: University of Alaska Fairbanks; 1991.
Roby DD. Thesis: Behavioural patterns of barren-ground caribou of the Central Arctic herd adjacent to the Trans-Alaska Oil Pipeline. Fairbanks: University of Alaska; 1978.
Ferguson SH, Elkie PC. Seasonal movement patterns of woodland caribou (Rangifer tarandus caribou). J Zool. 2004;262(2):125–34. https://doi.org/10.1017/S0952836903004552.
Fancy SG, Pank LF, Whitten KR, Regelin WL. Seasonal movements of caribou in arctic Alaska as determined by satellite. Can J Zool. 1989;67(3):644–50. https://doi.org/10.1139/z89-093.
Person BT, Prichard AK, Carroll GM, Yokel DA, Suydam RS, George JC. Distribution and movements of the Teshekpuk caribou herd 1990-2005: prior to oil and gas development. Arctic. 2007;60(3):238–50. https://www.jstor.org/stable/40512892.
Turney L, Heard DC. The use of satellite images to estimate snow depth and distribution on the for-ested winter range of the Beverly caribou herd. Rangifer. 1991;11(4):164–6. https://doi.org/10.7557/22.214.171.1241.
Tucker B, Mahoney S, Greene B, Menchenton E, Russell L. The influence of snow depth and hardness on winter habitat selection by caribou on the southwest coast of Newfoundland. Rangifer. 1991;11(4):160–3. https://doi.org/10.7557/126.96.36.1990.
Joly K, Chapin FS III, Klein DR. Winter habitat selection by caribou in relation to lichen abundance, wildfires, grazing, and landscape characteristics in Northwest Alaska. Ecoscience. 2010;17(3):321–33. https://doi.org/10.2980/17-3-3337.
Joly K, Dale BW, Collins WB, Adams LG. Winter habitat use by female caribou in relation to wildland fires in interior Alaska. Can J Zool. 2003;81(7):1192–201. https://doi.org/10.1139/z03-109.
Joly K, Cameron MD. Early fall and late winter diets of migratory caribou in Northwest Alaska. Rangifer. 2018;38(1):27–38. https://doi.org/10.7557/188.8.131.5207.
Jandt R, Joly K, Randy Meyers C, Racine C. Slow recovery of lichen on burned caribou winter range in Alaska tundra: potential influences of climate warming and other disturbance factors. Arctic Antarct Alp Res. 2008;40(1):89–95. https://doi.org/10.1657/1523-0430(06-122)[JANDT]2.0.CO;2.
Llano GA. Utilization of lichens in the arctic and subarctic. Econ Bot. 1956;10(4):367–92. https://doi.org/10.1007/BF02859767.
Anderson TA, Johnson CJ. Distribution of barren-ground caribou during winter in response to fire. Ecosphere. 2014;5(10):140. https://doi.org/10.1890/ES14-00010.1.
Johnson CJ, Parker KL, Heard DC, Gillingham MP. A multiscale behavioral approach to understanding the movements of woodland caribou. Ecol Appl. 2002;12(6):1840–60. https://doi.org/10.1890/1051-0761(2002)012[1840:AMBATU]2.0.CO;2.
Thing H. Behavior, mechanics and energetics associated with winter cratering by caribou in northwestern Alaska. Biological Papers of the University of Alaska. Fairbanks: University of Alaska Fairbanks. 1977;18:1–31.
Collins WB, Smith TS. Effects of wind-hardened snow on foraging by reindeer (Rangifer tarandus). Arctic. 1991;44(3):217–22. http://www.jstor.org/stable/40511242.
Johnson CJ, Parker KL, Heard DC. Feeding site selection by woodland caribou in north-central British Columbia. Rangifer. 2000;20(5):158–72. https://doi.org/10.7557/184.108.40.2062.
Bergerud AT. The role of the environment in the aggregation, movement, and disturbance behavior of caribou. International Symposium Papers. The behaviour of ungulates and its relation to management. Morges: International Union for Conservation of Nature and Natural Resources. 1974;2:552–85.
Johnson CJ, Parker KL, Heard DC. Foraging across a variable landscape: behavioral decisions made by woodland caribou at multiple spatial scales. Oecologia. 2001;127(4):590–602. https://doi.org/10.1007/s004420000573.
Parker KL, Barboza PS, Gillingham MP. Nutrition integrates environmental responses of ungulates. Funct Ecol. 2009;23(1):57–69. https://doi.org/10.1111/j.1365-2435.2009.01528.x.
Jakimchuk RD, Ferguson SH, Sopuck LG. Differential habitat use and sexual segregation in the Central Arctic caribou herd. Can J Zool. 1987;65(3):534–41. https://doi.org/10.1139/z87-083.
Latombe G, Fortin D, Parrott L. Spatio-temporal dynamics in the response of woodland caribou and moose to the passage of grey wolf. J Anim Ecol. 2014;83(1):185–98. https://doi.org/10.1111/1365-2656.12108.
Nicholson KL, Arthur SM, Horne JS, Garton EO, Vecchio D. Modeling caribou movements: seasonal ranges and migration routes of the Central Arctic Herd. PLoS One. 2016;11(4):150333. https://doi.org/10.1371/journal.pone.0150333.
Arthur SM, Del VPA. Effects of oil field development on calf production and Survival in the Central Arctic Herd federal aid in wildlife restoration grants. Juneau: Alaska Department of Fish and Game; 2009.
Cameron RD, Whitten KR. Seasonal movements and sexual segregation of caribou determined by aerial survey. J Wildl Manag. 1979;43(3):626–33. https://doi.org/10.2307/3808740.
Cameron RD, Smith WT, White RG, Griffith B. Central Arctic caribou and petroleum development: distributional, nutritional, and reproductive implications. Arctic. 2005;58(1):1–9. https://www.jstor.org/stable/40512660.
Prichard AK, Parrett LS, Lenart EA, Caikoski JR, Joly K, Person BT. Interchange and overlap among four adjacent Arctic caribou herds. J Wildl Manag. 2020;84(8):1500–14. https://doi.org/10.1002/jwmg.21934.
Johnson HE, Golden TS, Adams LG, Gustine DD, Lenart EA. Caribou use of habitat near energy development in Arctic Alaska. J Wildl Manag. 2020;84(3):401–12. https://doi.org/10.1002/jwmg.21809.
White RG, Thomson BR, Skogland T, Person SJ, Russell DE, Holleman DF, et al. Ecology of caribou at Prudhoe Bay, Alaska. In: Brown J, editor. Ecological Investigations of the Tundra Biome in the Prudhoe Bay Region, Alaska. Fairbanks: Biological Papers, University of Alaska, Special Report No. 2; 1975:215.
Griffith B, Douglas DC, Walsh NE, Young DD, McCabe TR, Russell DE, et al. The Porcupine caribou herd. In: Douglas DC, Reynolds PE, Rhode EB, editors. Arctic Refuge coastal plain terrestrial wildlife research summaries. Reston: U. S. Geological Survey, Biological Science Report USGS/BRD/BSR-2002-0001, Biological Resources Division; 2002:8–37.
Cameron RD, Lenart EA, Reed DJ, Whitten KR, Smith WT. Abundance and movements of caribou in the oilfield complex near Prudhoe Bay, Alaska. Rangifer. 1995;15(1):3–7. https://doi.org/10.7557/220.127.116.110.
Lenart EA. Central Arctic caribou management report and plan, game management unit 26B: report period 1 July 2012–30 June 2017, and plan period 1 July 2017–30 June 2022. Juneau: Alaska Department of Fish and Game; 2021.
Lenart EA. Units 26B and 26C caribou. In: Harpe P, LA MC, editors. Caribou management report of survey and inventory activities, 1 July 2012–30 June 2014. Species Management Report and Plan. Juneau: Alaska Department of Fish and Game; 2015:38.
Lenart EA. 2019 Central Arctic Caribou winter distribution, parturition, postcalving, and captures Memorandum, 01/30/2019. Fairbanks: Alaska Department of Fish and Game; 2019.
Lenart EA. 2018 Central Arctic Caribou winter distribution, parturition, postcalving, and captures Memorandum, 01/15/2019. Fairbanks: Alaska Department of Fish and Game; 2019.
Lenart EA. 2020 Central Arctic Caribou winter distribution, parturition, postcalving, and fall composition Memorandum, 10/28/2020. Fairbanks: Alaska Department of Fish and Game; 2020.
Lenart EA. Central Arctic Caribou Herd - Units 26B and 26C. In: Harper P, editor. Caribou management report of survey-inventory activities 1 July 2010–30 June 2012. Species Management Report and Plan. Juneau: Alaska Department of Fish and Game; 2013:356–89.
Lenart EA. Units 26b and 26c caribou. In: Harper P, editor. Caribou management report of survey-inventory activities 1 July 2008–30 June 2010. Project 3. Juneau: Alaska Department of Fish and Game; 2011:315–45.
Lenart EA. Units 26b and 26c caribou. In: Harper P, editor. Caribou management report of survey-inventory activities 1 July 2006–30 June 2008. Project 3. Juneau: Alaska Department of Fish and Game; 2009:299–325.
Boggs K, Flagstad L, Boucher T, Kuo T, Fehringer D, Guyer S, et al. Vegetation map and classification: Northern, Western, and interior Alaska. Alaska Natural Heritage Program, University of Alaska Anchorage, Alaska Center for Conservation Science, Anchorage, AK, USA; 2016.
Walker DA. Vegetation and environmental gradients of the Prudhoe Bay region, Alaska. Hanover: US Army Cold Regions Research and Engineering Laboratory; 1985.
Nowacki GJ, Spencer P, Fleming M, Brock T, Jorgenson T. Unified Ecoregions of Alaska: 2001. Open-File Report 2002–297; 2003. https://doi.org/10.3133/ofr2002297.
Walker DA, Binnian E, Evans BM, Lederer ND, Nordstrand E, Webber PJ. Terrain, vegetation and landscape evolution of the R4D research site, Brooks Range foothills, Alaska. Ecography. 1989;12(3):238–61. https://doi.org/10.1111/j.1600-0587.1989.tb00844.x.
Pedersen SH, Liston GE, Welker JM. Snow depth and snow density measured in Arctic Alaska for caribou winter applications in 2018 and 2019. Arctic Data Center. https://arcticdata.io/catalog/view/doi:10.18739/A2MS3K24P.
Sturm M, Holmgren J, Liston GE. A seasonal snow cover classification system for local to global applications. J Clim. 1995;8(5):1261–83. https://doi.org/10.1175/1520-0442(1995)008<1261:ASSCCS>2.0.CO;2.
Benson CS. The seasonal snow cover of Arctic Alaska. Res Pap No 51. Arct Inst North Am. 1969;51:1–80.
Zhang T, Osterkamp TE, Stamnes K. Some characteristics of the climate in northern Alaska, USA. Arct Alp Res. 1996;28(4):509–18. https://doi.org/10.2307/1551862.
Overland JE, Wang M, Ballinger TJ. Recent increased warming of the Alaskan marine Arctic due to midlatitude linkages. Adv Atmos Sci. 2018;35(1):75–84. https://doi.org/10.1007/s00376-017-7026-1.
Bieniek PA, Bhatt US, Thoman RL, Angeloff H, Partain J, Papineau J, et al. Climate divisions for Alaska based on objective methods. J Appl Meteorol Climatol. 2012;51(7):1276–89. https://doi.org/10.1175/JAMC-D-11-0168.1.
Shulski M, Wendler G. The climate of Alaska. Fairbanks: University of Alaska Press; 2007.
Olsson PQ, Hinzman LD, Sturm M, Liston GE, Kane DL. Surface climate and snow–weather relationships of the Kuparuk Basin on Alaska’s Arctic Slope, Technical Report ERDC/CRRELTR-02-10. Hanover: US Army Corps of Engineers Engineer Research and Development Center; 2002.
Sturm M, Holmgren J. An automatic snow depth probe for field validation campaigns. Water Resour Res. 2018;54(11):9695–701. https://doi.org/10.1029/2018WR023559.
Liston GE, Itkin P, Stroeve J, Tschudi M, Stewart JS, Pedersen SH, et al. A lagrangian snow-evolution system for sea-ice applications (snowmodel-LG): part I–model description. J Geophys Res Ocean. 2020;125(10):e2019JC015913. https://doi.org/10.1029/2019JC015913.
Stuefer S, Kane DL, Liston GE. In situ snow water equivalent observations in the US arctic. Hydrol Res. 2013;44(1):21–34. https://doi.org/10.2166/nh.2012.177.
Liston GE, Hiemstra CA. A simple data assimilation system for complex snow distributions (SnowAssim). J Hydrometeorol. 2008;9(5):989–1004. https://doi.org/10.1175/2008JHM871.1.
Pedersen SH, Liston GE, Tamstorf MP, Abermann J, Lund M, Schmidt NM. Quantifying snow controls on vegetation greenness. Ecosphere. 2018;9(6):e02309. https://doi.org/10.1002/ecs2.2309.
Liston GE, Elder K. A meteorological distribution system for high-resolution terrestrial modeling (MicroMet). J Hydrometeorol. 2006;7(2):217–34. https://doi.org/10.1175/JHM486.1.
Hill K, Sousanes P. WC_G ARCN RAWS corrected data. 2020 (Accessed 1 May 2020); https://irma.nps.gov/DataStore/Reference/Profile/2239061.
NCEI. Land-based weather station data - Local Climatological Data (LCD). 2020 (Accessed 1 May 2020); https://www.ncdc.noaa.gov/data-access/land-based-station-data/land-based-datasets
Gelaro R, McCarty W, Suárez MJ, Todling R, Molod A, Takacs L, et al. The modern-era retrospective analysis for research and applications, version 2 (MERRA-2). J Clim. 2017;30(14):5419–54. https://doi.org/10.1175/JCLI-D-16-0758.1.
USGS. USGS National Elevation Dataset (NED), 1 arc-second. 2020. https://catalog.data.gov/dataset/1-arc-second-digital-elevation-models-dems-usgs-national-map-3dep-downloadable-data-collection. Accessed 1 Jan 2020.
Wickham J, Homer C, Vogelmann J, McKerrow A, Mueller R, Herold N, et al. The multi-resolution land characteristics (MRLC) consortium—20 years of development and integration of USA national land cover data. Remote Sens. 2014;6(8):7424–41. https://doi.org/10.3390/rs6087424.
Liston GE, Sturm M. A snow-transport model for complex terrain. J Glaciol. 1998;44(148):498–516. https://doi.org/10.1017/S0022143000002021.
Liston GE, Haehnel RB, Sturm M, Hiemstra CA, Berezovskaya S, Tabler RD. Instruments and methods simulating complex snow distributions in windy environments using SnowTran-3D. J Glaciol. 2007;53(181):241–56. https://doi.org/10.3189/172756507782202865.
Sturm M, Taras B, Liston GE, Derksen C, Jonas T, Lea J. Estimating snow water equivalent using snow depth data and climate classes. J Hydrometeorol. 2010;11(6):1380–94. https://doi.org/10.1175/2010JHM1202.1.
Stuefer SL, Kane DL, Dean KM. Snow water equivalent measurements in remote arctic Alaska watersheds. Water Resour Res. 2020;56(4):e2019WR025621. https://doi.org/10.1029/2019WR025621.
Sturm M, Parr C, Pedersen SH, Arp C, Urban F, Welker JM, et al. A report on the snow conditions of the north slope and Brooks Range of Alaska during the winter of 2018. Fairbanks: University of Alaska Fairbanks; 2018.
Macander MJ, Swingley CS, Joly K, Raynolds MK. Landsat-based snow persistence map for Northwest Alaska. Remote Sens Environ. 2015;163:23–31. https://doi.org/10.1016/j.rse.2015.02.028.
Liston GE, Sturm M. The role of winter sublimation in the Arctic moisture budget. Nord Hydrol. 2004;35(4–5):325–34. https://doi.org/10.2166/nh.2004.0024.
Moser TJ, Nash TH III, Thomson JW. Lichens of Anaktuvuk Pass, Alaska, with emphasis on the impact of caribou grazing. Bryologist. 1979;82(3):393–408. https://doi.org/10.2307/3242215.
Nelson PR, Roland C, Macander MJ, McCune B. Detecting continuous lichen abundance for mapping winter caribou forage at landscape spatial scales. Remote Sens Environ. 2013;137:43–54. https://doi.org/10.1016/j.rse.2013.05.026.
Kennedy B, Pouliot D, Manseau M, Fraser R, Duffe J, Pasher J, et al. Assessment of Landsat-based terricolous macrolichen cover retrieval and change analysis over caribou ranges in northern Canada and Alaska. Remote Sens Environ. 2020;240:111694. https://doi.org/10.1016/j.rse.2020.111694.
Macander MJ, Frost GV, Nelson PR, Swingley CS. Regional quantitative cover mapping of tundra plant functional types in Arctic Alaska. Remote Sens. 2017;9(10):1024. https://doi.org/10.3390/rs9101024.
Macander MJ, Palm EC, Frost GV, Herriges JD, Nelson PR, Roland C, et al. Lichen cover mapping for caribou ranges in interior Alaska and Yukon. Environ Res Lett. 2020;15(5):55001. https://doi.org/10.1088/1748-9326/ab6d38.
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. https://doi.org/10.1002/ece3.4823.
Thurfjell H, Ciuti S, Boyce MS. Applications of step-selection functions in ecology and conservation. Mov Ecol. 2014;2(4):1–12. https://doi.org/10.1186/2051-3933-2-4.
Olsson PQ, Sturm M, Racine CH, Romanovsky V, Liston GE. Five stages of the Alaskan Arctic cold season with ecosystem implications. Arctic Antarct Alp Res. 2003;35(1):74–81. https://doi.org/10.1657/1523-0430(2003)035[0074:FSOTAA]2.0.CO;2.
Crawley MJ. Statistics: an introduction using R. 1st ed. Chichester: Wiley; 2011:342.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48. https://doi.org/10.18637/jss.v067.i01.
Calenge C. Analysis of animal movements in R: the adehabitatLT Package. adehabitatLT.pdf; 2019:85. https://cran.r-project.org/web/packages/adehabitatLT/adehabitatLT.pdf. Accessed 1 Jan 2020.
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. https://doi.org/10.1016/j.ecolmodel.2006.03.017.
Lavielle M. Using penalized contrasts for the change-point problem. Signal Process. 2005;85(8):1501–10. https://doi.org/10.1016/j.sigpro.2005.01.012.
Lavielle M. Detection of multiple changes in a sequence of dependent variables. Stoch Process their Appl. 1999;83(1):79–102. https://doi.org/10.1016/S0304-4149(99)00023-X.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2019.
Avgar T, Potts JR, Lewis MA, Boyce MS. Integrated step selection analysis: bridging the gap between resource selection and animal movement. Methods Ecol Evol. 2016;7(5):619–30. https://doi.org/10.1111/2041-210X.12528.
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. https://doi.org/10.1890/04-0953.
Matthiopoulos J. The use of space by animals as a function of accessibility and preference. Ecol Model. 2003;159(2):239–68. https://doi.org/10.1016/S0304-3800(02)00293-4.
Forester JD, Im HK, Rathouz PJ. Accounting for animal movement in estimation of resource selection functions: sampling and data analysis. Ecology. 2009;90(12):3554–65. https://doi.org/10.1890/08-0874.1.
Fieberg J, Johnson DH. MMI: multimodel inference or models with management implications? J Wildl Manag. 2015;79(5):708–18. https://doi.org/10.1002/jwmg.894.
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. https://doi.org/10.1002/ece3.3122.
Cramer EM. Multicollinearity. In: S K, NL J, editors. Encyclopedia of Statistical Sciences. New York: Wiley; 1985:639–43.
Dormann CF, Elith J, Bacher S, Buchmann C, Carl G, Carré G, et al. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography. 2013;36(1):27–46. https://doi.org/10.1111/j.1600-0587.2012.07348.x.
Sturm M, Parr C, Larsen C. A report on the snow cover of the 1002 area of the Arctic National Wildlife Refuge, 2014–2019. Fairbanks: University of Alaska Fairbanks; 2019.
Baltensperger AP, Joly K. Using seasonal landscape models to predict space use and migratory patterns of an arctic ungulate. Mov Ecol. 2019;7(18):1–19. https://doi.org/10.1186/s40462-019-0162-8.
Skoog RO. Ph.D. Thesis: Ecology of the caribou (Rangifer tarandus granti) in Alaska. Berkeley: University of California; 1968.
Kohler J, Aanes R. Effect of winter snow and ground-icing on a Svalbard reindeer population: results of a simple snowpack model. Arctic Antarct Alp Res. 2004;36(3):333–41. https://doi.org/10.1657/1523-0430(2004)036[0333:EOWSAG]2.0.CO;2.
Carroll G. Unit 26A caribou management report. In: Brown C, editor. Caribou management report of survey and inventory activities 1 July 2002–30 June 2004. Juneau: Alaska Department of Fish and Game; 2005;246:18–68.
Adams L, Dale BW, Roffler GH. Extraordinary movements of the Denali caribou herd following the perfect storm. Rangifer. 2005;25(4):19–25. https://doi.org/10.7557/18.104.22.1686.
Le Corre M, Dussault C, Côté SD. Weather conditions and variation in timing of spring and fall migrations of migratory caribou. J Mammal. 2016;98(1):260–71. https://doi.org/10.1093/jmammal/gyw177.
Collins WB, Dale BW, Adams LG, McElwain DE, Joly K. Fire, grazing history, lichen abundance, and winter distribution of caribou in Alaska’s taiga. J Wildl Manag. 2011;75(2):369–77. https://doi.org/10.1002/jwmg.39.
Le Corre M, Dussault C, Côté SD. Where to spend the winter? The role of intraspecific competition and climate in determining the selection of wintering areas by migratory caribou. Oikos. 2020;129(4):512–25. https://doi.org/10.1111/oik.06668.
Cooper EJ, Wookey PA. Field measurements of the growth rates of forage lichens, and the implications of grazing by Svalbard reindeer. Symbiosis. 2001;31(1–3):173–86. http://hdl.handle.net/10222/77837.
Adams LG, Stephenson RO, Dale BW, Ahgook RT, Demma DJ. Population dynamics and harvest characteristics of wolves in the Central Brooks Range, Alaska. Wildl Monogr. 2008;170(1):1–25. https://doi.org/10.2193/2008-012.
Whitten KR, Garner GW, Mauer FJ, Harris RB. Productivity and early calf survival in the Porcupine caribou herd. J Wildl Manag. 1992;56(2):201–12. https://doi.org/10.2307/3808814.
Klimstra R. Wolf management report and plan, Game Management Unit 26A: Report period 1 July 2010–30 June 2015, and plan period 1 July 2015–30 June 2020. Species Management Report and Plan ADF&G/DWC/SMR&P-2020-12. Juneau: Alaska Department of Fish and Game; 2020.
Joly K, Bente P, Dau J. Response of overwintering caribou to burned habitat in Northwest Alaska. Arctic. 2007;60(4):401–10. https://www.jstor.org/stable/40512963.
Gustine DD, Parker KL, Lay RJ, Gillingham MP, Heard DC. Interpreting resource selection at different scales for woodland caribou in winter. J Wildl Manag. 2006;70(6):1601–14. https://doi.org/10.2193/0022-541X(2006)70[1601:IRSADS]2.0.CO;2.
Prichard AK, Yokel DA, Rea CL, Person BT, Parrett LS. The effect of frequency of telemetry locations on movement-rate calculations in arctic caribou. Wildl Soc Bull. 2014;38(1):78–88. https://doi.org/10.1002/wsb.357.
Henshaw J. The activities of the wintering caribou in northwestern Alaska in relation to weather and snow conditions. Int J Biometeorol. 1968;12(1):21–7. https://doi.org/10.1007/BF01552974.
Klein DR, Meldgaard M, Fancy SG. Factors determining leg length in Rangifer tarandus. J Mammal. 1987;68(3):642–55. https://doi.org/10.2307/1381597.
Mellor M. A review of basic snow mechanics. Hanover: US Army Cold Regions Research and Engineering Laboratory; 1974:42.
Pomeroy JW, Brun E. Physical properties of snow. In: Jones HG, et al., editors. Snow ecology: an interdisciplinary examination of snow-covered ecosystems. Cambridge: Cambridge University Press; 2001:45–118.
Nawrocki TW, Carlson ML, Osnas JLD, Trammell EJ, Witmer FDW. Regional mapping of species-level continuous foliar cover: beyond categorical vegetation mapping. Ecol Appl. 2020;30(4):e02081. https://doi.org/10.1002/eap.2081.
Beumer LT, Varpea Ø, Hansen BB. Cratering behaviour and faecal C: N ratio in relation to seasonal snowpack characteristics in a high-arctic ungulate. Polar Res. 2017;36(1):1–12. https://doi.org/10.1080/17518369.2017.1286121.
Beumer LT, Pohle J, Schmidt NM, Chimienti M, Desforges J-P, Hansen LH, et al. An application of upscaled optimal foraging theory using hidden Markov modelling: year-round behavioural variation in a large arctic herbivore. Mov Ecol. 2020;8(1):25. https://doi.org/10.1186/s40462-020-00213-x.
Ellington EH, Muntz EM, Gehrt SD. Seasonal and daily shifts in behavior and resource selection: how a carnivore navigates costly landscapes. Oecologia. 2020;194(1):87–100. https://doi.org/10.1007/s00442-020-04754-1.
Pedersen SH, Liston GE, Tamstorf MP, Westergaard-Nielsen A, Schmidt NM. Quantifying episodic snowmelt events in arctic ecosystems. Ecosystems. 2015;18(5):839–56. https://doi.org/10.1007/s10021-015-9867-8.
Bieniek PA, Bhatt US, Walsh JE, Lader R, Griffith B, Roach JK, et al. Assessment of Alaska rain-on-snow events using dynamical downscaling. J Appl Meteorol Climatol. 2018;57(8):1847–63. https://doi.org/10.1175/JAMC-D-17-0276.1.
Hansen BB, Lorentzen JR, Welker JM, Varpe Ø, Aanes R, Beumer LT, et al. Reindeer turning maritime: Ice-locked tundra triggers changes in dietary niche utilization. Ecosphere. 2019;10(4):e02672. https://doi.org/10.1002/ecs2.2672.
Hansen BB, Aanes R, Sæther B-E. Feeding-crater selection by high-arctic reindeer facing ice-blocked pastures. Can J Zool. 2010;88(2):170–7. https://doi.org/10.1139/Z09-130.
Paoli A, Weladji RB, Holand Ø, Kumpula J. Winter and spring climatic conditions influence timing and synchrony of calving in reindeer. PLoS One. 2018;13(4):e0195603. https://doi.org/10.1371/journal.pone.0195603.
Gurarie E, Hebblewhite M, Joly K, Kelly AP, Adamczewski J, Davidson SC, et al. Tactical departures and strategic arrivals: divergent effects of climate and weather on caribou spring migrations. Ecosphere. 2019;10(12):e02971. https://doi.org/10.1002/ecs2.2971.
Cameron RD, Smith WT, Fancy SG, Gerhart KL, White RG. Calving success of female caribou in relation to body weight. Can J Zool. 1993;71(3):480–6. https://doi.org/10.1139/z93-069.
Adams LG, Dale BW. Timing and synchrony of parturition in Alaskan caribou. J Mammal. 1998;79(1):287–94. https://doi.org/10.2307/1382865.
Skogland T. The effects of density dependent resource limitation on size of wild reindeer. Oecologia. 1983;60(2):156–68. https://doi.org/10.1007/BF00379517.
Adams LG, Dale BW. Reproductive performance of female Alaskan caribou. J Wildl Manag. 1998;62(4):1184–95. https://doi.org/10.2307/3801982.
Albon SD, Irvine RJ, Halvorsen O, Langvatn R, Loe LE, Ropstad E, et al. Contrasting effects of summer and winter warming on body mass explain population dynamics in a food-limited Arctic herbivore. Glob Chang Biol. 2017;23(4):1374–89. https://doi.org/10.1111/gcb.13435.
Dolant C, Montpetit B, Langlois A, Brucker L, Zolina O, Johnson CA, et al. Assessment of the barren ground caribou die-off during winter 2015–2016 using passive microwave observations. Geophys Res Lett. 2018;45(10):4908–16. https://doi.org/10.1029/2017GL076752.
Loe LE, Liston GE, Pigeon G, Barker K, Horvitz N, Stien A, et al. The neglected season: warmer autumns counteract harsher winters and promote population growth in Arctic reindeer. Glob Chang Biol. 2021;27(5):993–1002. https://doi.org/10.1111/gcb.15458.
Mech LD, Meier TJ, Burch JW, Adams LG. Patterns of prey selection by wolves in Denali National Park, Alaska. Ecol Conserv wolves a Chang world Proc Second North Am Symp wolves Occas Publ Can Circumpolar Institute, Edmonton, Canada. 1995;35:231–43.
Ferguson MAD, Messier F. Collection and analysis of traditional ecological knowledge about a population of Arctic tundra caribou. Arctic. 1997;50(1):17–28. https://www.jstor.org/stable/40512038.
Kruse JA, White RG, Epstein HE, Archie B, Berman M, Braund SR, et al. Modeling sustainability of Arctic communities: an interdisciplinary collaboration of researchers and local knowledge holders. Ecosystems. 2004;7(8):815–28. https://doi.org/10.1007/s10021-004-0008-z.
We thank Kaj Lynöe, Carrie Vuyovich, Luke Metherell, Pertti Ala-aho, and Dirk Nickisch for their assistance and companionship during fieldwork 2018 and 2019; Matthew Sturm and Charlie Parr for companionship and discussions during fieldwork 2018 and 2019, and generously sharing field snow observations and teaching Stine Højlund Pedersen to navigate and lead expeditions in this massive landscape. We are grateful to Coyote Air, staff at Toolik Field Station of University of Alaska Fairbanks, and Polar Fields staff for providing logistic support. Support for and access to caribou GPS location data was provided by Alaska Department of Fish and Game; Chrissy Pohl and Robyn McGhee at ConocoPhillips Alaska, Inc.; and Emily Sappington at ExxonMobil Alaska Production. Anaktuvuk Pass city council, Gates of Arctic National Park and Preserve, Arctic National Wildlife Refuge, State of Alaska Department of Natural Resources, and U.S. Department of the Interior Bureau of Land Management granted us research and land access permits. Reviews by Robert G. White, Donald E. Spalinger, Cleve Cowles, John W. Merickel, and two anonymous reviewers improved this manuscript. We recognize and appreciate access to traditional lands of Iñupiat, Nunamiut, Gwich’in, and Koyukon Athabascans to conduct this research.
This research was funded by United States National Science Foundation grants 1604249, 1602898, and 1839198; National Aeronautics and Space Administration grant 80NSSC18K0571; and JMW’s UArctic Research Chair position.
Ethics approval and consent to participate
Procedures for handling live animals adhered to the guidelines of the American Society of Mammalogists and were approved by the Alaska Department of Fish and Game Animal Care and Use Committee (Permits held by ADFG Division of wildlife conservation 04-012, 09-03, 11-013, 2013-17, 2015-06, 2016-30, 0019-2017-19, 0019-2018-49, 0019-2019-44).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Animation showing daily, modeled snow depth (background map shading) produced using SnowModel from 1 September 2018–31 May 2019 (dates will iterate at the top of the animation) overlaid by daily Central Arctic Herd caribou GPS-locations (white circles).
Additional file 2:
R scripts to reproduce Fig. 4, Fig. 6, and Fig. 7 and the movement modeling using integrated Step Selection Analysis.
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.
About this article
Cite this article
Pedersen, S.H., Bentzen, T.W., Reinking, A.K. et al. Quantifying effects of snow depth on caribou winter range selection and movement in Arctic Alaska. Mov Ecol 9, 48 (2021). https://doi.org/10.1186/s40462-021-00276-4
- Arctic Alaska
- Barren-ground caribou
- Movement ecology
- Rangifer tarandus
- Integrated step selection analysis
- Snow depth
- Winter range