Six personas to adopt when framing theoretical research questions in biology

Theory is a critical component of the biological research process, and complements observational and experimental approaches. However, most biologists receive little training on how to frame a theoretical question and, thus, how to evaluate when theory has successfully answered the research question. Here, we develop a guide with six verbal framings for theoretical models in biology. These correspond to different personas one might adopt as a theorist: 'Advocate', 'Explainer', 'Instigator', 'Mediator', 'Semantician' and 'Tinkerer'. These personas are drawn from combinations of two starting points (pattern or mechanism) and three foci (novelty, robustness or conflict). We illustrate each of these framings with examples of specific theoretical questions, by drawing on recent theoretical papers in the fields of ecology and evolutionary biology. We show how the same research topic can be approached from slightly different perspectives, using different framings. We show how clarifying a model's framing can debunk common misconceptions of theory: that simplifying assumptions are bad, more detail is always better, models show anything you want and modelling requires substantial maths knowledge. Finally, we provide a roadmap that researchers new to theoretical research can use to identify a framing to serve as a blueprint for their own theoretical research projects.

Dispersal synchronizes giant kelp forests

Spatial synchrony is the tendency for population fluctuations to be correlated among different locations. This phenomenon is a ubiquitous feature of population dynamics and is important for ecosystem stability, but several aspects of synchrony remain unresolved. In particular, the extent to which any particular mechanism, such as dispersal, contributes to observed synchrony in natural populations has been difficult to determine. To address this gap, we leveraged recent methodological improvements to determine how dispersal structures synchrony in giant kelp (Macrocystis pyrifera), a global marine foundation species that has served as a useful system for understanding synchrony. We quantified population synchrony and fecundity with satellite imagery across 11 years and 880 km of coastline in southern California, USA, and estimated propagule dispersal probabilities using a high-resolution ocean circulation model. Using matrix regression models that control for the influence of geographic distance, resources (seawater nitrate), and disturbance (destructive waves), we discovered that dispersal was an important driver of synchrony. Our findings were robust to assumptions about propagule mortality during dispersal and consistent between two metrics of dispersal: (1) the individual probability of dispersal and (2) estimates of demographic connectivity that incorporate fecundity (the number of propagules dispersing). We also found that dispersal and environmental conditions resulted in geographic clusters with distinct patterns of synchrony. This study is among the few to statistically associate synchrony with dispersal in a natural population and the first to do so in a marine organism. The synchronizing effects of dispersal and environmental conditions on foundation species, such as giant kelp, likely have cascading effects on the spatial stability of biodiversity and ecosystem function.

Stabilizing effects of spatially heterogeneous disturbance via reduced spatial synchrony on a rocky shore community

Understanding how synchronous species fluctuations affect community stability is a main research topic in ecology. Yet experimental studies evaluating how changes in disturbance regimes affect the synchrony and stability of populations and communities remain rare. We hypothesized that spatially heterogeneous disturbances of moderate intensity would promote metacommunity stability by decreasing the spatial synchrony of species fluctuations. To test this hypothesis, we exposed rocky shore communities of algae and invertebrates to homogeneous and gradient-like spatial patterns of disturbance at two levels of intensity for 4 years and used synchrony networks to characterize community responses to these disturbances. The gradient-like disturbance at low intensity enhanced spatial β diversity compared to the other treatments and produced the most heterogeneous and least synchronized network, which was also the most stable in terms of population and community fluctuations. In contrast, homogeneous disturbance destabilized the community, enhancing spatial synchronization. Intense disturbances always reduced spatial β diversity, indicating that strong perturbations could destabilize communities via biotic homogenization regardless of their spatial structure. Our findings corroborated theoretical predictions, emphasizing the importance of spatially heterogeneous disturbances in promoting stability by amplifying asynchronous spatial and temporal fluctuations in population and community abundance. In contrast to other networks, synchrony networks are vulnerable to the removal of most peripheral nodes, which are less synchronized, but may contribute more to stability than other nodes by dampening large fluctuations in species abundance. Our findings suggest that climate change and direct anthropogenic disturbance can compromise the stability of ecological communities through combined effects on diversity and synchrony, as well as further affecting ecosystems through habitat loss.

Seasonality in Environment and Population Processes Alters Population Spatial Synchrony

AbstractPopulation spatial synchrony-the tendency for temporal population fluctuations to be correlated across locations-is common and important to metapopulation stability and persistence. One common cause of spatial synchrony, termed the Moran effect, occurs when populations respond to environmental fluctuations, such as weather, that are correlated over space. Although the degree of spatial synchrony in environmental fluctuations can differ between seasons and different population processes occur in different seasons, the impact on population spatial synchrony is uncertain because prior work has largely assumed that the spatial synchrony of environmental fluctuations and their effect on populations are consistent over annual sampling intervals. We used theoretical models to examine how seasonality in population processes and the spatial synchrony of environmental drivers affect population spatial synchrony. We found that population spatial synchrony can depend not only on the spatial synchrony of environmental drivers but also on the degree to which environmental fluctuations are correlated across seasons, locally, and across space. Moreover, measurements of synchrony from "snapshot" population censuses may not accurately reflect synchrony during other parts of the year. Together, these results show that neglecting seasonality in environmental conditions and population processes is consequential for understanding population spatial synchrony and its driving mechanisms.

Forest Landscape Effects on Dispersal of Spruce Budworm Choristoneura fumiferana (Clemens, 1865) (Lepidoptera, Tortricidae) and Forest Tent Caterpillar Malacosoma disstria Hübner, 1820 (Lepidoptera, Lasiocampidae) Female Moths in Alberta, Canada

Leaf-rollers and tent caterpillars, the families Torticidae and Lasiocampidae, represent a significant component of the Lepidoptera, and are well-represented in the forest insect pest literature of North America. Two species in particular-spruce budworm (Choristoneura fumiferana (Clem.)) and forest tent caterpillar (Malacosoma disstria Hbn.)-are the most significant pests of the Pinaceae and Salicacae, respectively, in the boreal forest of Canada, each exhibiting periodic outbreaks of tremendous extent. Dispersal is thought to play a critical role in the triggering of population eruptions and in the synchronization of outbreak cycling, but formal studies of dispersal, in particular studies of long-range dispersal by egg-bearing adult females, are rare. Here, it is shown in two independent studies that adult females of both species tend to disperse away from sparse or defoliated forest, and toward intact or undefoliated forest, suggesting that long-range dispersal during an outbreak peak is adaptive to the species and an important factor in their population dynamics, and hence their evolutionary biology.

