Population fluctuations of the wood lemming<Emphasis Type="Italic">Myopus schisticolor</Emphasis> in eastern and western Finland

Olavi Eskelinen, Pertti Sulkava and Risto Sulkava From 1982 to 2003 we studied fluctuations in populations of the wood lemmingMyopus schisticolor (Liljeborg, 1844) in the Heinävesi (eastern Finland) and Keuruu regions (western Finland) by counting field signs and dead animals in standardized field surveys. We compared the population fluctuations of lemmings to those of other voles, owls and small mustelids in these regions. The lemming population in Heinävesi fluctuated regularly in 3-year cycles and in synchrony with the field vole population. Populations of owls also fluctuated in synchrony with wood lemming and field vole populations. In the Keuruu region, oscillations in the wood lemming population were irregular, and neither lemming and vole populations nor lemming and owl populations were correlated. Although direct mechanistic evidence is lacking, specialist predators such as weasels and owls probably cause the cyclicity in the field vole and wood lemming populations in the Heinävesi area. On the other hand, scarcity of high-quality habitats, unfavourable winter weather conditions and generalist predators may prevent the development of cyclicity in the wood lemming population of Keuruu.     

The breeding productivity of dark-bellied brent geese and curlew sandpipers in relation to changes in the numbers of arctic foxes and lemmings on the Taimyr Peninsula, Siberia —

There were about three-year cycles in the populations of arctic foxes, and the breeding productivities of brent geese and curlew sandpipers on the Taimyr Peninsula, Russia, The populations of arctic foxes and lemmings changed in synchrony. The breeding productivities of the birds tended to be good when the arctic foxes were increasing in numbers and poor when the arctic foxes were decreasing. There was a negative relationship between arctic fox numbers (or occupied lairs) and the breeding productivity of brent geese in the following year. Although there was evidence of wide-spread synchrony In the lemming cycle across the Taimyr Peninsula, some localities showed differences, However, such sites would still have been influenced by the general pattern of fox abundance in the typical tundra zone of the Taimyr Peninsula, where most of the arctic foxes breed and from which extensive movements of foxes occur after a decline in lemming numbers. The results support a prey-switching hypothesis (also known as the alternative prey hypothesis) whereby arctic foxes, and other predators, feed largely on lemmings when these are abundant or increasing, but switch to birds when the lemming population is small or declining. The relationships between arctic foxes, lemmings and brent geese may be further influenced by snowny owls which create fox-exclusion zones around their nests, thus providing safe nesting areas for the geese.     

Asynchronous population dynamics of Siberian lemmings across the Palaearctic tundra

The synchrony of Siberian lemming (Lemmus sibiricus L.) population dynamics was investigated during a ship-borne expedition along the Palaearctic tundra coast in the summer of 1994. On 12 sites along the coast from the Kola Peninsula to Wrangel Island, relative densities of lemmings were recorded using a standardised snap-trapping programme. The phase position of the lemming cycle in each of the studied populations was determined based on current density estimates, signs of previous density and the age profile of each population (ageing based on eye lens mass). In addition, dendrochronological methods were used to determine when the last peak in the density of microtine populations occurred at each site. The examined lemming populations were in different phases of the lemming cycle. Some populations were in the peak phase, as indicated by high current densities, an age profile in which older individuals were well represented, and signs of high previous density (abundant old lemming faeces). Other populations were in the decline phase, as reflected in a moderate current density, a predominance of older individuals and signs of high previous density. Populations in the low phase had an extremely low current density and showed signs of high previous density, while populations in the increase phase had a moderate current density, a predominance of younger individuals and showed signs of low previous density. The results of phase determinations based on dendrochronological methods support the findings based on lemming demography. Recent Russian studies carried out on some of the sites also agreed with our phase determination results. Thus, on a regional scale (across the whole Palaearctic tundra), the population dynamics of Siberian lemmings can be considered asynchronous. However, sites situated adjacent to each other were often phase synchronous, suggesting a more fine-grained pattern of dynamics with synchrony over distances as long as 1000 km or so, e.g. the Yamal and Taymyr Peninsulas. Received: 19 August 1998 / Accepted: 1 March 1999     

Response of an arctic predator guild to collapsing lemming cycles

