Iterative model predictions for wildlife populations impacted by rapid climate change

To improve understanding and management of the consequences of current rapid environmental change, ecologists advocate using long-term monitoring data series to generate iterative near-term predictions of ecosystem responses. This approach allows scientific evidence to increase rapidly and management strategies to be tailored simultaneously. Iterative near-term forecasting may therefore be particularly useful for adaptive monitoring of ecosystems subjected to rapid climate change. Here, we show how to implement near-term forecasting in the case of a harvested population of rock ptarmigan in high-arctic Svalbard, a region subjected to the largest and most rapid climate change on Earth. We fitted state-space models to ptarmigan counts from point transect distance sampling during 2005-2019 and developed two types of predictions: (1) explanatory predictions to quantify the effect of potential drivers of ptarmigan population dynamics, and (2) anticipatory predictions to assess the ability of candidate models of increasing complexity to forecast next-year population density. Based on the explanatory predictions, we found that a recent increasing trend in the Svalbard rock ptarmigan population can be attributed to major changes in winter climate. Currently, a strong positive effect of increasing average winter temperature on ptarmigan population growth outweighs the negative impacts of other manifestations of climate change such as rain-on-snow events. Moreover, the ptarmigan population may compensate for current harvest levels. Based on the anticipatory predictions, the near-term forecasting ability of the models improved nonlinearly with the length of the time series, but yielded good forecasts even based on a short time series. The inclusion of ecological predictors improved forecasts of sharp changes in next-year population density, demonstrating the value of ecosystem-based monitoring. Overall, our study illustrates the power of integrating near-term forecasting in monitoring systems to aid understanding and management of wildlife populations exposed to rapid climate change. We provide recommendations for how to improve this approach.

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.

Circumpolar status of Arctic ptarmigan: Population dynamics and trends

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

Seasonal difference in temporal transferability of an ecological model: near-term predictions of lemming outbreak abundances

Ecological models have been criticized for a lack of validation of their temporal transferability. Here we answer this call by investigating the temporal transferability of a dynamic state-space model developed to estimate season-dependent biotic and climatic predictors of spatial variability in outbreak abundance of the Norwegian lemming. Modelled summer and winter dynamics parametrized by spatial trapping data from one cyclic outbreak were validated with data from a subsequent outbreak. There was a distinct difference in model transferability between seasons. Summer dynamics had good temporal transferability, displaying ecological models' potential to be temporally transferable. However, the winter dynamics transferred poorly. This discrepancy is likely due to a temporal inconsistency in the ability of the climate predictor (i.e. elevation) to reflect the winter conditions affecting lemmings both directly and indirectly. We conclude that there is an urgent need for data and models that yield better predictions of winter processes, in particular in face of the expected rapid climate change in the Arctic.

The PREDICTS database: a global database of how local terrestrial biodiversity responds to human impacts