Disturbance and nutrients synchronise kelp forests across scales through interacting Moran effects

Spatial synchrony is a ubiquitous and important feature of population dynamics, but many aspects of this phenomenon are not well understood. In particular, it is largely unknown how multiple environmental drivers interact to determine synchrony via Moran effects, and how these impacts vary across spatial and temporal scales. Using new wavelet statistical techniques, we characterised synchrony in populations of giant kelp Macrocystis pyrifera, a widely distributed marine foundation species, and related synchrony to variation in oceanographic conditions across 33 years (1987-2019) and >900 km of coastline in California, USA. We discovered that disturbance (storm-driven waves) and resources (seawater nutrients)-underpinned by climatic variability-act individually and interactively to produce synchrony in giant kelp across geography and timescales. Our findings demonstrate that understanding and predicting synchrony, and thus the regional stability of populations, relies on resolving the synergistic and antagonistic Moran effects of multiple environmental drivers acting on different timescales.

Quantifying and predicting population connectivity of an outbreaking forest insect pest

Context

Dispersal has a key role in the population dynamics of outbreaking species such as the spruce budworm (Choristoneura fumiferana) as it can synchronize the demography of distant populations and favor the transition from endemic to epidemic states. However, we know very little about how landscape structure influences dispersal in such systems while such knowledge is essential for better forecasting of spatially synchronous population dynamics and to guide management strategies.

Objectives

We aimed to characterize the spatial environmental determinants of spruce budworm dispersal to determine how these features affect outbreak spread in Quebec (Canada). We then apply our findings to predict expected future landscape connectivity and explore its potential consequences on future outbreaks.

Methods

We used a machine-learning landscape genetics approach on 447 larvae covering most of the outbreak area and genotyped at 3562 SNP loci to identify the main variables affecting connectivity.

Results

We found that the connectivity between outbreak populations was driven by the combination of precipitation and host cover. Our forecasting suggests that between the current and next outbreaks, connectivity may increase between Ontario and Quebec, and might decrease in the eastern part, which could have the effect of limiting outbreak spread from Ontario and Quebec to the eastern provinces.

Conclusions

Although we did not identify any discrete barriers, low connectivity areas might constrain dispersal in the current and future outbreaks and should in turn, be intensively monitored. However, continued sampling as the outbreak progresses is needed to confirm the temporal stability of the observed patterns.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10980-021-01382-9.

Climate-induced outbreaks in high-elevation pines are driven primarily by immigration of bark beetles from historical hosts

Warming temperatures are allowing native insect herbivores to expand into regions that previously exceeded their thermal tolerance, encounter new host species, and pose significant threats to native communities. However, the dynamics of these expansions remain poorly understood, particularly in the extent to which outbreaks remain reliant on emigration from historical hosts or are driven by local reproduction within novel hosts in the expanded range. We tested these non-mutually exclusive hypotheses using spatially explicit data on mountain pine beetle (Dendroctonus ponderosae), which historically undergoes intermittent outbreaks in low-elevation lodgepole pine (Pinus contorta), but is now causing severe mortality in a high-elevation endangered species, whitebark pine (Pinus albicaulis). We compiled data from 2000 to 2019 across British Columbia, Canada, at 1-km² resolution, and analyzed spatiotemporal patterns of beetle infestations, lodgepole pine distributions, expansion into habitats dominated by whitebark pine, and the likelihood of future outbreaks in all pine communities under simulated conditions. Overall, we found strong support for the hypothesis of emigration from the historical host species continuing to be a major driver of outbreaks in the more recently accessed host. First, beetle population pressure was consistently the best predictor of infestation severity in both lodgepole and whitebark pine, and appeared to be mostly unidirectional from lodgepole to whitebark pine. Second, infestations in lodgepole pine were of a longer duration than those in whitebark pine, which appeared too brief to sustain transitions from endemic to eruptive dynamics. Furthermore, resource depletion appears to drive emigration from lodgepole pine, whereas in whitebark pine drought appears to favor establishment of immigrants although bioclimatic factors and stand structure preclude self-sustaining outbreaks. Finally, we project that most pine in British Columbia will be at risk in the event of a new major outbreak. We describe implications for conserving and protecting whitebark pine and to other climate-driven range expansions.

Landscape-scale population connectivity in two parasitoid species associated with the spruce budworm: Testing the birdfeeder effect using genetic data

Periodic and spatially synchronous outbreaks of insect pests have dramatic consequences for boreal and sub-boreal forests. Within these multitrophic systems, parasitoids can be stabilizing agents by dispersing toward patches containing higher host density (the so-called birdfeeder effect). However, we know little about the dispersal abilities of parasitoids in continuous forested landscapes, limiting our understanding of the spatiotemporal dynamics of host-parasitoid systems, and constraining our ability to predict forest resilience in the context of global changes. In this study, we investigate the spatial genetic structure and spatial variation in genetic diversity of two important species of spruce budworm larval parasitoids during outbreaks: Apanteles fumiferanae Viereck (Braconidae) and Glypta fumiferanae (Viereck) (Ichneumonidae). Using parasitoids sampled in 2014 from 26 and 29 locations across a study area of 350,000 km² , we identified 1,012 and 992 neutral SNP loci for A. fumiferanae (N = 279 individuals) and G. fumiferanae (N = 382), respectively. Using DAPC, PCA, AMOVA, and IBD analyses, we found evidence for panmixia and high genetic connectivity for both species, matching the previously described genetic structure of the spruce budworm within the same context, suggesting similar effective dispersal during outbreaks and high parasitoid population densities between outbreaks. We also found a significant negative relationship between genetic diversity and latitude for A. fumiferanae but not for G. fumiferanae, suggesting that northern range limits may vary by species within the spruce budworm parasitoid community. These spatial dynamics should be considered when predicting future insect outbreak severities in boreal landscapes.

The Vision of Managing for Pest-Resistant Landscapes: Realistic or Utopic?

PURPOSE OF REVIEW

