Spatial dynamics of Microtus vole populations in continuous and fragmented agricultural landscapes

Small mammal populations often exhibit large-scale spatial synchrony, which is purportedly caused by stochastic weather-related environmental perturbations, predation or dispersal. To elucidate the relative synchronizing effects of environmental perturbations from those of dispersal movements of small mammalian prey or their predators, we investigated the spatial dynamics of Microtus vole populations in two differently structured landscapes which experience similar patterns of weather and climatic conditions. Vole and predator abundances were monitored for three years on 28 agricultural field sites arranged into two 120-km-long transect lines in western Finland. Sites on one transect were interconnected by continuous agricultural farmland (continuous landscape), while sites on the other were isolated from one another to a varying degree by mainly forests (fragmented landscape). Vole populations exhibited large-scale (>120 km) spatial synchrony in fluctuations, which did not differ in degree between the landscapes or decline with increasing distance between trapping sites. However, spatial variation in vole population growth rates was higher in the fragmented than in the continuous landscape. Although vole-eating predators were more numerous in the continuous agricultural landscape than in the fragmented, our results suggest that predators do not exert a great influence on the degree of spatial synchrony of vole population fluctuations, but they may contribute to bringing out-of-phase prey patches towards a regional density level. The spatial dynamics of vole populations were similar in both fragmented and continuous landscapes despite inter-landscape differences in both predator abundance and possibilities of vole dispersal. This implies that the primary source of synchronization lies in a common weather-related environment.     

Population synchrony and stability in environmentally forced metacommunities

A general prediction from simple metapopulation models is that spatially synchronized forcing can spatially synchronize population dynamics and destabilize metapopulations. In contrast, spatially asynchronous forcing is predicted to decrease population synchrony and promote temporal stability and population persistence, especially in the presence of dispersal. Only recently have studies begun to experimentally address these predictions. Moreover, few studies have experimentally examined how such processes operate in the context of competition communities. Stabilizing processes may continue to operate when placed within a metacommunity context with multiple competing consumers but only at low to intermediate levels of dispersal. High dispersal rates can reverse these predictions and lead to destabilization. We tested this under controlled conditions using an experimental aquatic system composed of three competing species of zooplankton. Metacommunities experienced different levels of dispersal and environmental forcing in the form of spatially synchronous or asynchronous pH perturbations. We found support that dispersal can have contrasting effects on population stability depending on the degree to which population dynamics were synchronized in space. Dispersal under synchronous forcing or no forcing had either neutral of positive effects on spatial population synchrony of all three zooplankton species. In these treatments, dispersal reduced population stability at the local and metapopulation levels for two of three species. In contrast, asynchronously varying environments reduced population synchrony relative to unforced systems, regardless of dispersal level. In these treatments, dispersal enhanced temporal stability and persistence of populations not by reducing population synchrony but by enhancing population minima and spatial averaging of abundances. High dispersal rates under asynchronous forcing reduced the abundance of one species, consistent with increasing regional competition and general metacommunity theory. However, no effects on its stability or persistence were observed. Our work highlights the context‐dependent effects of dispersal on population dynamics in varying environments.     

Geographic variation in density-dependent dynamics impacts the synchronizing effect of dispersal and regional stochasticity

Explanations for the ubiquitous presence of spatially synchronous population dynamics have assumed that density-dependent processes governing the dynamics of local populations are identical among disjunct populations, and low levels of dispersal or small amounts of regionalized stochasticity (Moran effect) can act to synchronize populations. In this study we used historical spatially referenced data on gypsy moth (Lymantria dispar) outbreaks to document that density-dependent processes can vary substantially across geographical landscapes. This variation may be due in part to geographical variation in habitat (e.g., variation in forest composition). We then used a second-order log-linear stochastic model to explore how inter-population variation in density-dependent processes affects synchronization via either synchronous stochastic forcing or dispersal. We found that geographical variation in direct density-dependence (first order) greatly diminishes synchrony caused by stochasticity but only slightly decreases synchronization via dispersal. Variation in delayed density-dependence (second order) diluted synchrony caused by regional stochasticity to a lesser extent than first-order variation, but it did not have any influence on synchrony caused by dispersal. In general, synchronization caused by dispersal was primarily dependent upon the instability of populations and only weakly, if at all, affected by similarities in density-dependence among populations. We conclude that studies of synchrony should carefully consider both the nature of the synchronizing agents and the pattern of local density-dependent processes, including how these vary geographically.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Effects of environmental synchrony and density-dependent dispersal on temporal and spatial slopes of Taylor's law

