Evolution of body condition‐dependent dispersal in metapopulations

Body condition‐dependent dispersal strategies are common in nature. Although it is obvious that environmental constraints may induce a positive relationship between body condition and dispersal, it is not clear whether positive body conditional dispersal strategies may evolve as a strategy in metapopulations. We have developed an individual‐based simulation model to investigate how body condition–dispersal reaction norms evolve in metapopulations that are characterized by different levels of environmental stochasticity and dispersal mortality. In the model, body condition is related to fecundity and determined either by environmental conditions during juvenile development (adult dispersal) or by those experienced by the mother (natal dispersal). Evolutionarily stable reaction norms strongly depend on metapopulation conditions: positive body condition dependency of dispersal evolved in metapopulation conditions with low levels of dispersal mortality and high levels of environmental stochasticity. Negative body condition‐dependent dispersal evolved in metapopulations with high dispersal mortality and low environmental stochasticity. The latter strategy is responsible for higher dispersal rates under kin competition when dispersal decisions are based on body condition reached at the adult life stage. The evolution of both positive and negative body condition‐dependent dispersal strategies is consequently likely in metapopulations and depends on the prevalent environmental conditions.     

Spatial connectedness imposes local‐ and metapopulation‐level selection on life history through feedbacks on demography

Dispersal evolution impacts the fluxes of individuals and hence, connectivity in metapopulations. Connectivity is therefore decoupled from the structural connectedness of the patches within the spatial network. Because of demographic feedbacks, local selection also drives the evolution of other life history traits. We investigated how different levels of connectedness affect trait evolution in experimental metapopulations of the two‐spotted spider mite. We separated local‐ and metapopulation‐level selection and linked trait divergence to population dynamics. With lower connectedness, an increased starvation resistance and delayed dispersal evolved. Reproductive performance evolved locally by transgenerational plasticity or epigenetic processes. Costs of dispersal, but also changes in local densities and temporal fluctuations herein are found to be putative drivers. In addition to dispersal, demographic traits are able to evolve in response to metapopulation connectedness at both the local and metapopulation level by genetic and/or non‐genetic inheritance. These trait changes impact the persistence of spatially structured populations.     

Classical metapopulation dynamics and eco‐evolutionary feedbacks in dendritic networks

Eco‐evolutionary dynamics are now recognized to be highly relevant for population and community dynamics. However, the impact of evolutionary dynamics on spatial patterns, such as the occurrence of classical metapopulation dynamics, is less well appreciated. Here, we analyse the evolutionary consequences of spatial network connectivity and topology for dispersal strategies and quantify the eco‐evolutionary feedback in terms of altered classical metapopulation dynamics. We find that network properties, such as topology and connectivity, lead to predictable spatio‐temporal correlations in fitness expectations. These spatio‐temporally stable fitness patterns heavily impact evolutionarily stable dispersal strategies and lead to eco‐evolutionary feedbacks on landscape level metrics, such as the number of occupied patches, the number of extinctions and recolonizations as well as metapopulation extinction risk and genetic structure. Our model predicts that classical metapopulation dynamics are more likely to occur in dendritic networks, and especially in riverine systems, compared to other types of landscape configurations. As it remains debated whether classical metapopulation dynamics are likely to occur in nature at all, our work provides an important conceptual advance for understanding the occurrence of classical metapopulation dynamics which has implications for conservation and management of spatially structured populations.     

Tests of two methods for identifying founder effects in metapopulations reveal substantial type II error

