ORIGINAL RESEARCH article

Front. Ecol. Evol., 02 January 2026

Sec. Population, Community, and Ecosystem Dynamics

Volume 13 - 2025 | https://doi.org/10.3389/fevo.2025.1595890

The triangulation between Arctic birds, Arctic foxes and Arctic lemmings

  • 1. Department of Zoology, Stockholm University, Stockholm, Sweden

  • 2. Department of Ecology, Lund University, Lund, Sweden

Abstract

Many studies have suggested that Arctic predators preying on lemmings, especially the Arctic fox, take advantage of the high prey densities during peak years, but have to switch to alternative prey during years of low lemming density, thereby forcing other prey species into the basic lemming cycle. Data on the number of breeding geese and waders in Arctic areas show three-year cycles, supporting this hypothesis. There are two alternative explanations for such a prey switch. The switch to the alternative prey could be due to a functional prey response, which assumes that the proportion of birds in fox diet should be related to lemming population phase due to a type III functional response between foxes and lemmings. The numerical switch explanation on the other hand, assumes that bird numbers should be related to lemming population phase because of a time lag in the numerical response between foxes and lemmings. With data on the diet and population numbers of Arctic foxes, Arctic lemmings and Arctic birds from eight different sites in northern Siberia, we tested these hypotheses. The use of birds in fox diet was not related to lemming population phase. Instead, we found a time lag in numerical response to lemmings, thus supporting the numerical switch explanation for the alternative prey hypothesis.

Introduction

Predation is a powerful process in terrestrial ecological communities (). Most attention has so far been paid to the larger predators (e.g. ; ) but there is an increasing interest in how meso-carnivores function in mammal communities (; ; ). Studies have shown a strong top-down effect of predation control on a preferred prey (e.g. ), but also a cascading effect on secondary prey populations ().

Most of these studies have been performed in complex mammal communities including many species. However, in ecological systems with fewer links, the community effects can be more pronounced and easier to study. Arctic mammal communities have very few but dominating species where Arctic lemmings (Lemmus spp, Dicrostonyx spp.) for example show dramatic fluctuations in their population dynamics (; ). During some years densities can be very high, followed by a sharp decline and crash in numbers. These density changes are known to follow a more or less cyclic pattern of three to four years between peak years and have attracted scientists for more than a century (; ; ; ; ). These cyclic dynamics of lemming numbers are of substantial interest in analyses of Arctic ecosystems (; ; ; ). Similar cyclic pattern can be seen for many birds of prey and mammalian carnivores (; ; ; ). Further, breeding success of Arctic birds such as waders and ducks have also been shown to follow the lemming cyclicity, often with one-year time lag (; ; ; ). Many studies have therefore suggested that Arctic predators preying on lemmings reach high densities during peak years, but switch to alternative prey during years of low lemming density, and thereby might force other prey species into the basic lemming cycle, the so called “Roselaar - Summers hypothesis” or “the alternative prey hypothesis” (; ; ; ; ). Recent warming in Arctic areas have severe consequences () and have raised concerns for disruption to the population dynamics both of lemmings (; ); ; ; ; but see ) and of Arctic birds (). This can have profound effects on several ecological processes ().

Data on the number of breeding geese and waders in Arctic areas typically indicate three-year cycles, supporting “the alternative prey hypothesis” (). Further, found a strong cross-correlation between Arctic fox (Vulpes lagopus) numbers and first-year brent geese (Branta bernicla), indicating that the breeding success of geese was better in years when lemming numbers increased than during lemming peak years. It was suggested that predators, especially the Arctic fox, after a decline in lemming numbers would switch diet from lemmings to birds (; , , ; ). Most studies have noted how the phase of the lemming cycle affects different species of Arctic birds. But interestingly, the apex of the widely accepted “alternative prey hypothesis” in Arctic areas, i.e. the mechanism behind a proposed switch in predator diet, is not clear.

