Nonstationary spatio-temporal small rodent dynamics: evidence from long-term Norwegian fox bounty data

Taking the beat of the Arctic: are lemming population cycles changing due to winter climate?

Reports of fading vole and lemming population cycles and persisting low populations in some parts of the Arctic have raised concerns about the spread of these fundamental changes to tundra food web dynamics. By compiling 24 unique time series of lemming population fluctuations across the circumpolar region, we show that virtually all populations displayed alternating periods of cyclic/non-cyclic fluctuations over the past four decades. Cyclic patterns were detected 55% of the time (n = 649 years pooled across sites) with a median periodicity of 3.7 years, and non-cyclic periods were not more frequent in recent years. Overall, there was an indication for a negative effect of warm spells occurring during the snow onset period of the preceding year on lemming abundance. However, winter duration or early winter climatic conditions did not differ on average between cyclic and non-cyclic periods. Analysis of the time series shows that there is presently no Arctic-wide collapse of lemming cycles, even though cycles have been sporadic at most sites during the last decades. Although non-stationary dynamics appears a common feature of lemming populations also in the past, continued warming in early winter may decrease the frequency of periodic irruptions with negative consequences for tundra ecosystems.

Spatiotemporal patterns of red fox scavenging in forest and tundra: the influence of prey fluctuations and winter conditions

Concern has been raised regarding red fox (Vulpes Vulpes) population increase and range expansion into alpine tundra, directly and indirectly enhanced by human activities, including carrion supply, and its negative impact on native fauna. In this study, we used cameras on bait stations and hunting remains to investigate how spatiotemporal patterns of red fox scavenging were influenced by abundance and accessibility of live prey, i.e., small rodent population cycles, snow depth, and primary productivity. We found contrasting patterns of scavenging between habitats during winter. In alpine areas, use of baits was highest post rodent peaks and when snow depth was low. This probably reflected relatively higher red fox abundance due to increased reproduction or migration of individuals from neighboring areas, possibly also enhanced by a diet shift. Contrastingly, red fox use of baits in the forest was highest during rodent low phase, and when snow was deep, indicating a higher dependency of carrion under these conditions. Scavenging patterns by red fox on the pulsed but predictable food resource from hunting remains in the autumn revealed no patterns throughout the rodent cycle. In this study, we showed that small rodent dynamics influenced red fox scavenging, at least in winter, but with contrasting patterns depending on environmental conditions. In marginal alpine areas, a numerical response to higher availability of rodents possible lead to the increase in bait visitation the proceeding winter, while in more productive forest areas, low availability of rodents induced a functional diet shift towards scavenging.

The delayed effect of cooling reinforced the NAO-plague connection in pre-industrial Europe

Previous studies on the connection between climate and plague were mostly conducted without considering the influence of large-scale atmospheric circulations and long-term historical observations. The current study seeks to reveal the sophisticated role of climatic control on plague by investigating the combined effect of North Atlantic Oscillation (NAO) and temperature on plague outbreaks in Europe from 1347 to 1760 CE. Moving correlation analysis is applied to explore the non-linear relationship between NAO and plague transmission over time. Also, we apply the cross-correlation function to identify the role of temperature in mediating the NAO-plague connection and the lead-lag relationship in between. Our statistical results show that the pathway from climate change to plague incidence is distinctive in its spatial, temporal, and non-linear patterns. The multi-decadal temperature change exerted a 15-22 years lagged impact on the NAO-plague correlation in different European regions. The NAO-plague correlation in Atlantic-Central Europe primarily remained positive, while the correlation in Mediterranean Europe switched between positive and negative alternately. The modulating effect of temperature over the NAO-plague correlation increases exponentially with the magnitude of the temperature anomaly, but the effect is negligible between 0.3 and -0.3 °C anomaly. Our findings show that a lagged influence from the temperature extremes dominantly controls the correlation between NAO and plague incidence. A forecast from our study suggests that large-scale plague outbreaks are unlikely to happen in Europe if NAO remains at its current positive phase during the earth's future warming.

Population cycles and outbreaks of small rodents: ten essential questions we still need to solve

Most small rodent populations in the world have fascinating population dynamics. In the northern hemisphere, voles and lemmings tend to show population cycles with regular fluctuations in numbers. In the southern hemisphere, small rodents tend to have large amplitude outbreaks with less regular intervals. In the light of vast research and debate over almost a century, we here discuss the driving forces of these different rodent population dynamics. We highlight ten questions directly related to the various characteristics of relevant populations and ecosystems that still need to be answered. This overview is not intended as a complete list of questions but rather focuses on the most important issues that are essential for understanding the generality of small rodent population dynamics.