Alpine and arctic lemming populations appear to be highly sensitive to climate change, and when faced with warmer and shorter winters, their well-known high-amplitude population cycles may collapse. Being keystone species in tundra ecosystems, changed lemming dynamics may convey significant knock-on effects on trophically linked species. Here, we analyse long-term (1988–2010), community-wide monitoring data from two sites in high-arctic Greenland and document how a collapse in collared lemming cyclicity affects the population dynamics of the predator guild. Dramatic changes were observed in two highly specialized lemming predators: snowy owl and stoat. Following the lemming cycle collapse, snowy owl fledgling production declined by 98 per cent, and there was indication of a severe population decline of stoats at one site. The less specialized long-tailed skua and the generalist arctic fox were more loosely coupled to the lemming dynamics. Still, the lemming collapse had noticeable effects on their reproductive performance. Predator responses differed somewhat between sites in all species and could arise from site-specific differences in lemming dynamics, intra-guild interactions or subsidies from other resources. Nevertheless, population extinctions and community restructuring of this arctic endemic predator guild are likely if the lemming dynamics are maintained at the current non-cyclic, low-density state.     

Spatial dynamics of Norwegian tetraonid populations

Different species in a given site or population of a given species in different sites may fluctuate in synchrony if they are affected similarly by factors such as spatially autocorrelated climate, predation, or by dispersal between populations of one species. We used county wise time series of hunting bag records of four Norwegian tetraonid species covering 24 years to examine patterns of interspecific and intraspecific synchrony. We estimated synchrony at three spatial scales; national, regional (consisting of counties with similar climate), and county level. Ecologically related species with overlapping distributions exhibited strong synchrony across Norway, but there was much variation between the different regions and counties. Regions with a long coastline to both the North Sea and the Norwegian Ocean exhibited an overall stronger synchrony than those consisting of more continental areas. Intraspecific synchrony was generally low across all counties, but stronger synchrony between counties within regions defined by climatic conditions. Synchrony was negatively related to distance between populations in three of four species. Only the synchrony in willow ptarmigan showed a clear negative relationship with distance, while the other species had both strong positive and negative correlations at short distances. Strong interspecific synchrony between some species pairs within regions and weak intraspecific synchrony across counties within regions suggest a stronger synchronizing effect from environmental factors such as weather or predation and less effect from dispersal. Our results suggest that the complete tetraonid community is structured by environmental factors affecting the different species similarly and causes widespread interspecific synchrony. Local factors affecting the population dynamics nevertheless frequently forces neighbouring populations out of phase.     

Lack of concordance between mtDNA gene flow and population density fluctuations in the bank vole

The genetic structure of bank voles Clethrionomys glareolus was determined from analyses of mitochondrial DNA (mtDNA) sequences, and compared with previous data on geographical synchrony in population density fluctuations. From 31 sample sites evenly spaced out along a 256-km transect in SE Norway a total of 39 distinct mtDNA haplotypes were found. The geographical distribution of the haplotypes was significantly nonrandom, and a cladistic analysis of the evolutionary relationship among haplotypes shows that descendant types were typically limited to a single site, whereas the ancestral types were more widely distributed geographically. This geographical distribution pattern of mtDNA haplotypes strongly indicates that the range and amount of female dispersal is severely restricted and insufficient to account for the previously observed synchrony in population density fluctuations. We conclude that geographical synchrony in this species must be caused by factors that are external to the local population, such as e.g. mobile predators.     

The role of lemming herbivory in the sex ratio and shoot demography of willow populations

We investigated the role of lemming herbivory on the age structure and physical form of a stand of willows ( Salix lanata ) using lemming scars on stems as an indication of past herbivory. Salix lanata had a female-biased sex ratio comparable to that found in other willow species and it has been proposed that such a sex ratio may be due to selective herbivory on male genets. There was, however, no difference in the degree of scarring in male and female plants. Scarring did not appear to change the overall physical structure of the willows including the number of stems, the proportion of dead stems and the degree of branching. However, 72% of scarred 7 mm stems were dead, compared to only 22% of non-scarred stems, suggesting that lemming herbivory may play a role in ramet death. There was a significant positive correlation between stem ages and the lemming cycle as determined by the frequency of scars. With experimental scarring, S. lanata increased its production of new shoots in relation to the degree of scarring. We suggest that this demonstrates a pattern of compensatory growth in willows, whereby they take advantage of a period of relatively low herbivory following the crash of the lemming populations in order to compensate for damage that occurred during the lemming peak.     