Biodiversity continues to decline in the face of increasing anthropogenic pressures such as habitat destruction, exploitation, pollution and introduction of alien species. Existing global databases of species’ threat status or population time series are dominated by charismatic species. The collation of datasets with broad taxonomic and biogeographic extents, and that support computation of a range of biodiversity indicators, is necessary to enable better understanding of historical declines and to project – and avert – future declines. We describe and assess a new database of more than 1.6 million samples from 78 countries representing over 28,000 species, collated from existing spatial comparisons of local-scale biodiversity exposed to different intensities and types of anthropogenic pressures, from terrestrial sites around the world. The database contains measurements taken in 208 (of 814) ecoregions, 13 (of 14) biomes, 25 (of 35) biodiversity hotspots and 16 (of 17) megadiverse countries. The database contains more than 1% of the total number of all species described, and more than 1% of the described species within many taxonomic groups – including flowering plants, gymnosperms, birds, mammals, reptiles, amphibians, beetles, lepidopterans and hymenopterans. The dataset, which is still being added to, is therefore already considerably larger and more representative than those used by previous quantitative models of biodiversity trends and responses. The database is being assembled as part of the PREDICTS project (Projecting Responses of Ecological Diversity In Changing Terrestrial Systems – http://www.predicts.org.uk). We make site-level summary data available alongside this article. The full database will be publicly available in 2015.

Indirect food web interactions mediated by predator-rodent dynamics: relative roles of lemmings and voles

Production cycles in birds are proposed as prime cases of indirect interactions in food webs. They are thought to be driven by predators switching from rodents to bird nests in the crash phase of rodent population cycles. Although rodent cycles are geographically widespread and found in different rodent taxa, bird production cycles appear to be most profound in the high Arctic where lemmings dominate. We hypothesized that this may be due to arctic lemmings inducing stronger predator responses than boreal voles. We tested this hypothesis by estimating predation rates in dummy bird nests during a rodent cycle in low-Arctic tundra. Here, the rodent community consists of a spatially variable mix of one lemming (Lemmus lemmus) and two vole species (Myodes rufocanus and Microtus oeconomus) with similar abundances. In consistence with our hypothesis, lemming peak abundances predicted well crash-phase nest predation rates, whereas the vole abundances had no predictive ability. Corvids were found to be the most important nest predators. Lemmings appear to be accessible to the whole predator community which makes them particularly powerful drivers of food web dynamics.

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.

How spatial variation in areal extent and configuration of labile vegetation states affect the riparian bird community in Arctic tundra

The Arctic tundra is currently experiencing an unprecedented combination of climate change, change in grazing pressure by large herbivores and growing human activity. Thickets of tall shrubs represent a conspicuous vegetation state in northern and temperate ecosystems, where it serves important ecological functions, including habitat for wildlife. Thickets are however labile, as tall shrubs respond rapidly to both abiotic and biotic environmental drivers. Our aim was to assess how large-scale spatial variation in willow thicket areal extent, configuration and habitat structure affected bird abundance, occupancy rates and species richness so as to provide an empirical basis for predicting the outcome of environmental change for riparian tundra bird communities. Based on a 4-year count data series, obtained through a large-scale study design in low arctic tundra in northern Norway, statistical hierarchical community models were deployed to assess relations between habitat configuration and bird species occupancy and community richness. We found that species abundance, occupancy and richness were greatly affected by willow areal extent and configuration, habitat features likely to be affected by intense ungulate browsing as well as climate warming. In sum, total species richness was maximized in large and tall willow patches of small to intermediate degree of fragmentation. These community effects were mainly driven by responses in the occupancy rates of species depending on tall willows for foraging and breeding, while species favouring other vegetation states were not affected. In light of the predicted climate driven willow shrub encroachment in riparian tundra habitats, our study predicts that many bird species would increase in abundance, and that the bird community as a whole could become enriched. Conversely, in tundra regions where overabundance of large herbivores leads to decreased areal extent, reduced height and increased fragmentation of willow thickets, bird community richness and species-specific abundance are likely to be significantly reduced.

The importance of willow thickets for ptarmigan and hares in shrub tundra: the more the better?

In patchy habitats, the relationship between animal abundance and cover of a preferred habitat may change with the availability of that habitat, resulting in a functional response in habitat use. Here, we investigate the relationship of two specialized herbivores, willow ptarmigan (Lagopus lagopus) and mountain hare (Lepus timidus), to willows (Salix spp.) in three regions of the shrub tundra zone-northern Norway, northern European Russia and western Siberia. Shrub tundra is a naturally patchy habitat where willow thickets represent a major structural element and are important for herbivores both as food and shelter. Habitat use was quantified using feces counts in a hierarchical spatial design and related to several measures of willow thicket configuration. We document a functional response in the use of willow thickets by ptarmigan, but not by hares. For hares, whose range extends into forested regions, occurrence increased overall with willow cover. The occurrence of willow ptarmigan showed a strong positive relationship to willow cover and a negative relationship to thicket fragmentation in the region with lowest willow cover at landscape scale, where willow growth may be limited by reindeer browsing. In regions with higher cover, in contrast, such relationships were not observed. Differences in predator communities among the regions may contribute to the observed pattern, enhancing the need for cover where willow thickets are scarce. Such region-specific relationships reflecting regional characteristics of the ecosystem highlight the importance of large-scale investigations to understand the relationships of habitat availability and use, which is a critical issue considering that habitat availability changes quickly with climate change and human impact.

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

1. The geographical pattern in Fennoscandian small rodent population dynamics with a southern noncyclic and a northern cyclic region, and with latitudinal gradients in density-dependent structure, cycle period length and spatial synchrony within the northern cyclic region, has been widely publicized and interpreted in the ecological literature. However, the time-series data on which these inferences have been established are relatively short and originate from a specific time period (mostly around 1970-90). Hence, it can be questioned whether the geographical population dynamics patterns are consistent over time (i.e. whether they are stationary). 2. Here we analyse an almost century long (1880-1976) panel of fox bounty time series including 18 counties of Norway, thus spanning the whole range of latitudes of Fennoscandia (i.e. 15 latitudinal degrees). These fox time series mirror the dynamics of their dominant small rodent prey, in particular, with respect to cycle period length and spatial synchrony. 3. While we found some evidence consistent with previous analyses showing a clearly patterned dynamics according to latitude, such patterns were not stationary on a longer time-scale. In particular, we observed a shift from an extensively synchronous (i.e. regionalized) 4-year cycle north of 60 degrees N just after the 'Little Ice Age' (1880-1910) to a diversification of cycle period length (3-5 years) and eventually, partial loss of cyclicity and synchronicity in later periods. Incidents of loss of cyclicity appeared to be preceded by changes in cycle period (i.e. period lengthening and shortening). 4. These results show that the dynamics of Fennoscandian small rodents, and their associated guild of predators, are more prone to change than previously acknowledged.

Collapsing population cycles

During the past two decades population cycles in voles, grouse and insects have been fading out in Europe. Here, we discuss the cause and implication of these changes. Several lines of evidence now point to climate forcing as the general underlying cause. However, how climate interacts with demography to induce regime shifts in population dynamics is likely to differ among species and ecosystems. Herbivores with high-amplitude population cycles, such as voles, lemmings, snowshoe hares and forest Lepidoptera, form the heart of terrestrial food web dynamics. Thus, collapses of these cycles are also expected to imply collapses of important ecosystem functions, such as the pulsed flows of resources and disturbances.