The relationship between predators and their prey is governed by two major characteristics of the predator: their behavioral (functional) response to changes in prey numbers, and their numerical response to such changes. There are some scenarios in which Arctic foxes could affect the prey fluctuations observed by . One possibility is that Arctic foxes towards lemmings show a type III sigmoid response function typical of generalist predators (), where they readily switch to birds when lemming numbers decline. If so, the proportion of birds in fox diet should be related to the phase of the lemming cycle but independent of bird numbers. The switch to the alternative prey would in this case be due to a functional prey response. We call this the functional response hypothesis. If Arctic foxes instead show a type II functional response to lemmings () the switching to birds would be both weak and late.

An alternative explanation is that the cross-correlation pattern between foxes and birds is an effect of a delayed numerical response of fox numbers to lemming numbers, whereby fox numbers would be highest when the lemming population declines (which indeed was found by ). The switch to the alternative prey in this case would then be due to a numerical explanation. We call this the numerical response hypothesis. The observed variation in breeding success of for example geese would thus be an effect of a variation in predator numbers only, but not the result of a behavioral switch in diet preferences. In this case the number of birds should be related to the phase of the lemming cycle. The functional response to birds in this case could be of several types. A specialist (type II) functional response would emphasize fluctuations more than would a proportional (type I) response, as the bird population at intervals would be driven to very low levels by the predator. A type III response function to birds would instead tend to stabilize the predator-prey system ().

In a study of Arctic foxes and geese, pointed out that both numerical and functional response should be taken into account in analyses of prey switch and cyclic lemming populations. They demonstrated a similar numerical pattern as in other studies (e.g. ), but they had no quantitative data on functional response. In this study we attempt to test the two alternative mechanisms using a set of data from 17 study sites across Northern Siberia (). Thus, instead of temporal data that would be more natural to use in these analyses, we have used a geographical data set where natural variation within one season has given us contrasts in both predator and prey numbers and thereby generating comparable data as the temporal data set would. This assumes that there are similar ecological functions through Northern Siberia.

Materials and methods

Our study was performed during a ship-based expedition along the north coast of Siberia in the summer of 1994, where we visited 17 study sites from the Kola Peninsula in the west to Wrangel Island in the east (Figure 1; ). At each site, colleagues censused lemming populations, focusing on the Siberian lemming Lemmus sibiricus () and the collared lemming (Dicrostonyx torquatus) (; ), while we surveyed Arctic fox dens and collected scats for diet analysis (; ). All these inventories were about 10 kilometer apart (Hedberg 1985) and in similar habitats. From eight sites, 1 to 10, we also performed a quantitative survey of the abundance of birds during the first leg of the expedition (). There are thus two data sets, one for the first eight sites including bird inventories (Table 1) and the total set of 17 sites including estimates of abundance of lemmings and Arctic foxes (Table 2).

Figure 1

Table 1

SiteObserv.
Length km
Observ.
Time (h)
No. Birds/
10 h
No. Birds observedDominating species, no. obs.
DucksDivers, merganserGeese, swansGrouseDucksDivers, merganserGeese, swans
115.55.71032575127Clangula hyemalis 16Gavia stellate 7Anser fabilis 39
212.29.413634138220Anas crecca 11Mergus meganser 5Anser fabilis 44
314.38.640024216182Anas acuta 2Gavia arctica 4Anser albifrons 189
415.87.8265991013036Clangula hyemalis 58Gavia arctica 5Anser fabilis 63
542.717.2691014019Clangula hyemalis 8Gavia arctica 1Anser albifrons 23
830.813.874231456Polystica stelleri 15Gavia arctica 1Anser albifrons 30
927.19.8116901070Somateria spectabilis 90Branta bernicla 107
1032.812.922100102Somateria spectabilis 80Branta bernicla 10

Number of birds observed along line transects at the eight different sites in Siberia including the dominating species in each bird group on each site.

Table 2

