Spatial synchrony cascades across ecosystem boundaries and up food webs via resource subsidies.

Cross-ecosystem subsidies are critical to ecosystem structure and function, especially in recipient ecosystems where they are the primary source of organic matter to the food web. Subsidies are indicative of processes connecting ecosystems and can couple ecological dynamics across system boundaries. However, the degree to which such flows can induce cross-ecosystem cascades of spatial synchrony, the tendency for system fluctuations to be correlated across locations, is not well understood. Synchrony has destabilizing effects on ecosystems, adding to the importance of understanding spatiotemporal patterns of synchrony transmission. In order to understand whether and how spatial synchrony cascades across the marine-terrestrial boundary via resource subsidies, we studied the relationship between giant kelp forests on rocky nearshore reefs and sandy beach ecosystems that receive resource subsidies in the form of kelp wrack (detritus). We found that synchrony cascades from rocky reefs to sandy beaches, with spatiotemporal patterns mediated by fluctuations in live kelp biomass, wave action, and beach width. Moreover, wrack deposition synchronized local abundances of shorebirds that move among beaches seeking to forage on wrack-associated invertebrates, demonstrating that synchrony due to subsidies propagates across trophic levels in the recipient ecosystem. Synchronizing resource subsidies likely play an underappreciated role in the spatiotemporal structure, functioning, and stability of ecosystems.

Tail-dependence of masting synchrony results in continent-wide seed scarcity.

Spatial synchrony may be tail-dependent, meaning it is stronger for peaks rather than troughs, or vice versa. High interannual variation in seed production in perennial plants, called masting, can be synchronized at subcontinental scales, triggering extensive resource pulses or famines. We used data from 99 populations of European beech (Fagus sylvatica) to examine whether masting synchrony differs between mast peaks and years of seed scarcity. Our results revealed that seed scarcity occurs simultaneously across the majority of the species range, extending to populations separated by distances up to 1800 km. Mast peaks were spatially synchronized at distances up to 1000 km and synchrony was geographically concentrated in northeastern Europe. Extensive synchrony in the masting lower tail means that famines caused by beech seed scarcity are amplified by their extensive spatial synchrony, with diverse consequences for food web functioning and climate change biology.

Seasonality in Environment and Population Processes Alters Population Spatial Synchrony.

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

Demographic Basis of Spatially Structured Fluctuations in a Threespine Stickleback Metapopulation.

AbstractUncovering the demographic basis of population fluctuations is a central goal of population biology. This is particularly challenging for spatially structured populations, which require disentangling synchrony in demographic rates from coupling via movement between locations. In this study, we fit a stage-structured metapopulation model to a 29-year time series of threespine stickleback abundance in the heterogeneous and productive Lake Mývatn, Iceland. The lake comprises two basins (North and South) connected by a channel through which the stickleback disperse. The model includes time-varying demographic rates, allowing us to assess the potential contributions of recruitment and survival, spatial coupling via movement, and demographic transience to the population's large fluctuations in abundance. Our analyses indicate that recruitment was only modestly synchronized between the two basins, whereas survival probabilities of adults were more strongly synchronized, contributing to cyclic fluctuations in the lake-wide population size with a period of approximately 6 years. The analyses further show that the two basins were coupled through movement, with the North Basin subsidizing the South Basin and playing a dominant role in driving the lake-wide dynamics. Our results show that cyclic fluctuations of a metapopulation can be explained in terms of the combined effects of synchronized demographic rates and spatial coupling.

Cold winters drive consistent and spatially synchronous 8-year population cycles of cabbage stem flea beetle.