Genetic analysis has been promoted as a way to reconstruct recent historical dynamics (“historical demography”) by screening for signatures of events, such as bottlenecks, that disrupt equilibrium patterns of variation. Such analyses might also identify “metapopulation” processes like extinction and recolonization or source-sink dynamics, but this potential remains largely unrealized. Here we use simulations to test the ability of two currently used strategies to distinguish between a set of interconnected subpopulations (demes) that have undergone bottlenecks or extinction and recolonization events (metapopulation dynamics) from a set of static demes. The first strategy, decomposed pairwise regression, provides a holistic test for heterogeneity among demes in their patterns of isolation-by-distance. This method suffered from a type II error rate of 59–100 %, depending on parameter conditions. The second strategy tests for deviations from mutation-drift equilibrium on a deme-by-deme basis to identify sites likely to have experienced recent bottlenecks or founder effects. Although bottleneck tests have good statistical power for single populations with recent population declines, their validity in structured populations has been called into question, and they have not been tested in a metapopulation context with immigration (or colonization) and population recovery. Our simulations of hypothetical metapopulations show that population recovery can rapidly eliminate the statistical signature of a bottleneck, and that moderate levels of gene flow can generate a false signal of recent population growth for demes in equilibrium. Although we did not cover all possible metapopulation scenarios, the performance of the tests was disappointing. Our results indicate that these methods might often fail to identify population bottlenecks and founder effects if population recovery and/or gene flow are influential demographic features of the study system.     

Spatial Heterogeneity in Habitat Quality and Cross-Scale Interactions in Metapopulations

Abstract Integration of habitat heterogeneity into spatially realistic metapopulation approaches reveals the potential for key cross-scale interactions. Broad-scale environmental gradients and land-use practices can create autocorrelation of habitat quality of suitable patches at intermediate spatial scales. Patch occupancy then depends not only on habitat quality at the patch scale but also on feedbacks from surrounding neighborhoods of autocorrelated patches. Metapopulation dynamics emerge from how demographic and dispersal processes interact with relevant habitat heterogeneity. We provide an empirical example from a metapopulation of round-tailed muskrats (Neofiber alleni) in which habitat quality of suitable patches was spatially autocorrelated most strongly within 1,000 m, which was within the expected dispersal range of the species. After controlling for factors typically considered in metapopulation studies—patch size, local patch quality, patch connectivity—we use a cross-variogram analysis to demonstrate that patch occupancy by muskrats was correlated with habitat quality across scales ≤1,171 m. We also discuss general consequences of spatial heterogeneity of habitat quality for metapopulations related to potential cross-scale interactions. We focus on spatially correlated extinctions and metapopulation persistence, hierarchical scaling of source–sink dynamics, and dispersal decisions by individuals in relation to information constraints.     

Condition‐dependent movement and dispersal in experimental metacommunities

Dispersal and the underlying movement behaviour are processes of pivotal importance for understanding and predicting metapopulation and metacommunity dynamics. Generally, dispersal decisions are condition‐dependent and rely on information in the broad sense, like the presence of conspecifics. However, studies on metacommunities that include interspecific interactions generally disregard condition‐dependence. Therefore, it remains unclear whether and how dispersal in metacommunities is condition‐dependent and whether rules derived from single‐species contexts can be scaled up to (meta)communities. Using experimental protist metacommunities, we show how dispersal and movement depend on and are adjusted by the strength of interspecific interactions. We found that the predicting movement and dispersal in metacommunities requires knowledge on behavioural responses to intra‐ and interspecific interaction strengths. Consequently, metacommunity dynamics inferred directly from single‐species metapopulations without taking interspecific interactions into account are likely flawed. Our work identifies the significance of condition‐dependence for understanding metacommunity dynamics, stability and the coexistence and distribution of species.     

Does colonization asymmetry matter in metapopulations?

Despite the considerable evidence showing that dispersal between habitat patches is often asymmetric, most of the metapopulation models assume symmetric dispersal. In this paper, we develop a Monte Carlo simulation model to quantify the effect of asymmetric dispersal on metapopulation persistence. Our results suggest that metapopulation extinctions are more likely when dispersal is asymmetric. Metapopulation viability in systems with symmetric dispersal mirrors results from a mean field approximation, where the system persists if the expected per patch colonization probability exceeds the expected per patch local extinction rate. For asymmetric cases, the mean field approximation underestimates the number of patches necessary for maintaining population persistence. If we use a model assuming symmetric dispersal when dispersal is actually asymmetric, the estimation of metapopulation persistence is wrong in more than 50% of the cases. Metapopulation viability depends on patch connectivity in symmetric systems, whereas in the asymmetric case the number of patches is more important. These results have important implications for managing spatially structured populations, when asymmetric dispersal may occur. Future metapopulation models should account for asymmetric dispersal, while empirical work is needed to quantify the patterns and the consequences of asymmetric dispersal in natural metapopulations.     