Taylor's law (TL) is an empirical rule that describes an approximate relationship between the variance and mean of population density: log₁₀(variance) ≈ log₁₀(a) + b × log₁₀(mean). Population synchrony is another prevailing feature observed in empirical populations. This study investigated the effects of environmental synchrony and density-dependent dispersal on the temporal (b_T) and spatial (b_S) slopes of TL, using an empirical dataset of gray-sided vole populations and simulation analyses based on the second-order autoregressive (AR) model. Eighty-five empirical populations satisfied the temporal and spatial TLs with b_T = 1.943 (±SE 0.143) and b_S = 1.579 (±SE 0.136). The pairwise synchrony of population was 0.377 ± 0.199 (mean ± SD). Most simulated populations that obeyed the AR model satisfied the form of the temporal and spatial TLs without being affected by the environmental synchrony and density-dependent dispersal; however, those simulated slopes were too steep. The incorporation of environmental synchrony resulted in reduced simulated slopes, but those slopes, too, were still unrealistically steep. By incorporating density-dependent dispersal, simulated slopes decreased and fell within a realistic range. However, the simulated populations without environmental synchrony did not exhibit an adequate degree of density synchrony. In simulations that included both environmental synchrony and density-dependent dispersal, 92.7% of the simulated datasets provided realistic values for b_T, b_S and population synchrony. Because the two slopes were more sensitive to the variation of density-dependent dispersal than that of environmental synchrony, density-dependent dispersal may be the key to the determination of b_T and b_S.     

Spatial synchrony through density-independent versus density-dependent dispersal

Many theoretical studies support the notion that strong dispersal fosters spatial synchrony. Nonetheless, the effect of conditional vs. unconditional dispersal has remained a matter of controversy. We scrutinize recent findings on a desynchronizing effect of negative density-dependent dispersal based on spatially explicit simulation models. Keeping net emigration rates equivalent, we compared density-independent and density-dependent dispersal for different types of intraspecific density regulation, ranging from under-compensation to over-compensation. In general, density-independent dispersal possessed a slightly higher synchronizing potential but this effect was very small and sensitive compared to the influence of the type of local density regulation. Notably, consistent outcomes for the comparison of conditional dispersal strategies strongly relied on the control of equivalent emigration rates. We conclude that the strength of dispersal is more important for spatial synchrony than its density dependence. Most important is the mode of intraspecific density regulation.     

Spatial synchrony through density-independent versus density-dependent dispersal

Many theoretical studies support the notion that strong dispersal fosters spatial synchrony. Nonetheless, the effect of conditional vs. unconditional dispersal has remained a matter of controversy. We scrutinize recent findings on a desynchronizing effect of negative density-dependent dispersal based on spatially explicit simulation models. Keeping net emigration rates equivalent, we compared density-independent and density-dependent dispersal for different types of intraspecific density regulation, ranging from under-compensation to over-compensation. In general, density-independent dispersal possessed a slightly higher synchronizing potential but this effect was very small and sensitive compared to the influence of the type of local density regulation. Notably, consistent outcomes for the comparison of conditional dispersal strategies strongly relied on the control of equivalent emigration rates. We conclude that the strength of dispersal is more important for spatial synchrony than its density dependence. Most important is the mode of intraspecific density regulation.     

Dispersal and spatial scale affect synchrony in spatial population dynamics

A large body of theoretical studies has shown that synchrony among populations is critical for the long-term persistence of species in fragmented habitats. Although the effects of dispersal and environmental factors on synchrony have been investigated theoretically, empirical studies of these relationships have been lacking. We explored the interplay between environmental and demographic factors (fecundity, survival, dispersal) on population synchrony for 53 species of birds. We show that the interspecific differences in mean synchrony were determined by global environmental factors whose action was probably mediated by the abundance of each species. After removing the effect of these global factors on synchrony, the residual synchrony was strongly correlated with dispersal distance. The relationship between dispersal and synchrony was stronger for the species nesting in wet habitats than for those nesting in dry habitats. Our results indicate that different factors synchronize bird populations at different spatial scales, thus highlighting the role of scale in understanding spatial population dynamics and extinction risks.     

Spatially autocorrelated disturbances and patterns in population synchrony

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

Testing the interaction between environmental variation and dispersal strategy on population dynamics using a soil mite experimental system

