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.     

Synchronous dynamics and rates of extinction in spatially structured populations

We explore extinction rates using a spatially arranged set of subpopulations obeying Ricker dynamics. The population system is subjected to dispersal of individuals among the subpopulations as well as to local and global disturbances. We observe a tight positive correlation between global extinction rate and the level of synchrony in dynamics among thesubpopulations. Global disturbances and to a lesser extent, migration, are capable of synchronizing the temporal dynamics of the subpopulations over a rather wide span of the population growth rate r. Local noise decreases synchrony, as does increasing distance among the subpopulations. Synchrony also levels off with increasing r: in the chaotic region, subpopulations almost invariably behave asynchronously. We conclude that it is asynchrony that reduces the probability of global extinctions, not chaos as such: chaos is a special case only. The relationship between global extinction rate, synchronous dynamics and population growth rate is robust to changes in dispersal rates and ranges.     

Canonical functions for dispersal-induced synchrony

Two processes are universally recognized for inducing spatial synchrony in abundance: dispersal and correlated environmental stochasticity. In the present study we seek the expected relationship between synchrony and distance in populations that are synchronized by density-independent dispersal. In the absence of dispersal, synchrony among populations with simple dynamics has been shown to echo the correlation in the environment. We ask what functional form we may expect between synchrony and distance when dispersal is the synchronizing agent. We formulate a continuous-space, continuous-time model that explicitly represents the time evolution of the spatial covariance as a function of spatial distance. Solving this model gives us two simple canonical functions for dispersal-induced covariance in spatially extended populations. If dispersal is rare relative to birth and death, then covariances between nearby points will follow the dispersal distance distribution. At long distances, however, the covariance tails off according to exponential or Bessel functions (depending on whether the population moves in one or two dimensions). If dispersal is common, then the covariances will follow the mixture distribution that is approximately Gaussian around the origin and with an exponential or Bessel tail. The latter mixture results regardless of the original dispersal distance distribution. There are hence two canonical functions for dispersal-induced synchrony     

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

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

Synchronization in ecological systems by weak dispersal coupling with time delay

One of the most salient spatiotemporal patterns in population ecology is the synchronization of fluctuating local populations across vast spatial extent. Synchronization of abundance has been widely observed across a range of spatial scales in relation to the rate of dispersal among discrete populations. However, the dependence of synchrony on patterns of among-patch movement across heterogeneous landscapes has been largely ignored. Here, we consider the duration of movement between two predator–prey communities connected by weak dispersal and its effect on population synchrony. More specifically, we introduce time-delayed dispersal to incorporate the finite transmission time between discrete populations across a continuous landscape. Reducing the system to a phase model using weakly connected network theory, it is found that the time delay is an important factor determining the nature and stability of phase-locked states. Our analysis predicts enhanced convergence to stable synchronous fluctuations in general and a decreased ability of systems to produce in-phase synchronization dynamics in the presence of delayed dispersal. These results introduce delayed dispersal as a tool for understanding the importance of dispersal time across a landscape matrix in affecting metacommunity dynamics. They further highlight the importance of landscape and dispersal patterns for predicting the onset of synchrony between weakly coupled populations.     

The evolution of dispersal distance in spatially-structured populations

Most evolutionary models of dispersal have concentrated on dispersal rate, with emigration being either global or restricted to nearest neighbours. Yet most organisms fall into an intermediate region where most dispersal is local but there is a wide range of dispersal distances. We use an individual-based model with 2500 patches each with identical local dynamics and show that the dispersal distance is under selection pressure. The dispersal distance that evolves is critically dependent on the ecological dynamics. When the cost of dispersal increases linearly with distance, selection is for short-distance dispersal under stable and damped local dynamics but longer distance dispersal is favoured as local dynamics become more complex. For the cases of stable, damped and periodic patch dynamics global patch synchrony occurs even with very short-distance dispersal. Increasing the scale of dispersal for chaotic local dynamics increases the scale of synchrony but global synchrony does not neccesarily occur. We discuss these results in the light of other possible causes of dispersal and argue for the importance of incorporating non-equilibrium population dynamics into evolutionary models of dispersal distance.     

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.     

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.     

Synchrony in brown trout, Salmo trutta, population dynamics: a 'Moran effect' on early-life stages

