Is the diet cyclic phase-dependent in boreal vole populations?

Herbivorous rodents in boreal, alpine and arctic ecosystems are renowned for their multi-annual population cycles. Researchers have hypothesised that these cycles may result from herbivore-plant interactions in various ways. For instance, if the biomass of preferred food plants is reduced after a peak phase of a cycle, rodent diets can be expected to become dominated by less preferred food plants, leading the population to a crash. It could also be expected that the taxonomic diversity of rodent diets increases from the peak to the crash phase of a cycle. The present study is the first to use DNA metabarcoding to quantify the diets of two functionally important boreal rodent species (bank vole and tundra vole) to assess whether their diet changed systematically in the expected cyclic phase-dependent manner. We found the taxonomic diet spectrum broad in both vole species but with little interspecific overlap. There was no evidence of systematic shifts in diet diversity metrics between the phases of the population cycle in either species. While both species' diet composition changed moderately between cycle phases and seasons, these changes were small compared to other sources of diet variation-especially differences between individuals. Thus, the variation in diet that could be attributed to cyclic phases is marginal relative to the overall diet flexibility. Based on general consumer-resource theory, we suggest that the broad diets with little interspecific overlap render it unlikely that herbivore-plant interactions generate their synchronous population cycles. We propose that determining dietary niche width should be the first step in scientific inquiries about the role of herbivore-plant interactions in cyclic vole populations.

Long-term monitoring of cycles in Clethrionomys rutilus in the Yukon boreal forest

Baseline studies of small rodent populations in undisturbed ecosystems are rare. We report here 50 years of monitoring and experimentation in Yukon of a dominant rodent species in the North American boreal forest, the red-backed vole Clethrionomys rutilus. These voles breed in summer, weigh 20-25 g, and reach a maximum density of 20 to 25 per ha. Their populations have shown consistent 3-4-year cycles for the last 50 years with the only change being that peak densities averaged 8/ha until 2000 and 18/ha since that year. During the last 25 years, we have measured food resources, predator numbers, and winter weather, and for 1-year social interactions, to estimate their contribution to changes in the rate of summer increase and the rate of overwinter decline. All these potential limiting factors could contribute to changes in density, and we measured their relative contributions statistically with multiple regressions. The rate of winter decline in density was related to both food supply and winter severity. The rate of summer increase was related to summer berry crops and white spruce cone production. No measure of predator numbers was related to winter or summer changes in vole abundance. There was a large signal of climate change effects in these populations. There is no density dependence in summer population growth and only a weak one in winter population declines. None of our results provide a clear understanding of what generates 3-4-year cycles in these voles, and the major missing piece may be an understanding of social interactions at high density.

Are population cycles recovering?

The dynamics of populations of small mammals of Central Siberia was analyzed. The studies were carried out at the Yenisei ecological station "Mirnoye" of the A.N. Severtsov Institute of Ecology and Evolution, Russian Academy of Sciences. The time series analysis was performed by the wavelet transform using the statistical data processing language R. In the 20th century, the dynamics of the population of the community and some of its constituent species (Sorex araneus; S. caecutiens; S. isodon; S. tundrensis; S. minutus; Craseomys rufocanus; Clethrionomys rutilus; Microtus oeconomus; M. agrestis) were characterized by a 4-year periodicity. The type of dynamics changed to noncyclic by the nineties, but by 2022, four species (S. araneus, S. isodon, C. rutilus, and M. oeconomus) and the community as a whole showed a tendency toward recovery of population cycles. The remaining species were characterized by consistently low numbers with irregular low amplitude fluctuations.

Phenological mismatch affects individual fitness and population growth in the winter moth

Climate change can severely impact species that depend on temporary resources by inducing phenological mismatches between consumer and resource seasonal timing. In the winter moth, warmer winters caused eggs to hatch before their food source, young oak leaves, became available. This phenological mismatch changed the selection on the temperature sensitivity of egg development rate. However, we know little about the fine-scale fitness consequences of phenological mismatch at the individual level and how this mismatch affects population dynamics in the winter moth. To determine the fitness consequences of mistimed egg hatching relative to timing of oak budburst, we quantified survival and pupation weight in a feeding experiment. We found that mismatch greatly increased mortality rates of freshly hatched caterpillars, as well as affecting caterpillar growth and development time. We then investigated whether these individual fitness consequences have population-level impacts by estimating the effect of phenological mismatch on population dynamics, using our long-term data (1994-2021) on relative winter moth population densities at four locations in The Netherlands. We found a significant effect of mismatch on population density with higher population growth rates in years with a smaller phenological mismatch. Our results indicate that climate change-induced phenological mismatch can incur severe individual fitness consequences that can impact population density in the wild.

Cold winters drive consistent and spatially synchronous 8-year population cycles of cabbage stem flea beetle

Population cycles have been observed in mammals as well as insects, but consistent population cycling has rarely been documented in agroecosystems and never for a beetle. We analysed the long-term population patterns of the cabbage stem flea beetle Psylliodes chrysocephala in winter oilseed rape over 50 years. Psylliodes chrysocephala larval density from 3045 winter oilseed rape fields in southern Sweden showed strong 8-year population cycles in regional mean density. Fluctuations in larval density were synchronous over time across five subregional populations. Subregional mean environmental variables explained 90.6% of the synchrony in P. chrysocephala populations at the 7-11 year time-scale. The number of days below -10°C showed strong anti-phase coherence with larval densities in the 7-11 year time-scale, such that more cold days resulted in low larval densities. High levels of the North Atlantic Oscillation weather system are coherent and anti-phase with cold weather in Scania, Sweden. At the field-scale, later crop planting date and more cold winter days were associated with decreased overwintering larval density. Warmer autumn temperatures, resulting in greater larval accumulated degree days early in the season, increased overwintering larval density. Despite variation in environmental conditions and crop management, 8-year cycles persisted for cabbage stem flea beetle throughout the 50 years of data collection. Moran effects, influenced by the North Atlantic Oscillation weather patterns, are the primary drivers of this cycle and synchronicity. Insect pest data collected in commercial agriculture fields is an abundant source of long-term data. We show that an agricultural pest can have the same periodic population cycles observed in perennial and unmanaged ecosystems. This unexpected finding has implications for sustainable pest management in agriculture and shows the value of long-term pest monitoring projects as an additional source of time-series data to untangle the drivers of population cycles.

Demographic Basis of Spatially Structured Fluctuations in a Threespine Stickleback Metapopulation