Connectivity dominates larval replenishment in a coastal reef fish metapopulation

Direct estimates of larval retention and connectivity are essential to understand the structure and dynamics of marine metapopulations, and optimize the size and spacing of reserves within networks of marine-protected areas (MPAs). For coral reef fishes, while there are some empirical estimates of self-recruitment at isolated populations, exchange among sub-populations has been rarely quantified. Here, we used microsatellite DNA markers and a likelihood-based parentage analysis to assess the relative magnitude of self-recruitment and exchange among eight geographically distinct sub-populations of the panda clownfish Amphiprion polymnus along 30 km of coastline near Port Moresby, Papua New Guinea. In addition, we used an assignment/exclusion test to identify immigrants arriving from genetically distinct sources. Overall, 82 per cent of the juveniles were immigrants while 18 per cent were progeny of parents genotyped in our focal metapopulation. Of the immigrants, only 6 per cent were likely to be genetically distinct from the focal metapopulation, suggesting most of the connectivity is among sub-populations from a rather homogeneous genetic pool. Of the 18 per cent that were progeny of known adults, two-thirds dispersed among the eight sub-populations and only one-third settled back into natal sub-populations. Comparison of our data with previous studies suggested that variation in dispersal distances is likely to be influenced by the geographical setting and spacing of sub-populations.     

Non-random dispersal in the butterfly Maniola jurtina: implications for metapopulation models

The dispersal patterns of animals are important in metapopulation ecology because they affect the dynamics and survival of populations. Theoretical models assume random dispersal but little is known in practice about the dispersal behaviour of individual animals or the strategy by which dispersers locate distant habitat patches. In the present study, we released individual meadow brown butterflies (Maniola jurtina) in a non-habitat and investigated their ability to return to a suitable habitat. The results provided three reasons for supposing that meadow brown butterflies do not seek habitat by means of random flight. First, when released within the range of their normal dispersal distances, the butterflies orientated towards suitable habitat at a higher rate than expected at random. Second, when released at larger distances from their habitat, they used a non-random, systematic, search strategy in which they flew in loops around the release point and returned periodically to it. Third, butterflies returned to a familiar habitat patch rather than a non-familiar one when given a choice. If dispersers actively orientate towards or search systematically for distant habitat, this may be problematic for existing metapopulation models, including models of the evolution of dispersal rates in metapopulations.     

Dispersal and the metapopulation paradigm in amphibian ecology and conservation: are all amphibian populations metapopulations?

Amphibians are frequently characterized as having limited dispersal abilities, strong site fidelity and spatially disjunct breeding habitat. As such, pond‐breeding species are often alleged to form metapopulations. Amphibian species worldwide appear to be suffering population level declines caused, at least in part, by the degradation and fragmentation of habitat and the intervening areas between habitat patches. If the simplification of amphibians occupying metapopulations is accurate, then a regionally based conservation strategy, informed by metapopulation theory, is a powerful tool to estimate the isolation and extinction risk of ponds or populations. However, to date no attempt to assess the class‐wide generalization of amphibian populations as metapopulations has been made. We reviewed the literature on amphibians as metapopulations (53 journal articles or theses) and amphibian dispersal (166 journal articles or theses for 53 anuran species and 37 salamander species) to evaluate whether the conditions for metapopulation structure had been tested, whether pond isolation was based only on the assumption of limited dispersal, and whether amphibian dispersal was uniformly limited. We found that in the majority of cases (74%) the assumptions of the metapopulation paradigm were not tested. Breeding patch isolation via limited dispersal and/or strong site fidelity was the most frequently implicated or tested metapopulation condition, however we found strong evidence that amphibian dispersal is not as uniformly limited as is often thought. The frequency distribution of maximum movements for anurans and salamanders was well described by an inverse power law. This relationship predicts that distances beneath 11–13 and 8–9 km, respectively, are in a range that they may receive one emigrating individual. Populations isolated by distances approaching this range are perhaps more likely to exhibit metapopulation structure than less isolated populations. Those studies that covered larger areas also tended to report longer maximum movement distances – a pattern with implications for the design of mark‐recapture studies. Caution should be exercised in the application of the metapopulation approach to amphibian population conservation. Some amphibian populations are structured as metapopulations – but not all.     