Dispersal can play an important role in both the local and regional dynamics of populations. Empirical studies have shown that the proportion of individuals dispersing is often density dependent, which may have implications for the effect of dispersal on populations. In this study, we manipulate the dispersal strategy of adults within two-patch laboratory populations of soil mites and compare the consequences of fixed (density-independent) and density-dependent dispersal in environments of constant and temporally varying resource availability. Effects of dispersal on population dynamics were dependent on the presence of environmental variation. Both dispersal strategies tended to spatially homogenize the population abundance of adults in a variable environment. However, the effect of environmental variation on mean adult abundance was greater with density-dependent dispersal than with fixed dispersal. Adult dispersal did not affect juvenile or egg abundance. This study demonstrates the potential significance of density-dependent dispersal for population dynamics, but emphasizes the role of the environmental context.     

Effects of patch number and dispersal patterns on population dynamics and synchrony

In this paper, we examine the effects of patch number and different dispersal patterns on dynamics of local populations and on the level of synchrony between them. Local population renewal is governed by the Ricker model and we also consider asymmetrical dispersal as well as the presence of environmental heterogeneity. Our results show that both population dynamics and the level of synchrony differ markedly between two and a larger number of local populations. For two patches different dispersal rules give very versatile dynamics. However, for a larger number of local populations the dynamics are similar irrespective of the dispersal rule. For example, for the parameter values yielding stable or periodic dynamics in a single population, the dynamics do not change when the patches are coupled with dispersal. High intensity of dispersal does not guarantee synchrony between local populations. The level of synchrony depends also on dispersal rule, the number of local populations, and the intrinsic rate of increase. In our study, the effects of density-independent and density-dependent dispersal rules do not show any consistent difference. The results call for caution when drawing general conclusions from models of only two interacting populations and question the applicability of a large number of theoretical papers dealing with two local populations.     

Density-dependent dispersal complicates spatial synchrony in tri-trophic food chains

Spatial synchrony can increase extinction risk and undermines metapopulation persistence. Both dispersal and biotic interactions can strongly affect spatial synchrony. Here, we explore the spatial synchrony of a tri-trophic food chain in two patches connected by density-dependent dispersal, namely the strategies of prey evasion (PE) and predator pursuit (PP). The dynamics of the food chain are depicted by both the Hastings–Powell model and the chemostat model, with synchrony measured by the Pearson correlation coefficient. We use the density-independent dispersal in the system as a baseline for comparison. Results show that the density-independent dispersal of a species in the system can promote its dynamic synchrony. Dispersal of intermediate species in the tri-trophic food chain is the strongest synchronizer. In contrast, the density-dependent PP and PE of intermediate species can desynchronize the system. Highly synchronized dynamics emerged when the basal species has a strong PE strategy or when the top species has a moderate PP strategy. Our results reveal the complex relationship between density-dependent dispersal and spatial synchrony in tri-trophic systems.     

Searching for mechanisms of synchrony in spatially structured gamebird populations

1.  Time series data on five species of gamebird from the Dolomitic Alps were used to examine the relative importance of dispersal and common stochastic events in causing synchrony between spatially structured populations.
2.  Cross-correlation analysis of detrended time series was used to describe the spatial pattern of fluctuations in abundance, while standardized time series were used to describe both fluctuations and the trend in abundance. There were large variations in synchrony both within and between species and only weak negative relationships with distance.
3.  Species in neighbouring habitats were more likely to be in synchrony than species separated by several habitats. Species with similar density-dependent structure were more likely to be in synchrony.
4.  In order to estimate the relative importance of dispersal and environmental stochasticity, we modelled the spatial dynamics of each species using two different approaches. First, we used estimating functions and bootstrapping of time series data to calculate the relative importance of dispersal and stochastic effects for each species. Second, we estimated the intensity of environmental stochasticity from climatic records during the breeding season and then modelled the dispersal rate and dispersal distance for each species. The two models exhibited similar results for rock ptarmigan, black grouse, hazel grouse and rock partridge, while contrasting patterns were observed for capercaillie.
5.  The results suggest that environmental stochasticity plays the dominant role in synchronizing the fluctuations of these galliform species, although there will also be some dispersal between populations.     

Source–sink dynamics: how sinks affect demography of sources

Models of source–sink population dynamics have to make assumptions about whether, and eventually how, demographic parameters in source habitats are dependent on the demography in sink habitats. However, the empirical basis for making such assumptions has been weak. Here we report a study on experimental root vole populations, where estimates of demographic parameters were contrasted between source patches in source–sink (treatment) and source–source systems (control). In the presence of a sink patch (simulated by a pulsed removal of immigrants), source‐patch populations failed to increase over the breeding season, mainly due to a high spatially density‐dependent dispersal rate from source to sink patches. The per capita recruitment rate was almost two times higher in source–sink than in the source–source systems, but this did not compensate for the loss rate due to dispersal from source to sink patches. Sex ratio in the source–sink systems became less female biased, probably as a result of an enhanced frequency of dispersal movements in females. Good knowledge of the degree of density‐and habitat‐dependent dispersal is critical for predicting the dynamics of source–sink populations.     