SiteDateNo. trapped/24hArea
km2
Fox/
100km2
% in fox dietPopulation
LemmusDicroLem DicroLemmingsBirdsFresh Fox scatsPhase
Kola Peninsula19.60.890.001093.6726.120.0026.121.848Low
Kanin Peninsula211.60.000.00792.530.0035.002*Low
Kolguyev Island314.60.000.00705.710.098.219*n/a
Pechora Bay416.60.004.80847.1424.9713.8938.941.262Low
W Yamal Penins518.612.7116.479310.7544.8142.6587.54.16166Extended peak
N Yamal Penins621.60.008.73395.1380.946.7687.72.0050Peak/decline
NW Taymyr Peninsula824.636.264.607910.1351.0212.8463.929.335Extended peak
Chelyuskin Peninssula927.629.5448.92656.1539.2560.5599.80.0025Peak/decline
NE Taymyr Peninsula1030.623.522.401049.6264.4930.4695.00.5735Extended peak
Olenëkskiy Bay116.72.040.51902.2288.628.7297.30.6730Low
Yana Delta1212.714.889.126726.8779.8316.3796.23.2050Increase
N.S.I. Faadeyevskiy13a10.751.760.0010028.0098.800.0098.80.5375Peak
N.S.I. Kotel’nyy13b12.740.260.003528.5760.4528.7889.20.0040Peak/decline
Indigirka/Lopatka1415.714.290.5713024.6289.669.8499.50.2050Decline
Kolyma Delta1518.736.0011.001105.4552.9140.0993.07.0010Increase
Ayon Island1620.71.200.00508.0050.4735.0385.50.0030Low
Wrangel Island1725.733.0024.0016012.5035.3364.2399.60.4445Peak

Summary of lemming and arctic fox inventories.

Site numbers, dates (day.month) and names follow . Asterisk (*) denotes data that were excluded from calculations; at site 3 due to absence of rodents, at site 2 due to small sample size. N.S.I., New Siberian Islands. Data on arctic foxes are from and population phase of Lemmus from . For bird data see Table 1.

Arctic fox dens are usually situated in characteristic landforms and have lush vegetation due to digging and fertilization, making them relatively easy to locate at a distance (e.g. ; ; ). When scanning through the landscape we could identify dens from several kilometer distance and sometimes even at 10–20 kilometers. By walking at ridges we could thus scan through the landscape and cover large areas. A single visit at a den was sufficient to detect if it was occupied or not, and further inspections clarified whatever if it was a breeding den with a litter of Arctic foxes. When single foxes visit dens during the summer, there are no or very few food remains and no clear run ways (). Given the open and relatively homogenous tundra landscape, and well-established den monitoring routines, we assumed that we found a similar proportion of dens in all inventoried areas (). The number of occupied dens multiplied by two was used as an index of density of breeding Arctic foxes (). The area inventoried at each site varied from 35 to 160 km², mostly depending on the time spent at each site (Table 1). During the expedition, we covered 1464 km2 and inspected 142 Arctic fox dens ().

We collected Arctic fox scats (feces) at all occupied dens. Fresh scats, from the summer of 1994, were separated from older scats by appearance. We ignored scats that were two years or older, as determined by extensive weathering, generally being white and brittle, or overgrowth of vegetation from previous seasons. Fresh and old scats contained similar proportions of migrating birds (), indicating that scats on dens were from summers only, making age separation easier. Scats were dried at 90°C and prey remains were identified using reference material. In the analysis of scats we identified bird groups, rodent species, reindeer Rangifer tarandus, mountain hare Lepus timidus, insects and plant material, as far as possible. We could distinguish bird remains from Anseriformes and Galliformes (; ), but excluded waders since they are difficult to identify in fox scats (; ; ), and passerines because of their minor importance for Arctic foxes (; ). Egg shells are underestimated in scats and we only found less than 1% at the different sites. Molecular methods with DNA coding () or stable isotopes () would be needed to reveal the amount of eggs consumed. Thus. given with the problems of egg identification and bird inventories, our estimates regarding birds should be considered underestimations. However, these shortcomings would be similar in all areas and therefore not cause a systematic bias.