Circumpolar status of Arctic ptarmigan: Population dynamics and trends

Rock ptarmigan (Lagopus muta) and willow ptarmigan (L. lagopus) are Arctic birds with a circumpolar distribution but there is limited knowledge about their status and trends across their circumpolar distribution. Here, we compiled information from 90 ptarmigan study sites from 7 Arctic countries, where almost half of the sites are still monitored. Rock ptarmigan showed an overall negative trend on Iceland and Greenland, while Svalbard and Newfoundland had positive trends, and no significant trends in Alaska. For willow ptarmigan, there was a negative trend in mid-Sweden and eastern Russia, while northern Fennoscandia, North America and Newfoundland had no significant trends. Both species displayed some periods with population cycles (short 3-6 years and long 9-12 years), but cyclicity changed through time for both species. We propose that simple, cost-efficient systematic surveys that capture the main feature of ptarmigan population dynamics can form the basis for citizen science efforts in order to fill knowledge gaps for the many regions that lack systematic ptarmigan monitoring programs.

Relating the 4-year lemming (Lemmus spp. and Dicrostonyx spp.) population cycle to a 3.8-year lunar cycle and ENSO

Reported peak years of lemming (Lemmus spp. and Dicrostonyx spp.) and Arctic fox (Vulpes lagopus (Linnaeus, 1758)) abundance were compiled from the literature for 12 locations spanning 127 years. The mean period of the 34 reported lemming and Arctic fox cycles from 1868 to 1994 was 3.8 years, suggesting that the period of the 4-year cycle is actually 3.8 years. Peak population years were predicted using a simple model based on a 3.8-year lunar cycle. For nearly 130 years, reported years of peak abundance of lemmings and Arctic foxes were significantly correlated with and have persistently stayed in phase with predicted peak years of abundance. Over the same period, predicted peak years of lemming abundance have been closely aligned with peak (i.e., La Niña) years of the January–March Southern Oscillation Index (SOI). From 1952 to 1995, peak flowering in Norway tended to occur close to trough June–August SOI (El Niño) years. The hypothesis proposed is that the 3.8-year lunar cycle governs the timing of the lemming cycle, but it does not cause the population cycling itself. If this hypothesis is true, then the heretofore unexplained source of the persistent periodicity and quasi-metronomic regularity of the lemming cycle is identified.

Population spatial synchrony enhanced by periodicity and low detuning with environmental forcing

Explaining why fluctuations in abundances of spatially disjunct populations often are correlated through time is a major goal of population ecologists. We address two hypotheses receiving little to no testing in wild populations: (i) that population cycling facilitates synchronization given weak coupling among populations, and (ii) that the ability of periodic external forces to synchronize oscillating populations is a function of the mismatch in timescales (detuning) between the force and the population. Here, we apply new analytical methods to field survey data on gypsy moth outbreaks. We report that at timescales associated with gypsy moth outbreaks, spatial synchrony increased with population periodicity via phase locking. The extent to which synchrony in temperature and precipitation influenced population synchrony was associated with the degree of mismatch in dominant timescales of oscillation. Our study provides new empirical methods and rare empirical evidence that population cycling and low detuning can promote population spatial synchrony.

Rodent population cycle as a determinant of gastrointestinal nematode abundance in a low-arctic population of the red fox

We analyzed an 11-year time series (2005-2015) of parasite abundance for three intestinal nematode species in the red fox (Vulpes vulpes) as a function of the multi-annual rodent population cycle in low-arctic Norway, while correcting for other potential covariates that could influence prevalence and abundance. Rodents are paratenic and facultative intermediate hosts for the two Ascarididae species Toxascaris leonina and Toxocara canis, respectively and key prey for the red fox. Still the relative importance of indirect transmission through rodents and direct transmission through free-living stages is unclear. Abundance of these Ascarididae species in individual red foxes (N = 612) exhibited strongly cyclic dynamics that closely mirrored the 4-year rodent cycle. Negative binomial models provided evidence for a direct proportional increase in Ascarididae abundance with rodent density suggesting that predator functional response to rodent prey is the key transmission mechanism. In contrast, no cycles and constantly very low abundance were apparent for Uncinaria stenocephala - a third nematode species recorded without paratenic or intermediate stages.

