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.     

Environment,but not migration rate,influences extinction risk in experimental metapopulations

Ecological theory suggests that several demographic factors influence metapopulation extinction risk, including synchrony in population size between subpopulations, metapopulation size and the magnitude of fluctuations in population size. Theoretically, each of these is influenced by the rate of migration between subpopulations. Here we report on an experiment where we manipulated migration rate within metapopulations of the freshwater zooplankton Daphnia magna to examine how migration influenced each of these demographic variables, and subsequent effects on metapopulation extinction. In addition, our experimental procedures introduced unplanned but controlled differences between metapopulations in light intensity, enabling us to examine the relative influences of environmental and demographic factors. We found that increasing migration rate increased subpopulation synchrony. We failed to detect effects of migration on population size and fluctuations in population size at the metapopulation or subpopulation level, however. In contrast, light intensity did not influence synchrony, but was positively correlated with population size and negatively correlated with population fluctuation. Finally, synchrony did not influence time to extinction, while population size and the magnitude of fluctuations did. We conclude that environmental factors had a greater influence on extinction risk than demographic factors, and that metapopulation size and fluctuation were more important to extinction risk than metapopulation 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.     

Evolution of dispersal in metapopulations with local density dependence and demographic stochasticity

In this paper, we predict the outcome of dispersal evolution in metapopulations based on the following assumptions: (i) population dynamics within patches are density-regulated by realistic growth functions; (ii) demographic stochasticity resulting from finite population sizes within patches is accounted for; and (iii) the transition of individuals between patches is explicitly modelled by a disperser pool. We show, first, that evolutionarily stable dispersal rates do not necessarily increase with rates for the local extinction of populations due to external disturbances in habitable patches. Second, we describe how demographic stochasticity affects the evolution of dispersal rates: evolutionarily stable dispersal rates remain high even when disturbance-related rates of local extinction are low, and a variety of qualitatively different responses of adapted dispersal rates to varied levels of disturbance become possible. This paper shows, for the first time, that evolution of dispersal rates may give rise to monotonically increasing or decreasing responses, as well as to intermediate maxima or minima.     

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.     

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.     

The effects of small dispersal rates on extinction times in structured metapopulation models

Habitat destruction is a critical factor that affects persistence in several taxa, including Pacific salmon. Salmon are noted for their ability to home to their natal streams for reproduction. Since straying (i.e., spawners reproducing in nonnatal streams) is typically low in salmon, its effects have not been appreciated. In this article, we develop both a general analytical model and a simple simulation model describing structured metapopulations to study how weak connections between subpopulations affect the ability of a species to tolerate habitat destruction and/or declines in habitat quality. Our goals are to develop general principles and to relate these principles to salmon population dynamics. The analytical model describes the dynamics of two density-dependent subpopulations, connected by dispersal, whose growth rates fluctuate in response to environmental and demographic stochasticity. We find that, for moderate levels of environmental variability, small dispersal rates can significantly increase mean extinction times. This effect declines with increasing habitat quality, increasing temporal correlation, and increasing spatial correlation, but it is still significant for realistic parameter values. The simulation model shows there is a threshold rate of dispersal that minimizes extinction probabilities. These results cannot be seen in classical metapopulation models and provide new insights into the rescue effect.     

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.     

Coherence,conservation and patch‐occupancy analysis

Spatial coherence (synchrony) among subpopulations poses a danger to the metacommunity, as it increases the risk of regional extinction. When this effect is significant, the use of inference techniques based on the stochastic patch occupancy model (SPOM) may be inadequate, since SPOMs assume that each habitat patch is either occupied or empty, thereby neglecting the intra‐patch dynamics. Here we suggest a general classification of the dynamics that allows the identification, in a model‐independent manner, of the regimes where coherence effects are strong. We also present a new technique, based on patch occupancy (presence/absence) data, for identifying the role of spatial coherence in the stabilization of a metapopulation. If the chance of a local extinction grows with the connectivity, this implies that spatial synchronization is too strong and that regional‐scale extinction becomes possible. When this scenario occurs, a decrease in the movement of individuals (habitat fragmentation, reduced dispersal rates) has a positive effect on the sustainability of the spatially distributed population. The results of individual based simulations of a spatially structured population are analyzed with SPOM and the regime where the two‐state approximation fails is identified.     

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.     

Nonlinear Effect of Dispersal Rate on Spatial Synchrony of Predator-Prey Cycles