We used a modified frequency of occurrence measure to estimate the amount of each prey category (; ). When there were remains from more than one species in a single scat, we took into account the proportion of each prey species by dry volume. For example, one scat with 40% Lemmus and 60% Dicrostonyx plus another scat with 60% Lemmus and 40% Dicrostonyx, were considered to be equivalent to one scat with 100% Lemmus and one with 100% Dicrostonyx. We call this semi-quantitative measure “percent whole scat equivalents” (% WSE) (; ). Sample sizes remain the same as for frequency of occurrence. The advantage of this measure is that the relative amount of each prey category in the feces is taken into account. This is especially important for rare prey items that occur in small quantities in each scat. With a strict frequency of occurrence measure, these would be over-estimated.

We noted all birds observed along systematic line transects at the eight sites in June (Table 1; ). We chose routes to cover representative samples of the most important habitats. All transects were on low elevation tundra up to 15 km from the coast. During slow walks (the same speed on all transects), all birds seen or heard were recorded. The length of the transects differed between sites (Table 1), and we have thus used number of birds per 10 hours transect walk as an index of bird densities. We had to work during different hours a day, but light conditions do not change much between day and night in the Arctic summer. The main inventories were made between 09:00 and 21:00 (86%) and a few in the early mornings (06:00 – 09:00: 7%) or late evening (21:00 – 24:00: 7%). The weather conditions were stable with little rain (). In this paper we have used an estimate of the number of birds of interest for Arctic foxes, i.e. the sum of ducks, geese, swans and ptarmigan (Anseriformes and Galliformes). For a more detailed description of the bird surveys, see . Since we stayed at each site only 1–3 days, the breeding season for the birds was different during the expedition. However, we moved successively to higher latitudes through June, and therefore visited the different sites at approximately the same phenological time. The data emanate from the nest building, laying or early incubation period at sites 1-10. The statistical analyses are linear and non-linear regressions performed in excel.

Population densities of lemmings were estimated by snap trapping. However, different methods were used when trapping the two lemming species present. Relative density estimates of the Siberian lemming, Lemmus sibiricus, were obtained according to a standardized grid snap-trapping program with 300 snap-traps (). With data on lemming demography and dendrochronological methods, determined the phase of the population cycle for the Lemmus (Table 2). Traps for the collared lemming, Dicrostonyx torquatus, on the other hand, were set on selected sites with 200 Sherman live traps and 50 snap-traps (; ), and the results are not directly comparable to those from the grid trapping for Lemmus. We therefore calculated a “total lemming index” for each site (; ), based on total number of trapped lemmings (both Lemmus and Dicrostonyx) per 24 hours with a constant number of 550 traps. This is the lemming density index used throughout this study.

Results

The basic patterns in abundance of birds, lemmings and arctic foxes showed large variation over this vast area with clear population peaks of lemmings on especially the Arctic islands in eastern Siberia (Table 2) but also on Taymyr and Yamal Peninsulas (Table 2). The Arctic fox density followed this with similar high densities. For the birds the highest abundance was also on an island but in western Siberia, Kolguyev Island, followed by the close Pechora Bay (Table 1) with the dominating species White-fronted and Bean goose. However, on the Yamal Peninsula, Brant goose showed highest abundance (Table 1). Among the ducks, Long-tailed duck and King eiders were common (Table 1).

Arctic fox diet was dominated by lemmings of both Lemmus and Dicrostonyx genera. In populations with high lemming densities, there were between 60 and 99% lemmings in Arctic fox scats (Table 2). However, foxes also utilized birds to some extent and on Kolguyev Island (Table 2), where rodents are absent, birds constituted 98% of their diet. The highest bird density was also found on Kolguyev Island. Although we found rodent remains in one of the Arctic fox scats (n = 41) from Kolguyev Island, it is likely that this was from a migrating fox. We have therefore excluded data from that site in analyses concerning lemmings.