Density-dependent dispersal and spatial population dynamics

The synchronization of the dynamics of spatially subdivided populations is of both fundamental and applied interest in population biology. Based on theoretical studies, dispersal movements have been inferred to be one of the most general causes of population synchrony, yet no empirical study has mapped distance-dependent estimates of movement rates on the actual pattern of synchrony in species that are known to exhibit population synchrony. Northern vole and lemming species are particularly well-known for their spatially synchronized population dynamics. Here, we use results from an experimental study to demonstrate that tundra vole dispersal movements did not act to synchronize population dynamics in fragmented habitats. In contrast to the constant dispersal rate assumed in earlier theoretical studies, the tundra vole, and many other species, exhibit negative density-dependent dispersal. Simulations of a simple mathematical model, parametrized on the basis of our experimental data, verify the empirical results, namely that the observed negative density-dependent dispersal did not have a significant synchronizing effect.     

Does lemming winter grazing impact vegetation in the Canadian Arctic?

In low-productivity environments such as the tundra, it has been proposed that regular, multi-annual population cycles of lemmings could be driven by winter food depletion in years of peak abundance. If lemming population dynamics is controlled by food resources, we predict that (1) winter grazing should negatively impact the abundance of food plants, (2) this impact should be proportional to lemming density and (3) high lemming winter grazing pressure should result in reduced plant growth during the following summer. We tested these predictions on Bylot Island, Nunavut, Canada, where two species of lemmings are present: the brown (Lemmus trimucronatus) and collared lemming (Dicrostonyx groenlandicus). We installed 16 exclosures in their preferred wintering habitat (snowbeds) and annually sampled plant biomass inside and outside exclosures at snow melt and at peak growth during the summers of 2009–2012, covering a full population cycle. Winter grazing had no impact on total vascular plant or moss biomass at snow melt in all years. Among plant families, only Caryophyllaceae, which was uncommon, showed a decline. In moss taxa, a negative effect was found on Polytrichum in only 1 year out of three. Overall, plant regrowth during the subsequent summer showed annual variation and tended to be reduced in the 2 years of high lemming abundance. However, this could be a consequence of summer grazing. Overall, the impact of lemming winter grazing on plants was weak and short-lived, even in years of high lemming abundance. Therefore, our results are not consistent with the hypothesis that food depletion during winter was the cause of the lemming decline following peak abundance at our study site. Other factors may limit lemming populations and prevent them from reaching densities high enough to exhaust their food resources.     

Indirect effects of lemming cycles on sandpiper dynamics: 50 years of counts from southern Sweden

The bird-lemming hypothesis postulates that breeding success of tundra-nesting geese and waders in Siberia follows the cyclic pattern of lemming populations, as a result of predators switching from lemmings to birds when the lemming population crashes. We present 50 years of data on constant-effort catches of red knot Calidris canutus and curlew sandpiper C. ferruginea at an autumn migratory stopover site (Ottenby) at the Baltic Sea, supplemented with literature data on winter censuses of dark-bellied brent goose Branta b. bernicla and white-fronted goose Anser albifrons in northwestern Europe, and waders in Germany and Southern Africa. Number and proportion of juveniles in these bird populations (both our own and literature data) were compared with an index of predation pressure (calculated from the abundance of lemmings on the Taimyr peninsula), and climate indices for the North Eurasia and the North Atlantic regions. The index of predation pressure correlated significantly with the number of juveniles of red knot and curlew sandpiper, but not with number of adults. Also, this index correlated with the reproductive performance of geese and waders reported in the literature. Fourier analysis revealed a significant deviation from random noise with the maximum spectral density at the period length of 3 years for number of juvenile red knots and curlew sandpipers captured at Ottenby, abundance of lemmings, reproduction in arctic fox Alopex lagopus, and reproductive performance in geese on the Siberian tundra. Also, the date of passage at Ottenby for adult red knot and curlew sandpiper showed a spectral density peak at a period length of 3 years, the latter species also showing a peak at a period length of 5-6 years. Passage dates for adult red knot and curlew sandpiper were earlier in years of high predation pressure compared with years of low predation pressure. The fluctuations in reproductive success of the studied Siberian goose and wader species appear to be primarily influenced by biotic factors in the breeding area, rather than by abiotic factors, such as climate oscillations. Annual variations in migratory arctic bird populations may have far reaching effects in habitats along their migration routes and in their wintering areas. We suggest a link between lemming cyclicity in the Northern Hemisphere and predation pressure on Southern Hemisphere benthos, in which the signal is carried between continents by long distance migrating waders.     