The alternative prey hypothesis revisited: Still valid for willow ptarmigan population dynamics

The alternative prey hypothesis predicts that the interaction between generalist predators and their main prey is a major driver of population dynamics of alternative prey species. In Fennoscandia, changes in climate and human land use are assumed to alter the dynamics of cyclic small rodents (main prey) and lead to increased densities and range expansion of an important generalist predator, the red fox Vulpes vulpes. In order to better understand the role of these potential changes in community structure on an alternative prey species, willow ptarmigan Lagopus lagopus, we analyzed nine years of population census data from SE Norway to investigate how community interactions affected their population dynamics. The ptarmigan populations showed no declining trend during the study period, and annual variations corresponded with marked periodic small rodent peaks and declines. Population growth and breeding success were highly correlated, and both demographic variables were influenced by an interaction between red fox and small rodents. Red foxes affected ptarmigan negatively only when small rodent abundance was low, which is in accordance with the alternative prey hypothesis. Our results confirm the important role of red fox predation in ptarmigan dynamics, and indicate that if small rodent cycles are disrupted, this may lead to decline in ptarmigan and other alternative prey species due to elevated predation pressure.

Spatial and temporal variation in the range-wide cyclic dynamics of greater sage-grouse

Periodic changes in abundance, or population cycles, are common in a variety of species and is one of the most widely studied ecological phenomena. The strength of, and synchrony between population cycles can vary across time and space and understanding these patterns can provide insight into the mechanisms generating population cycles and their variability within and among species. Here, we used wavelet and spectral analysis on a range-wide dataset of abundance for the greater sage-grouse (Centrocercus urophasianus) to test for regional differences in temporal cyclicity. Overall, we found that most populations (11 of 15) were cyclic at some point in a 50-year time series (1965-2015), but the patterns varied over both time and space. Several peripheral populations demonstrated amplitude dampening or loss of cyclicity following population lows in the mid-1990s. Populations through the core of the range in the Great and Wyoming Basins had more consistent cyclic dynamics, but period length appeared to shorten from 10-12 to 6-8 years. In one time period, where cyclicity was greatest overall, increased pairwise population synchrony was correlated with cycle intensity. Our work represents a comprehensive range-wide assessment of cyclic dynamics and revealed substantial variation in temporal and spatial trends of cyclic dynamics across populations.

Why Do the Boreal Forest Ecosystems of Northwestern Europe Differ from Those of Western North America?

The boreal forest is one of the largest terrestrial biomes on Earth. Conifers normally dominate the tree layer across the biome, but other aspects of ecosystem structure and dynamics vary geographically. The cause of the conspicuous differences in the understory vegetation and the herbivore-predator cycles between northwestern Europe and western North America presents an enigma. Ericaceous dwarf shrubs and 3- to 4-year vole-mustelid cycles characterize the European boreal forests, whereas tall deciduous shrubs and 10-year snowshoe hare-lynx cycles characterize the North American ones. We discuss plausible explanations for this difference and conclude that it is bottom-up: Winter climate is the key determinant of the dominant understory vegetation that then determines the herbivore-predator food-web interactions. The crucial unknown for the twenty-first century is how climate change and increasing instability will affect these forests, both with respect to the dynamics of individual plant and animal species and to their community interactions.

Predator-vole interactions in Northern Europe: the role of small mustelids revised

The cyclic population dynamics of vole and predator communities is a key phenomenon in northern ecosystems, and it appears to be influenced by climate change. Reports of collapsing rodent cycles have attributed the changes to warmer winters, which weaken the interaction between voles and their specialist subnivean predators. Using population data collected throughout Finland during 1986-2011, we analyse the spatio-temporal variation in the interactions between populations of voles and specialist, generalist and avian predators, and investigate by simulations the roles of the different predators in the vole cycle. We test the hypothesis that vole population cyclicity is dependent on predator-prey interactions during winter. Our results support the importance of the small mustelids for the vole cycle. However, weakening specialist predation during winters, or an increase in generalist predation, was not associated with the loss of cyclicity. Strengthening of delayed density dependence coincided with strengthening small mustelid influence on the summer population growth rates of voles. In conclusion, a strong impact of small mustelids during summers appears highly influential to vole population dynamics, and deteriorating winter conditions are not a viable explanation for collapsing small mammal population cycles.

Dampening prey cycle overrides the impact of climate change on predator population dynamics: a long-term demographic study on tawny owls