The numerical response of breeding Arctic foxes (per 100 km2) to density of birds (observed per 10 h) showed a negative trend and thus no numerical relationship (y = 8.2 - 0.009*x; r2 = 0.13; P > 0.05) (Figure 2). However, the numerical response of breeding Arctic foxes to lemmings (no. trapped per 24 h) showed signs of a numerical time lag (Figure 3, see further ). The bold arrows in Figure 3 show the population trend for the lemmings (Table 2) and the thin arrows indicate the direction of a limit population cycle (e.g. ). There would thus be a high density of foxes the year after a crash as predicted by the numerical switch explanation. Further, there was a tendency that bird density was higher at sites where the lemming population dynamics was in a low phase compared to where it was in peak phase (Figure 2; ANOVA: F (1,4) = 4.70, P = 0.09).

Figure 2

), from eight sites in northern Siberia.

Figure 3

and .

When we compared the amount of birds and lemmings in fox scats, the percentage of birds in fox diet was not related to the lemming population phase (ANOVA, F (3,19) = 0.76, P = 0.54), contradicting the functional switch explanation. On the other hand, there was a clear functional response for foxes preying upon birds, i.e. the percentage of birds in fox scats was related to bird density (Table 2, Figure 3). The response function for birds (Figure 4) was a linear type I function (y = 0.0032*x – 0.0493, r2 = 0.60, P<0.01). This indicates that Arctic foxes prey on birds with little behavioral specificity but rather in direct relation to their availability. On the other hand, Arctic foxes show a strong functional response to lemming species of both genera with a type II functional response (y = 94.9 * x/(x + 0.57); ).

Figure 4

Discussion

Our results indicates that the picture is somewhat more complicated than originally proposed by and , . The variation in reproductive success observed in Arctic geese (e.g. ; ) cannot be explained only by a diet switch in the Arctic fox, i.e. a switch due to the functional response hypothesis. We found no relationship between percentage of birds in fox diet and population phase for lemmings, and no evidence for a specific behavioral response other than a linear increase of birds in diet with increased bird density. Instead, there was a tendency that bird density was related to the phase of the lemming cycle, as many other studies have found (; ; ). However, for the numerical hypothesis we plot the number of breeding Arctic foxes against the lemming index and we get a pattern suggesting a time lag based on the population phase for the rodents (). This also indicates a pattern of a classic limit cycle for a predator-prey system that can be generated with a functional response (). It would be very interesting to analyze the empirical predator-prey data with the theoretical framework in e.g. . Although not formally tested, this indicates a time lag in numerical response for foxes preying on lemmings (, see also ), thus supports the numerical response hypothesis. This would result in a high density of Arctic foxes the year after a lemming population peak. The functional response of Arctic foxes preying on their major prey, the lemmings, was a steep type II, which would postpone the switch to birds and therefore strengthen the effect on bird numbers. showed that the functional response to Lemmus was of type II but for Dicrostonyx they found a type III response, although Lemmus is considered as the dominating species in Siberian tundra areas (; ). This could actually indicate a limit cycle (). Altogether, our results indicates that the population dynamics of lemmings influence breeding Arctic birds through a time lag in the numerical response that Arctic foxes show towards especially Lemmus. There might be another relationship between Dicrostonyx, Arctic foxes and birds, which could explain some anomalies in the breeding pattern of birds (). There might be similar difference in how Arctic foxes hunt different geese species or waders, but we were unable to detect any patterns of this in our data. More detailed studies, including longer visits during different phases of the lemming cycle, would be needed to validate this conclusion.

Arctic foxes seem to be specialized on Lemmus with a steep functional II response but show a type III response function to Dicrostonyx suggesting a stabilizing effect on the latter (). The predation pattern further indicates that Arctic foxes exploit birds in an opportunistic way in relation to availability (; , ). There are thus similarities in how Arctic foxes utilize birds and Dicrostonyx. However, in areas where Lemmus are absent or rare, Dicrostonyx might assume the role of the dominating, and thus preferred, prey species as for example on Greenland (). However, found a type III functional response for Arctic foxes and Dicrostonyx lemmings in Greenland. An interesting question is why a predator like the Arctic fox is not perfectly opportunistic, but shows preference for certain prey types. As shown on islands without rodents, or in contacts with man, Arctic foxes can quickly adapt to new food types and hunting techniques. Yet, in Siberia and Fennoscandia they are specialized on Lemmus (; ; ). It might be an effect of specific search images and hunting techniques, but this question remains unresolved.