Forest managers have long suggested that forests can be made more resilient to insect pests by reducing the abundance of hosts, yet this has rarely been done. The goal of our paper is to review whether recent scientific evidence supports forest manipulation to decrease vulnerability. To achieve this goal, we first ask if outbreaks of forest insect pests have been more severe in recent decades. Next, we assess the relative importance of climate change and forest management-induced changes in forest composition/structure in driving these changes in severity.

RECENT FINDINGS

Forest structure and composition continue to be implicated in pest outbreak severity. Mechanisms, however, remain elusive. Recent research elucidates how forest compositional and structural diversity at neighbourhood, stand, and landscape scales can increase forest resistance to outbreaks. Many recent outbreaks of herbivorous forest insects have been unprecedented in terms of duration and spatial extent. Climate change may be a contributing factor, but forest structure and composition have been clearly identified as contributing to these unprecedented outbreaks.

SUMMARY

Current research supports using silviculture to create pest-resistant forest landscapes. However, the precise mechanisms by which silviculture can increase resistance remains uncertain. Further, humans tend to more often create pest-prone forests due to political, economic, and human resistance to change and a short-sighted risk management perspective that focuses on reactive rather than proactive responses to insect outbreak threats. Future research efforts need to focus on social, political, cultural, and educational mechanisms to motivate implementation of proven ecological solutions if pest-resistant forests are to be favoured by management.

Pine Looper Bupalus piniaria (L.) Outbreaks Reconstruction: A Case Study for Southern Siberia

Simple Summary

The pine looper damages Scots pine forests over vast areas. However, the study of its population fluctuations is hampered by the lack of long-term observation series. The dendrochronological reconstruction is often used to study the history of its outbreaks. In some cases, such reconstructions require methods that work without comparison with other tree species. We have proposed such a technique based on the analysis of the early- and latewood growth. The technique makes it possible to separate the effect of defoliation on tree rings from weather influence. Besides, it is more sensitive than previously developed methods when reconstructing outbreaks of the pine looper. The history of outbreaks reconstructed by the technique for the West Siberian Plain’s forest-steppe includes 11 defoliation events from 1914 to 2017. The results obtained using the developed method are useful to better understand the patterns of population dynamics of the pine looper and other phyllophagous pests.

Abstract

The pine looper Bupalus piniaria is one of the most widespread phyllophagous insect species across Northern Eurasia, defoliating Scots pine forests over vast territories. Since there are not enough long-term documented observations on a series of outbreaks, there is a need for methods allowing them to be reconstructed to study their dynamics patterns. Previously, dendrochronological methods were successfully used to solve such issues. However, the most common approach is not applicable for the Western Siberian forest-steppe since it requires comparison with a non-damaged tree species close to pine in terms of longevity and resistance to rot. In the pine forests of the steppe and forest-steppe zones of Western Siberia, there are no species that are not damaged by the pine looper that meets these requirements. Methods allowing not using control species are also not free from disadvantages (e.g., weak specificity). Therefore, we have developed a new method based on the analysis, not of the tree-ring width but the early- and latewood width to reconstruct past defoliation events. The past defoliation by the pine looper is indicated by the presence of a negative pointer year for latewood, followed by a negative pointer year for earlywood in a subsequent year among the majority of individuals. Linear modeling showed a difference between the climate impact on radial growth and the defoliation one. The obtained reconstruction was compared with the results of other methods (mowing window, OUTBREAK, independent component analysis), literature, and Forest Service data. The developed new method (pointer year method; PYM) showed high efficiency confirmed by results of the tree-ring series analysis (11 revealed outbreaks in the past). Compared with other reconstruction techniques under the given conditions (a favorable combination of heat and humidity; probably low-intense and short defoliation), the proposed method provided more precise results than those proposed earlier. Due to high accuracy, the PYM can be useful for detecting late-summer and autumn past defoliations of tree species with clear difference between early- and latewood even though the damage was weak.

A general mathematical method for predicting spatio-temporal correlations emerging from agent-based models

Agent-based models are used to study complex phenomena in many fields of science. While simulating agent-based models is often straightforward, predicting their behaviour mathematically has remained a key challenge. Recently developed mathematical methods allow the prediction of the emerging spatial patterns for a general class of agent-based models, whereas the prediction of spatio-temporal pattern has been thus far achieved only for special cases. We present a general and mathematically rigorous methodology that allows deriving the spatio-temporal correlation structure for a general class of individual-based models. To do so, we define an auxiliary model, in which each agent type of the primary model expands to three types, called the original, the past and the new agents. In this way, the auxiliary model keeps track of both the initial and current state of the primary model, and hence the spatio-temporal correlations of the primary model can be derived from the spatial correlations of the auxiliary model. We illustrate the agreement between analytical predictions and agent-based simulations using two example models from theoretical ecology. In particular, we show that the methodology is able to correctly predict the dynamical behaviour of a host–parasite model that shows spatially localized oscillations.

Dynamical Ising model of spatially coupled ecological oscillators

Long-range synchrony from short-range interactions is a familiar pattern in biological and physical systems, many of which share a common set of 'universal' properties at the point of synchronization. Common biological systems of coupled oscillators have been shown to be members of the Ising universality class, meaning that the very simple Ising model replicates certain spatial statistics of these systems at stationarity. This observation is useful because it reveals which aspects of spatial pattern arise independently of the details governing local dynamics, resulting in both deeper understanding of and a simpler baseline model for biological synchrony. However, in many situations a system's dynamics are of greater interest than their static spatial properties. Here, we ask whether a dynamical Ising model can replicate universal and non-universal features of ecological systems, using noisy coupled metapopulation models with two-cycle dynamics as a case study. The standard Ising model makes unrealistic dynamical predictions, but the Ising model with memory corrects this by using an additional parameter to reflect the tendency for local dynamics to maintain their phase of oscillation. By fitting the two parameters of the Ising model with memory to simulated ecological dynamics, we assess the correspondence between the Ising and ecological models in several of their features (location of the critical boundary in parameter space between synchronous and asynchronous dynamics, probability of local phase changes and ability to predict future dynamics). We find that the Ising model with memory is reasonably good at representing these properties of ecological metapopulations. The correspondence between these models creates the potential for the simple and well-known Ising class of models to become a valuable tool for understanding complex biological systems.

Drivers of amphibian population dynamics and asynchrony at local and regional scales