Consequences of variable larval dispersal pathways and resulting phenotypic mixtures to the dynamics of marine metapopulations

Larval dispersal can connect distant subpopulations, with important implications for marine population dynamics and persistence, biodiversity conservation and fisheries management. However, different dispersal pathways may affect the final phenotypes, and thus the performance and fitness of individuals that settle into subpopulations. Using otolith microchemical signatures that are indicative of ‘dispersive’ larvae (oceanic signatures) and ‘non-dispersive’ larvae (coastal signatures), we explore the population-level consequences of dispersal-induced variability in phenotypic mixtures for the common triplefin (a small reef fish). We evaluate lipid concentration and otolith microstructure and find that ‘non-dispersive’ larvae (i) have greater and less variable lipid reserves at settlement (and this variability attenuates at a slower rate), (ii) grow faster after settlement, and (iii) experience similar carry-over benefits of lipid reserves on post-settlement growth relative to ‘dispersive’ larvae. We then explore the consequences of phenotypic mixtures in a metapopulation model with two identical subpopulations replenished by variable contributions of ‘dispersive’ and ‘non-dispersive’ larvae and find that the resulting phenotypic mixtures can have profound effects on the size of the metapopulation. We show that, depending upon the patterns of connectivity, phenotypic mixtures can lead to larger metapopulations, suggesting dispersal-induced demographic heterogeneity may facilitate metapopulation persistence.     

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.     

How does temporal variation in habitat connectivity influence metapopulation dynamics?

Metapopulation persistence depends on connectivity between habitat patches. While emphasis has been placed on the spatial dynamics of connectivity, much less has been placed on its short‐term temporal dynamics. In many terrestrial and aquatic ecosystems, however, transient (short‐term) changes in connectivity occur as habitat patches are connected and disconnected due, for example, to climatic or hydrological variability. We evaluated the implications of transient connectivity using a network‐based metapopulation model and a series of scenarios representing temporal changes in connectivity. The transient loss of connectivity can influence metapopulation persistence, and more strongly autocorrelated temporal dynamics affect metapopulation persistence more severely. Given that many ecosystems experience short‐term and temporary loss of habitat connectivity, it is important that these dynamics are adequately represented in metapopulation models; failing to do so may yield overly optimistic‐estimates of metapopulation persistence in fragmented landscapes.     

Consequences of forest fragmentation for the dynamics of bird populations: conceptual issues and the evidence

This paper reviews the consequences of forest fragmentation for the dynamics of bird populations. Owing to high mobility and large home ranges, birds usually perceive fragmented forests in a finegrained manner, i.e. embrace several forest fragments in functional home ranges. On a regional scale, however, coarse-grained clusters of fine-grained fragments (hierarchical fragmentation may sub-divide bird populations into isolated demes, which enter a domain of metapopulation dynamics. Distinctions are made between pure distance-area or population-level effects and more indirect community-level effects due to changes in landscape composition. Distance-area effects, such as insularization and decreasing fragment size, directly prevent dispersal and reduce population size. Landscape effects, such as reduced fragment-matrix and interior-edge ratios, increase the pressure from surrounding predators, competitors, parasites and disease. In short, forest fragmentation can be viewed as a two-step process. Initially, fine-grained fragmentation triggers distance-area and landscape effects on a local scale, which in turn, results in a range retraction of a population to non-fragmented or less fragmented parts of a region. At a certain point, non-fragmented areas become so widely spaced out that regional distance-area effects come into operation, giving rise to metapopulation dynamics. Although few bird metapopulations have yet been documented, metapopulation dynamics probably is a common characteristic of bird populations confined to 'hierarchical' fragmented forests.     

Movements and source–sink dynamics of a Masai giraffe metapopulation

Spatial variation in habitat quality and anthropogenic factors, as well as social structure, can lead to spatially structured populations of animals. Demographic approaches can be used to improve our understanding of the dynamics of spatially structured populations and help identify subpopulations critical for the long-term persistence of regional metapopulations. We provide a regional metapopulation analysis to inform conservation management for Masai giraffes (Giraffa camelopardalis tippelskirchi) in five subpopulations defined by land management designations. We used data from an individual-based mark–recapture study to estimate subpopulation sizes, subpopulation growth rates, and movement probabilities among subpopulations. We assessed the source–sink structure of the study population by calculating source–sink statistics, and we created a female-based matrix metapopulation model composed of all subpopulations to examine how variation in demographic components of survival, reproduction, and movement affected metapopulation growth rate. Movement data indicated no subpopulation was completely isolated, but movement probabilities varied among subpopulations. Source–sink statistics and net flow of individuals indicated three subpopulations were sources, while two subpopulations were sinks. We found areas with higher wildlife protection efforts and fewer anthropogenic impacts were sources, and less-protected areas were identified as sinks. Our results highlight the importance of identifying source–sink dynamics among subpopulations for effective conservation planning and emphasize how protected areas can play an important role in sustaining metapopulations.     

Maintenance of the sexes and persistence of a clonal organism in spatially complex metapopulations

Clonal organisms persist at a range of population sex ratios, from equal numbers of males and females to single-sex systems. When intersexual competition is strong enough to drive one sex locally extinct, the maintenance of the sexes is facilitated by the semi-independent dynamics of populations within a metapopulation. These semi-independent dynamics are influenced by dispersal and recolonization rates, which are affected by the spatial arrangement of populations. To establish the quantitative relationship between spatially complex metapopulations and the maintenance of the sexes, we used a mathematical model of the liverwort Marchantia inflexa. This clonal organism is found in discrete patches on rocks and along the banks of streams, which form single-sex and two-sex metapopulations. In this system, asexual propagules mainly disperse short distances. Long-distance between-patch dispersal and recolonization mainly occurs via sexual propagules, which require both sexes to be present. Dispersal of these two types of propagules could interact with the spatial arrangement of populations to affect the maintenance of the sexes. With our mathematical model, we found that at intermediate distances between populations, metapopulations maintained both sexes, and the spatial arrangement of populations changed the threshold at which one sex was lost. On the other hand, when populations were close to one another, one sex was lost and the single-sex metapopulation persisted through dispersal of asexual propagules. When populations were far apart, one sex was lost, and the metapopulation either went extinct due to lack of recolonization by asexual propagules or persisted because clumped populations facilitated recolonization. These idealized spatial arrangements help clarify the effects of the spatial arrangement on the maintenance of the sexes and the persistence of metapopulations of clonal organisms, which can help explain geographic parthenogenesis and the distribution of asexual populations, the persistence of asexual species, and inform the conservation of clonal organisms.     

Two‐species metacommunity dynamics mediated by habitat preference

The Levins model is a simple and widely used metapopulation model that describes temporal changes in the regional abundance of a single species and has increasingly been applied to metacommunity contexts including multiple species. Although a fundamental assumption commonly made when using the model is that species randomly move between habitat patches, most organisms exhibit habitat preference in reality. A method of incorporating habitat preference (directed dispersal) into the Levins metapopulation model was developed in a previous study. In the current study, we extended the approach to explore two‐species metacommunity dynamics (i.e. competition and predation) mediated by habitat preference. Our results theoretically revealed that coexistence of competing metapopulations requires conspecific aggregation and heterospecific segregation whereas the conspecific segregation of prey and effective avoidance of unsuitable prey‐free patches are crucial for persistence of predator metapopulations. In addition, we qualitatively and quantitatively demonstrated the effect of habitat preference on the outcomes of interspecific interactions. The present study opens a new research avenue in metacommunity ecology in complex nature and contributes to improved landscape management for the conservation of species (e.g. territorial and group‐living animals) and biodiversity.     

Consequences of metapopulation collapse: comparison of genetic attributes between two Allegheny woodrat metapopulations

Disruptions in metapopulation connectivity due to demographic pressure can leave genetically isolated subpopulations susceptible to genetic drift, accumulation of deleterious alleles, and inbreeding depression. Such a scenario may be playing out within Allegheny woodrat (Neotoma magister) metapopulations as a series of synergistic extrinsic pressures have contributed to the rangewide decline of the species over the last 40 years. Our goal was to elucidate the effects of demographic collapse on metapopulation function by using 11 microsatellites markers to quantify differences in patterns of connectivity and genetic diversity between a demographically stable metapopulation and one in severe demographic decline. The demographically diminished metapopulation had lower levels of genetic diversity than the stable metapopulation at all levels evaluated (metapopulation-, subpopulation-, and individual-scales). In contrast to patterns of connectivity observed within the stable metapopulation, peripheral subpopulations in the diminished metapopulation had become completely isolated and were drifting toward genetic fixation, likely as a result of the extirpation of stepping-stone subpopulations. The declining genetic parameters observed within these isolated peripheral subpopulations suggest that inbreeding depression may be contributing significantly to their demographic decline. Allegheny woodrats readily express the genetic consequences of metapopulation decline due to the low effective population sizes of subpopulations and the species’ limited dispersal capacity. Differences in genetic parameters observed between demographically stable and diminished Allegheny woodrat metapopulations emphasize the risks posed to metapopulation function and associated genetic processes introduced with demographic decline.     

Observational evidence that maladaptive gene flow reduces patch occupancy in a wild insect metapopulation

Theory predicts that dispersal throughout metapopulations has a variety of consequences for the abundance and distribution of species. Immigration is predicted to increase abundance and habitat patch occupancy, but gene flow can have both positive and negative demographic consequences. Here, we address the eco‐evolutionary effects of dispersal in a wild metapopulation of the stick insect Timema cristinae, which exhibits variable degrees of local adaptation throughout a heterogeneous habitat patch network of two host‐plant species. To disentangle the ecological and evolutionary contributions of dispersal to habitat patch occupancy and abundance, we contrasted the effects of connectivity to populations inhabiting conspecific host plants and those inhabiting the alternate host plant. Both types of connectivity should increase patch occupancy and abundance through increased immigration and sharing of beneficial alleles through gene flow. However, connectivity to populations inhabiting the alternate host‐plant species may uniquely cause maladaptive gene flow that counters the positive demographic effects of immigration. Supporting these predictions, we find the relationship between patch occupancy and alternate‐host connectivity to be significantly smaller in slope than the relationship between patch occupancy and conspecific‐host connectivity. Our findings illustrate the ecological and evolutionary roles of dispersal in driving the distribution and abundance of species.     

Demography and genetics suggest reversal of dolphin source‐sink dynamics,with implications for conservation

The forecast for the viability of populations depends upon metapopulation dynamics: the combination of reproduction and mortality within populations, as well as dispersal between populations. This study focuses on an Indo‐Pacific bottlenose dolphin (Tursiops aduncus) population in coastal waters near Bunbury, Western Australia. Demographic modeling of this population suggested that recent reproductive output was not sufficient to offset mortality. Migrants from adjacent populations might make up this deficit, so that Bunbury would act as a “sink,” or net recipient population. We investigated historical dispersal in and out of Bunbury, using microsatellites and mitochondrial DNA of 193 dolphins across five study locations along the southwestern Australian coastline. Our results indicated limited gene flow between Bunbury and adjacent populations. The data also revealed a net‐dispersal from Bunbury to neighboring populations, with microsatellites showing that more than twice as many individuals per generation dispersed out of Bunbury than into Bunbury. Therefore, in historic times, Bunbury appears to have acted as a source population, supporting nearby populations. In combination with the prior finding that Bunbury is currently not producing surplus offspring to support adjacent populations, this potential reversal of source‐sink dynamics may have serious conservation implications for Bunbury and other populations nearby.