Population cycles have been observed in mammals as well as insects, but consistent population cycling has rarely been documented in agroecosystems and never for a beetle. We analysed the long-term population patterns of the cabbage stem flea beetle Psylliodes chrysocephala in winter oilseed rape over 50 years. Psylliodes chrysocephala larval density from 3045 winter oilseed rape fields in southern Sweden showed strong 8-year population cycles in regional mean density. Fluctuations in larval density were synchronous over time across five subregional populations. Subregional mean environmental variables explained 90.6% of the synchrony in P. chrysocephala populations at the 7-11 year time-scale. The number of days below -10°C showed strong anti-phase coherence with larval densities in the 7-11 year time-scale, such that more cold days resulted in low larval densities. High levels of the North Atlantic Oscillation weather system are coherent and anti-phase with cold weather in Scania, Sweden. At the field-scale, later crop planting date and more cold winter days were associated with decreased overwintering larval density. Warmer autumn temperatures, resulting in greater larval accumulated degree days early in the season, increased overwintering larval density. Despite variation in environmental conditions and crop management, 8-year cycles persisted for cabbage stem flea beetle throughout the 50 years of data collection. Moran effects, influenced by the North Atlantic Oscillation weather patterns, are the primary drivers of this cycle and synchronicity. Insect pest data collected in commercial agriculture fields is an abundant source of long-term data. We show that an agricultural pest can have the same periodic population cycles observed in perennial and unmanaged ecosystems. This unexpected finding has implications for sustainable pest management in agriculture and shows the value of long-term pest monitoring projects as an additional source of time-series data to untangle the drivers of population cycles.

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

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

Dispersal Increases Spatial Synchrony of Populations but Has Weak Effects on Population Variability: A Meta-analysis.

AbstractThe effects of dispersal on spatial synchrony and population variability have been well documented in theoretical research, and a growing number of empirical tests have been performed. Yet a synthesis is still lacking. Here, we conducted a meta-analysis of relevant experiments and examined how dispersal affected spatial synchrony and temporal population variability across scales. Our analyses showed that dispersal generally promoted spatial synchrony, and such effects increased with dispersal rate and decreased with environmental correlation among patches. The synchronizing effect of dispersal, however, was detected only when spatial synchrony was measured using the correlation-based index, not when the covariance-based index was used. In contrast to theoretical predictions, the effect of dispersal on local population variability was generally nonsignificant, except when environmental correlation among patches was negative and/or the experimental period was long. At the regional scale, while low dispersal stabilized metapopulation dynamics, high dispersal led to destabilization. Overall, the sign and strength of dispersal effects on spatial synchrony and population variability were modulated by taxa, environmental heterogeneity, type of perturbations, patch number, and experimental length. Our synthesis demonstrates that dispersal can affect the dynamics of spatially distributed populations, but its effects are context dependent on abiotic and biotic factors.

The geography of metapopulation synchrony in dendritic river networks.

Dendritic habitats, such as river ecosystems, promote the persistence of species by favouring spatial asynchronous dynamics among branches. Yet, our understanding of how network topology influences metapopulation synchrony in these ecosystems remains limited. Here, we introduce the concept of fluvial synchrogram to formulate and test expectations regarding the geography of metapopulation synchrony across watersheds. By combining theoretical simulations and an extensive fish population time-series dataset across Europe, we provide evidence that fish metapopulations can be buffered against synchronous dynamics as a direct consequence of network connectivity and branching complexity. Synchrony was higher between populations connected by direct water flow and decayed faster with distance over the Euclidean than the watercourse dimension. Likewise, synchrony decayed faster with distance in headwater than mainstem populations of the same basin. As network topology and flow directionality generate fundamental spatial patterns of synchrony in fish metapopulations, empirical synchrograms can aid knowledge advancement and inform conservation strategies in complex habitats.

Drought effects on invertebrate metapopulation dynamics and quasi-extinction risk in an intermittent river network.

Ecological communities can remain stable in the face of disturbance if their constituent species have different resistance and resilience strategies. In turn, local stability scales up regionally if heterogeneous landscapes maintain spatial asynchrony across discrete populations-but not if large-scale stressors synchronize environmental conditions and biological responses. Here, we hypothesized that droughts could drastically decrease the stability of invertebrate metapopulations both by filtering out poorly adapted species locally, and by synchronizing their dynamics across a river network. We tested this hypothesis via multivariate autoregressive state-space (MARSS) models on spatially replicated, long-term data describing aquatic invertebrate communities and hydrological conditions in a set of temperate, lowland streams subject to seasonal and supraseasonal drying events. This quantitative approach allowed us to assess the influence of local (flow magnitude) and network-scale (hydrological connectivity) drivers on invertebrate long-term trajectories, and to simulate near-future responses to a range of drought scenarios. We found that fluctuations in species abundances were heterogeneous across communities and driven by a combination of hydrological and stochastic drivers. Among metapopulations, increasing extent of dry reaches reduced the abundance of functional groups with low resistance or resilience capacities (i.e. low ability to persist in situ or recolonize from elsewhere, respectively). Our simulations revealed that metapopulation quasi-extinction risk for taxa vulnerable to drought increased exponentially as flowing habitats contracted within the river network, whereas the risk for taxa with resistance and resilience traits remained stable. Our results suggest that drought can be a synchronizing agent in riverscapes, potentially leading to regional quasi-extinction of species with lower resistance and resilience abilities. Better recognition of drought-driven synchronization may increase realism in species extinction forecasts as hydroclimatic extremes continue to intensify worldwide.

Synergistic effects of harvest and climate drive synchronous somatic growth within key New Zealand fisheries.

Fisheries harvest has pervasive impacts on wild fish populations, including the truncation of size and age structures, altered population dynamics and density, and modified habitat and assemblage composition. Understanding the degree to which harvest-induced impacts increase the sensitivity of individuals, populations and ultimately species to environmental change is essential to ensuring sustainable fisheries management in a rapidly changing world. Here we generated multiple long-term (44-62 years), annually resolved, somatic growth chronologies of four commercially important fishes from New Zealand's coastal and shelf waters. We used these novel data to investigate how regional- and basin-scale environmental variability, in concert with fishing activity, affected individual somatic growth rates and the magnitude of spatial synchrony among stocks. Changes in somatic growth can affect individual fitness and a range of population and fishery metrics such as recruitment success, maturation schedules and stock biomass. Across all species, individual growth benefited from a fishing-induced release of density controls. For nearshore snapper and tarakihi, regional-scale wind and temperature also additively affected growth, indicating that future climate change-induced warming and potentially strengthened winds will initially promote the productivity of more poleward populations. Fishing increased the sensitivity of deep-water hoki and ling growth to the Interdecadal Pacific Oscillation (IPO). A forecast shift to a positive IPO phase, in concert with current harvest strategies, will likely promote individual hoki and ling growth. At the species level, historical fishing practices and IPO synergized to strengthen spatial synchrony in average growth between stocks separated by 400-600 nm of ocean. Increased spatial synchrony can, however, increase the vulnerability of stocks to deleterious stochastic events. Together, our individual- and species-level results show how fishing and environmental factors can conflate to initially promote individual growth but then possibly heighten the sensitivity of stocks to environmental change.

Spatial and temporal variation of ecosystem properties at macroscales.