Predicting the dynamics of animal populations with different life histories requires careful understanding of demographic responses to multifaceted aspects of global changes, such as climate and trophic interactions. Continent-scale dampening of vole population cycles, keystone herbivores in many ecosystems, has been recently documented across Europe. However, its impact on guilds of vole-eating predators remains unknown. To quantify this impact, we used a 27-year study of an avian predator (tawny owl) and its main prey (field vole) collected in Kielder Forest (UK) where vole dynamics shifted from a high- to a low-amplitude fluctuation regime in the mid-1990s. We measured the functional responses of four demographic rates to changes in prey dynamics and winter climate, characterized by wintertime North Atlantic Oscillation (wNAO). First-year and adult survival were positively affected by vole density in autumn but relatively insensitive to wNAO. The probability of breeding and number of fledglings were higher in years with high spring vole densities and negative wNAO (i.e. colder and drier winters). These functional responses were incorporated into a stochastic population model. The size of the predator population was projected under scenarios combining prey dynamics and winter climate to test whether climate buffers or alternatively magnifies the impact of changes in prey dynamics. We found the observed dampening vole cycles, characterized by low spring densities, drastically reduced the breeding probability of predators. Our results illustrate that (i) change in trophic interactions can override direct climate change effect; and (ii) the demographic resilience entailed by longevity and the occurrence of a floater stage may be insufficient to buffer hypothesized environmental changes. Ultimately, dampened prey cycles would drive our owl local population towards extinction, with winter climate regimes only altering persistence time. These results suggest that other vole-eating predators are likely to be threatened by dampening vole cycles throughout Europe.

Developing, testing and application of rodent population dynamics and capture models based on an adjusted Leslie matrix-based population approach

Small rodents in general and the multimammate rat Apodemus agrarius in particular, damage crops and cause major economic losses in China. Therefore, accurate predictions of the population size of A. agrarius and an efficient control strategy are urgently needed. We developed a population dynamics model by applying a Leslie matrix method, and a capture model based on optimal harvesting theory for A. agrarius. Our models were parametrized using demographic estimates from a capture–mark–recapture (CMR) study conducted on the Qinshui Forest Farm in Northwestern China. The population dynamics model incorporated 12 equally balanced age groups and included immigration and emigration parameters. The model was evaluated by assessing the predictions for four years based on the known starting population in 2004 from the 2004–2007 CMR data. The capture model incorporated two functional age categories (juvenile and adult) and used density-dependent and density-independent factors. The models were used to assess the effect of rodent control measures between 2004 and 2023 on population dynamics and the resulting numbers of rats. Three control measures affecting survival rates were considered. We found that the predicted population dynamics of A. agrarius between 2004 and 2007 compared favorably with the observed population dynamics. The models predicted that the population sizes of A. agrarius in the period between 2004 and 2023 under the control measure applied in August 2004 were very similar to the optimal population sizes, and no significant difference was found between the two population sizes. We recommend using the population dynamics and capture models based on CMR-estimated demographic schedules for rodent, provided these data are available. The models that we have developed have the potential to play an important role in predicting the effects of rodent management and in evaluating different control strategies.

The relationship between wolverine and larger predators, lynx and wolf, in a historical ecosystem context

Apex predators play an important role in shaping ecosystem structure. They may suppress smaller predators (mesopredators) but also subsidize scavengers via carrion provisioning. However, the importance of these interactions can change with ecosystem context. The wolverine (Gulo gulo) is a cold-adapted carnivore and facultative scavenger. It has a circumboreal distribution, where it could be either suppressed or subsidized by larger predators. In Scandinavia, the wolverine might interact with two larger predators, wolf (Canis lupus) and lynx (Lynx lynx), but human persecution decimated the populations in the nineteenth and early twentieth century. We investigated potential relationships between wolverine and the larger predators using hunting bag statistics from 15 Norwegian and Swedish counties in 1846-1922. Our best models showed a positive association between wolverine and lynx trends, taking ecological and human factors into account. There was also a positive association between year-to-year fluctuations in wolverine and wolf in the latter part of the study period. We suggest these associations could result from positive lynx-wolverine interactions through carrion provisioning, while wolves might both suppress wolverine and provide carrion with the net effect becoming positive when wolf density drops below a threshold. Wolverines could thus benefit from lynx presence and low-to-intermediate wolf densities.

Demographic responses of a site-faithful and territorial predator to its fluctuating prey: long-tailed skuas and arctic lemmings