Causes of cyclicity of Epirrita autumnata (Lepidoptera, Geometridae): grandiose theory and tedious practice

Creating multiyear cycles in population density demands, in traditional models, causal factors that operate on local populations in a density-dependent way with time lags. However, cycles of the geometrid Epirrita autumnata in northern Europe may be regional, not local; i.e., successive outbreaks occur in different localities. We review possible causes of cycles of E. autumnata under both local and regional scenarios, including large-scale synchrony. Assuming cyclicity is a local phenomenon, individual populations of E. autumnata display peaks but populations all over the outbreak range fluctuate in synchrony. This concept assumes that the peaks at most localities are so low that they do not lead to visible defoliation and easily remain unnoticed. In this scenario, populations are able to start recovery a few years after the crash, i.e., at the time of the mitigation of detrimental delayed density-dependent factors, such as delayed inducible resistance of the host plant or parasitism. In that case, the same factors that lead to crashes also explain the periodicity of cyclic fluctuations. According to the regional cyclicity scenario, different factors can be important in different phases of the cycle. The key is to identify the factors that tend to produce outbreaks with a periodicity of about 10 years. Initiation of the increase phase seems to coincide with maxima in sunspot activity, but causal connections remain unclear. Climatic factor(s) associated with the solar cycle could contribute to the large-scale geographic synchrony.     

Are goose nesting success and lemming cycles linked? Interplay between nest density and predators

The suggested link between lemming cycles and reproductive success of arctic birds is caused by potential effects of varying predation pressure (the Alternative Prey Hypothesis, APH) and protective association with birds of prey (the Nesting Association Hypothesis, NAH). We used data collected over two complete lemming cycles to investigate how fluctuations in lemming density were associated with nesting success of greater snow geese ( Anser caerulescens atlanticus ) in the Canadian High Arctic. We tested predictions of the APH and NAH for geese breeding at low and high densities. Goose nesting success varied from 22% to 91% between years and the main egg predator was the arctic fox ( Alopex lagopus ). Nesting associations with snowy owls ( Nyctea scandiaca ) were observed but only during peak lemming years for geese nesting at low density. Goose nesting success declined as distance from owls increased and reached a plateau at 550 m. Artificial nest experiments indicated that owls can exclude predators from the vicinity of their nests and thus reduce goose egg predation rate. Annual nest failure rate was negatively associated with rodent abundance and was generally highest in low lemming years. This relationship was present even after excluding goose nests under the protective influence of owls. However, nest failure was inversely density-dependent at high breeding density. Thus, annual variations in nest density influenced the synchrony between lemming cycles and oscillations in nesting success. Our results suggest that APH is the main mechanism linking lemming cycles and goose nesting success and that nesting associations during peak lemming years (NAH) can enhance this positive link at the local level. The study also shows that breeding strategies used by birds (the alternative prey) could affect the synchrony between oscillations in avian reproductive success and rodent cycles.     

Predator–prey relationships: arctic foxes and lemmings

1. The number of breeding dens and litter sizes of arctic foxes Alopex lagopus were recorded and the diet of the foxes was analysed during a ship-based expedition to 17 sites along the Siberian north coast. At the same time the cyclic dynamics of co-existing lemming species were examined.
2. The diet of arctic foxes was dominated by the Siberian lemming Lemmus sibiricus (on one site the Norwegian lemming L. lemmus ), followed by the collared lemming Dicrostonyx torquatus .
3. The examined Lemmus sibiricus populations were in different phases of the lemming cycle as determined by age profiles and population densities.
4. The numerical response of arctic foxes to varying densities of Lemmus had a time lag of 1 year, producing a pattern of limit cycles in lemming–arctic fox interactions. Arctic fox litter sizes showed no time lag, but a linear relation to Lemmus densities. We found no evidence for a numerical response to population density changes in Dicrostonyx .
5. The functional or dietary response of arctic foxes followed a type II curve for Lemmus , but a type III response curve for Dicrostonyx .
6. Arctic foxes act as resident specialist for Lemmus and may increase the amplitude and period of their population cycles. For Dicrostonyx , on the other hand, arctic foxes act as generalists which suggests a capacity to dampen oscillations.     