Spatially-separated populations often exhibit positively correlated fluctuations in abundance and other population variables, a phenomenon known as spatial synchrony. Generation and maintenance of synchrony requires forces that rapidly restore synchrony in the face of desynchronizing forces such as demographic and environmental stochasticity. One such force is dispersal, which couples local populations together, thereby synchronizing them. Theory predicts that average spatial synchrony can be a nonlinear function of dispersal rate, but the form of the dispersal rate-synchrony relationship has never been quantified for any system. Theory also predicts that in the presence of demographic and environmental stochasticity, realized levels of synchrony can exhibit high variability around the average, so that ecologically-identical metapopulations might exhibit very different levels of synchrony. We quantified the dispersal rate-synchrony relationship using a model system of protist predator-prey cycles in pairs of laboratory microcosms linked by different rates of dispersal. Paired predator-prey cycles initially were anti-synchronous, and were subject to demographic stochasticity and spatially-uncorrelated temperature fluctuations, challenging the ability of dispersal to rapidly synchronize them. Mean synchrony of prey cycles was a nonlinear, saturating function of dispersal rate. Even extremely low rates of dispersal (<0.4% per prey generation) were capable of rapidly bringing initially anti-synchronous cycles into synchrony. Consistent with theory, ecologically-identical replicates exhibited very different levels of prey synchrony, especially at low to intermediate dispersal rates. Our results suggest that even the very low rates of dispersal observed in many natural systems are sufficient to generate and maintain synchrony of cyclic population dynamics, at least when environments are not too spatially heterogeneous.     

Neighbourhood size,dispersal distance and the complex dynamics of the spatial Ricker model

Amongst the most frequently made assumptions in simple population models are that individuals interact equally with every other individual and that dispersal occurs with equal likelihood to any location. This is especially true for models of a single population (as opposed to a patchy population or metapopulation). For many species of animals and probably for all plant species these assumptions are unlikely to hold true. Here one much-studied population model—the Ricker model—is reformulated such that interactions occur only between individuals located within a certain distance of each other and dispersal distance is finite. Two alternative reformulations are presented. Results demonstrate that both limiting the interaction neighbourhood and reducing dispersal distance tend to stabilise the global population dynamics, although the extent to which this occurs depends upon the reformulation used. Spatial pattern formation is a feature of the simulated population. At lower intrinsic rates of growth (r) these patterns tend to be static, while for higher r, they are dynamic. Both the stabilisation of global dynamics and spatial pattern formation are well-described features of metapopulation models. Here, similar effects are shown to occur on a single contiguous patch of habitat.     

荒漠破碎化生境中长爪沙鼠集合种群野外验证研究

近年来,人类活动和自然干扰,导致内蒙古阿拉善荒漠区生境的破碎化,出现了长爪沙鼠在不同斑块间的不连续分布,每一斑块内可能存在一个局域种群,而集合种群建立的前提条件,是局域种群斑块状分布在离散的栖息地环境中。2002~2012年每年的4~10月,在阿拉善荒漠区禁牧、轮牧、过牧和开垦4种人为不同利用方式形成的生境斑块中,采用标志重捕法对长爪沙鼠（Meriones unguiculatus）种群进行定点监测。通过分析长爪沙鼠种群动态,计算各局域种群的灭绝风险,利用Spearman秩相关系数检验种群动态的空间同步性,同时以种群周转率对长爪沙鼠扩散能力进行评估,以检验阿拉善荒漠区长爪沙鼠种群空间结构是否具有经典集合种群的功能。结果表明：（1） 不同生境斑块可被长爪沙鼠局域种群占据,11年间捕获长爪沙鼠2~7次不等;（2） 长爪沙鼠所有局域种群均具有灭绝风险,在轮牧区和禁牧区灭绝率高达1.000 0,开垦区灭绝率最低,也达到0.333 4,而本研究期间最大局域种群（2008年过牧区,26只/hm²）,在2010年发生了局域灭绝;（3） 不同生境斑块间没有明显的空间隔离而阻碍局域种群的重新建立,长爪沙鼠扩散能力较强,绝大部分月份的种群周转率在50.0%以上,特别是周转率达到100.0%的月份较多;（4） 不同生境斑块间仅轮牧区和禁牧区中长爪沙鼠种群密度显著正相关（P<0.05）,而其他生境斑块间相关性均不显著（P >0.05）,长爪沙鼠局域种群整体显示出明显的非同步空间动态。阿拉善荒漠区长爪沙鼠种群满足作为经典集合种群物种区域续存的4个条件,具有作为研究小哺乳动物集合种群的潜在价值。     

Quantifying spatial correlations in extinction risk for an aphid metapopulation