Identifying the drivers of population fluctuations in spatially distinct populations remains a significant challenge for ecologists. Whereas regional climatic factors may generate population synchrony (i.e. the Moran effect), local factors including the level of density dependence may reduce the level of synchrony. Although divergences in the scaling of population synchrony and spatial environmental variation have been observed, the regulatory factors that underlie such mismatches are poorly understood. Few previous studies have investigated how density-dependent processes and population-specific responses to weather variation influence spatial synchrony at both local and regional scales. We addressed this issue in a pond-breeding amphibian, the great crested newt Triturus cristatus. We used capture-recapture data collected through long-term surveys in five T. cristatus populations in Western Europe. In all populations-and subpopulations within metapopulations-population size, annual survival and recruitment fluctuated over time. Likewise, there was considerable variation in these demographic rates between populations and within metapopulations. These fluctuations and variations appear to be context-dependent and more related to site-specific characteristics than local or regional climatic drivers. We found a low level of demographic synchrony at both local and regional levels. Weather has weak and spatially variable effects on survival, recruitment and population growth rate. In contrast, density dependence was a common phenomenon (at least for population growth) in almost all populations and subpopulations. Our findings support the idea that the Moran effect is low in species where the population dynamics more closely depends on local factors (e.g. population density and habitat characteristics) than on large-scale environmental fluctuation (e.g. regional climatic variation). Such responses may have far-reaching consequences for the long-term viability of spatially structured populations and their ability to respond to large-scale climatic anomalies.

Hidden similarities in the dynamics of a weakly synchronous marine metapopulation

Populations of many marine species are only weakly synchronous, despite coupling through larval dispersal and exposure to synchronous environmental drivers. Although this is often attributed to observation noise, factors including local environmental differences, spatially variable dynamics, and chaos might also reduce or eliminate metapopulation synchrony. To differentiate spatially variable dynamics from similar dynamics driven by spatially variable environments, we applied hierarchical delay embedding. A unique output of this approach, the "dynamic correlation," quantifies similarity in intrinsic dynamics of populations, independently of whether their abundance is correlated through time. We applied these methods to 17 populations of blue crab (Callinectes sapidus) along the US Atlantic coast and found that their intrinsic dynamics were broadly similar despite largely independent fluctuations in abundance. The weight of evidence suggests that the latitudinal gradient in temperature, filtered through a unimodal response curve, is sufficient to decouple crab populations. As unimodal thermal performance is ubiquitous in ectotherms, we suggest that this may be a general explanation for the weak synchrony observed at large distances in many marine species, although additional studies are needed to test this hypothesis.

Temporal variation in spatial genetic structure during population outbreaks: Distinguishing among different potential drivers of spatial synchrony

Spatial synchrony is a common characteristic of spatio-temporal population dynamics across many taxa. While it is known that both dispersal and spatially autocorrelated environmental variation (i.e., the Moran effect) can synchronize populations, the relative contributions of each, and how they interact, are generally unknown. Distinguishing these mechanisms and their effects on synchrony can help us to better understand spatial population dynamics, design conservation and management strategies, and predict climate change impacts. Population genetic data can be used to tease apart these two processes as the spatio-temporal genetic patterns they create are expected to be different. A challenge, however, is that genetic data are often collected at a single point in time, which may introduce context-specific bias. Spatio-temporal sampling strategies can be used to reduce bias and to improve our characterization of the drivers of spatial synchrony. Using spatio-temporal analyses of genotypic data, our objective was to identify the relative support for these two mechanisms to the spatial synchrony in population dynamics of the irruptive forest insect pest, the spruce budworm (Choristoneura fumiferana), in Quebec (Canada). AMOVA, cluster analysis, isolation by distance, and sPCA were used to characterize spatio-temporal genomic variation using 1,370 SBW larvae sampled over four years (2012-2015) and genotyped at 3,562 SNP loci. We found evidence of overall weak spatial genetic structure that decreased from 2012 to 2015 and a genetic diversity homogenization among the sites. We also found genetic evidence of a long-distance dispersal event over >140 km. These results indicate that dispersal is the key mechanism involved in driving population synchrony of the outbreak. Early intervention management strategies that aim to control source populations have the potential to be effective through limiting dispersal. However, the timing of such interventions relative to outbreak progression is likely to influence their probability of success.

Long-range dispersal promotes species persistence in climate extremes

Anthropogenic global warming in this century can act as a leading factor for large scale species extinctions in the near future. Species, in order to survive, need to develop dispersal strategies depending upon their environmental niche. Based on empirical evidence only a few previous studies have addressed how dispersal can evolve with changing temperature. However, for the analytical tractability, there is a need to develop an explicit model to ask how the temperature-dependent dispersal alters ecological dynamics. We investigate the persistence of species in a spatial ecological model, where dispersal is considered as a function of temperature. Spatial persistence is of major concern and dispersal is reasonably an important factor for extinction risk in the context of promoting synchrony. Our study yields how the temperature influences species decision of dispersal, resulting in either short-range or long-range dispersal. We examine synchronous or asynchronous behavior of species under their thermal dependence of dispersal. Moreover, we also analyze the transients to study the collective behavior of species away from their final or asymptotic dynamics. One of the key findings is at the most unfavorable environmental conditions long-range dispersal works out as the driving force for the persistence of species.

Density Dependence of Egg Recruitment and Moth Dispersal in Spruce Budworms

Egg recruitment quantifies the relative importance of realized fecundity and migration rates in the population dynamics of highly mobile insects. We develop here a formal context upon which to base the measurement and interpretation of egg recruitment in population dynamics of eastern and western spruce budworms, two geographically separated species that share a very similar ecology. Under most circumstances, per capita egg recruitment rates in these budworms are higher in low-density populations and lower in high-density populations, relative to the regional mean: Low-density populations are nearly always migration sinks for gravid moths, and dense populations nearly always sources. The slope of this relationship, measured on a log scale, is negatively correlated with migration rate, and ranges between 0 and −1. The steeper the slope, the more marked net migration. Using our western spruce budworm observations, we found strong evidence of density-dependent emigration in budworms, so migration is not simply a random perturbation in the lagged, density-dependent stochastic process leading to budworm outbreaks. It is itself statistically and biologically density-dependent. Therefore, moth migration is a synchronizing factor and a spread mechanism that is essential to understanding the development and expansion of spruce budworm outbreaks at regional scales in the boreal forests of North America.

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.