Similar prey switching cases have been studied in a few other mammalian communities. , . found evidence that the mountain hare (Lepus timidus) population was forced into a red fox-vole cycle of four years through an alternative prey switching in the boreal ecosystem in Fennoscandia. Red foxes showed a time lag in numerical response to their main prey, Microtus agrestis and Myodes glareolus, and a type II functional response to both vole species (). There is no information on to what extent the secondary prey is consumed in relation to population phase of the primary prey species. However, all other information indicates that the cyclic population fluctuations for the alternative prey populations can be explained by a delayed numerical response. In boreal areas in North America, the dominating mammalian herbivore is the snowshoe hare (Lepus americanus) with its 10-year cycle. Predators like the Canadian lynx (Lynx candensis), great horned owl (Bubo virigianus) and coyote (Canis latrans) follow these cycles, but through prey switching different grouse populations and squirrels are also forced into the same cycle (; ; , , ). The predators showed a delayed numerical response to their main prey, the snowshoe hare, both in the Alberta and Yukon studies (; ; ). The functional response for lynx and great horned owl to snowshoe hares was of type II in both studies, and for coyotes a type II in the Yukon but type III in Alberta (; ).

In another red fox-vole system in southern Sweden, instead found that generalist predators stabilized vole fluctuations by switching between rabbits and voles due to a numerical response. came to the same conclusion when red foxes could switch between Arctic ground squirrels, Spermophillus paryii, and lemmings, and thereby stabilizing lemming fluctuations. Red foxes in Arctic areas do not prey on lemmings to the same extent as Arctic foxes do (; ; ). Thus, in these cases red foxes functioned as generalist predators and could stabilize prey populations through a high population density. It is possible that cases where cyclic predator population cause secondary prey populations to cycle are confined to animal communities with few species and low productivity (), where even “generalist” predators behave as specialists ().

This study gives support to the hypothesis that lemmings are key species in Arctic ecosystems and that changes in lemming abundance can cause cascade effects in several other organism groups. Lemming population fluctuations affect Arctic foxes, which in turn affect breeding success and population numbers of Brent geese. These observed cascade effects were caused by changes in predator numbers with a time lag, not by changes in predator behavior, thus supporting the numerical switch explanation for the alternative prey hypothesis.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.

Ethics statement

The animal study was approved by Försöksdjursorganisation vid Stockholms universitet. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

AA: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Writing – original draft, Writing – review & editing. MT: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Writing – review & editing. SS: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Writing – review & editing.

Funding

The author(s) declared that financial support was not received for this work and/or its publication. Funding was also received from the Carl Trygger foundation, Hierta-Retzius’ foundation, Ebba and Sven Schwartz foundation, the Ymer foundation, and foundations at the Royal Swedish Academy of Science.

Acknowledgments

We are most grateful to the Swedish Polar Research Secretariat for organizing the Tundra -94 expedition. We had very valuable help from D. Fucik with scat analysis and from N. Kjellén and E. Isakson with field work.

Conflict of interest

The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Summary

Keywords

predation, lemming cycle, geese, waders, Vulpes lagopus

Citation

Angerbjörn A, Tannerfeldt M and Svensson S (2026) The triangulation between Arctic birds, Arctic foxes and Arctic lemmings. Front. Ecol. Evol. 13:1595890. doi: 10.3389/fevo.2025.1595890

Received

18 March 2025

Revised

21 November 2025

Accepted

30 November 2025

Published

02 January 2026

Volume

13 - 2025

Edited by

Aurelio F. Malo, University of Alcalá, Spain

Reviewed by

Justin Bastow, Eastern Washington University, United States

Lech Stempniewicz, University of Gdańsk, Poland

Updates

Copyright

*Correspondence: Anders Angerbjörn,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics