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.     

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.     

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.     

How does intraspecific density regulation influence metapopulation synchrony and persistence?

Intraspecific density regulation influences the synchronization of local population dynamics through dispersal. Spatial synchrony in turn may jeopardize metapopulation persistence. Joining results from previous studies suggests that spatial synchrony is highest at moderate over-compensation and is low at compensating and at very strong over-compensating density regulation. We scrutinize this supposition of a unimodal relationship using a process-based metapopulation model with explicit local population dynamics. We extend the usually studied range of density regulation to under-compensation and analyse resulting metapopulation persistence. We find peaks of spatial synchrony not only at over-compensatory but also under-compensatory density regulation and show that effects of local density compensation on synchrony follow a bimodal rather than unimodal relationship. Persistence of metapopulations however, shows a unimodal relationship with a broad plateau of high persistence from compensatory to over-compensatory density regulation. This range of high persistence comprises both levels of low and high spatial synchrony. Thus, not synchrony alone jeopardizes metapopulation persistence, but only in interplay with high local extinction risk. The functional forms of the relations of density compensation with spatial synchrony and persistence are robust to increases in dispersal mortality, landscape dynamics, or density dependence of dispersal. However, with each of these increases the maxima of spatial synchrony and persistence shift to higher over-compensation and levels of synchrony are reduced. Overall, for over-compensation high landscape connectivity has negative effects while for under-compensation connectivity affects persistence positively. This emphasizes the importance of species life-history traits for management decisions with regard to landscape connectivity: while dispersal corridors are essential for species with under-compensatory density regulation, they may have detrimental effects for endangered species with over-compensation.     

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.     

The influence of multiple dispersal mechanisms and landscape structure on population clustering and connectivity in fragmented artesian spring snail populations

Many organisms occupy heterogeneous landscapes that contain both barriers to movement as well as corridors that facilitate dispersal. The extent to which such features determine population connectivity will depend on the mechanisms utilized by organisms to disperse. Here we examined the interaction between landscape structure and dispersal in the endemic aquatic snail, Fonscochlea accepta , in the fragmented artesian spring ecosystem of arid central Australia. We used frequentist and Bayesian analyses of microsatellite data to identify population structure and immigration for 1130 snails sampled from 50 springs across an entire spring complex. We introduce a modified isolation-by-distance analysis to test hypotheses about how populations are clustered and to distinguish the most likely dispersal pathways within and between those clusters. Highly significant differences in F _ST values and significant isolation-by-distance patterns were detected among springs across the entire complex, while Bayesian assignment tests revealed the presence of two hierarchical levels of spring clustering. Clusters were defined by the spatial aggregation of springs, dynamic aquatic habitat connections between springs and the ecology of the snails. Bayesian immigrant identification and our modified isolation-by-distance analysis revealed that dispersal occurs at two geographical scales via two very different mechanisms. Short range dispersal (usually ≤ 300 m) occurs via active movement facilitated by aquatic connections among springs while long-range dispersal (≥ 3 km) is likely facilitated by an animal vector (phoresy). These results underline the importance of both dispersal mode and landscape structure in influencing connectivity rates and patterns among 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.     

Spatial synchrony of wader populations in inland lakes of the Iberian Peninsula

Spatial synchronization refers to similarity in temporal variations between spatially separated populations. Three mechanisms have been associated with the spatial synchrony of populations: Moran effect, dispersal and trophic interactions. In this study, we explored the degree of spatial synchrony of three wader species populations (Pied Avocet, Black-winged Stilt and Kentish Plover) using monthly estimates of their abundance in inland lakes of the Iberian Peninsula. The effect of several types of wetland variables (structural, hydroperiod and landscape) on spatial synchronization was explored. Groups of lakes with significant synchronization were identified for all three species. The lakes with wastewater input presented longer hydroperiods than those that did not receive these effluents, and this factor was positively related to the spatial synchrony of the Pied Avocet and Kentish Plover populations. The distance between lakes (used as an indicator of the dispersal effect on synchronization) was significant only in Pied Avocet. No structural or landscape variables were related to spatial synchronization in any species. It was impossible to identify any variable related to the spatial synchronization of Black-winged Stilt abundance as a possible result of the high ecological plasticity of this species. Our data provides the first evidence for mechanisms that act on the spatial synchronizing of wader populations in temporary continental lakes in central Spain, and show that the hydroperiod of lakes acts as an important factor in the spatial synchronization of aquatic species and that its effect is mediated by the reception of urban wastewater.     

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.     

Density-dependent dispersal and spatial population dynamics

The synchronization of the dynamics of spatially subdivided populations is of both fundamental and applied interest in population biology. Based on theoretical studies, dispersal movements have been inferred to be one of the most general causes of population synchrony, yet no empirical study has mapped distance-dependent estimates of movement rates on the actual pattern of synchrony in species that are known to exhibit population synchrony. Northern vole and lemming species are particularly well-known for their spatially synchronized population dynamics. Here, we use results from an experimental study to demonstrate that tundra vole dispersal movements did not act to synchronize population dynamics in fragmented habitats. In contrast to the constant dispersal rate assumed in earlier theoretical studies, the tundra vole, and many other species, exhibit negative density-dependent dispersal. Simulations of a simple mathematical model, parametrized on the basis of our experimental data, verify the empirical results, namely that the observed negative density-dependent dispersal did not have a significant synchronizing effect.     

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.     

Synchrony of woodland bird populations: the effect of landscape structure

The influence of environmental stochasticity and dispersal in producing patterns in population synchrony was examined for 53 woods censused annually from 1990 to 1999 for nine resident bird species (wren Troglodytes troglodytes , dunnock Prunella modularis , robin Erithacus rubecula , blackbird Turdus merula , song thrush Turdus philomelos , long-tailed tit Aegithalos caudatus , blue tit Parus caeruleus , great tit Parus major , and chaffinch Fringilla coelebs ) and four migrant bird species (garden warbler Sylvia borin , blackcap Sylvia atricapilla , chiffchaff Phylloscopus collybita and willow warbler Phylloscopus trochilus ). Twelve species showed global synchrony of population counts due to regional population trends and widespread annual population fluctuations. There was a clear link between population fluctuations and winter weather for wren, and three other species showed their largest population declines after the coldest winters. Eight species showed a decline in synchrony with distance between woods, and there was evidence for dispersal causing this pattern in three species. Landscape structure affected patterns of synchrony in several species, with lower synchrony in landscapes with less woodland. For three species, this difference in synchrony across a landscape gradient of decreasing woodland cover accounted for the decline in synchrony over distance. Three species showed greater synchrony between woods with similar amounts of hedgerow in the surrounding landscape, suggesting that the surroundings of a wood influence the population dynamics of some species breeding in the wood. Habitat fragmentation can alter the processes contributing to population synchrony. Loss of woodland reduces the relative abundance of woodland bird species. The remaining patches of habitat are smaller, more isolated and are set in a more hostile landscape, all of which may disrupt dispersal between patches and alter the population dynamics within woods.     

The walk is never random: subtle landscape effects shape gene flow in a continuous white-tailed deer population in the Midwestern United States

One of the pervasive challenges in landscape genetics is detecting gene flow patterns within continuous populations of highly mobile wildlife. Understanding population genetic structure within a continuous population can give insights into social structure, movement across the landscape and contact between populations, which influence ecological interactions, reproductive dynamics or pathogen transmission. We investigated the genetic structure of a large population of deer spanning the area of Wisconsin and Illinois, USA, affected by chronic wasting disease. We combined multiscale investigation, landscape genetic techniques and spatial statistical modelling to address the complex questions of landscape factors influencing population structure. We sampled over 2000 deer and used spatial autocorrelation and a spatial principal components analysis to describe the population genetic structure. We evaluated landscape effects on this pattern using a spatial autoregressive model within a model selection framework to test alternative hypotheses about gene flow. We found high levels of genetic connectivity, with gradients of variation across the large continuous population of white-tailed deer. At the fine scale, spatial clustering of related animals was correlated with the amount and arrangement of forested habitat. At the broader scale, impediments to dispersal were important to shaping genetic connectivity within the population. We found significant barrier effects of individual state and interstate highways and rivers. Our results offer an important understanding of deer biology and movement that will help inform the management of this species in an area where overabundance and disease spread are primary concerns.     

Spatial dynamics of Norwegian tetraonid populations

Different species in a given site or population of a given species in different sites may fluctuate in synchrony if they are affected similarly by factors such as spatially autocorrelated climate, predation, or by dispersal between populations of one species. We used county wise time series of hunting bag records of four Norwegian tetraonid species covering 24 years to examine patterns of interspecific and intraspecific synchrony. We estimated synchrony at three spatial scales; national, regional (consisting of counties with similar climate), and county level. Ecologically related species with overlapping distributions exhibited strong synchrony across Norway, but there was much variation between the different regions and counties. Regions with a long coastline to both the North Sea and the Norwegian Ocean exhibited an overall stronger synchrony than those consisting of more continental areas. Intraspecific synchrony was generally low across all counties, but stronger synchrony between counties within regions defined by climatic conditions. Synchrony was negatively related to distance between populations in three of four species. Only the synchrony in willow ptarmigan showed a clear negative relationship with distance, while the other species had both strong positive and negative correlations at short distances. Strong interspecific synchrony between some species pairs within regions and weak intraspecific synchrony across counties within regions suggest a stronger synchronizing effect from environmental factors such as weather or predation and less effect from dispersal. Our results suggest that the complete tetraonid community is structured by environmental factors affecting the different species similarly and causes widespread interspecific synchrony. Local factors affecting the population dynamics nevertheless frequently forces neighbouring populations out of phase.     

Spatial synchrony in red grouse population dynamics