Spatially autocorrelated disturbances and patterns in population synchrony

Spatially synchronous population dynamics have been documented in many taxa. The prevailing view is that the most plausible candidates to explain this pattern are extrinsic disturbances (the Moran effect) and dispersal. In most cases disentangling these factors is difficult. Theoretical studies have shown that dispersal between subpopulations is more likely to produce a negative relationship between population synchrony and distance between the patches than perturbations. As analyses of empirical data frequently show this negative relationship between the level of synchrony and distance between populations, this has emphasized the importance of dispersal as a synchronizing agent. However, several weather patterns show spatial autocorrelation, which could potentially produce patterns in population synchrony similar to those caused by dispersal. By using spatially extended versions of several population dynamic models, we show that this is indeed the case. Our results show that, especially when both factors (spatially autocorrelated perturbations and distance-dependent dispersal) act together, there may exist groups of local populations in synchrony together but fluctuating asynchronously with some other groups of local populations. We also show, by analysing 56 long-term population data sets, that patterns of population synchrony similar to those found in our simulations are found in natural populations as well. This finding highlights the subtlety in the interactions of dispersal and noise in organizing spatial patterns in population fluctuations.     

The effect of snow cover on lemming population cycles in the Canadian High Arctic

Rising temperatures and changes in the precipitation regime will have a strong impact on the quality of the snow cover in the Arctic. A snow cover of good quality protecting lemmings from cold temperatures and predators is thought to be an important factor for maintaining the cyclic dynamic of their populations in the tundra. We examined if the characteristics of annual fluctuations (amplitude and shape of phases) in brown lemming (Lemmus trimucronatus) density could be determined by snow depth, snow density, sub-nivean temperature and persistence of snow. Using an 18-year time series of brown lemming abundance on Bylot Island in the Canadian Arctic, we tested if snow variables could explain the residual variation between the observed lemming density and the one predicted by models where cyclicity had been accounted for. Our analysis provides support for the hypothesis that snow cover can affect the amplitude and possibly also the periodicity of lemming population cycles in the High Arctic. Summer abundance of brown lemmings was higher following winters with a deep snow cover and a low-density snow pack near the ground but was unaffected by the date of establishment or melting and duration of the snow cover. Two snow variables showed a temporal trend; mean winter snow depth tended to increase and date of establishment of the hiemal threshold occurred earlier over time. These temporal trends, which should be favourable to lemmings, may explain why healthy population cycles have apparently been maintained at our study site contrary to other Arctic sites.     

Spatial (a)synchrony in population fluctuations of five plant species in fragmented habitats

The spatial scale at which populations show synchronous temporal fluctuations in abundance, relative to the spatial scale over which they can disperse, may influence the persistence of local and regional populations. There have been frequent demonstrations of spatial synchrony in population dynamics of animal populations. But few studies have investigated the degree of spatial synchrony in less mobile taxa, e.g. plants, where life history, dispersal and interaction with the environment would be different due to a sessile phase. This study has during three years investigated the synchrony in local population size changes in four short-lived species, and during a nine-year period for one long-lived species, in a semi-natural grassland landscape in southern Sweden. The spatial scale of this study was less than 15 km, which is quite small in comparison with other studies, but the temporal scale was of similar magnitude as the few studies on plant abundances and synchrony. When using detrended estimates of population size change, a significant pattern of decreasing synchrony with increasing distance was found for the two short-lived species that were most confined to manage semi-natural grasslands. Spatial synchrony was detected up to a few km. However, the species displayed synchrony in different years. The degree of synchrony can thus vary considerably across years and among species. Spatially autocorrelated weather conditions could partly explain the spatial scale of synchrony found during certain time intervals. However, the prevailing asynchrony suggests that local factors dominate the dynamics of the populations at the investigated scale.     

Dispersal, Environmental Correlation, and Spatial Synchrony in Population Dynamics