Synchrony among populations (i.e. spatial covariation in temporal fluctuations of population size or growth rate) is a common feature to many animals. Both large-scale autocorrelated climatic factors (the 'Moran effect') and dispersal between populations are candidates to explain synchrony, although their relative influence is difficult to assess. Only a few investigations have reported patterns of synchrony among freshwater populations, and even fewer directly related these patterns to an environmental variable. In the present study, we analysed the spatio-temporal patterns of fluctuation of 57 brown trout populations widespread across France, each sampled continuously during 5 years. We compared the respective influence of connectivity and stream distance within basins (i.e. that potentially allow a basin-scale dispersal) and environmental factors (hydrological and air temperature variables, available for 37 sites) on the synchrony of brown trout cohort densities (0+, 1+ and adults). A series of Mantel tests revealed that the degree of synchrony was not related to connectivity or stream distance between sites, suggesting no effect of dispersal at the basin-scale. The degree of synchrony among sites for the 0+ fish was significantly related to the degree of hydrological synchrony (based on high flows during the emergence period). For all three age classes, the synchrony in the temperature patterns did not explain synchrony in trout dynamics. Our results allow us to discuss the respective influence of dispersal and climatic factors on the spatio-temporal patterns of trout dynamics at the basin scale.     

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.     

Self-organized dynamics in spatially structured populations

Self-organization and pattern formation represent the emergence of order in temporal and spatial processes. Self-organization in population ecology is gaining attention due to the recent advances concerning temporal fluctuations in the population size of dispersal-linked subunits. We shall report that spatially structured models of population renewal promote the emergence of a complex power law order in spatial population dynamics. We analyse a variety of population models showing that self-organization can be identified as a temporal match in population dynamics among local units, and how the synchrony changes in time. Our theoretical results are concordant with analyses of population data on the Canada lynx.     

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.     

When can dispersal synchronize populations?

While spatial synchrony of oscillating populations has been observed in many ecological systems, the causes of this phenomenon are still not well understood. The most common explanations have been the Moran effect (synchronous external stochastic influences) and the effect of dispersal among populations. Since ecological systems are typically subject to large spatially varying perturbations which destroy synchrony, a plausible mechanism explaining synchrony must produce rapid convergence to synchrony. We analyze the dynamics through time of the synchronizing effects of dispersal and, consequently, determine whether dispersal can be the mechanism which produces synchrony. Specifically, using methods new to ecology, we analyze a two patch predator-prey model, with identical weak dispersal between the patches. We find that a difference in time scales (i.e. one population has dynamics occurring much faster than the other) between the predator and prey species is the most important requirement for fast convergence to synchrony.     

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.     

The genetic underpinnings of population cyclicity: establishing expectations for the genetic anatomy of cycling populations

Despite extensive research into the mechanisms underlying population cyclicity, we have little understanding of the impacts of numerical fluctuations on the genetic variation of cycling populations. Thus, the potential implications of natural and anthropogenically‐driven variation in population cycle dynamics on the diversity and evolutionary potential of cyclic populations is unclear. Here, we use Canada lynx Lynx canadensis matrix population models, set up in a linear stepping‐stone, to generate demographic replicates of biologically realistic cycling populations. Overall, increasing cycle amplitude predictably reduced genetic diversity and increased genetic differentiation, with cyclic effects increased by population synchrony. Modest dispersal rates (1–3% of the population) between high and low amplitude cyclic populations did not diminish these effects suggesting that spatial variation in cyclic amplitude should be reflected in patterns of genetic diversity and differentiation at these rates. At high dispersal rates (6%) groups containing only high amplitude cyclic populations had higher diversity and lower differentiation than those mixed with low amplitude cyclic populations. Negative density‐dependent dispersal did not impact genetic diversity, but did homogenize populations by reducing differentiation and patterns of isolation by distance. Surprisingly, temporal changes in diversity and differentiation throughout a cycle were not always consistent with population size. In particular, negative density‐dependent dispersal simultaneously decreased differences in genetic diversity while increasing differences in genetic differentiation between numerical peaks and nadirs. Combined, our findings suggest demographic changes at fine temporal scales can impact genetic variation of interacting populations and provide testable predictions relating population cyclicty to genetic variation. Further, our results suggest that including realistic demographic and dispersal parameters in population genetic models and using information from temporal changes in genetic variation could help to discern complex demographic scenarios and illuminate population dynamics at fine temporal scales.     

Global patterns of environmental synchrony and the Moran effect