Local adaptation, dispersal evolution, and the spatial eco-evolutionary dynamics of invasion

Local adaptation and dispersal evolution are key evolutionary processes shaping the invasion dynamics of populations colonizing new environments. Yet their interaction is largely unresolved. Using a single-species population model along a one-dimensional environmental gradient, we show how local competition and dispersal jointly shape the eco-evolutionary dynamics and speed of invasion. From a focal introduction site, the generic pattern predicted by our model features a temporal transition from wave-like to pulsed invasion. Each regime is driven primarily by local adaptation, while the transition is caused by eco-evolutionary feedbacks mediated by dispersal. The interaction range and cost of dispersal arise as key factors of the duration and speed of each phase. Our results demonstrate that spatial eco-evolutionary feedbacks along environmental gradients can drive strong temporal variation in the rate and structure of population spread, and must be considered to better understand and forecast invasion rates and range dynamics.

Spatial synchrony in sub-arctic geometrid moth outbreaks reflects dispersal in larval and adult life cycle stages

Spatial synchrony in population dynamics can be caused by dispersal or spatially correlated variation in environmental factors like weather (Moran effect). Distinguishing between these mechanisms is challenging for natural populations, and the study of dispersal-induced synchrony in particular has been dominated by theoretical modelling and laboratory experiments. The goal of the present study was to evaluate the evidence for dispersal as a cause of meso-scale (distances of tens of kilometres) spatial synchrony in natural populations of the two cyclic geometrid moths Epirrita autumnata and Operophtera brumata in sub-arctic mountain birch forest in northern Norway. To infer the role of dispersal in geometrid synchrony, we applied three complementary approaches, namely estimating the effect of design-based dispersal barriers (open sea) on synchrony, comparing the strength of synchrony between E. autumnata (winged adults) and the less dispersive O. brumata (wingless adult females), and relating the directionality (anisotropy) of synchrony to the predominant wind directions during spring, when geometrid larvae engage in windborne dispersal (ballooning). The estimated effect of dispersal barriers on synchrony was almost three times stronger for the less dispersive O. brumata than E. autumnata. Inter-site synchrony was also weakest for O. brumata at all spatial lags. Both observations argue for adult dispersal as an important synchronizing mechanism at the spatial scales considered. Further, synchrony in both moth species showed distinct anisotropy and was most spatially extensive parallel to the east-west axis, coinciding closely to the overall dominant wind direction. This argues for a synchronizing effect of windborne larval dispersal. Congruent with most extensive dispersal along the east-west axis, E. autumnata also showed evidence for a travelling wave moving southwards at a speed of 50-80 km/year. Our results suggest that dispersal processes can leave clear signatures in both the strength and directionality of synchrony in field populations, and highlight wind-driven dispersal as promising avenue for further research on spatial synchrony in natural insect populations.

Proximate determinants of Taylor's law slopes

Taylor's law (TL), a commonly observed and applied pattern in ecology, describes variances of population densities as related to mean densities via log(variance) = log(a) + b*log(mean). Variations among datasets in the slope, b, have been associated with multiple factors of central importance in ecology, including strength of competitive interactions and demographic rates. But these associations are not transparent, and the relative importance of these and other factors for TL slope variation is poorly studied. TL is thus a ubiquitously used indicator in ecology, the understanding of which is still opaque. The goal of this study was to provide tools to help fill this gap in understanding by providing proximate determinants of TL slopes, statistical quantities that are correlated to TL slopes but are simpler than the slope itself and are more readily linked to ecological factors. Using numeric simulations and 82 multi-decadal population datasets, we here propose, test and apply two proximate statistical determinants of TL slopes which we argue can become key tools for understanding the nature and ecological causes of TL slope variation. We find that measures based on population skewness, coefficient of variation and synchrony are effective proximate determinants. We demonstrate their potential for application by using them to help explain covariation in slopes of spatial and temporal TL (two common types of TL). This study provides tools for understanding TL, and demonstrates their usefulness.

Occasional long distance dispersal increases spatial synchrony of population cycles

Spatially separated populations of the same species often exhibit correlated fluctuations in abundance, a phenomenon known as spatial synchrony. Dispersal can generate spatial synchrony. In nature, most individuals disperse short distances with a minority dispersing long distances. The effect of occasional long distance dispersal on synchrony is untested, and theoretical predictions are contradictory. Occasional long distance dispersal might either increase both overall synchrony and the spatial scale of synchrony, or reduce them. We conducted a protist microcosm experiment to test whether occasional long distance dispersal increases or decreases overall synchrony and the spatial scale of synchrony. We assembled replicate 15-patch ring metapopulations of the protist predator Euplotes patella and its protist prey Tetrahymena pyriformis. All metapopulations experienced the same dispersal rate, but differed in dispersal distance. Some metapopulations experienced strictly short distance (nearest neighbour) dispersal, others experienced a mixture of short- and long distance dispersal. Occasional long distance dispersal increased overall spatial synchrony and the spatial scale of synchrony for both prey and predators, though the effects were not statistically significant for predators. As predicted by theory, dispersal generated spatial synchrony by entraining the phases of the predator-prey cycles in different patches, a phenomenon known as phase locking. Our results are consistent with theoretical models predicting that occasional long distance dispersal increases spatial synchrony. However, our results also illustrate that the spatial scale of synchrony need not match the spatial scale of the processes generating synchrony. Even strictly short distance dispersal maintained high spatial synchrony for many generations at spatial scales much longer than the dispersal distance, thanks to phase locking.

Metapopulation dynamics in a changing climate: Increasing spatial synchrony in weather conditions drives metapopulation synchrony of a butterfly inhabiting a fragmented landscape

