Long-term patterns of muskox (Ovibos moschatus) demographics in high arctic Greenland

Population genomics of the muskox' resilience in the near absence of genetic variation

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Geochemical landscapes as drivers of wildlife reproductive success: Insights from a high-Arctic ecosystem

The bioavailability of essential and non-essential elements in vegetation is expected to influence the performance of free-ranging terrestrial herbivores. However, attempts to relate the use of geochemical landscapes by animal populations directly to reproductive output are currently lacking. Here we measured concentrations of 14 essential and non-essential elements in soil and vegetation samples collected in the Zackenberg valley, northeast Greenland, and linked these to environmental conditions to spatially predict and map geochemical landscapes. We then used long-term (1996-2021) survey data of muskoxen (Ovibos moschatus) to quantify annual variation in the relative use of essential and non-essential elements in vegetated sites and their relationship to calf recruitment the following year. Results showed that the relative use of the geochemical landscape by muskoxen varied substantially between years and differed among elements. Selection for vegetated sites with higher levels of the essential elements N, Cu, Se, and Mo was positively linked to annual calf recruitment. In contrast, selection for vegetated sites with higher concentrations of the non-essential elements As and Pb was negatively correlated to annual calf recruitment. Based on the concentrations measured in our study, we found no apparent associations between annual calf recruitment and levels of C, Mn, Co, Zn, Cd, Ba, Hg, and C:N ratio in the vegetation. We conclude that the spatial distribution and access to essential and non-essential elements are important drivers of reproductive output in muskoxen, which may also apply to other wildlife populations. The value of geochemical landscapes to assess habitat-performance relationships is likely to increase under future environmental change.

Sequential analysis of δ15 N in guard hair suggests late gestation is the most critical period for muskox calf recruitment

RATIONALE

Analysis of stable isotopes in tissue and excreta may provide information about animal diets and their nutritional state. As body condition may have a major influence on reproduction, linking stable isotope values to animal demographic rates may help unravel the drivers behind animal population dynamics.

METHODS

We performed sequential analysis of δ¹⁵ N values in guard hair from 21 muskoxen (Ovibos moschatus) from Zackenberg in high arctic Greenland. We were able to reconstruct the dietary history for the population over a 5-year period with contrasting environmental conditions. We examined the linkage between guard hair δ¹⁵ N values in 12 three-month periods and muskox calf recruitment to detect critical periods for muskox reproduction. Finally, we conducted similar analyses of the correlation between environmental conditions (snow depth and air temperature) and calf recruitment.

RESULTS

δ¹⁵ N values exhibited a clear seasonal pattern with high levels in summer and low levels in winter. However, large inter-annual variation was found in winter values, suggesting varying levels of catabolism depending on snow conditions. In particular δ¹⁵ N values during January-March were linked to muskox recruitment rates, with higher values coinciding with lower calf recruitment. δ¹⁵ N values were a better predictor of muskox recruitment rates than environmental conditions.

CONCLUSIONS

Although environmental conditions may ultimately determine the dietary δ¹⁵ N signal in muskox guard hairs, muskox calf recruitment was more strongly correlated with δ¹⁵ N values than ambient snow and temperature. The period January-March, corresponding to late gestation, appears particularly critical for muskox reproduction.

Herbivores in Arctic ecosystems: Effects of climate change and implications for carbon and nutrient cycling

Arctic terrestrial herbivores influence tundra carbon and nutrient dynamics through their consumption of resources, waste production, and habitat-modifying behaviors. The strength of these effects is likely to change spatially and temporally as climate change drives shifts in herbivore abundance, distribution, and activity timing. Here, we review how herbivores influence tundra carbon and nutrient dynamics through their consumptive and nonconsumptive effects. We also present evidence for herbivore responses to climate change and discuss how these responses may alter the spatial and temporal distribution of herbivore impacts. Several current knowledge gaps limit our understanding of the changing functional roles of herbivores; these include limited characterization of the spatial and temporal variability in herbivore impacts and of how herbivore activities influence the cycling of elements beyond carbon. We conclude by highlighting approaches that will promote better understanding of herbivore effects on tundra ecosystems, including their integration into existing biogeochemical models, new applications of remote sensing techniques, and the continued use of distributed experiments.