AbstractUncovering the demographic basis of population fluctuations is a central goal of population biology. This is particularly challenging for spatially structured populations, which require disentangling synchrony in demographic rates from coupling via movement between locations. In this study, we fit a stage-structured metapopulation model to a 29-year time series of threespine stickleback abundance in the heterogeneous and productive Lake Mývatn, Iceland. The lake comprises two basins (North and South) connected by a channel through which the stickleback disperse. The model includes time-varying demographic rates, allowing us to assess the potential contributions of recruitment and survival, spatial coupling via movement, and demographic transience to the population's large fluctuations in abundance. Our analyses indicate that recruitment was only modestly synchronized between the two basins, whereas survival probabilities of adults were more strongly synchronized, contributing to cyclic fluctuations in the lake-wide population size with a period of approximately 6 years. The analyses further show that the two basins were coupled through movement, with the North Basin subsidizing the South Basin and playing a dominant role in driving the lake-wide dynamics. Our results show that cyclic fluctuations of a metapopulation can be explained in terms of the combined effects of synchronized demographic rates and spatial coupling.

An experimental analysis of density dependence in meadow voles: Within-season and delayed effects

Wild mammal populations exhibit a variety of dynamics, ranging from fairly stable with little change in population size over time to high-amplitude cyclic or erratic fluctuations. A persistent question in population ecology is why populations fluctuate as they do. Answering this seemingly simple question has proven to be challenging. Broadly, density-dependent feedback mechanisms should allow populations to grow at low density and slow or halt growth at high density. However, experimental tests of what demographic processes result in density-dependent feedback and on what timescale have proven elusive. Here, we used replicated density perturbation experiments and capture-mark-recapture analyses to test density-dependent population growth in populations of meadow voles (Microtus pennsylvanicus) during the summer breeding season by manipulating founding population density and observing the pattern of survival, reproduction, and population growth. High population density had no consistent effect on survival rates but generally negatively influenced recruitment and population growth rates. However, these density-dependent effects varied within the breeding season and across years. Our study provides evidence that density-dependent feedback mechanisms operate at finer time scales than previously believed and that process, additively with delayed year effects, is key to understanding multiyear population demography.

Integrating eco-evolutionary dynamics and modern coexistence theory

Community ecology typically assumes that competitive exclusion and species coexistence are unaffected by evolution on the time scale of ecological dynamics. However, recent studies suggest that rapid evolution operating concurrently with competition may enable species coexistence. Such findings necessitate general theory that incorporates the coexistence contributions of eco-evolutionary processes in parallel with purely ecological mechanisms and provides metrics for quantifying the role of evolution in shaping competitive outcomes in both modelling and empirical contexts. To foster the development of such theory, here we extend the interpretation of the two principal metrics of modern coexistence theory-niche and competitive ability differences-to systems where competitors evolve. We define eco-evolutionary versions of these metrics by considering how invading and resident species adapt to conspecific and heterospecific competitors. We show that the eco-evolutionary niche and competitive ability differences are sums of ecological and evolutionary processes, and that they accurately predict the potential for stable coexistence in previous theoretical studies of eco-evolutionary dynamics. Finally, we show how this theory frames recent empirical assessments of rapid evolution effects on species coexistence, and how empirical work and theory on species coexistence and eco-evolutionary dynamics can be further integrated.

From pattern to process? Dual travelling waves, with contrasting propagation speeds, best describe a self-organised spatio-temporal pattern in population growth of a cyclic rodent

The dynamics of cyclic populations distributed in space result from the relative strength of synchronising influences and the limited dispersal of destabilising factors (activators and inhibitors), known to cause multi‐annual population cycles. However, while each of these have been well studied in isolation, there is limited empirical evidence of how the processes of synchronisation and activation–inhibition act together, largely owing to the scarcity of datasets with sufficient spatial and temporal scale and resolution. We assessed a variety of models that could be underlying the spatio‐temporal pattern, designed to capture both theoretical and empirical understandings of travelling waves using large‐scale (>35,000 km²), multi‐year (2011–2017) field monitoring data on abundances of common vole (Microtus arvalis), a cyclic agricultural rodent pest. We found most support for a pattern formed from the summation of two radial travelling waves with contrasting speeds that together describe population growth rates across the region.

Model responses to CO2 and warming are underestimated without explicit representation of Arctic small-mammal grazing

We use a simple model of coupled carbon and nitrogen cycles in terrestrial ecosystems to examine how "explicitly representing grazers" vs. "having grazer effects implicitly aggregated in with other biogeochemical processes in the model" alters predicted responses to elevated carbon dioxide and warming. The aggregated approach can affect model predictions because grazer-mediated processes can respond differently to changes in climate compared with the processes with which they are typically aggregated. We use small-mammal grazers in a tundra as an example and find that the typical three-to-four-year cycling frequency is too fast for the effects of cycle peaks and troughs to be fully manifested in the ecosystem biogeochemistry. We conclude that implicitly aggregating the effects of small-mammal grazers with other processes results in an underestimation of ecosystem response to climate change, relative to estimations in which the grazer effects are explicitly represented. The magnitude of this underestimation increases with grazer density. We therefore recommend that grazing effects be incorporated explicitly when applying models of ecosystem response to global change.

Regime shift tipping point in hare population collapse associated with climatic and agricultural change during the very early 20th century

Animal populations at northern latitudes may have cyclical dynamics that are degraded by climate change leading to trophic cascade. Hare populations at more southerly latitudes are characterized by dramatic declines in abundance associated with agricultural intensification. We focus on the impact of historical climatic and agricultural change on a mid-latitude population of mountain hares, Lepus timidus hibernicus. Using game bag records from multiple sites throughout Ireland, the hare population index exhibited a distinct regime shift. Contrary to expectations, there was a dynamical structure typical of northern latitude hare populations from 1853 to 1908, during which numbers were stable but cyclic with a periodicity of 8 years. This regime was replaced by dynamics more typical of southern latitude hare populations from 1909 to 1970, in which cycles were lost and numbers declined dramatically. Destabilization of the autumn North Atlantic Oscillation (NAO) led to the collapse of similar cycles in the hare population, coincident with the onset of agricultural intensification (a shift from small-to-large farms) in the first half of the 20th century. Similar, but more recent regime shifts have been observed in Arctic ecosystems and attributed to anthropogenic climate change. The present study suggests such shifts may have occurred at lower latitudes more than a century ago during the very early 20th century. It seems likely that similar tipping points in the population collapse of other farmland species may have occurred similarly early but went undocumented. As northern systems are increasingly impacted by climate change and probable expansion of agriculture, the interaction of these processes is likely to disrupt the pulsed flow of resources from cyclic populations impacting ecosystem function.

A mechanism for population self-regulation: Social density suppresses GnRH expression and reduces reproductivity in voles