Habitat fragmentation and climate change are both prominent manifestations of global change, but there is little knowledge on the specific mechanisms of how climate change may modify the effects of habitat fragmentation, for example, by altering dynamics of spatially structured populations. The long-term viability of metapopulations is dependent on independent dynamics of local populations, because it mitigates fluctuations in the size of the metapopulation as a whole. Metapopulation viability will be compromised if climate change increases spatial synchrony in weather conditions associated with population growth rates. We studied a recently reported increase in metapopulation synchrony of the Glanville fritillary butterfly (Melitaea cinxia) in the Finnish archipelago, to see if it could be explained by an increase in synchrony of weather conditions. For this, we used 23 years of butterfly survey data together with monthly weather records for the same period. We first examined the associations between population growth rates within different regions of the metapopulation and weather conditions during different life-history stages of the butterfly. We then examined the association between the trends in the synchrony of the weather conditions and the synchrony of the butterfly metapopulation dynamics. We found that precipitation from spring to late summer are associated with the M. cinxia per capita growth rate, with early summer conditions being most important. We further found that the increase in metapopulation synchrony is paralleled by an increase in the synchrony of weather conditions. Alternative explanations for spatial synchrony, such as increased dispersal or trophic interactions with a specialist parasitoid, did not show paralleled trends and are not supported. The climate driven increase in M. cinxia metapopulation synchrony suggests that climate change can increase extinction risk of spatially structured populations living in fragmented landscapes by altering their dynamics.

Do natural enemies explain fluctuations in low-density spruce budworm populations?

Understanding the causal pathways through which forest insect outbreaks are triggered is important for resource managers. However, detailed population dynamics studies are hard to conduct in low-density, pre-outbreak populations because the insects are difficult to sample in sufficient numbers. Using laboratory-raised larvae installed in the field across a 1,000 km east-west gradient in Québec (Canada) over an 11-yr period, we examined if parasitism and predation were likely to explain fluctuations in low-density spruce budworm (Choristoneura fumiferana; SBW) populations. Parasitism rates by the two main larval parasitoid species, Elachertus cacoeciae and Tranosema rostrale, peaked during different years. This suggests that temporal fluctuations in overall parasitism were partly buffered by compensatory dynamics among parasitoid species. Still, spatial covariance analyses indicate that the residual interannual variation in parasitism rates was substantial and correlated over large distances (up to 700 km). On the other hand, interannual variation in predation rates was not spatially correlated. Piecewise structural equation models indicate that temporal variation in parasitism and predation does not influence temporal variation in wild SBW abundance. Spatially, however, SBWs installed in warmer locations tended to show higher parasitism rates, and these higher rates correlated with lower wild SBW population levels. Overall, the results indicate that large-scale drops in parasitism occur and could potentially contribute to SBW population increases. However, during the period covered by this study, other factors such as direct effects of weather on SBW larval development or indirect effects through host tree physiology or phenology were more likely to explain large-scale variation in wild SBW populations.

Synchrony, Weather, and Cycles in Southern Pine Beetle (Coleoptera: Curculionidae)

Spatial synchrony and cycles are common features of forest insect pests, but are often studied as separate phenomenon. Using time series of timber damage caused by Dendroctonus frontalis Zimmermann (Coleoptera: Curculionidae) (southern pine beetle) in 10 states within the southern United States, this study examines synchrony in D. frontalis abundance, the synchronizing effects of temperature extremes, and the evidence for shared cycles among state populations. Cross-correlation and cluster analyses are used to quantify synchrony across a range of geographic distances and to identify groups of states with synchronous dynamics. Similar techniques are used to quantify spatial synchrony in temperature extremes and to examine their relationship to D. frontalis fluctuations. Cross-wavelet analysis is then used to examine pairs of time series for shared cycles. These analyses suggest there is substantial synchrony among states in D. frontalis fluctuations, and there are regional groups of states with similar dynamics. Synchrony in D. frontalis fluctuations also appears related to spatial synchrony in summer and winter temperature extremes. The cross-wavelet results suggest that D. frontalis dynamics may differ among regions and are not stationary. Significant oscillations were present in some states over certain time intervals, suggesting an endogenous feedback mechanism. Management of D. frontalis outbreaks could potentially benefit from a multistate regional approach because populations are synchronous on this level. Extreme summer temperatures are likely to become the most important synchronizing agent due to climate change.

Catching ghosts with a coarse net: use and abuse of spatial sampling data in detecting synchronization

Synchronization of population dynamics in different habitats is a frequently observed phenomenon. A common mathematical tool to reveal synchronization is the (cross)correlation coefficient between time courses of values of the population size of a given species where the population size is evaluated from spatial sampling data. The corresponding sampling net or grid is often coarse, i.e. it does not resolve all details of the spatial configuration, and the evaluation error-i.e. the difference between the true value of the population size and its estimated value-can be considerable. We show that this estimation error can make the value of the correlation coefficient very inaccurate or even irrelevant. We consider several population models to show that the value of the correlation coefficient calculated on a coarse sampling grid rarely exceeds 0.5, even if the true value is close to 1, so that the synchronization is effectively lost. We also observe 'ghost synchronization' when the correlation coefficient calculated on a coarse sampling grid is close to 1 but in reality the dynamics are not correlated. Finally, we suggest a simple test to check the sampling grid coarseness and hence to distinguish between the true and artifactual values of the correlation coefficient.

Drivers of synchrony of acorn production in the valley oak (Quercus lobata) at two spatial scales

We investigated spatial synchrony of acorn production by valley oaks (Quercus lobata) among individual trees at the within-population, local level and at the among-population, statewide level spanning the geographic range of the species. At the local level, the main drivers of spatial synchrony were water availability and flowering phenology of individual trees, while proximity, temperature differences between trees, and genetic similarity failed to explain a significant proportion of variance in spatial synchrony. At the statewide level, annual rainfall was the primary driver, while proximity was significant by itself but not when controlling for rainfall; genetic similarity was again not significant. These results support the hypothesis that environmental factors, the Moran effect, are key drivers of spatial synchrony in acorn production at both small and large geographic scales. The specific environmental factors differed depending on the geographic scale, but were in both cases related to water availability. In addition, flowering phenology, potentially affecting either density-independent pollination failure (the pollination Moran effect) or density-dependent pollination efficiency (pollen coupling), plays a key role in driving spatial synchrony at the local geographic scale.

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

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

Novel insights on population and range edge dynamics using an unparalleled spatiotemporal record of species invasion