Many species exhibit widespread spatial synchrony in population fluctuations. This pattern is of great ecological interest and can be a source of concern when the species is rare or endangered. Both dispersal and spatial correlations in the environment have been implicated as possible causes of this pattern, but these two factors have rarely been studied in combination. We develop a spatially structured population model, simple enough to obtain analytic solutions for the population correlation, that incorporates both dispersal and environmental correlation. We ask whether these two synchronizing factors contribute additively to the total spatial population covariance. We find that there is always an interaction between these two factors and that this interaction is small only when one or both of the environmental correlation and the dispersal rate are small. The interaction is opposite in sign to the environmental correlation; so, in the normal case of positive environmental correlation across sites, the population synchrony will be lower than predicted by simply adding the effects of dispersal and environmental correlation. We also find that population synchrony declines as the strength of population regulation increases. These results indicate that dispersal and environmental correlation need to be considered in combination as explanations for observed patterns of population synchrony.     

Cycles and synchrony in the Collared Lemming (Dicrostonyx groenlandicus) in Arctic North America

Lemming populations are generally characterised by their cyclic nature, yet empirical data to support this are lacking for most species, largely because of the time and expense necessary to collect long-term population data. In this study we use the relative frequency of yearly willow scarring by lemmings as an index of lemming abundance, allowing us to plot population changes over a 34-year period. Scars were collected from 18 sites in Arctic North America separated by 2-1,647 km to investigate local synchrony among separate populations. Over the period studied, populations at all 18 sites showed large fluctuations but there was no regular periodicity to the patterns of population change. Over all possible combinations of pairs of sites, only sites that were geographically connected and close (<6 km) showed significant synchrony in fluctuations. The populations studied may not even be cyclic, at least for the time period 1960 to 1994, and although fluctuating, randomisation tests could not reject the null hypothesis of random fluctuations. These data have implications for the testing of hypotheses regarding lemming cycles and highlight the need for long-term trapping data to characterise the lemming cycle.     

Environment,but not migration rate,influences extinction risk in experimental metapopulations

Ecological theory suggests that several demographic factors influence metapopulation extinction risk, including synchrony in population size between subpopulations, metapopulation size and the magnitude of fluctuations in population size. Theoretically, each of these is influenced by the rate of migration between subpopulations. Here we report on an experiment where we manipulated migration rate within metapopulations of the freshwater zooplankton Daphnia magna to examine how migration influenced each of these demographic variables, and subsequent effects on metapopulation extinction. In addition, our experimental procedures introduced unplanned but controlled differences between metapopulations in light intensity, enabling us to examine the relative influences of environmental and demographic factors. We found that increasing migration rate increased subpopulation synchrony. We failed to detect effects of migration on population size and fluctuations in population size at the metapopulation or subpopulation level, however. In contrast, light intensity did not influence synchrony, but was positively correlated with population size and negatively correlated with population fluctuation. Finally, synchrony did not influence time to extinction, while population size and the magnitude of fluctuations did. We conclude that environmental factors had a greater influence on extinction risk than demographic factors, and that metapopulation size and fluctuation were more important to extinction risk than metapopulation synchrony.     

Geographical variation in the spatial synchrony of a forest-defoliating insect: isolation of environmental and spatial drivers

Despite the pervasiveness of spatial synchrony of population fluctuations in virtually every taxon, it remains difficult to disentangle its underlying mechanisms, such as environmental perturbations and dispersal. We used multiple regression of distance matrices (MRMs) to statistically partition the importance of several factors potentially synchronizing the dynamics of the gypsy moth, an invasive species in North America, exhibiting outbreaks that are partially synchronized over long distances (approx. 900 km). The factors considered in the MRM were synchrony in weather conditions, spatial proximity and forest-type similarity. We found that the most likely driver of outbreak synchrony is synchronous precipitation. Proximity played no apparent role in influencing outbreak synchrony after accounting for precipitation, suggesting dispersal does not drive outbreak synchrony. Because a previous modelling study indicated weather might indirectly synchronize outbreaks through synchronization of oak masting and generalist predators that feed upon acorns, we also examined the influence of weather and proximity on synchrony of acorn production. As we found for outbreak synchrony, synchrony in oak masting increased with synchrony in precipitation, though it also increased with proximity. We conclude that precipitation could synchronize gypsy moth populations directly, as in a Moran effect, or indirectly, through effects on oak masting, generalist predators or diseases.