Nearly 100 years ago, Charles Elton described lemming and vole population cycles as ecological models for understanding population regulation in nature. Yet, the mechanisms driving these cycles are still not fully understood. These rodent populations can continue to cycle in the absence of predation and with food supplementation, and represent a major unsolved problem in population ecology. It has been hypothesized that the social environment at high population density can drive selection for a low-reproduction phenotype, resulting in population self-regulation as an intrinsic mechanism driving the cycles. However, a physiological mechanism for this self-regulation has not been demonstrated. We manipulated population density in wild meadow voles Microtus pennsylvanicus using large-scale field enclosures over 3 years and examined reproductive performance and physiology. Within the field enclosures, we assessed the proportion of breeding animals, mass at sexual maturation, and faecal androgen and oestrogen metabolites. We then collected brain tissue from juvenile voles born at high or low density, quantified mRNA expression of gonadotropin-releasing hormone (GnRH) and oestrogen receptor alpha (ERα) and measured DNA methylation at six CpG sites in a region that was highly conserved with the mouse GnRH promoter. At high density, there was a lower proportion of reproductive animals. Juvenile voles born at high densities had reduced expression of GnRH in the hypothalamus, accompanied by marginally lower faecal sex hormone metabolites. Female juvenile voles born at high density also had higher methylation levels at two CpG sites while males did not, aligning with prior observations that females (but not males) from high-density environments retain reduced reproduction long term. Our results support a physiological basis for population self-regulation in vole cycles, as altering population density alone induced reproductive downregulation at the hypothalamic level. Our results demonstrate that altering the early-life social environment can fundamentally impact reproductive function in the brain. This, in turn, can drive population demography changes in wild animals.

Impacts of predator-mediated interactions along a climatic gradient on the population dynamics of an alpine bird

According to classic theory, species' population dynamics and distributions are less influenced by species interactions under harsh climatic conditions compared to under more benign climatic conditions. In alpine and boreal ecosystems in Fennoscandia, the cyclic dynamics of rodents strongly affect many other species, including ground-nesting birds such as ptarmigan. According to the ‘alternative prey hypothesis’ (APH), the densities of ground-nesting birds and rodents are positively associated due to predator–prey dynamics and prey-switching. However, it remains unclear how the strength of these predator-mediated interactions change along a climatic harshness gradient in comparison with the effects of climatic variation. We built a hierarchical Bayesian model to estimate the sensitivity of ptarmigan populations to interannual variation in climate and rodent occurrence across Norway during 2007–2017. Ptarmigan abundance was positively linked with rodent occurrence, consistent with the APH. Moreover, we found that the link between ptarmigan abundance and rodent dynamics was strongest in colder regions. Our study highlights how species interactions play an important role in population dynamics of species at high latitudes and suggests that they can become even more important in the most climatically harsh regions.

Age structured discrete-time disease models with demographic population cycles

