Winter snow and spring temperature have differential effects on vegetation phenology and productivity across Arctic plant communities.

Tundra dominates two-thirds of the unglaciated, terrestrial Arctic. Although this region has experienced rapid and widespread changes in vegetation phenology and productivity over the last several decades, the specific climatic drivers responsible for this change remain poorly understood. Here we quantified the effect of winter snowpack and early spring temperature conditions on growing season vegetation phenology (timing of the start, peak, and end of the growing season) and productivity of the dominant tundra vegetation communities of Arctic Alaska. We used daily remotely sensed normalized difference vegetation index (NDVI), and daily snowpack and temperature variables produced by SnowModel and MicroMet, coupled physically based snow and meteorological modeling tools, to (1) determine the most important snowpack and thermal controls on tundra vegetation phenology and productivity and (2) describe the direction of these relationships within each vegetation community. Our results show that soil temperature under the snowpack, snowmelt timing, and air temperature following snowmelt are the most important drivers of growing season timing and productivity among Arctic vegetation communities. Air temperature after snowmelt was the most important control on timing of season start and end, with warmer conditions contributing to earlier phenology in all vegetation communities. In contrast, the controls on the timing of peak season and productivity also included snowmelt timing and soil temperature under the snowpack, dictated in part by the snow insulating capacity. The results of this novel analysis suggest that while future warming effects on phenology may be consistent across communities of the tundra biome, warming may result in divergent, community-specific productivity responses if coupled with reduced snow insulating capacity lowers winter soil temperature and potential nutrient cycling in the soil.

Spatiotemporally variable snow properties drive habitat use of an Arctic mesopredator.

Climate change is rapidly altering the composition and availability of snow, with implications for snow-affected ecological processes, including reproduction, predation, habitat selection, and migration. How snowpack changes influence these ecological processes is mediated by physical snowpack properties, such as depth, density, hardness, and strength, each of which is in turn affected by climate change. Despite this, it remains difficult to obtain meaningful snow information relevant to the ecological processes of interest, precluding a mechanistic understanding of these effects. This problem is acute for species that rely on particular attributes of the subnivean space, for example depth, thermal resistance, and structural stability, for key life-history processes like reproduction, thermoregulation, and predation avoidance. We used a spatially explicit snow evolution model to investigate how habitat selection of a species that uses the subnivean space, the wolverine, is related to snow depth, snow density, and snow melt on Arctic tundra. We modeled these snow properties at a 10 m spatial and a daily temporal resolution for 3 years, and used integrated step selection analyses of GPS collar data from 21 wolverines to determine how these snow properties influenced habitat selection and movement. We found that wolverines selected deeper, denser snow, but only when it was not undergoing melt, bolstering the evidence that these snow properties are important to species that use the Arctic snowpack for subnivean resting sites and dens. We discuss the implications of these findings in the context of climate change impacts on subnivean species.

Navigating snowscapes: scale-dependent responses of mountain sheep to snowpack properties.

Winters are limiting for many terrestrial animals due to energy deficits brought on by resource scarcity and the increased metabolic costs of thermoregulation and traveling through snow. A better understanding of how animals respond to snow conditions is needed to predict the impacts of climate change on wildlife. We compared the performance of remotely sensed and modeled snow products as predictors of winter movements at multiple spatial and temporal scales using a data set of 20,544 locations from 30 GPS-collared Dall sheep (Ovis dalli dalli) in Lake Clark National Park and Preserve, Alaska, USA from 2005 to 2008. We used daily 500-m MODIS normalized difference snow index (NDSI), and multi-resolution snow depth and density outputs from a snowpack evolution model (SnowModel), as covariates in step selection functions. We predicted that modeled snow depth would perform best across all scales of selection due to more informative spatiotemporal variation and relevance to animal movement. Our results indicated that adding any of the evaluated snow metrics substantially improved model performance and helped characterize winter Dall sheep movements. As expected, SnowModel-simulated snow depth outperformed NDSI at fine-to-moderate scales of selection (step scales < 112 h). At the finest scale, Dall sheep selected for snow depths below mean chest height (<54 cm) when in low-density snows (100 kg/m3 ), which may have facilitated access to ground forage and reduced energy expenditure while traveling. However, sheep selected for higher snow densities (>300 kg/m3 ) at snow depths above chest height, which likely further reduced energy expenditure by limiting hoof penetration in deeper snows. At moderate-to-coarse scales (112-896 h step scales), however, NDSI was the best-performing snow covariate. Thus, the use of publicly available, remotely sensed, snow cover products can substantially improve models of animal movement, particularly in cases where movement distances exceed the MODIS 500-m grid threshold. However, remote sensing products may require substantial data thinning due to cloud cover, potentially limiting its power in cases where complex models are necessary. Snowpack evolution models such as SnowModel offer users increased flexibility at the expense of added complexity, but can provide critical insights into fine-scale responses to rapidly changing snow properties.

Quantifying effects of snow depth on caribou winter range selection and movement in Arctic Alaska.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

A Lagrangian Snow-Evolution System for Sea-Ice Applications (SnowModel-LG): Part I-Model Description.

A Lagrangian snow-evolution model (SnowModel-LG) was used to produce daily, pan-Arctic, snow-on-sea-ice, snow property distributions on a 25 × 25-km grid, from 1 August 1980 through 31 July 2018 (38 years). The model was forced with NASA's Modern Era Retrospective-Analysis for Research and Applications-Version 2 (MERRA-2) and European Centre for Medium-Range Weather Forecasts (ECMWF) ReAnalysis-5th Generation (ERA5) atmospheric reanalyses, and National Snow and Ice Data Center (NSIDC) sea ice parcel concentration and trajectory data sets (approximately 61,000, 14 × 14-km parcels). The simulations performed full surface and internal energy and mass balances within a multilayer snowpack evolution system. Processes and features accounted for included rainfall, snowfall, sublimation from static-surfaces and blowing-snow, snow melt, snow density evolution, snow temperature profiles, energy and mass transfers within the snowpack, superimposed ice, and ice dynamics. The simulations produced horizontal snow spatial structures that likely exist in the natural system but have not been revealed in previous studies spanning these spatial and temporal domains. Blowing-snow sublimation made a significant contribution to the snowpack mass budget. The superimposed ice layer was minimal and decreased over the last four decades. Snow carryover to the next accumulation season was minimal and sensitive to the melt-season atmospheric forcing (e.g., the average summer melt period was 3 weeks or 50% longer with ERA5 forcing than MERRA-2 forcing). Observed ice dynamics controlled the ice parcel age (in days), and ice age exerted a first-order control on snow property evolution.