Muskoxen homogenise soil microbial communities and affect the abundance of methanogens and methanotrophs

Grazing herbivores may affect soil microbial communities indirectly by impacting soil structure and vegetation composition. In high arctic environments, this impact is poorly elucidated, while having potentially wide-reaching effects on the ecosystem. This study examines how a key arctic herbivore, the muskox Ovibos moschatus, affects the soil microbial community in a high arctic fen. Environmental DNA was extracted from soil samples taken from grazed control plots and from muskox exclosures established 5 years prior. We sequenced amplicons of the 16S rRNA gene to provide insight into the microbial communities. We found that in the grazed control plots, microbial communities exhibited high evenness and displayed highly similar overall diversity. In plots where muskoxen had been excluded, microbial diversity was significantly reduced, and had more uneven intra-sample populations and overall lower ecological richness and evenness. We observed that the composition of microbial communities in grazed soils were significantly affected by the presence of muskoxen, as seen by elevated relative abundances of Bacteroides and Firmicutes, two major phyla found in muskox faeces. Furthermore, an increase in relative abundance of bacteria involved in degradation of recalcitrant carbohydrates and cycling of nitrogen was observed in grazed soil. Ungrazed soils displayed increased abundances of bacteria potentially involved in anaerobic oxidation of methane, whereas some methanogens were more abundant in grazed soils. This corroborates a previous finding that methane emissions are higher in arctic fens under muskox grazing. Our results show that the presence of large herbivores stimulates soil microbial diversity and has a homogenizing influence on the inter-species dynamics in soil microbial communities. The findings of this study, thus, improve our understanding of the effect of herbivore grazing on arctic ecosystems and the derived methane cycling.

Environment and physiology shape Arctic ungulate population dynamics

Species conservation in a rapidly changing world requires an improved understanding of how individuals and populations respond to changes in their environment across temporal scales. Increased warming in the Arctic puts this region at particular risk for rapid environmental change, with potentially devastating impacts on resident populations. Here, we make use of a parameterized full life cycle, individual-based energy budget model for wild muskoxen, coupling year-round environmental data with detailed ontogenic metabolic physiology. We show how winter food accessibility, summer food availability, and density dependence drive seasonal dynamics of energy storage and thus life history and population dynamics. Winter forage accessibility defined by snow depth, more than summer forage availability, was the primary determinant of muskox population dynamics through impacts on calf recruitment and longer term carryover effects of maternal investment. Simulations of various seasonal snow depth and plant biomass and quality profiles revealed that timing of and improved/limited winter forage accessibility had marked influence on calf recruitment (±10-80%). Impacts on recruitment were the cumulative result of condition-driven reproductive performance at multiple time points across the reproductive period (ovulation to calf weaning) as a trade-off between survival and reproduction. Seasonal and generational condition effects of snow-rich winters interacted with age structure and density to cause pronounced long-term consequences on population growth and structure, with predicted population recovery times from even moderate disturbances of 10 years or more. Our results show how alteration in winter forage accessibility, mediated by snow depth, impacts the dynamics of northern herbivore populations. Further, we present here a mechanistic and state-based model framework to assess future scenarios of environmental change, such as increased or decreased snowfall or plant biomass and quality to impact winter and summer forage availability across the Arctic.

Quantifying energetic and fitness consequences of seasonal heterothermy in an Arctic ungulate