The Moran effect revisited: spatial population synchrony under global warming

The world is spatially autocorrelated. Both abiotic and biotic properties are more similar among neighboring than distant locations, and their temporal co-fluctuations also decrease with distance. P. A. P. Moran realized the ecological importance of such ‘spatial synchrony’ when he predicted that isolated populations subject to identical log-linear density-dependent processes should have the same correlation in fluctuations of abundance as the correlation in environmental noise. The contribution from correlated weather to synchrony of populations has later been coined the ‘Moran effect’. Here, we investigate the potential role of the Moran effect in large-scale ecological outcomes of global warming. Although difficult to disentangle from dispersal and species interaction effects, there is compelling evidence from across taxa and ecosystems that spatial environmental synchrony causes population synchrony. Given this, and the accelerating number of studies reporting climate change effects on local population dynamics, surprisingly little attention has been paid to the implications of global warming for spatial population synchrony. However, a handful of studies of insects, birds, plants, mammals and marine plankton indicate decadal-scale changes in population synchrony due to trends in environmental synchrony. We combine a literature review with modeling to outline potential pathways for how global warming, through changes in the mean, variability and spatial autocorrelation of weather, can impact population synchrony over time. This is particularly likely under a ‘generalized Moran effect’, i.e. when relaxing Moran's strict assumption of identical log-linear density-dependence, which is highly unrealistic in the wild. Furthermore, climate change can influence spatial population synchrony indirectly, through its effects on dispersal and species interactions. Because changes in population synchrony may cascade through food-webs, we argue that the (generalized) Moran effect is key to understanding and predicting impacts of global warming on large-scale ecological dynamics, with implications for extinctions, conservation and management.     

Landscape effects on temporal and spatial properties of vole population fluctuations

Populations of northern small rodents have previously been observed to fluctuate in spatial synchrony over distances ranging from tens to hundreds of kilometers between sites. It has been suggested that this phenomenon is caused by common environmental perturbations, mobile predators or dispersal movements. However, very little focus has been given to how the physical properties of the geographic area over which synchrony occurs, such as landscape composition and climate, affect spatial population dynamics. This study reports on the spatial and temporal properties of vole population fluctuations in two areas of western Finland: one composed of large interconnected areas of agricultural farmland interspersed by forests and the other highly dominated by forest areas, containing more isolated patches of agricultural land. Furthermore, the more agricultural area exhibits somewhat milder winters with less snow than the forested area. We found the amplitude of vole cycles to be essentially the same in the two areas, suggesting that the relative amount of predation on small rodents by generalist versus specialist predators is similar in both areas. No seasonal differences in the timing of synchronization were observable for Microtus voles, whereas bank vole populations in field habitats appeared to become synchronized primarily during winter. Microtus populations in field habitats exhibited smaller spatial variation and a higher degree of synchrony in the more continuous agricultural landscape than in the forest-dominated landscape. We suggest that this inter-areal difference is due to differences in the degree of inter-patch connectivity, with predators and dispersal acting as the primary synchronizing agents. Bank vole populations in field habitats were more synchronized within the forest-dominated landscape, most likely reflecting the suitability of the inter-patch matrix and the possibility of dispersal. Our study clearly indicates that landscape composition needs to be taken into account when describing the spatial properties of small rodent population dynamics.     

Travelling wave dynamics and self-organization in a spatio-temporally structured population

We analyse spatial population dynamics showing that periodic or period-like chaotic dynamics produce self-organization structures, such as travelling waves. We suggest that self-organized patterns are associated with spatial synchrony patterns that often depend on geographical distance between subpopulations. The population dynamics also show statistical spatial autocorrelation patterns. We contrast our theoretical simulations with empirical data on annual damages in young sapling stands caused by voles. We conclude, on the basis of the periodicity, synchrony, and spatial autocorrelation patterns, and our simulation results, that vole dynamics represent travelling waves in population dynamics. We suggest that because such synchrony patterns are frequently observed in natural populations, spatial self-organization may be more common in population dynamics than reported in the literature.     

New views on how population‐intrinsic and community‐extrinsic processes interact during the vole population cycles