Environmental variability, through interannual variation in food availability or climatic variables, is usually detrimental to population growth. It can even select for constancy in key life-history traits, though some exceptions are known. Changes in the level of environmental variability are therefore important to predict population growth or life-history evolution. Recently, several cyclic vole and lemming populations have shown large dynamical changes that might affect the demography or life-histories of rodent predators. Skuas constitute an important case study among rodent predators, because of their strongly saturating breeding productivity (they lay only two eggs) and high degree of site fidelity, in which they differ from nomadic predators raising large broods in good rodent years. This suggests that they cannot capitalize on lemming peaks to the same extent as nomadic predators and might be more vulnerable to collapses of rodent cycles. We develop a model for the population dynamics of long-tailed skuas feeding on lemmings to assess the demographic consequences of such variable and non-stationary prey dynamics, based on data collected in NE Greenland. The model shows that populations of long-tailed skua sustain well changes in lemming dynamics, including temporary collapses (e.g. 10 years). A high floater-to-breeder ratio emerges from rigid territorial behaviour and a long-life expectancy, which buffers the impact of adult abundance's decrease on the population reproductive output. The size of the floater compartment is affected by changes in both mean and coefficient of variation of lemming densities (but not cycle amplitude and periodicity per se). In Greenland, the average lemming density is below the threshold density required for successful breeding (including during normally cyclic periods). Due to Jensen's inequality, skuas therefore benefit from lemming variability; a positive effect of environmental variation. Long-tailed skua populations are strongly adapted to fluctuating lemming populations, an instance of demographic lability in the reproduction rate. They are also little affected by poor lemming periods, if there are enough floaters, or juveniles disperse to neighbouring populations. The status of Greenland skua populations therefore strongly depends upon floater numbers and juvenile movements, which are not known. This reveals a need to intensify colour-ringing efforts on the long-tailed skua at a circumpolar scale.

Synchrony in the snowshoe hare (Lepus americanus) cycle in northwestern North America, 1970–2012

Snowshoe hares (Lepus americanus Erxleben, 1777) fluctuate in 9–10 year cycles throughout much of their North American range. Regional synchrony has been assumed to be the rule for these cycles, so that hare populations in virtually all of northwestern North America have been assumed to be in phase. We gathered qualitative and quantitative data on hare numbers and fur returns of Canada lynx (Lynx canadensis Kerr, 1792) in the boreal forest regions of Alaska, Yukon, Northwest Territories, and northern British Columbia to describe synchrony in the time window of 1970–2012. Broad-scale synchrony in lynx fur returns was strong from 1970 to about 1995 but then seemed to break down in different parts of this region. Hare populations at 20 sites in Alaska, the Yukon, and Northwest Territories showed peak populations that lagged by 1–4 years during the 1990s and 2000s cycles. The simplest hypothesis to explain these patterns of asynchrony in hare cycles is the movement of predators from British Columbia north into the Yukon and then east into the Northwest Territories and west into Alaska. A traveling wave of these cycles is clearly seen in the lynx fur returns from western Canada and Alaska from 1970 to 2009. One consequence of a failure of synchrony is that hare predators like Canada lynx and Great-horned Owls (Bubo virginianus (Gmelin, 1788)) can move from one adjacent area to the next within this region and survive long enough to prolong low densities in hare populations that have declined earlier.

Nonlinear effects of climate on boreal rodent dynamics: mild winters do not negate high-amplitude cycles

Small rodents are key species in many ecosystems. In boreal and subarctic environments, their importance is heightened by pronounced multiannual population cycles. Alarmingly, the previously regular rodent cycles appear to be collapsing simultaneously in many areas. Climate change, particularly decreasing snow quality or quantity in winter, is hypothesized as a causal factor, but the evidence is contradictory. Reliable analysis of population dynamics and the influence of climate thereon necessitate spatially and temporally extensive data. We combined data on vole abundances and climate, collected at 33 locations throughout Finland from 1970 to 2011, to test the hypothesis that warming winters are causing a disappearance of multiannual vole cycles. We predicted that vole population dynamics exhibit geographic and temporal variation associated with variation in climate; reduced cyclicity should be observed when and where winter weather has become milder. We found that the temporal patterns in cyclicity varied between climatically different regions: a transient reduction in cycle amplitude in the coldest region, low-amplitude cycles or irregular dynamics in the climatically intermediate regions, and strengthening cyclicity in the warmest region. Our results did not support the hypothesis that mild winters are uniformly leading to irregular dynamics in boreal vole populations. Long and cold winters were neither a prerequisite for high-amplitude multiannual cycles, nor were mild winters with reduced snow cover associated with reduced winter growth rates. Population dynamics correlated more strongly with growing season than with winter conditions. Cyclicity was weakened by increasing growing season temperatures in the cold, but strengthened in the warm regions. High-amplitude multiannual vole cycles emerge in two climatic regimes: a winter-driven cycle in cold, and a summer-driven cycle in warm climates. Finally, we show that geographic climatic gradients alone may not reliably predict biological responses to climate change.