We use juvenile-adult discrete-time infectious disease models with intrinsically generated demographic population cycles to study the effects of age structure on the persistence or extinction of disease and the basic reproduction number, [Formula: see text]. Our juvenile-adult Susceptible-Infectious-Recovered (SIR) and Infectious-Salmon Anemia-Virus (ISA[Formula: see text] models share a common disease-free system that exhibits equilibrium dynamics for the Beverton-Holt recruitment function. However, when the recruitment function is the Ricker model, a juvenile-adult disease-free system exhibits a range of dynamic behaviours from stable equilibria to deterministic period k population cycles to Neimark-Sacker bifurcations and deterministic chaos. For these two models, we use an extension of the next generation matrix approach for calculating [Formula: see text] to account for populations with locally asymptotically stable period k cycles in the juvenile-adult disease-free system. When [Formula: see text] and the juvenile-adult demographic system (in the absence of the disease) has a locally asymptotically stable period k population cycle, we prove that the juvenile-adult disease goes extinct whenever [Formula: see text]. Under the same period k juvenile-adult demographic assumption but with [Formula: see text], we prove that the juvenile-adult disease-free period k population cycle is unstable and the disease persists. When [Formula: see text], our simulations show that the juvenile-adult disease-free period k cycle dynamics drives the juvenile-adult SIR disease dynamics, but not the juvenile-adult ISAv disease dynamics.

End-user involvement to improve predictions and management of populations with complex dynamics and multiple drivers

Sustainable management of wildlife populations can be aided by building models that both identify current drivers of natural dynamics and provide near-term predictions of future states. We employed a Strategic Foresight Protocol (SFP) involving stakeholders to decide the purpose and structure of a dynamic state-space model for the population dynamics of the Willow Ptarmigan, a popular game species in Norway. Based on local knowledge of stakeholders, it was decided that the model should include food web interactions and climatic drivers to provide explanatory predictions. Modeling confirmed observations from stakeholders that climate change impacts Ptarmigan populations negatively through intensified outbreaks of insect defoliators and later onset of winter. Stakeholders also decided that the model should provide anticipatory predictions. The ability to forecast population density ahead of the harvest season was valued by the stakeholders as it provides the management extra time to consider appropriate harvest regulations and communicate with hunters prior to the hunting season. Overall, exploring potential drivers and predicting short-term future states, facilitate collaborative learning and refined data collection, monitoring designs, and management priorities. Our experience from adapting a SFP to a management target with inherently complex dynamics and drivers of environmental change, is that an open, flexible, and iterative process, rather than a rigid step-wise protocol, facilitates rapid learning, trust, and legitimacy.

Impact of climate change on the small mammal community of the Yukon boreal forest

Long‐term monitoring is critical to determine the stability and sustainability of wildlife populations, and if change has occurred, why. We have followed population density changes in the small mammal community in the boreal forest of the southern Yukon for 46 years with density estimates by live trapping on 3–5 unmanipulated grids in spring and autumn. This community consists of 10 species and was responsible for 9% of the energy flow in the herbivore component of this ecosystem from 1986 to 1996, but this increased to 38% from 2003 to 2014. Small mammals, although small in size, are large in the transfer of energy from plants to predators and decomposers. Four species form the bulk of the biomass. There was a shift in the dominant species from the 1970s to the 2000s, with Myodes rutilus increasing in relative abundance by 22% and Peromyscus maniculatus decreasing by 22%. From 2007 to 2018, Myodes comprised 63% of the catch, Peromyscus 20%, and Microtus species 17%. Possible causes of these changes involve climate change, which is increasing primary production in this boreal forest, and an associated increase in the abundance of 3 rodent predators, marten (Martes americana), ermine (Mustela ermine) and coyotes (Canis latrans). Following and understanding these and potential future changes will require long‐term monitoring studies on a large scale to measure metapopulation dynamics. The small mammal community in northern Canada is being affected by climate change and cannot remain stable. Changes will be critically dependent on food–web interactions that are species‐specific.

Voles and weasels in the boreal Fennoscandian small mammal community: what happens if the least weasel disappears due to climate change?

Climate change, habitat loss and fragmentation are major threats for populations and a challenge for individual behavior, interactions and survival. Predator–prey interactions are modified by climate processes. In the northern latitudes, strong seasonality is changing and the main predicted feature is shortening and instability of winter. Vole populations in the boreal Fennoscandia exhibit multiannual cycles. High amplitude peak numbers of voles and dramatic population lows alternate in 3–5‐year cycles shortening from North to South. One key factor, or driver, promoting the population crash and causing extreme extended lows, is suggested to be predation by the least weasel. We review the arms race between prey voles and weasels through the multiannual density fluctuation, affected by climate change, and especially the changes in the duration and stability of snow cover. For ground‐dwelling small mammals, snow provides thermoregulation and shelter for nest sites, and helps them hide from predators. Predicted increases in the instability of winter forms a major challenge for species with coat color change between brown summer camouflage and white winter coat. One of these is the least weasel, Mustela nivalis nivalis. Increased vulnerability of wrong‐colored weasels to predation affects vole populations and may have dramatic effects on vole dynamics. It may have cascading effects on other small rodent–predator interactions and even on plant–animal interactions and forest dynamics.

Numerical response of a mammalian specialist predator to multiple prey dynamics in Mediterranean farmlands

The study of rodent population cycles has greatly contributed, both theoretically and empirically, to our understanding of the circumstances under which predator-prey interactions destabilize populations. According to the specialist predator hypothesis, reciprocal interactions between voles and small predators that specialize on voles, such as weasels, can cause multiannual cycles. A fundamental feature of classical weasel-vole models is a long time-lag in the numerical response of the predator to variations in prey abundance: weasel abundance increases with that of voles and peaks approximately 1 yr later. We investigated the numerical response of the common weasel (Mustela nivalis) to fluctuating abundances of common voles (Microtus arvalis) in recently colonized agrosteppes of Castilla-y-Léon, northwestern Spain, at the southern limit of the species' range. Populations of both weasels and voles exhibited multiannual cycles with a 3-yr period. Weasels responded quickly and numerically to changes in common-vole abundance, with a time lag between prey and weasel abundance that did not exceed 4 months and occurred during the breeding season, reflecting the quick conversion of prey into predator offspring and/or immigration to sites with high vole populations. We found no evidence of a sustained, high weasel abundance following vole abundance peaks. Weasel population growth rates showed spatial synchrony across study sites approximately 60 km apart. Weasel dynamics were more synchronized with that of common voles than with other prey species (mice or shrews). However, asynchrony within, as well as among sites, in the abundance of voles and alternative prey suggests that weasel mobility could allow them to avoid starvation during low-vole phases, precluding the emergence of prolonged time lag in the numerical response to voles. Our observations are inconsistent with the specialist predator hypothesis as currently formulated, and suggest that weasels might follow rather than cause the vole cycles in northwestern Spain. The reliance of a specialized predator on a functional group of prey such as small rodents does not necessarily lead to a long delay in the numerical response by the predator, depending on the spatial and interspecific synchrony in prey dynamics.

Environmental degradation amplifies species' responses to temperature variation in a trophic interaction

Land‐use and climate change are two of the primary drivers of the current biodiversity crisis. However, we lack understanding of how single‐species and multispecies associations are affected by interactions between multiple environmental stressors.

We address this gap by examining how environmental degradation interacts with daily stochastic temperature variation to affect individual life history and population dynamics in a host–parasitoid trophic interaction, using the Indian meal moth, Plodia interpunctella, and its parasitoid wasp Venturia canescens.

We carried out a single‐generation individual life‐history experiment and a multigeneration microcosm experiment during which individuals and microcosms were maintained at a mean temperature of 26°C that was either kept constant or varied stochastically, at four levels of host resource degradation, in the presence or absence of parasitoids.

At the individual level, resource degradation increased juvenile development time and decreased adult body size in both species. Parasitoids were more sensitive to temperature variation than their hosts, with a shorter juvenile stage duration than in constant temperatures and a longer adult life span in moderately degraded environments. Resource degradation also altered the host's response to temperature variation, leading to a longer juvenile development time at high resource degradation. At the population level, moderate resource degradation amplified the effects of temperature variation on host and parasitoid populations compared with no or high resource degradation and parasitoid overall abundance was lower in fluctuating temperatures. Top‐down regulation by the parasitoid and bottom‐up regulation driven by resource degradation contributed to more than 50% of host and parasitoid population responses to temperature variation.

Our results demonstrate that environmental degradation can strongly affect how species in a trophic interaction respond to short‐term temperature fluctuations through direct and indirect trait‐mediated effects. These effects are driven by species differences in sensitivity to environmental conditions and modulate top‐down (parasitism) and bottom‐up (resource) regulation. This study highlights the need to account for differences in the sensitivity of species’ traits to environmental stressors to understand how interacting species will respond to simultaneous anthropogenic changes.

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.

Predator-prey feedback in a gyrfalcon-ptarmigan system?

Specialist predators with oscillating dynamics are often strongly affected by the population dynamics of their prey, yet they are not always the cause of prey cycling. Only those that exert strong (delayed) regulation of their prey can be. Inferring predator-prey coupling from time series therefore requires contrasting models with top-down versus bottom-up predator-prey dynamics. We study here the joint dynamics of population densities of the Icelandic gyrfalcon Falco rusticolus, and its prey, the rock ptarmigan Lagopus muta. The dynamics of both species are likely not only linked to each other but also to stochastic weather variables acting as confounding factors. We infer the degree of coupling between populations, as well as forcing by abiotic variables, using multivariate autoregressive models MAR(p), with p = 1 and 2 time lags. MAR(2) models, allowing for species to cycle independently from each other, further suggest alternative scenarios where a cyclic prey influences its predator but not the other way around (i.e., bottom-up scenarios). The classical MAR(1) model predicts that the time series exhibit predator-prey feedback (i.e., reciprocal dynamic influence between prey and predator), and that weather effects are weak and only affecting the gyrfalcon population. Bottom-up MAR(2) models produced a better fit but less realistic cross-correlation patterns. Simulations of MAR(1) and MAR(2) models further demonstrate that the top-down MAR(1) models are more likely to be misidentified as bottom-up dynamics than vice versa. We therefore conclude that predator-prey feedback in the gyrfalcon-ptarmigan system is likely the main cause of observed oscillations, though bottom-up dynamics cannot yet be excluded with certainty. Overall, we showed how to make more out of ecological time series by using simulations to gauge the quality of model identification, and paved the way for more mechanistic modeling of this system by narrowing the set of important biotic and abiotic drivers.

The relative strengths of rapid and delayed density dependence acting on a terrestrial herbivore change along a pollution gradient

Animal populations vary in response to a combination of density-dependent and density-independent forces, which interact to drive their population dynamics. Understanding how abiotic forces mediate the form and strength of density-dependent processes remains a central goal of ecology, and is of increasing urgency in a rapidly changing world. Here, we report for the first time that industrial pollution determines the relative strength of rapid and delayed density dependence operating on an animal population. We explored the impacts of pollution and climate on the population dynamics of an eruptive leafmining moth, Phyllonorycter strigulatella, around a coal-fired power plant near Apatity, north-western Russia. Populations were monitored at 14 sites over 26 years. The relative strengths of rapid and delayed density dependence varied with distance from the power plant. Specifically, the strength of rapid density dependence increased while the strength of delayed density dependence decreased with increasing distance from the pollution source. Paralleling the increasing strength of rapid density dependence, we observed declines in the densities of P. strigulatella, increases in predation pressure from birds and ants, and declines in an unknown source of mortality (perhaps plant antibiosis) with increasing distance from the power plant. In contrast to the associations with pollution, associations between climate change and leafminer population densities were negligible. Our results may help to explain the outbreaks of insect herbivores that are frequently observed in polluted environments. We show that they can result from the weakening of rapid (stabilizing) density dependence relative to the effects of destabilizing delayed density dependence. Moreover, our results may explain some of the variation reported in published studies of animal populations in polluted habitats. Variable results may emerge in part because of the location of the study sites on different parts of pollution gradients. Finally, in a rapidly changing world, effects of anthropogenic pollution may be as, or more, important than are effects of climate change on the future dynamics of animal populations.

Aquatic versus Terrestrial Insects: Real or Presumed Differences in Population Dynamics?

The study of insect populations is dominated by research on terrestrial insects. Are aquatic insect populations different or are they just presumed to be different? We explore the evidence across several topics. (1) Populations of terrestrial herbivorous insects are constrained most often by enemies, whereas aquatic herbivorous insects are constrained more by food supplies, a real difference related to the different plants that dominate in each ecosystem. (2) Population outbreaks are presumed not to occur in aquatic insects. We report three examples of cyclical patterns; there may be more. (3) Aquatic insects, like terrestrial insects, show strong oviposition site selection even though they oviposit on surfaces that are not necessarily food for their larvae. A novel outcome is that density of oviposition habitat can determine larval densities. (4) Aquatic habitats are often largely 1-dimensional shapes and this is presumed to influence dispersal. In rivers, drift by insects is presumed to create downstream dispersal that has to be countered by upstream flight by adults. This idea has persisted for decades but supporting evidence is scarce. Few researchers are currently working on the dynamics of aquatic insect populations; there is scope for many more studies and potentially enlightening contrasts with terrestrial insects.

Population-level manipulations of field vole densities induce subsequent changes in plant quality but no impacts on vole demography

Grazing-induced changes in plant quality have been suggested to drive the negative delayed density dependence exhibited by many herbivore species, but little field evidence exists to support this hypothesis. We tested a key premise of the hypothesis that reciprocal feedback between vole grazing pressure and the induction of anti-herbivore silicon defenses in grasses drives observed population cycles in a large-scale field experiment in northern England. We repeatedly reduced population densities of field voles (Microtus agrestis) on replicated 1-ha grassland plots at Kielder Forest, northern England, over a period of 1 year. Subsequently, we tested for the impact of past density on vole life history traits in spring, and whether these effects were driven by induced silicon defenses in the voles' major over-winter food, the grass Deschampsia caespitosa. After several months of density manipulation, leaf silicon concentrations diverged and averaged 22% lower on sites where vole density had been reduced, but this difference did not persist beyond the period of the density manipulations. There were no significant effects of our density manipulations on vole body mass, spring population growth rate, or mean date for the onset of spring reproduction the following year. These findings show that grazing by field voles does induce increased silicon defenses in grasses at a landscape scale. However, at the vole densities encountered, levels of plant damage appear to be below those needed to induce changes in silicon levels large and persistent enough to affect vole performance, confirming the threshold effects we have previously observed in laboratory-based studies. Our findings do not support the plant quality hypothesis for observed vole population cycles in northern England, at least over the range of vole densities that now prevail here.

Mechanism matters: the cause of fluctuations in boom-bust populations governs optimal habitat restoration strategy

Many populations exhibit boom-bust dynamics in which abundance fluctuates dramatically over time. Past research has focused on identifying whether the cause of fluctuations is primarily exogenous, e.g., environmental stochasticity coupled with weak density dependence, or endogenous, e.g., over-compensatory density dependence. Far fewer studies have addressed whether the mechanism responsible for boom-bust dynamics matters with respect to at-risk species management. Here, we ask whether the best strategy for restoring habitat across a landscape differs under exogenously vs. endogenously driven boom-bust dynamics. We used spatially explicit individual-based models to assess how butterfly populations governed by the two mechanisms would respond to habitat restoration strategies that varied in the level of resource patchiness, from a single large patch to multiple patches spaced at different distances. Our models showed that the restoration strategy that minimized extinction risk and boom-bust dynamics would be markedly different depending on the governing mechanism. Exogenously governed populations fared best in a single large habitat patch, whereas for endogenously driven populations, boom-bust dynamics were dampened and extinction risk declined when the total restored area was split into multiple patches with low to moderate inter-patch spacing. Adding environmental stochasticity to the endogenous model did not alter this result. Habitat fragmentation lowered extinction risk in the endogenously driven populations by reducing their growth rate, precluding both "boom" phases and, more importantly, "bust" phases. Our findings suggest that (1) successful restoration will depend on understanding the causes of fluctuations in at-risk populations, (2) the level and pattern of spatiotemporal environmental heterogeneity will also affect the ideal management approach, and (3) counterintuitively, for at-risk species with endogenously governed boom-bust dynamics, lowering the intrinsic population growth rate may decrease extinction risk.

Impact of rewilding, species introductions and climate change on the structure and function of the Yukon boreal forest ecosystem

Community and ecosystem changes are happening in the pristine boreal forest ecosystem of the Yukon for 2 reasons. First, climate change is affecting the abiotic environment (temperature, rainfall and growing season) and driving changes in plant productivity and predator–prey interactions. Second, simultaneously change is occurring because of mammal species reintroductions and rewilding. The key ecological question is the impact these faunal changes will have on trophic dynamics. Primary productivity in the boreal forest is increasing because of climatic warming, but plant species composition is unlikely to change significantly during the next 50–100 years. The 9–10‐year population cycle of snowshoe hares will persist but could be reduced in amplitude if winter weather increases predator hunting efficiency. Small rodents have increased in abundance because of increased vegetation growth. Arctic ground squirrels have disappeared from the forest because of increased predator hunting efficiency associated with shrub growth. Reintroductions have occurred for 2 reasons: human reintroductions of large ungulates and natural recolonization of mammals and birds extending their geographic ranges. The deliberate rewilding of wood bison (Bison bison) and elk (Cervus canadensis) has changed the trophic structure of this boreal ecosystem very little. The natural range expansion of mountain lions (Puma concolor), mule deer (Odocoileus hemionus) and American marten (Martes americana) should have few ecosystem effects. Understanding potential changes will require long‐term monitoring studies and experiments on a scale we rarely deem possible. Ecosystems affected by climate change, species reintroductions and human alteration of habitats cannot remain stable and changes will be critically dependent on food web interactions.

Modification of a Pollen Trap Design To Capture Airborne Conidia of Entomophaga maimaiga and Detection of Conidia by Quantitative PCR

The goal of this study was to develop effective and practical field sampling methods for quantification of aerial deposition of airborne conidia of Entomophaga maimaiga over space and time. This important fungal pathogen is a major cause of larval death in invasive gypsy moth (Lymantria dispar) populations in the United States. Airborne conidia of this pathogen are relatively large (similar in size to pollen), with unusual characteristics, and require specialized methods for collection and quantification. Initially, dry sampling (settling of spores from the air onto a dry surface) was used to confirm the detectability of E. maimaiga at field sites with L. dispar deaths caused by E. maimaiga, using quantitative PCR (qPCR) methods. We then measured the signal degradation of conidial DNA on dry surfaces under field conditions, ultimately rejecting dry sampling as a reliable method due to rapid DNA degradation. We modified a chamber-style trap commonly used in palynology to capture settling spores in buffer. We tested this wet-trapping method in a large-scale (137-km) spore-trapping survey across gypsy moth outbreak regions in Pennsylvania undergoing epizootics, in the summer of 2016. Using 4-day collection periods during the period of late instar and pupal development, we detected variable amounts of target DNA settling from the air. The amounts declined over the season and with distance from the nearest defoliated area, indicating airborne spore dispersal from outbreak areas.IMPORTANCE We report on a method for trapping and quantifying airborne spores of Entomophaga maimaiga, an important fungal pathogen affecting gypsy moth (Lymantria dispar) populations. This method can be used to track dispersal of E. maimaiga from epizootic areas and ultimately to provide critical understanding of the spatial dynamics of gypsy moth-pathogen interactions.

Coast-wide recruitment dynamics of Olympia oysters reveal limited synchrony and multiple predictors of failure

Recruitment of new propagules into a population can be a critical determinant of adult density. We examined recruitment dynamics in the Olympia oyster (Ostrea lurida), a species occurring almost entirely in estuaries. We investigated spatial scales of interannual synchrony across 37 sites in eight estuaries along 2,500 km of Pacific North American coastline, predicting that high vs. low recruitment years would coincide among neighboring estuaries due to shared exposure to regional oceanographic factors. Such synchrony in recruitment has been found for many marine species and some migratory estuarine species, but has never been examined across estuaries in a species that can complete its entire life cycle within the same estuary. To inform ongoing restoration efforts for Olympia oysters, which have declined in abundance in many estuaries, we also investigated predictors of recruitment failure. We found striking contrasts in absolute recruitment rate and frequency of recruitment failure among sites, estuaries, and years. Although we found a positive relationship between upwelling and recruitment, there was little evidence of synchrony in recruitment among estuaries along the coast, and only limited synchrony of sites within estuaries, suggesting recruitment rates are affected more strongly by local dynamics within estuaries than by regional oceanographic factors operating at scales encompassing multiple estuaries. This highlights the importance of local wetland and watershed management for the demography of oysters, and perhaps other species that can complete their entire life cycle within estuaries. Estuaries with more homogeneous environmental conditions had greater synchrony among sites, and this led to the potential for estuary-wide failure when all sites had no recruitment in the same year. Environmental heterogeneity within estuaries may thus buffer against estuary-wide recruitment failure, analogous to the portfolio effect for diversity. Recruitment failure was correlated with lower summer water temperature, higher winter salinity, and shorter residence time: all indicators of stronger marine influence on estuaries. Recruitment failure was also more common in estuaries with limited networks of nearby adult oysters. Large existing oyster networks are thus of high conservation value, while estuaries that lack them would benefit from restoration efforts to increase the extent and connectivity of sites supporting oysters.

Negative effects of density on space use of small mammals differ with the phase of the masting-induced population cycle

Home range size generally decreases with increasing population density, but testing how this relationship is influenced by other factors (e.g., food availability, kin structure) is a difficult task. We used spatially explicit capture–recapture models to examine how home range size varies with population density in the yellow‐necked mouse (Apodemus flavicollis). The relationship between population density and home range size was studied at two distinct phases of population fluctuations induced by beech (Fagus sylvatica) masting: post‐mast peak in abundance (first summer after mast, n = 2) and subsequent crash (second summer after mast, n = 2). We live‐trapped mice from June to September to avoid the confounding effects of autumn seedfall on home range size. In accordance with general predictions, we found that home range size was negatively associated with population density. However, after controlling for the effect of density, home ranges of mice were larger in post‐mast years than during the crash phase. This indicates a higher spatial overlap among neighbors in post‐mast years. We suggest that the increased spatial overlap is caused by negative density‐dependent dispersal that leads to high relatedness of individuals within population in the peak phase of the cycle.

Demographic differences in diet breadth of Canada lynx during a fluctuation in prey availability

Population dynamics of specialist carnivores are closely linked to prey availability, but the extent of variability in diet breadth of individual carnivores relative to natural variability in the abundance of their primary prey is not well understood. Canada lynx (Lynx canadensis) specialize on snowshoe hares (Lepus americanus) and exhibit cyclic fluctuations in abundance that lag 1–2 years behind those of snowshoe hares. Declining hare densities spur demographic changes in lynx, but it is unclear whether a corresponding increase in diet breadth occurs: (1) broadly across a lynx population; (2) only among individuals who are able to effectively switch to alternative prey; or (3) only among individuals who cannot capture sufficient primary prey. We measured stable isotope ratios of lynx muscle tissue spanning a cyclic increase and decline in hare density (1998–2001) in Fort Providence, NT, Canada. We found that lynx cohorts responded differently to hare population change, with yearling animals having broader diets at low hare densities, while adults and dependent juveniles maintained a constant diet through the initial decline in hare density. This result was consistent irrespective of lynx sex and indicates that yearling lynx likely are forced to adopt a broader diet when primary prey densities decline. Our results imply that select cohorts of specialist carnivores can exhibit high dietary plasticity in response to changes in primary prey abundance, prompting the need to determine whether increased diet breadth in young lynx is a successful strategy for surviving through periods of snowshoe hare scarcity. In this way, cohort‐specific niche expansion could strongly affect the dynamics of organisms exhibiting population cycles.

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.

De novo lipogenesis is suppressed during fasting but upregulated at population decline in cyclic voles

Arvicolines are susceptible to the development of fatty liver during short-term fasting. We examined the potential role of de novo lipogenesis (DNL) (i) in the development of fasting-induced fatty liver and (ii) during a population cycle by measuring the mRNA expression of acetyl-CoA carboxylase-1 (ACC1) and fatty acid synthase (FAS). Laboratory voles (Microtus oeconomus and Microtus arvalis) were fed or fasted for 12 or 18 h and their liver mRNA levels were determined. Both species showed decreased mRNA expression of ACC1 and FAS during fasting. This suggests that DNL does not participate in the development of fatty liver in voles, different from human non-alcoholic fatty liver disease. In wild bank voles (Myodes glareolus), the mRNA levels of the genes of interest were higher during the population decline compared to the increase phase. In conclusion, DNL was suppressed during acute fasting but upregulated during a long-term population decline-a period of purported scarcity of high-quality food.

Interactions between pupae of the pine processionary moth (Thaumetopoea pityocampa) and parasitoids in a Pinus forest

Parasitoids are significant enemies of many economically important insects and there is some evidence to suggest that their actions have a role in terminating the outbreaks of forest Lepidoptera populations. In this study, we examined the impact of parasitoids on the pupae of the pine processionary moth, and highlighted the presence of several parasitoid species for this developmental stage. A higher rate of parasitism was found when the pupal density in the soil was reduced, but the rate of parasitism was not influenced by pupal morphological traits or by the presence or absence of a cocoon around a pupa. Of the external factors examined, a delay in the time of descent of larvae from the trees had a positive effect on the level of parasitism. Observational data indicated that dipteran and hymenopteran were the most abundant parasitoids to emerge from moth pupae. Our study highlights the complexity of the parasitoid-host dynamics, and stresses the importance of carefully determining environmental effects on host-parasitoid relations.

Snakebites are associated with poverty, weather fluctuations, and El Niño

Snakebites are environmental and occupational health hazards that mainly affect rural populations worldwide. The ectothermic nature of snakes raises the issue of how climate change's impact on snake ecology could influence the incidence of snakebites in humans in ways that echo the increased predation pressure of snakes on their prey. We thus ask whether snakebites reported in Costa Rica from 2005 to 2013 were associated with meteorological fluctuations. We emphasize El Niño Southern Oscillation (ENSO), a climatic phenomenon associated with cycles of other neglected tropical diseases (NTDs) in the region and elsewhere. We ask how spatial heterogeneity in snakebites and poverty are associated, given the importance of the latter for NTDs. We found that periodicity in snakebites reflects snake reproductive phenology and is associated with ENSO. Snakebites are more likely to occur at high temperatures and may be significantly reduced after the rainy season. Nevertheless, snakebites cluster in Costa Rican areas with the heaviest rainfall, increase with poverty indicators, and decrease with altitude. Altogether, our results suggest that snakebites might vary as a result of climate change.

Seasonal demography of a cyclic lemming population in the Canadian Arctic

1. The causes of cyclical fluctuations in animal populations remain a controversial topic in ecology. Food limitation and predation are two leading hypotheses to explain small mammal population dynamics in northern environments. We documented the seasonal timing of the decline phases and demographic parameters (survival and reproduction) associated with population changes in lemmings, allowing us to evaluate some predictions from these two hypotheses. 2. We studied the demography of brown lemmings (Lemmus trimucronatus), a species showing 3- to 4-year population cycles in the Canadian Arctic, by combining capture-mark-recapture analysis of summer live-trapping with monitoring of winter nests over a 10-year period. We also examined the effects of some weather variables on survival. 3. We found that population declines after a peak occurred between the summer and winter period and not during the winter. During the summer, population growth was driven by change in survival, but not in fecundity or proportion of juveniles, whereas in winter population growth was driven by changes in late summer and winter reproduction. 4. We did not find evidence for direct density dependence on summer demographic parameters, though our analysis was constrained by the paucity of data during the low phase. Body mass, however, was highest in peak years. 5. Weather effects were detected only in early summer when lemming survival was positively related to snow depth at the onset of melt but negatively related to rainfall. 6. Our results show that high mortality causes population declines of lemmings during summer and fall, which suggests that predation is sufficient to cause population crashes, whereas high winter fecundity is the primary factor leading to population irruptions. The positive association between snow depth and early summer survival may be due to the protective cover offered by snow against predators. It is still unclear why reproduction remains low during the low phase.

Contrasting effects of large density changes on relative testes size in fluctuating populations of sympatric vole species

Across species, there is usually a positive relationship between sperm competition level and male reproductive effort on ejaculates, typically measured using relative testes size (RTS). Within populations, demographic and ecological processes may drastically alter the level of sperm competition and thus, potentially affect the evolution of testes size. Here, we use longitudinal records (across 38 years) from wild sympatric Fennoscandian populations of five species of voles to investigate whether RTS responds to natural fluctuations in population density, i.e. variation in sperm competition risk. We show that for some species RTS increases with density. However, our results also show that this relationship can be reversed in populations with large-scale between-year differences in density. Multiple mechanisms are suggested to explain the negative RTS-density relationship, including testes size response to density-dependent species interactions, an evolutionary response to sperm competition levels that is lagged when density fluctuations are over a certain threshold, or differing investment in pre- and post-copulatory competition at different densities. The results emphasize that our understanding of sperm competition in fluctuating environments is still very limited.

A boreal invasion in response to climate change? Range shifts and community effects in the borderland between forest and tundra

It has been hypothesized that climate warming will allow southern species to advance north and invade northern ecosystems. We review the changes in the Swedish mammal and bird community in boreal forest and alpine tundra since the nineteenth century, as well as suggested drivers of change. Observed changes include (1) range expansion and increased abundance in southern birds, ungulates, and carnivores; (2) range contraction and decline in northern birds and carnivores; and (3) abundance decline or periodically disrupted dynamics in cyclic populations of small and medium-sized mammals and birds. The first warm spell, 1930-1960, stands out as a period of substantial faunal change. However, in addition to climate warming, suggested drivers of change include land use and other anthropogenic factors. We hypothesize all these drivers interacted, primarily favoring southern generalists. Future research should aim to distinguish between effects of climate and land-use change in boreal and tundra ecosystems.

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.

The interplay between seasonality and density: consequences for female breeding decisions in a small cyclic herbivore

BACKGROUND

Cyclic rodent population dynamics are subjected to both intrinsic regulatory processes such as density-dependence and extrinsic environmental forcing. Among extrinsic factors, seasonal environmental variation is understood to facilitate cycles. In rodents, these processes have been studied mostly independently and their relative importance for population dynamics is poorly known.

RESULTS

We performed a detailed analysis of common vole (Microtus arvalis) reproduction in a cyclic population using a spatially extensive data set over 17 years in central-western France. Environmental seasonality was the main source of explained variation in common vole reproduction. Additionally, inter-annual variation in the environment explained a smaller part of the variance in reproduction in spring and summer than in winter, whereas the effect of density was only found in autumn and winter. In particular, we detected a strong impact of plant productivity on fecundity during the breeding season, with low vegetation productivity being able to bring vole reproduction nearly to a halt. In contrast, vole reproduction during autumn and winter was mainly shaped by intrinsic factors, with only the longer and heavier females being able to reproduce. The effect of population density on reproduction was negative, mediated by direct negative effects on the proportion of breeders in autumn and winter during outbreak years and by a delayed negative effect on litter size the following year.

CONCLUSIONS

During the main breeding season, variability of female vole reproduction is predominantly shaped by food resources, suggesting that only highly productive environment may induce vole outbreaks. During fall and winter, variability of female vole reproduction is mainly controlled by intrinsic factors, with high population density suppressing reproduction. This suggests, in this cyclic population, that negative direct density dependence on reproduction could explain winter declines after outbreaks.

Anatomy of a population cycle: the role of density dependence and demographic variability on numerical instability and periodicity

Determining the causes of cyclic fluctuations in population size is a central tenet in population ecology and provides insights into population regulatory mechanisms. We have a firm understanding of how direct and delayed density dependence affects population stability and cyclic dynamics, but there remains considerable uncertainty in the specific processes contributing to demographic variability and consequent change in cyclic propensity. Spatiotemporal variability in cyclic propensity, including recent attenuation or loss of cyclicity among several temperate populations and the implications of habitat fragmentation and climate change on this pattern, highlights the heightened need to understand processes underlying cyclic variation. Because these stressors can differentially impact survival and productivity and thereby impose variable time delays in density dependence, there is a specific need to elucidate how demographic vital rates interact with the type and action of density dependence to contribute to population stability and cyclic variation. Here, we address this knowledge gap by comparing the stability of time series derived from general and species-specific (Canada lynx: Lynx canadensis; small rodents: Microtus, Lemmus and Clethrionomys spp.) matrix population models, which vary in their demographic rates and the direct action of density dependence. Our results reveal that density dependence acting exclusively on survival as opposed to productivity is destabilizing, suggesting that a shift in the action of population regulation toward reproductive output may decrease cyclic propensity and cycle amplitude. This result was the same whether delayed density dependence was pulsatile and acted on a single time period (e.g. t-1, t-2 or t-3) vs. more constant by affecting a successive range of years (e.g. t-1,…, t-3). Consistent with our general models, reductions in reproductive potential in both the lynx and small rodent systems led to notably large drops in cyclic propensity and amplitude, suggesting that changes in this vital rate may contribute to the spatial or temporal variability observed in the cyclic dynamics of both systems. Collectively, our results reveal that the type of density dependence and its effect on different demographic parameters can profoundly influence numeric stability and cyclic propensity and therefore may shift populations across the cyclic-to-noncyclic boundary.

Linking climate change to population cycles of hares and lynx

The classic 10-year population cycle of snowshoe hares (Lepus americanus, Erxleben 1777) and Canada lynx (Lynx canadensis, Kerr 1792) in the boreal forests of North America has drawn much attention from both population and community ecologists worldwide; however, the ecological mechanisms driving the 10-year cyclic dynamic pattern are not fully revealed yet. In this study, by the use of historic fur harvest data, we constructed a series of generalized additive models to study the effects of density dependence, predation, and climate (both global climate indices of North Atlantic Oscillation index (NAO), Southern Oscillation index (SOI) and northern hemispheric temperature (NHT) and local weather data including temperature, rainfall, and snow). We identified several key pathways from global and local climate to lynx with various time lags: rainfall shows a negative, and snow shows a positive effect on lynx; NHT and NAO negatively affect lynx through their positive effect on rainfall and negative effect on snow; SOI positively affects lynx through its negative effect on rainfall. Direct or delayed density dependency effects, the prey effect of hare on lynx and a 2-year delayed negative effect of lynx on hare (defined as asymmetric predation) were found. The simulated population dynamics is well fitted to the observed long-term fluctuations of hare and lynx populations. Through simulation, we find density dependency and asymmetric predation, only producing damped oscillation, are necessary but not sufficient factors in causing the observed 10-year cycles; while extrinsic climate factors are important in producing and modifying the sustained cycles. Two recent population declines of lynx (1940-1955 and after 1980) were likely caused by ongoing climate warming indirectly. Our results provide an alternative explanation to the mechanism of the 10-year cycles, and there is a need for further investigation on links between disappearance of population cycles and global warming in hare-lynx system.

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.

Dynamic energy budget theory and population ecology: lessons from Daphnia

Dynamic energy budget (DEB) theory offers a perspective on population ecology whose starting point is energy utilization by, and homeostasis within, individual organisms. It is natural to ask what it adds to the existing large body of individual-based ecological theory. We approach this question pragmatically--through detailed study of the individual physiology and population dynamics of the zooplankter Daphnia and its algal food. Standard DEB theory uses several state variables to characterize the state of an individual organism, thereby making the transition to population dynamics technically challenging, while ecologists demand maximally simple models that can be used in multi-scale modelling. We demonstrate that simpler representations of individual bioenergetics with a single state variable (size), and two life stages (juveniles and adults), contain sufficient detail on mass and energy budgets to yield good fits to data on growth, maturation and reproduction of individual Daphnia in response to food availability. The same simple representations of bioenergetics describe some features of Daphnia mortality, including enhanced mortality at low food that is not explicitly incorporated in the standard DEB model. Size-structured, population models incorporating this additional mortality component resolve some long-standing questions on stability and population cycles in Daphnia. We conclude that a bioenergetic model serving solely as a 'regression' connecting organismal performance to the history of its environment can rest on simpler representations than those of standard DEB. But there are associated costs with such pragmatism, notably loss of connection to theory describing interspecific variation in physiological rates. The latter is an important issue, as the type of detailed study reported here can only be performed for a handful of species.