Quantifying the complex spatial dynamics taking place at range edges is critical for understanding future distributions of species, yet very few systems have sufficient data or the spatial resolution to empirically test these dynamics. This paper reviews how data from a large-scale pest management programme have provided important contributions to the fields of population dynamics and invasion biology. The invasion of gypsy moth (Lymantria dispar) is well-documented from its introduction near Boston, Massachusetts USA in 1869 to its current extent of over 900,000 km² in Eastern North America. Over the past two decades, the USDA Forest Service Slow the Spread (STS) programme for managing the future spread of gypsy moth has produced unrivalled spatiotemporal data across the invasion front. The STS programme annually deploys a grid of 60,000-100,000 pheromone-baited traps, currently extending from Minnesota to North Carolina. The data from this programme have provided the foundation for investigations of complex population dynamics and the ability to examine ecological hypotheses previously untestable outside of theoretical venues, particularly regarding invasive spread and Allee effects. This system provides empirical data on the importance of long-distance dispersal and time-lags on population establishment and spatial spread. Studies showing high rates of spatiotemporal variation of the range edge, from rapid spread to border stasis and even retraction, highlight future opportunities to test mechanisms that influence both invasive and native species ranges. The STS trap data have also created a unique opportunity to study low-density population dynamics and quantify Allee effects with empirical data. Notable contributions include evidence for spatiotemporal variation in Allee effects, demonstrating empirical links between Allee effects and spatial spread, and testing mechanisms of population persistence and growth rates at range edges. There remain several outstanding questions in spatial ecology and population biology that can be tested within this system, such as the scaling of local ecological processes to large-scale dynamics across landscapes. The gypsy moth is an ideal model of how important ecological questions can be answered by thinking more broadly about monitoring data.

The geography of spatial synchrony

Spatial synchrony, defined as correlated temporal fluctuations among populations, is a fundamental feature of population dynamics, but many aspects of synchrony remain poorly understood. Few studies have examined detailed geographical patterns of synchrony; instead most focus on how synchrony declines with increasing linear distance between locations, making the simplifying assumption that distance decay is isotropic. By synthesising and extending prior work, we show how geography of synchrony, a term which we use to refer to detailed spatial variation in patterns of synchrony, can be leveraged to understand ecological processes including identification of drivers of synchrony, a long-standing challenge. We focus on three main objectives: (1) showing conceptually and theoretically four mechanisms that can generate geographies of synchrony; (2) documenting complex and pronounced geographies of synchrony in two important study systems; and (3) demonstrating a variety of methods capable of revealing the geography of synchrony and, through it, underlying organism ecology. For example, we introduce a new type of network, the synchrony network, the structure of which provides ecological insight. By documenting the importance of geographies of synchrony, advancing conceptual frameworks, and demonstrating powerful methods, we aim to help elevate the geography of synchrony into a mainstream area of study and application.

Density dependence in demography and dispersal generates fluctuating invasion speeds

Density dependence plays an important role in population regulation and is known to generate temporal fluctuations in population density. However, the ways in which density dependence affects spatial population processes, such as species invasions, are less understood. Although classical ecological theory suggests that invasions should advance at a constant speed, empirical work is illuminating the highly variable nature of biological invasions, which often exhibit nonconstant spreading speeds, even in simple, controlled settings. Here, we explore endogenous density dependence as a mechanism for inducing variability in biological invasions with a set of population models that incorporate density dependence in demographic and dispersal parameters. We show that density dependence in demography at low population densities-i.e., an Allee effect-combined with spatiotemporal variability in population density behind the invasion front can produce fluctuations in spreading speed. The density fluctuations behind the front can arise from either overcompensatory population growth or density-dependent dispersal, both of which are common in nature. Our results show that simple rules can generate complex spread dynamics and highlight a source of variability in biological invasions that may aid in ecological forecasting.

Quantifying dispersal of a non-aggressive saprophytic bark beetle

Long distance dispersal to locate suitable breeding sites is recognized as a key trait influencing the population dynamics and distribution of bark beetles and other saprophytic insects. While dispersal behavior has been studied for a range of aggressive 'tree killing' bark beetles, few have considered the dispersal behaviour of non-aggressive saprophytic bark beetles that utilize kairomones (host volatiles). We present the results of a mark-recapture experiment that examined adult dispersal patterns of the saprophytic bark beetle Hylurgus ligniperda. Releases took place in summer and autumn 2014, in a clearcut pine forest in the central North Island, New Zealand. Both flight-experienced and flight-naïve adults were marked and released in the center of a circular trap grid that extended to 960 m with 170 or 200 panel traps baited with a kairomone blend of alpha-pinene and ethanol. Of the 18,464 released H. ligniperda, 9,209 (49.9%) of the beetles flew, and 96 (1.04%) of the beetles that flew were recaptured. Individuals were recaptured at all distances. The recapture of flight-experienced beetles declined with dispersal distance, and a diffusion model showed heterogeneous dispersal tendencies within the population. Our best model estimated that 46% of flight-experienced beetles disperse > 1 km, and 1.6% > 5 km. Conversely, no declining pattern was shown in the recapture of flight-naïve beetles, suggesting that emerging H. ligniperda may require a period of flight to initiate chemotropic orientation behavior and subsequent attraction to traps. We discuss the implications of these findings for the management of phytosanitary risks. For instance, combining landscape knowledge of source populations with dispersal processes facilitates estimation of pest pressure at economically sensitive areas such as harvest and timber storage sites. Quantitative dispersal estimates also inform pest risk assessments by predicting spread rates for H. ligniperda that has proven establishment capabilities in other countries.

Postglacial recolonization shaped the genetic diversity of the winter moth (Operophtera brumata) in Europe

Changes in climate conditions, particularly during the Quaternary climatic oscillations, have long been recognized to be important for shaping patterns of species diversity. For species residing in the western Palearctic, two commonly observed genetic patterns resulting from these cycles are as follows: (1) that the numbers and distributions of genetic lineages correspond with the use of geographically distinct glacial refugia and (2) that southern populations are generally more diverse than northern populations (the “southern richness, northern purity” paradigm). To determine whether these patterns hold true for the widespread pest species the winter moth (Operophtera brumata), we genotyped 699 individual winter moths collected from 15 Eurasian countries with 24 polymorphic microsatellite loci. We find strong evidence for the presence of two major genetic clusters that diverged ~18 to ~22 ka, with evidence that secondary contact (i.e., hybridization) resumed ~ 5 ka along a well‐established hybrid zone in Central Europe. This pattern supports the hypothesis that contemporary populations descend from populations that resided in distinct glacial refugia. However, unlike many previous studies of postglacial recolonization, we found no evidence for the “southern richness, northern purity” paradigm. We also find evidence for ongoing gene flow between populations in adjacent Eurasian countries, suggesting that long‐distance dispersal plays an important part in shaping winter moth genetic diversity. In addition, we find that this gene flow is predominantly in a west‐to‐east direction, suggesting that recently debated reports of cyclical outbreaks of winter moth spreading from east to west across Europe are not the result of dispersal.