Red grouse Lagopus lagopus scoticus populations exhibit unstable dynamics that are often characterised by regular periodic fluctuations in abundance. Time-series' of grouse harvesting records collected from 287 management units (moors) across Scotland, England and Wales were analysed to investigate the broad scale patterns of synchrony in these fluctuations. Estimation of the spatial autocorrelation of grouse population dynamics across moors indicates relatively high levels of synchrony between populations on adjacent moors, but that this synchrony declines sharply with increasing inter-moor distance. At distances of greater than 100  km, grouse population time-series exhibit only weakly positive cross-correlation coefficients. Twenty-eight geographical, environmental and other candidate variables were examined to construct a general linear model to explain variation in local synchrony. Grouse moor productivity (average size of shooting bag), distance from the Atlantic coast moving in a north-easterly direction, April and June temperatures, and June rainfall significantly increased the explanatory power of this model. An understanding of the processes underlying synchrony in red grouse population dynamics is a prerequisite to anticipating the effects of large-scale environmental change on regional patterns of grouse distribution and abundance.     

Implications of recurrent disturbance for genetic diversity

Exploring interactions between ecological disturbance, species’ abundances and community composition provides critical insights for ecological dynamics. While disturbance is also potentially an important driver of landscape genetic patterns, the mechanisms by which these patterns may arise by selective and neutral processes are not well‐understood. We used simulation to evaluate the relative importance of disturbance regime components, and their interaction with demographic and dispersal processes, on the distribution of genetic diversity across landscapes. We investigated genetic impacts of variation in key components of disturbance regimes and spatial patterns that are likely to respond to climate change and land management, including disturbance size, frequency, and severity. The influence of disturbance was mediated by dispersal distance and, to a limited extent, by birth rate. Nevertheless, all three disturbance regime components strongly influenced spatial and temporal patterns of genetic diversity within subpopulations, and were associated with changes in genetic structure. Furthermore, disturbance‐induced changes in temporal population dynamics and the spatial distribution of populations across the landscape resulted in disrupted isolation by distance patterns among populations. Our results show that forecast changes in disturbance regimes have the potential to cause major changes to the distribution of genetic diversity within and among populations. We highlight likely scenarios under which future changes to disturbance size, severity, or frequency will have the strongest impacts on population genetic patterns. In addition, our results have implications for the inference of biological processes from genetic data, because the effects of dispersal on genetic patterns were strongly mediated by disturbance regimes.     

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.     

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.     

A general modeling framework for describing spatially structured population dynamics

Variation in movement across time and space fundamentally shapes the abundance and distribution of populations. Although a variety of approaches model structured population dynamics, they are limited to specific types of spatially structured populations and lack a unifying framework. Here, we propose a unified network‐based framework sufficiently novel in its flexibility to capture a wide variety of spatiotemporal processes including metapopulations and a range of migratory patterns. It can accommodate different kinds of age structures, forms of population growth, dispersal, nomadism and migration, and alternative life‐history strategies. Our objective was to link three general elements common to all spatially structured populations (space, time and movement) under a single mathematical framework. To do this, we adopt a network modeling approach. The spatial structure of a population is represented by a weighted and directed network. Each node and each edge has a set of attributes which vary through time. The dynamics of our network‐based population is modeled with discrete time steps. Using both theoretical and real‐world examples, we show how common elements recur across species with disparate movement strategies and how they can be combined under a unified mathematical framework. We illustrate how metapopulations, various migratory patterns, and nomadism can be represented with this modeling approach. We also apply our network‐based framework to four organisms spanning a wide range of life histories, movement patterns, and carrying capacities. General computer code to implement our framework is provided, which can be applied to almost any spatially structured population. This framework contributes to our theoretical understanding of population dynamics and has practical management applications, including understanding the impact of perturbations on population size, distribution, and movement patterns. By working within a common framework, there is less chance that comparative analyses are colored by model details rather than general principles.     

Spatial synchrony of recruitment in mountain-dwelling woodland caribou

Spatial synchrony in population dynamics is a ubiquitous feature across a range of taxa. Understanding factors influencing this synchrony may shed light on important drivers of population dynamics. Three mechanisms influence the degree of spatial synchrony between populations: dispersal, shared predators, and spatial environmental covariance (the Moran effect). We assessed demographic spatial synchrony in recruitment (calf:cow ratio) of 10 northern mountain caribou herds in the Yukon Territory, Canada (1982–2008). Shared predators and dispersal were ruled out as causal mechanisms of spatial recruitment synchrony in these herds and therefore any spatial synchrony should be due to the Moran effect. We also assessed the degree of spatial synchrony in April snow depth to represent environmental variability. The regional average spatial synchrony in detrended residuals of April snow depth was 0.46 (95% CI 0.37 to 0.55). Spatial synchrony in caribou recruitment was weak at 0.13 (95% CI −0.06 to 0.32). The spatial scale of synchrony in April snow depth and caribou recruitment was 330.2 km (95% CI 236.3 to 370.0 km) and 170.0 km (95% CI 69.5 to 282.8 km), respectively. We also investigated how the similarity in terrain features between herds influenced the degree of spatial synchrony using exponential decay models. Only the difference in elevation variability between herds during calving was supported by the data. Herds with more similar elevation variability may track snowmelt ablation patterns in a more similar fashion, which would subsequently result in more synchronized predation rates on calves and/or nutritional effects impacting juvenile survival. Interspecific interactions with predators and alternate prey may also influence spatial synchrony of recruitment in these herds.