Animals have adapted behavioral and physiological strategies to conserve energy during periods of adverse conditions. Heterothermy is one such adaptation used by endotherms. While heterothermy-fluctuations in body temperature and metabolic rate-has been shown in large vertebrates, little is known of the costs and benefits of this strategy, both in terms of energy and in terms of fitness. Hence, our objective was to model the energetics of seasonal heterothermy in the largest Arctic ungulate, the muskox (Ovibos moschatus), using an individual-based energy budget model of metabolic physiology. We found that the empirically based drop in body temperature (winter max ~-0.8°C) overwinter in adult females resulted in substantial fitness benefits in terms of reduced daily energy expenditure and body mass loss. Body mass and energy reserves were 8.98% and 14.46% greater in modeled heterotherms compared to normotherms by end of winter. Based on environmental simulations, we show that seasonal heterothermy can, to some extent, buffer the negative consequences of poor prewinter body condition or reduced winter food accessibility, leading to greater winter survival (+20%-30%) and spring energy reserves (+10%-30%), and thus increased probability of future reproductive success. These results indicate substantial adaptive short-term benefits of seasonal heterothermy at the individual level, with potential implications for long-term population dynamics in highly seasonal environments.

Environmental conditions alter behavioural organization and rhythmicity of a large Arctic ruminant across the annual cycle

The existence and persistence of rhythmicity in animal activity during phases of environmental change is of interest in ecology, evolution and chronobiology. A wide diversity of biological rhythms in response to exogenous conditions and internal stimuli have been uncovered, especially for polar vertebrates. However, empirical data supporting circadian organization in behaviour of large ruminating herbivores remains inconclusive. Using year-round tracking data of the largest Arctic ruminant, the muskox (Ovibos moschatus), we modelled rhythmicity as a function of behaviour and environmental conditions. Behavioural states were classified based on patterns in hourly movements, and incorporated within a periodicity analyses framework. Although circadian rhythmicity in muskox behaviour was detected throughout the year, ultradian rhythmicity was most prevalent, especially when muskoxen were foraging and resting in mid-winter (continuous darkness). However, when combining circadian and ultradian rhythmicity together, the probability of behavioural rhythmicity declined with increasing photoperiod until largely disrupted in mid-summer (continuous light). Individuals that remained behaviourally rhythmic during mid-summer foraged in areas with lower plant productivity (NDVI) than individuals with arrhythmic behaviour. Based on our study, we conclude that muskoxen may use an interval timer to schedule their behavioural cycles when forage resources are low, but that the importance and duration of this timer are reduced once environmental conditions allow energetic reserves to be replenished ad libitum. We argue that alimentary function and metabolic requirements are critical determinants of biological rhythmicity in muskoxen, which probably applies to ruminating herbivores in general.

An application of upscaled optimal foraging theory using hidden Markov modelling: year-round behavioural variation in a large arctic herbivore

BACKGROUND

In highly seasonal environments, animals face critical decisions regarding time allocation, diet optimisation, and habitat use. In the Arctic, the short summers are crucial for replenishing body reserves, while low food availability and increased energetic demands characterise the long winters (9-10 months). Under such extreme seasonal variability, even small deviations from optimal time allocation can markedly impact individuals' condition, reproductive success and survival. We investigated which environmental conditions influenced daily, seasonal, and interannual variation in time allocation in high-arctic muskoxen (Ovibos moschatus) and evaluated whether results support qualitative predictions derived from upscaled optimal foraging theory.

METHODS

Using hidden Markov models (HMMs), we inferred behavioural states (foraging, resting, relocating) from hourly positions of GPS-collared females tracked in northeast Greenland (28 muskox-years). To relate behavioural variation to environmental conditions, we considered a wide range of spatially and/or temporally explicit covariates in the HMMs.

RESULTS

While we found little interannual variation, daily and seasonal time allocation varied markedly. Scheduling of daily activities was distinct throughout the year except for the period of continuous daylight. During summer, muskoxen spent about 69% of time foraging and 19% resting, without environmental constraints on foraging activity. During winter, time spent foraging decreased to 45%, whereas about 43% of time was spent resting, mediated by longer resting bouts than during summer.

CONCLUSIONS

Our results clearly indicate that female muskoxen follow an energy intake maximisation strategy during the arctic summer. During winter, our results were not easily reconcilable with just one dominant foraging strategy. The overall reduction in activity likely reflects higher time requirements for rumination in response to the reduction of forage quality (supporting an energy intake maximisation strategy). However, deep snow and low temperatures were apparent constraints to winter foraging, hence also suggesting attempts to conserve energy (net energy maximisation strategy). Our approach provides new insights into the year-round behavioural strategies of the largest Arctic herbivore and outlines a practical example of how to approximate qualitative predictions of upscaled optimal foraging theory using multi-year GPS tracking data.