Three Centuries of Synchronous Forest Defoliator Outbreaks in Western North America

Insect outbreaks often occur synchronously across large spatial scales, but the long-term temporal stability of the phenomenon and the mechanisms behind it are not well understood. In this study, I use a widespread lepidopteran defoliator native to western North America—the western spruce budworm—as a case study to explore patterns of and potential causes for synchronous population fluctuations. Analyses of synchrony are typically severely limited by the short historical records available for many species. To overcome this limitation, I compiled multi-century dendrochronological reconstructions of western spruce budworm outbreaks from across much of the species’ range. This allowed me to analyze synchrony at a sub-continental spatial scale over the last three centuries. I found statistically significant synchrony among regional outbreak records up to 2,000 km apart and identified numerous outbreak periods that occurred synchronously across much of the species’ range. I quantified spatial and temporal associations between climate and synchronous outbreak periods using paleoclimate reconstructions. The spatial patterns of outbreak histories and climate records were remarkably similar, with higher similarity in outbreak histories apparent between regions with more similar climate conditions. Synchronous outbreaks typically occurred during periods of average or above average moisture availability preceded by periods of low moisture availability. My results suggest that climatic variability has played a key role in synchronizing western spruce budworm population fluctuations in disjunct forests across western North America for at least the last three centuries. Widespread synchrony appears to be a natural part of this species’ population dynamics, though synchronous outbreaks have occurred more frequently during the 20^th century than during prior centuries. This study uses a novel combination of statistical methods and dendrochronological data to provide analyses of this species’ population dynamics with an unprecedented combination of spatial extent and temporal depth.

Climate change-related regime shifts have altered spatial synchrony of plankton dynamics in the North Sea

During the 1980s, the North Sea plankton community underwent a well-documented ecosystem regime shift, including both spatial changes (northward species range shifts) and temporal changes (increases in the total abundances of warmer water species). This regime shift has been attributed to climate change. Plankton provide a link between climate and higher trophic-level organisms, which can forage on large spatial and temporal scales. It is therefore important to understand not only whether climate change affects purely spatial or temporal aspects of plankton dynamics, but also whether it affects spatiotemporal aspects such as metapopulation synchrony. If plankton synchrony is altered, higher trophic-level feeding patterns may be modified. A second motivation for investigating changes in synchrony is that the possibility of such alterations has been examined for few organisms, in spite of the fact that synchrony is ubiquitous and of major importance in ecology. This study uses correlation coefficients and spectral analysis to investigate whether synchrony changed between the periods 1959-1980 and 1989-2010. Twenty-three plankton taxa, sea surface temperature (SST), and wind speed were examined. Results revealed that synchrony in SST and plankton was altered. Changes were idiosyncratic, and were not explained by changes in abundance. Changes in the synchrony of Calanus helgolandicus and Para-pseudocalanus spp appeared to be driven by changes in SST synchrony. This study is one of few to document alterations of synchrony and climate-change impacts on synchrony. We discuss why climate-change impacts on synchrony may well be more common and consequential than previously recognized.

Model-specification uncertainty in future forest pest outbreak

Climate change will modify forest pest outbreak characteristics, although there are disagreements regarding the specifics of these changes. A large part of this variability may be attributed to model specifications. As a case study, we developed a consensus model predicting spruce budworm (SBW, Choristoneura fumiferana [Clem.]) outbreak duration using two different predictor data sets and six different correlative methods. The model was used to project outbreak duration and the uncertainty associated with using different data sets and correlative methods (=model-specification uncertainty) for 2011-2040, 2041-2070 and 2071-2100, according to three forcing scenarios (RCP 2.6, RCP 4.5 and RCP 8.5). The consensus model showed very high explanatory power and low bias. The model projected a more important northward shift and decrease in outbreak duration under the RCP 8.5 scenario. However, variation in single-model projections increases with time, making future projections highly uncertain. Notably, the magnitude of the shifts in northward expansion, overall outbreak duration and the patterns of outbreaks duration at the southern edge were highly variable according to the predictor data set and correlative method used. We also demonstrated that variation in forcing scenarios contributed only slightly to the uncertainty of model projections compared with the two sources of model-specification uncertainty. Our approach helped to quantify model-specification uncertainty in future forest pest outbreak characteristics. It may contribute to sounder decision-making by acknowledging the limits of the projections and help to identify areas where model-specification uncertainty is high. As such, we further stress that this uncertainty should be strongly considered when making forest management plans, notably by adopting adaptive management strategies so as to reduce future risks.

Dissecting geographic variation in population synchrony using the common vole in central Europe as a test bed

Spatial synchrony of population fluctuations is ubiquitous in nature. Theoretical models suggest that correlated environmental stochasticity, dispersal, and trophic interactions are important promoters of synchrony in nature to leave characteristic signatures of distance‐dependent decays in synchrony. Recent refinements of this theory have clarified how distance‐decay curves may steepen if local dynamics are governed by different density‐dependent feedbacks and how synchrony should vary regionally if the importance and correlation of environmental stochasticity is location‐specific. We analysed spatiotemporal data for the common vole, Microtus arvalis from 49 districts in the Czech Republic to examine the pattern of population synchrony between 2000 and 2014. By extending the nonparametric covariation function, we develop a quantitative method that allows a dissection of the effects of distance and additional variables such as altitude on synchrony. To examine the pattern of local synchrony, we apply the noncentered local‐indicators of spatial association (ncLISA) which highlights areas with different degrees of synchrony than expected by the region‐wide average. Additionally, in order to understand the obtained pattern of local spatial correlations, we have regressed LISA results against the proportion of forest in each district. The common vole abundances fluctuated strongly and exhibited synchronous dynamics with the typical tendency for a decline of synchrony with increasing distance but, not with altitude. The correlation between the neighbor districts decreases as the proportion of forest increases. Forested areas are suboptimum habitats and are strongly avoided by common voles. The investigation of spatiotemporal dynamics in animal populations is a key issue in ecology. Although the majority of studies are focused on testing hypotheses about which mechanisms are involved in shaping this dynamics it is crucial to understand the sources of variation involved in order to understand the underlying processes.