Metapopulation persistence depends on dispersal offsetting local extinction risks. The mitigating power of dispersal to offset local extinction risks is substantially reduced if extinction risks are highly spatially correlated. Extinction risks are likely to be spatially correlated when mechanisms driving local extinction act on large spatial scales, such as climatic events or far-ranging predators. In this paper I develop a method to quantify spatial correlations in local extinction risks within a region. I then show that colonies of the aphid Aphis lugentis within clusters of their host plant, Senecio bigelovii, do not have strong spatial correlations in extinction risks, even though the primary mechanism for colony extinction is over-predation by far-ranging predators (ladybugs and syrphid flies). Consequently, dispersal potentially plays an important role offsetting local extinction. Aphids did experience a weak, but significantly positive correlation in extinction risk in one of the two years of this study as well as increasing extinction risks across all plant clusters in both years. These observations point to important gaps in our understanding of how population synchrony affects metapopulation persistence.     

Dynamics of Canadian lynx populations in space and time1

The records of Canadian lynx fur returns indicate rather regular self-repeating pattern with 9–11 yr cycles That the lynx populations over a vast geographical range fluctuate in synchrony appears to be well established The finding that the degree of synchrony levels off with increasing distance among provinces compared to increase again between the furthest–away provinces (O–shape) has largely escaped notice in previous analyses Moreover, the pattern of synchrony against distance varies depending on the period of time used m the analysis, sometimes there is O–shape, sometimes mere negative correlation, and sometimes no relationship whatsoever The entire rich pattern of lynx dynamics in space and time can be emulated by sets of local subpopulations acknowledging dispersal among populations This all suggest that local dynamics of populations cannot be fully understood unless the significance of spatial linkage of subpopulations via dispersal of individuals is acknowledged     

Recently introduced invasive geckos quickly reach population genetic equilibrium dynamics

Invasive species are spreading at high rates, yet fundamental processes allowing them to progress through the stages of invasion are unclear. The establishment stage is a critical point because this is when exotic species can survive, reproduce, and begin to spread. Unfortunately, inference of population dynamics during this stage may be impossible if historical and observational data are incomplete. Nonetheless, critical inferences on population dynamics during the establishment stage can be acquired indirectly by characterizing demographic history via the population genetics of recently introduced populations. Geckos have been introduced at a global scale and are one of the most successfully establishing families of alien reptile known. Here we conduct a series of population genetic analyses among five close subpopulations of the introduced Mediterranean gecko Hemidactylus turcicus. We tested for non-equilibrium genetic signatures, a pattern expected during early stages of invasion if there were few founders or repeated introductions led to population turnover. Genetic analyses showed no evidence of non-equilibrium dynamics such as genetic bottlenecks. Moreover, we found strong support for population genetic equilibrium dynamics. The observed results may have been generated via an introduction that involved high propagule pressure. However, given the life history of H. turcicus including generation time and dispersal potential, we favor the hypothesis that the invasive metapopulation has rapidly reached the establishment stage as indicated by relatively constant effective sizes and migration rates among introduced subpopulations. The ability to rapidly pass through the establishment stage may in part explain the invasion success of these geckos.     

Synchrony of butterfly populations across species' geographic ranges

Understanding the mechanisms by which global climate change and habitat loss impact upon biodiversity is essential in order to mitigate any negative impacts. One such impact may be changes to population synchrony (defined as correlated fluctuations in the density of separate populations). It is well established that synchrony depends on both dispersal ability and correlated environmental conditions, for example shared climate. However, what is not clear is whether differences in habitat or position within a species' range also mediate synchrony. Since synchronous metapopulations are thought to be more extinction‐prone, establishing the drivers of synchrony has clear conservation implications. Using three butterfly species (Maniola jurtina, Pyronia tithonus and Aphantopus hyperantus) we investigated the effects of habitat similarity and range position on population synchrony, after accounting for the effects of distance and climate. Range position was present in all minimum adequate models, though non‐significant using Mantel randomization tests in one case. We show that M. jurtina and P. tithonus synchrony is not consistent across species' ranges, with marginal populations showing more synchronous dynamics. Increased climatic constraints on marginal populations, leading to a narrower range of suitable microhabitats may be responsible for this, which is supported by the result that habitat similarity between sites was also positively correlated with population synchrony. As the landscape becomes increasingly homogeneous, overall population synchrony may be expected to rise. We conclude that habitat modification and climate change have the capacity to drive changes in population synchrony that could make species more vulnerable to extinction.     

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.