Muskox status, recent variation, and uncertain future

Muskoxen (Ovibos moschatus) are an integral component of Arctic biodiversity. Given low genetic diversity, their ability to respond to future and rapid Arctic change is unknown, although paleontological history demonstrates adaptability within limits. We discuss status and limitations of current monitoring, and summarize circumpolar status and recent variations, delineating all 55 endemic or translocated populations. Acknowledging uncertainties, global abundance is ca 170 000 muskoxen. Not all populations are thriving. Six populations are in decline, and as recently as the turn of the century, one of these was the largest population in the world, equaling ca 41% of today's total abundance. Climate, diseases, and anthropogenic changes are likely the principal drivers of muskox population change and result in multiple stressors that vary temporally and spatially. Impacts to muskoxen are precipitated by habitat loss/degradation, altered vegetation and species associations, pollution, and harvest. Which elements are relevant for a specific population will vary, as will their cumulative interactions. Our summaries highlight the importance of harmonizing existing data, intensifying long-term monitoring efforts including demographics and health assessments, standardizing and implementing monitoring protocols, and increasing stakeholder engagement/contributions.

On the interplay between hypothermia and reproduction in a high arctic ungulate

For free-ranging animals living in seasonal environments, hypometabolism (lowered metabolic rate) and hypothermia (lowered body temperature) can be effective physiological strategies to conserve energy when forage resources are low. To what extent such strategies are adopted by large mammals living under extreme conditions, as those encountered in the high Arctic, is largely unknown, especially for species where the gestation period overlaps with the period of lowest resource availability (i.e. winter). Here we investigated for the first time the level to which high arctic muskoxen (Ovibos moschatus) adopt hypothermia and tested the hypothesis that individual plasticity in the use of hypothermia depends on reproductive status. We measured core body temperature over most of the gestation period in both free-ranging muskox females in Greenland and captive female muskoxen in Alaska. We found divergent overwintering strategies according to reproductive status, where pregnant females maintained stable body temperatures during winter, while non-pregnant females exhibited a temporary decrease in their winter body temperature. These results show that muskox females use hypothermia during periods of resource scarcity, but also that the use of this strategy may be limited to non-reproducing females. Our findings suggest a trade-off between metabolically-driven energy conservation during winter and sustaining foetal growth, which may also apply to other large herbivores living in highly seasonal environments elsewhere.

Density, snow, and seasonality lead to variation in muskox (Ovibos moschatus) habitat selection during summer

Understanding how environmental conditions influence habitat selection and suitability of free-ranging animals is critical, as the outcome may have implications for individual fitness and population dynamics. Density and snow are among the most influential environmental conditions driving habitat-selection patterns of northern ungulates. We used two decades of census data from high Arctic Greenland to quantify inter- and intra-annual variations in muskox (Ovibos moschatus (Zimmermann, 1780)) habitat selection and suitability during the Arctic summer (July through October). Across years, habitat selection varied considerably, and the strength of habitat selection appeared negatively related to both muskox density and spring snow cover. In early summer, habitat suitability was high and spatially rather uniform. Towards the autumn, suitable habitats contracted to just the lower elevations, when muskoxen exhibited increasingly stronger habitat selection towards low elevations and dense vegetation. This selection strategy clearly reflects the need to build up fat reserves for the upcoming winter, highlighting the energetic importance of the Arctic summer. Extreme climatic events such as freezing rain in autumn are increasing in frequency in Greenland and limit muskox access to high-quality forage in fens. Such events may therefore negatively affect the energy acquisition process of muskox with potential cascading consequences on population dynamics.

An ecosystem-wide reproductive failure with more snow in the Arctic