Based on evidence from a series of recent studies linking behaviour to demography in experimental vole populations we propose how intrinsic and extrinsic factors interact through the various phases of the multi‐annual population cycles of voles and lemmings. We hypothesise that population growth in the increase phase of the cycle is enhanced by a high degree of space sharing (sociality) among reproductive females who share resource patches, especially during winter. These social females enjoy a high reproductive output due to good resource conditions, and facilitation provided by communal thermoregulation, breeding and defence of weanlings towards infanticidal conspecifics. We hypothesise on the other hand that the crash phase is initiated and enhanced by predation of adult males that leads to a series of cascading events involving infanticidal behaviour by immigrant males, increased mortality of adult social females, and inversely density‐dependent and/or disturbance‐induced dispersal. These events further enhance predation‐induced mortality and thus a negative feed‐back loop to the rate of the crash. In this framework we may explain how extrinsic factors such as predation and winter resource distribution contribute to transitions between docile and aggressive behaviours, and how this transition is spatially synchronised by inversely density‐dependent dispersal that may act to mediate a rapidly spreading wave throughout the population. We propose that innate differences among rodent species in their propensities for different social organizations also determine their propensity for exhibiting multi‐annual cycles as well as other cycle‐related phenomena such as shape of the population cycles and spatial synchrony. We provide a set of testable predictions for further empirical evaluation.     

The effect of habitat size and predation on the time to extinction of prairie vole populations: simulation studies via SERDYCA

We present a spatially explicit individual-based model of rodent dynamics, customized for the prairie vole species, Microtus ochrogaster. The model strives to represent the complexity of intertwining factors that determine the spatio-temporal dynamics of small rodents. It is based on trophic relationships and incorporates important features such as territorial competition, mating behavior, density-dependent predation and dispersal out of the modeled spatial region. Vegetation growth and vole fecundity are dependent on climatic components. The results of simulations show that the model correctly predicts the overall temporal dynamics of the population density. Time-series analysis shows a very good match between the periods corresponding to the peak population density frequencies predicted by the model and the ones reported in the literature. The model is used to study the relation between persistence, landscape area and predation. We use the notions of average time to extinction (ATE) and persistence frequency to quantify persistence. While the ATE decreases with decrease of area, it is a bell-shaped function of the predation level: increasing for “small” and decreasing for “large” predation levels.     

Population synchrony in small-world networks

Network topography ranges from regular graphs (linkage between nearest neighbours only) via small-world graphs (some random connections between nodes) to completely random graphs. Small-world linkage is seen as a revolutionary architecture for a wide range of social, physical and biological networks, and has been shown to increase synchrony between oscillating subunits. We study small-world topographies in a novel context: dispersal linkage between spatially structured populations across a range of population models. Regular dispersal between population patches interacting with density-dependent renewal provides one ecological explanation for the large-scale synchrony seen in the temporal fluctuations of many species, for example, lynx populations in North America, voles in Fennoscandia and grouse in the UK. Introducing a small-world dispersal kernel leads to a clear reduction in synchrony with both increasing dispersal rate and small-world dispersal probability across a variety of biological scenarios. Synchrony is also reduced when populations are affected by globally correlated noise. We discuss ecological implications of small-world dispersal in the frame of spatial synchrony in population fluctuations.     

Density-dependent dispersal and relative dispersal affect the stability of predator-prey metacommunities

Although density-dependent dispersal and relative dispersal (the difference in dispersal rates between species) have been documented in natural systems, their effects on the stability of metacommunities are poorly understood. Here we investigate the effects of intra- and interspecific density-dependent dispersal on the regional stability in a predator-prey metacommunity model. We show that, when the dynamics of the populations reach equilibrium, the stability of the metacommunity is not affected by density-dependent dispersal. However, the regional stability, measured as the regional variability or the persistence, can be modified by density-dependent dispersal when local populations fluctuate over time. Moreover these effects depend on the relative dispersal of the predator and the prey. Regional stability is modified through changes in spatial synchrony. Interspecific density-dependent dispersal always desynchronizses local dynamics, whereas intraspecific density-dependent dispersal may either synchronize or desynchronize it depending on dispersal rates. Moreover, intra- and interspecific density-dependent dispersal strengthen the top-down control of the prey by the predator at intermediate dispersal rates. As a consequence the regional stability of the metacommunity is increased at intermediate dispersal rates. Our results show that density-dependent dispersal and relative dispersal of species are keys to understanding the response of ecosystems to fragmentation.