There is considerable debate over the relative importance of dispersal and environmental disturbances (the Moran effect) as causes of spatial synchrony in fluctuations of animal populations. If environmental factors generally exhibit high levels of spatial autocorrelation, they may be playing a more important role in synchronizing animal populations than sometimes recognized. Here I examine this issue by analyzing spatial autocorrelation in annual rainfall and mean annual temperatures from sites throughout the world using the database maintained by the Global Historical Climatology Network. Both annual precipitation and mean annual temperatures exhibit high synchrony declining with distance and are statistically significant over large distance, often on a continental scale. In general, synchrony was slightly higher in annual precipitation at short distances, but greater in mean annual temperatures at long distances. No latitudinal gradient in synchrony of either variable was detected. The high overall synchrony observed in these environmental variables combined with a pattern of decline with distance similar to that observed in many animal populations suggest that the Moran effect can potentially play an important role in driving synchrony in a wide variety of ecological phenomena regardless of scale.     

Synchrony of spatial populations induced by colored environmental noise and dispersal

Spatial synchrony of oscillating populations has been observed in many ecological systems, and its influences and causes have attracted the interest of ecologists. Spatially correlated environmental noises, dispersal, and trophic interactions have been considered as the causes of spatial synchrony. In this study, we develop a spatially structured population model, which is described by coupled-map lattices and incorporates both dispersal and colored environmental noise. A method for generating time series with desired spatial correlation and color is introduced. Then, we use these generated time series to analyze the influence of noise color on synchrony in population dynamics. The noise color refers to the temporal correlation in the time series data of the noise, and is expressed as the degree of (first-order) autocorrelation for autoregressive noise. Patterns of spatial synchrony are considered for stable, periodic and chaotic population dynamics. Numerical simulations verify that environmental noise color has a major influence on the level of synchrony, which depends strongly on how noise is introduced into the model. Furthermore, the influence of noise color also depends on patterns of dispersal between local populations. In addition, the desynchronizing effect of reddened noise is always weaker than that of white noise. From our results, we notice that the role of reddened environmental noise on spatial synchrony should be treated carefully and cautiously, especially for the spatially structured populations linked by dispersal.     

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.     

Patterns and persistence of larval retention and connectivity in a marine fish metapopulation

Connectivity, the demographic linking of local populations through the dispersal of individuals, is one of the most poorly understood processes in population dynamics, yet has profound implications for conservation and harvest strategies. For marine species with pelagic larvae, direct estimation of connectivity remains logistically challenging and has mostly been limited to single snapshots in time. Here, we document seasonal and interannual patterns of larval dispersal in a metapopulation of the coral reef fish Amphiprion polymnus. A 3‐year record of larval trajectories within and among nine discrete local populations from an area of approximately 35 km was established by determining the natal origin of settled juveniles through DNA parentage analysis. We found that spatial patterns of both self‐recruitment and connectivity were remarkably consistent over time, with a low level of self‐recruitment at the scale of individual sites. Connectivity among sites was common and multidirectional in all years and was not significantly influenced by seasonal variability of predominant surface current directions. However, approximately 75% of the sampled juveniles could not be assigned to parents within the study area, indicating high levels of immigrations from sources outside the study area. The data support predictions that the magnitude and temporal stability of larval connectivity decreases significantly with increasing distance between subpopulations, but increases with the size of subpopulations. Given the considerable effort needed to directly measure larval exchange, the consistent patterns suggest snapshot parentage analyses can provide useful dispersal estimates to inform spatial management decisions.     

Using mechanistic models to understand synchrony in forest insect populations: the North American gypsy moth as a case study

In many forest insects, subpopulations fluctuate concurrently across large geographical areas, a phenomenon known as population synchrony. Because of the large spatial scales involved, empirical tests to identify the causes of synchrony are often impractical. Simple models are, therefore, a useful aid to understanding, but data often seem to contradict model predictions. For instance, chaotic population dynamics and limited dispersal are not uncommon among synchronous forest defoliators, yet both make it difficult to achieve synchrony in simple models. To test whether this discrepancy can be explained by more realistic models, we introduced dispersal and spatially correlated stochasticity into a mechanistic population model for the North American gypsy moth Lymantria dispar. The resulting model shows both chaotic dynamics and spatial synchrony, suggesting that chaos and synchrony can be reconciled by the incorporation of realistic dynamics and spatial structure. By relating alterations in model structure to changes in synchrony levels, we show that the synchrony is due to a combination of spatial covariance in environmental stochasticity and the origins of chaos in our multispecies model.