2018: Arctic researchers have just witnessed another extreme summer—but in a new sense of the word. Although public interest has long been focused on general warming trends and trends towards a lower sea ice cover in the Arctic Ocean, this summer saw the realization of another predicted trend: that of increasing precipitation during the winter months and of increased year-to-year variability. In a well-studied ecosystem in Northeast Greenland, this resulted in the most complete reproductive failure encountered in the terrestrial ecosystem during more than two decades of monitoring: only a few animals and plants were able to reproduce because of abundant and late melting snow. These observations, we suggest, should open our eyes to potentially drastic consequences of predicted changes in both the mean and the variability of arctic climate.

In summer 2018, many arctic regions experienced unprecedented large amounts of snow. This Perspective article reports the ecological impacts of this extreme event, and shows that the severe snow conditions resulted in an almost complete reproductive failure across the entire ecosystem. The 2018 may be a rare event, but it also offers a gloomy peep into the future.

High resistance towards herbivore-induced habitat change in a high Arctic arthropod community

Mammal herbivores may exert strong impacts on plant communities, and are often key drivers of vegetation composition and diversity. We tested whether such mammal-induced changes to a high Arctic plant community are reflected in the structure of other trophic levels. Specifically, we tested whether substantial vegetation changes following the experimental exclusion of muskoxen (Ovibos moschatus) altered the composition of the arthropod community and the predator-prey interactions therein. Overall, we found no impact of muskox exclusion on the arthropod community: the diversity and abundance of both arthropod predators (spiders) and of their prey were unaffected by muskox presence, and so was the qualitative and quantitative structure of predator-prey interactions. Hence, high Arctic arthropod communities seem highly resistant towards even large biotic changes in their habitat, which we attribute to the high connectance in the food web.

Interaction webs in arctic ecosystems: Determinants of arctic change?

How species interact modulate their dynamics, their response to environmental change, and ultimately the functioning and stability of entire communities. Work conducted at Zackenberg, Northeast Greenland, has changed our view on how networks of arctic biotic interactions are structured, how they vary in time, and how they are changing with current environmental change: firstly, the high arctic interaction webs are much more complex than previously envisaged, and with a structure mainly dictated by its arthropod component. Secondly, the dynamics of species within these webs reflect changes in environmental conditions. Thirdly, biotic interactions within a trophic level may affect other trophic levels, in some cases ultimately affecting land-atmosphere feedbacks. Finally, differential responses to environmental change may decouple interacting species. These insights form Zackenberg emphasize that the combination of long-term, ecosystem-based monitoring, and targeted research projects offers the most fruitful basis for understanding and predicting the future of arctic ecosystems.

Show Me Your Rump Hair and I Will Tell You What You Ate - The Dietary History of Muskoxen (Ovibos moschatus) Revealed by Sequential Stable Isotope Analysis of Guard Hairs

The nutritional state of animals is tightly linked to the ambient environment, and for northern ungulates the state strongly influences vital population demographics, such as pregnancy rates. Continuously growing tissues, such as hair, can be viewed as dietary records of animals over longer temporal scales. Using sequential data on nitrogen stable isotopes (δ15N) in muskox guard hairs from ten individuals in high arctic Northeast Greenland, we were able to reconstruct the dietary history of muskoxen over approximately 2.5 years with a high temporal resolution of app. 9 days. The dietary chronology included almost three full summer and winter periods. The diet showed strong intra- and inter-annual seasonality, and was significantly linked to changes in local environmental conditions (temperature and snow depth). The summer diets were highly similar across years, reflecting a graminoid-dominated diet. In contrast, winter diets were markedly different between years, a pattern apparently linked to snow conditions. Snow-rich winters had markedly higher δ15N values than snow-poor winters, indicating that muskoxen had limited access to forage, and relied more heavily on their body stores. Due to the close link between body stores and calf production in northern ungulates, the dietary winter signals could eventually serve as an indicator of calf production the following spring. Our study opens the field for further studies and longer chronologies to test such links. The method of sequential stable isotope analysis of guard hairs thus constitutes a promising candidate for population-level monitoring of animals in remote, arctic areas.