Evolution of predator dispersal in relation to spatio-temporal prey dynamics: how not to get stuck in the wrong place!

The eco-evolutionary dynamics of dispersal are recognised as key in determining the responses of populations to environmental changes. Here, by developing a novel modelling approach, we show that predators are likely to have evolved to emigrate more often and become more selective over their destination patch when their prey species exhibit spatio-temporally complex dynamics. We additionally demonstrate that the cost of dispersal can vary substantially across space and time. Perhaps as a consequence of current environmental change, many key prey species are currently exhibiting major shifts in their spatio-temporal dynamics. By exploring similar shifts in silico, we predict that predator populations will be most vulnerable when prey dynamics shift from stable to complex. The more sophisticated dispersal rules, and greater variance therein, that evolve under complex dynamics will enable persistence across a broader range of prey dynamics than the rules which evolve under relatively stable prey conditions.

The importance of marine vs. human-induced subsidies in the maintenance of an expanding mesocarnivore in the arctic tundra

1. Most studies addressing the causes of the recent increases and expansions of mesopredators in many ecosystems have focused on the top-down, releasing effect of extinctions of large apex predators. However, in the case of the northward expansion of the red fox into the arctic tundra, a bottom-up effect of increased resource availability has been proposed, an effect that can counteract prey shortage in the low phase of the multi-annual rodent cycle. Resource subsidies both with marine and with terrestrial origins could potentially be involved. 2. During different phases of a multi-annual rodent cycle, we investigated the seasonal dynamics and spatial pattern of resource use by red foxes across a coast to inland low arctic tundra gradient, Varanger Peninsula, Norway. We employed two complementary methods of diet analyses: stomach contents and stable isotope analysis. 3. We found that inland red foxes primarily subsisted on reindeer carrions during the low phase of a small rodent population cycle. Lemmings became the most important food item towards the peak phase of the rodent cycle, despite being less abundant than sympatric voles. Isotopic signatures of tissue from both predator and prey also revealed that red foxes near the coast used marine-derived subsidies in the winter, but these allochthonous resources did not spillover to adult foxes living beyond 20-25 km from the coast. 4. Although more needs to be learned about the link between increasing primary productivity due to climatic warming and trophic dynamics in tundra ecosystems, we suggest that changes in reindeer management through a bottom-up effect, at least regionally, may have paved the way towards the establishment of a new mesopredator in the tundra biome.

Determinants of lemming outbreaks

Population outbreaks in tundra rodents have intrigued scientists for a century as a result of their spectacular appearances and their general lessons in ecology. One outstanding question that has led to competing hypotheses is why sympatric lemmings and voles differ in regularity and shape of their outbreaks. Lemming outbreaks may be lost for decades while vole populations maintain regular population cycles. Moreover, when lemming populations eventually irrupt, they do so more steeply than the vole populations. Norwegian lemmings exhibited a large-scale outbreak synchronously with gray-sided voles in Finnmark, northern Fennoscandia, during 2006 to 2007 for the first time in two decades. Analyses of spatial variability of this outbreak across altitudinal gradients allowed us to identify determinants of the contrasting lemming and vole dynamics. The steeper lemming outbreak trajectories were caused by breeding and population growth during winter, when nonbreeding vole populations consistently declined. The differently shaped lemming and vole outbreaks appear to result from a particular demographic tactic of lemmings that evolved as an adaptation to the long and cold Arctic-Alpine winters. The lemming outbreak amplitude increased with altitude and vole density, indicating that lemming outbreaks are jointly facilitated by low temperatures and apparent mutualism with voles mediated by shared predators. High sensitivity to variation in climate and predation is likely to be the reasons why lemmings have more erratic population dynamics than sympatric voles. The combination of continued climatic warming and dampened vole cycles is expected to further decrease the frequency, amplitude, and geographic range of lemming outbreaks in tundra ecosystems.