Although spatial and temporal variation in ecological properties has been well-studied, crucial knowledge gaps remain for studies conducted at macroscales and for ecosystem properties related to material and energy. We test four propositions of spatial and temporal variation in ecosystem properties within a macroscale (1000 km's) extent. We fit Bayesian hierarchical models to thousands of observations from over two decades to quantify four components of variation - spatial (local and regional) and temporal (local and coherent); and to model their drivers. We found strong support for three propositions: (1) spatial variation at local and regional scales are large and roughly equal, (2) annual temporal variation is mostly local rather than coherent, and, (3) spatial variation exceeds temporal variation. Our findings imply that predicting ecosystem responses to environmental changes at macroscales requires consideration of the dominant spatial signals at both local and regional scales that may overwhelm temporal signals.

Population spatial synchrony enhanced by periodicity and low detuning with environmental forcing.

Explaining why fluctuations in abundances of spatially disjunct populations often are correlated through time is a major goal of population ecologists. We address two hypotheses receiving little to no testing in wild populations: (i) that population cycling facilitates synchronization given weak coupling among populations, and (ii) that the ability of periodic external forces to synchronize oscillating populations is a function of the mismatch in timescales (detuning) between the force and the population. Here, we apply new analytical methods to field survey data on gypsy moth outbreaks. We report that at timescales associated with gypsy moth outbreaks, spatial synchrony increased with population periodicity via phase locking. The extent to which synchrony in temperature and precipitation influenced population synchrony was associated with the degree of mismatch in dominant timescales of oscillation. Our study provides new empirical methods and rare empirical evidence that population cycling and low detuning can promote population spatial synchrony.

The potential role of intrinsic processes in generating abrupt and quasi-synchronous tree declines during the Holocene.

Multiple abrupt and sometimes near-synchronous declines in tree populations have been detected in the temperate forests of eastern North America and Europe during the Holocene. Traditional approaches to understanding these declines focus on searching for climatic or other broad-scale extrinsic drivers. These approaches include multi-proxy studies that match reconstructed changes in tree abundance to reconstructed changes in precipitation or temperature. Although these correlative approaches are informative, they neglect the potential role of intrinsic processes, such as competition and dispersal, in shaping tree community dynamics. We developed a simple process-based community model that includes competition among tree species, density-dependent survival, and dispersal to investigate how these processes might generate abrupt changes in tree abundances even when extrinsic climatic factors do not themselves change abruptly. Specifically, a self-reinforcing (i.e., positive) feedback between abundance and survival can produce abrupt changes in tree abundance in the absence of long-term climatic changes. Furthermore, spatially correlated, short-term environmental variation and seed dispersal can increase the synchrony of abrupt changes. Using the well-studied, late-Holocene crash of Tsuga canadensis (eastern hemlock) populations as an empirical case study, we find that our model generates abrupt and quasi-synchronized crashes qualitatively similar to the observed hemlock patterns. Other tree taxa vary in the frequency and clustering of abrupt change and the proportion of increases and decreases. This complexity argues for caution in interpreting abrupt changes in species abundances as indicative of abrupt climatic changes. Nonetheless, some taxa show patterns that the model cannot produce: observed abrupt declines in hemlock abundance are more synchronized than abrupt increases, whereas the degree of synchronization is the same for abrupt decreases and increases in the model. Our results show that intrinsic processes can be significant contributing factors in abrupt tree population changes and highlight the diagnostic value of analyzing entire time series rather than single events when testing hypotheses about abrupt changes. Thus, intrinsic processes should be considered along with extrinsic drivers when seeking to explain rapid changes in community composition.

Spatial heterogeneity in climate change effects decouples the long-term dynamics of wild reindeer populations in the high Arctic.

The 'Moran effect' predicts that dynamics of populations of a species are synchronized over similar distances as their environmental drivers. Strong population synchrony reduces species viability, but spatial heterogeneity in density dependence, the environment, or its ecological responses may decouple dynamics in space, preventing extinctions. How such heterogeneity buffers impacts of global change on large-scale population dynamics is not well studied. Here, we show that spatially autocorrelated fluctuations in annual winter weather synchronize wild reindeer dynamics across high-Arctic Svalbard, while, paradoxically, spatial variation in winter climate trends contribute to diverging local population trajectories. Warmer summers have improved the carrying capacity and apparently led to increased total reindeer abundance. However, fluctuations in population size seem mainly driven by negative effects of stochastic winter rain-on-snow (ROS) events causing icing, with strongest effects at high densities. Count data for 10 reindeer populations 8-324 km apart suggested that density-dependent ROS effects contributed to synchrony in population dynamics, mainly through spatially autocorrelated mortality. By comparing one coastal and one 'continental' reindeer population over four decades, we show that locally contrasting abundance trends can arise from spatial differences in climate change and responses to weather. The coastal population experienced a larger increase in ROS, and a stronger density-dependent ROS effect on population growth rates, than the continental population. In contrast, the latter experienced stronger summer warming and showed the strongest positive response to summer temperatures. Accordingly, contrasting net effects of a recent climate regime shift-with increased ROS and harsher winters, yet higher summer temperatures and improved carrying capacity-led to negative and positive abundance trends in the coastal and continental population respectively. Thus, synchronized population fluctuations by climatic drivers can be buffered by spatial heterogeneity in the same drivers, as well as in the ecological responses, averaging out climate change effects at larger spatial scales.

Occasional long distance dispersal increases spatial synchrony of population cycles.

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

Forests synchronize their growth in contrasting Eurasian regions in response to climate warming.

Forests play a key role in the carbon balance of terrestrial ecosystems. One of the main uncertainties in global change predictions lies in how the spatiotemporal dynamics of forest productivity will be affected by climate warming. Here we show an increasing influence of climate on the spatial variability of tree growth during the last 120 y, ultimately leading to unprecedented temporal coherence in ring-width records over wide geographical scales (spatial synchrony). Synchrony in growth patterns across cold-constrained (central Siberia) and drought-constrained (Spain) Eurasian conifer forests have peaked in the early 21st century at subcontinental scales (∼ 1,000 km). Such enhanced synchrony is similar to that observed in trees co-occurring within a stand. In boreal forests, the combined effects of recent warming and increasing intensity of climate extremes are enhancing synchrony through an earlier start of wood formation and a stronger impact of year-to-year fluctuations of growing-season temperatures on growth. In Mediterranean forests, the impact of warming on synchrony is related mainly to an advanced onset of growth and the strengthening of drought-induced growth limitations. Spatial patterns of enhanced synchrony represent early warning signals of climate change impacts on forest ecosystems at subcontinental scales.

Vertical and Seasonal Patterns Control Bacterioplankton Communities at Two Horizontally Coherent Coastal Upwelling Sites off Galicia (NW Spain).

Analysis of seasonal patterns of marine bacterial community structure along horizontal and vertical spatial scales can help to predict long-term responses to climate change. Several recent studies have shown predictable seasonal reoccurrence of bacterial assemblages. However, only a few have assessed temporal variability over both horizontal and vertical spatial scales. Here, we simultaneously studied the bacterial community structure at two different locations and depths in shelf waters of a coastal upwelling system during an annual cycle. The most noticeable biogeographic patterns observed were seasonality, horizontal homogeneity, and spatial synchrony in bacterial diversity and community structure related with regional upwelling-downwelling dynamics. Water column mixing eventually disrupted bacterial community structure vertical heterogeneity. Our results are consistent with previous temporal studies of marine bacterioplankton in other temperate regions and also suggest a marked influence of regional factors on the bacterial communities inhabiting this coastal upwelling system. Bacterial-mediated carbon fluxes in this productive region appear to be mainly controlled by community structure dynamics in surface waters, and local environmental factors at the base of the euphotic zone.

The geography of spatial synchrony.

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

Modification of a Pollen Trap Design To Capture Airborne Conidia of Entomophaga maimaiga and Detection of Conidia by Quantitative PCR.

The goal of this study was to develop effective and practical field sampling methods for quantification of aerial deposition of airborne conidia of Entomophaga maimaiga over space and time. This important fungal pathogen is a major cause of larval death in invasive gypsy moth (Lymantria dispar) populations in the United States. Airborne conidia of this pathogen are relatively large (similar in size to pollen), with unusual characteristics, and require specialized methods for collection and quantification. Initially, dry sampling (settling of spores from the air onto a dry surface) was used to confirm the detectability of E. maimaiga at field sites with L. dispar deaths caused by E. maimaiga, using quantitative PCR (qPCR) methods. We then measured the signal degradation of conidial DNA on dry surfaces under field conditions, ultimately rejecting dry sampling as a reliable method due to rapid DNA degradation. We modified a chamber-style trap commonly used in palynology to capture settling spores in buffer. We tested this wet-trapping method in a large-scale (137-km) spore-trapping survey across gypsy moth outbreak regions in Pennsylvania undergoing epizootics, in the summer of 2016. Using 4-day collection periods during the period of late instar and pupal development, we detected variable amounts of target DNA settling from the air. The amounts declined over the season and with distance from the nearest defoliated area, indicating airborne spore dispersal from outbreak areas.IMPORTANCE We report on a method for trapping and quantifying airborne spores of Entomophaga maimaiga, an important fungal pathogen affecting gypsy moth (Lymantria dispar) populations. This method can be used to track dispersal of E. maimaiga from epizootic areas and ultimately to provide critical understanding of the spatial dynamics of gypsy moth-pathogen interactions.

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

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