Colonization and extinction dynamics of an annual plant metapopulation in an urban environment

The metapopulation framework considers that the spatiotemporal distribution of organisms results from a balance between the colonization and extinction of populations in a suitable and discrete habitat network. Recent spatially realistic metapopulation models have allowed patch dynamics to be investigated in natural populations but such models have rarely been applied to plants. Using a simple urban fragmented population system in which favourable habitat can be easily mapped, we studied patch dynamics in the annual plant Crepis sancta (Asteraceae). Using stochastic patch occupancy models (SPOMs) and multi‐year occupancy data we dissected extinction and colonization patterns in our system. Overall, our data were consistent with two distinct metapopulation scenarios. A metapopulation (sensu stricto) dynamic in which colonization occurs over a short distance and extinction is lowered by nearby occupied patches (rescue effect) was found in a set of patches close to the city centre, while a propagule rain model in which colonization occurs from a large external population was most consistent with data from other networks. Overall, the study highlights the importance of external seed sources in urban patch dynamics. Our analysis emphasizes the fact that plant distributions are governed not only by habitat properties but also by the intrinsic properties of colonization and dispersal of species. The metapopulation approach provides a valuable tool for understanding how colonization and extinction shape occupancy patterns in highly fragmented plant populations. Finally, this study points to the potential utility of more complex plant metapopulation models than traditionally used for analysing ecological and evolutionary processes in natural metapopulations.     

Long‐term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia): survey methods,data management,and long‐term population trends

Long‐term observational studies conducted at large (regional) spatial scales contribute to better understanding of landscape effects on population and evolutionary dynamics, including the conditions that affect long‐term viability of species, but large‐scale studies are expensive and logistically challenging to keep running for a long time. Here, we describe the long‐term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia) that has been conducted since 1991 in a large network of 4000 habitat patches (dry meadows) within a study area of 50 by 70 km in the Åland Islands in Finland. We explain how the landscape structure has been described, including definition, delimitation, and mapping of the habitat patches; methods of field survey, including the logistics, cost, and reliability of the survey; and data management using the EarthCape biodiversity platform. We describe the long‐term metapopulation dynamics of the Glanville fritillary based on the survey. There has been no long‐term change in the overall size of the metapopulation, but the level of spatial synchrony and hence the amplitude of fluctuations in year‐to‐year metapopulation dynamics have increased over the years, possibly due to increasing frequency of exceptional weather conditions. We discuss the added value of large‐scale and long‐term population studies, but also emphasize the need to integrate more targeted experimental studies in the context of long‐term observational studies. For instance, in the case of the Glanville fritillary project, the long‐term study has produced an opportunity to sample individuals for experiments from local populations with a known demographic history. These studies have demonstrated striking differences in dispersal rate and other life‐history traits of individuals from newly established local populations (the offspring of colonizers) versus individuals from old, established local populations. The long‐term observational study has stimulated the development of metapopulation models and provided an opportunity to test model predictions. This combination of empirical studies and modeling has facilitated the study of key phenomena in spatial dynamics, such as extinction threshold and extinction debt.     

Metapopulation viability of an endangered holoparasitic plant in a dynamic landscape

By creating transient patch mosaics, disturbance can influence the dynamics of interacting populations in many ecosystems. In European heathland, traditional land use created such dynamic systems favourable for both early and later successional species. Little empirical evidence is, however, available on the impact of current management on metapopulations occurring in such landscapes. This paper looks at the metapopulation viability of the endangered holoparasite Cuscuta epithymum, a species that typically occurs in early successional stages of recently managed heathlands. We used both observational and experimental data from a 4‐yr study to parameterise a spatially explicit metapopulation model. This model explores the impact of demographic characteristics and spatiotemporal landscape patterns created by management events on metapopulation viability. Both occasional long‐distance dispersal and dormant seeds are shown to be critical for the long‐term survival of C. epithymum in a dynamic heathland landscape subjected to a fixed rotational mowing of 15 yr. A relatively high management frequency (<15 yr between two consecutive mowing events) appeared to be necessary to sustain a viable C. epithymum metapopulation. When there is a longer interval between management events, grazing can counterbalance the negative effects of vegetation succession. Our results indicate that small‐scale cyclical management events combined with extensive grazing are the most appropriate management strategy to maintain viable populations of C. epithymum instead of the current large‐scale management events. Our results further emphasise the importance of incorporating both spatiotemporal patch availability and key demographic characteristics, especially seed banks, for a realistic view of metapopulation dynamics in disturbed landscapes. This study clearly demonstrates the usefulness of metapopulation models to understand the impact of management events and to provide new ecological insights into processes acting at a landscape scale.     

Occupancy versus colonization–extinction models for projecting population trends at different spatial scales

Understanding spatiotemporal population trends and their drivers is a key aim in population ecology. We further need to be able to predict how the dynamics and sizes of populations are affected in the long term by changing landscapes and climate. However, predictions of future population trends are sensitive to a range of modeling assumptions. Deadwood‐dependent fungi are an excellent system for testing the performance of different predictive models of sessile species as these species have different rarity and spatial population dynamics, the populations are structured at different spatial scales, and they utilize distinct substrates. We tested how the projected large‐scale occupancies of species with differing landscape‐scale occupancies are affected over the coming century by different modeling assumptions. We compared projections based on occupancy models against colonization–extinction models, conducting the modeling at alternative spatial scales and using fine‐ or coarse‐resolution deadwood data. We also tested effects of key explanatory variables on species occurrence and colonization–extinction dynamics. The hierarchical Bayesian models applied were fitted to an extensive repeated survey of deadwood and fungi at 174 patches. We projected higher occurrence probabilities and more positive trends using the occupancy models compared to the colonization–extinction models, with greater difference for the species with lower occupancy, colonization rate, and colonization:extinction ratio than for the species with higher estimates of these statistics. The magnitude of future increase in occupancy depended strongly on the spatial modeling scale and resource resolution. We encourage using colonization–extinction models over occupancy models, modeling the process at the finest resource‐unit resolution that is utilizable by the species, and conducting projections for the same spatial scale and resource resolution at which the model fitting is conducted. Further, the models applied should include key variables driving the metapopulation dynamics, such as the availability of suitable resource units, habitat quality, and spatial connectivity.     

Metapopulation dynamics in a changing climate: Increasing spatial synchrony in weather conditions drives metapopulation synchrony of a butterfly inhabiting a fragmented landscape

Habitat fragmentation and climate change are both prominent manifestations of global change, but there is little knowledge on the specific mechanisms of how climate change may modify the effects of habitat fragmentation, for example, by altering dynamics of spatially structured populations. The long‐term viability of metapopulations is dependent on independent dynamics of local populations, because it mitigates fluctuations in the size of the metapopulation as a whole. Metapopulation viability will be compromised if climate change increases spatial synchrony in weather conditions associated with population growth rates. We studied a recently reported increase in metapopulation synchrony of the Glanville fritillary butterfly (Melitaea cinxia) in the Finnish archipelago, to see if it could be explained by an increase in synchrony of weather conditions. For this, we used 23 years of butterfly survey data together with monthly weather records for the same period. We first examined the associations between population growth rates within different regions of the metapopulation and weather conditions during different life‐history stages of the butterfly. We then examined the association between the trends in the synchrony of the weather conditions and the synchrony of the butterfly metapopulation dynamics. We found that precipitation from spring to late summer are associated with the M. cinxia per capita growth rate, with early summer conditions being most important. We further found that the increase in metapopulation synchrony is paralleled by an increase in the synchrony of weather conditions. Alternative explanations for spatial synchrony, such as increased dispersal or trophic interactions with a specialist parasitoid, did not show paralleled trends and are not supported. The climate driven increase in M. cinxia metapopulation synchrony suggests that climate change can increase extinction risk of spatially structured populations living in fragmented landscapes by altering their dynamics.     

Testing metapopulation dynamics using genetic, demographic and ecological data

The metapopulation concept is a cornerstone in the recent history of ecology and evolution. However, determining whether a natural system fits a metapopulation model is a complex issue. Extinction-colonization dynamics are indeed often difficult to quantify because species detectability is not always 100%, resulting in an imperfect record of extinctions. Here, we explore whether combining population genetics with demographic and ecological surveys can yield more realistic estimates of metapopulation dynamics. We apply this approach to the freshwater snail Drepanotrema depressissimum in a fragmented landscape of tropical ponds. In addition to studying correlations between genetic diversity and demographical or ecological characteristics, we undertake, for the first time, a detailed search for genetic signatures of extinction-recolonization events using temporal changes in allele frequencies within sites. Surprisingly, genetic data indicate that extinction is much rarer than suggested by demographic surveys. Consequently, this system is better described as a set of populations with different sizes and immigration rates than as a true metapopulation. We identify several cases of apparent extinction owing to nondetection of low-density populations, and of aestivating individuals in desiccated ponds. More generally, we observed a frequent mismatch between genetic and demographical/ecological information at small spatial and temporal scales. We discuss the causes of these discrepancies and show how these two types of data provide complementary information on population dynamics and history, especially when temporal genetic samples are available.     

Contrasting Patterns of Demography and Population Viability Among Gopher Tortoise Populations in Alabama

Population viability analyses are useful tools to predict abundance and extinction risk for imperiled species. In southeastern North America, the federally threatened gopher tortoise (Gopherus polyphemus) is a keystone species in the diverse and imperiled longleaf pine (Pinus palustris) ecosystem, and researchers have suggested that tortoise populations are declining and characterized by high extinction risk. We report results from a 30-year demographic study of gopher tortoises in southern Alabama (1991–2020), where 3 populations have been stable and 3 others have declined. To better understand the demographic vital rates associated with stable and declining tortoise populations, we used a multi-state hierarchical mark-recapture model to estimate sex- and stage-specific patterns of demographic vital rates at each population. We then built a predictive population model to project population dynamics and evaluate extinction risk in a population viability context. Population structure did not change significantly in stable populations, but juveniles became less abundant in declining populations over 30 years. Apparent survival varied by age, sex, and site; adults had higher survival than juveniles, but female survival was substantially lower in declining populations than in stable ones. Using simulations, we predicted that stable populations with high female survival would persist over the next 100 years but sites with lower female survival would decline, become male-biased, and be at high risk of extirpation. Stable populations were most sensitive to changes in apparent survival of adult females. Because local populations varied greatly in vital rates, our analysis improves upon previous demographic models for northern populations of gopher tortoises by accounting for population-level variation in demographic patterns and, counter to previous model predictions, suggests that small tortoise populations can persist when habitat is managed effectively. © 2021 The Wildlife Society.     

Demographic processes shaping genetic variation of the solitarious phase of the desert locust

Between plagues, the solitarious desert locust (Schistocerca gregaria) is generally thought to exist as small populations, which are particularly prone to extinction events in arid regions of Africa and Asia. Given the high genetic structuring observed in one geographical area (the Eritrean coast) by former authors, a metapopulation dynamics model involving repeated extinction and colonization events was favoured. In this study, we assessed the validity of a demographic scenario involving temporary populations of the solitarious phase of the desert locust by analysing large‐scale population genetic data. We scored 24 microsatellites in 23 solitarious population samples collected over most of the species range during remission. We found very little genetic structuring and little evidence of declining genetic diversity. A Bayesian clustering method distinguished four genetically differentiated units. Three groups were largely consistent with three population samples which had undergone recent bottleneck events. Nevertheless, the last genetically homogeneous unit included all individuals from the remaining 18 population samples and did not show evidence of demographic disequilibrium. An approximate Bayesian computation treatment indicated a large population size for this main genetic group, moderately reduced between plague and remission but still containing tens of thousands of individuals. Our results diverge from the hypothesis of a classical metapopulation dynamics model. They instead support the scenario in which large populations persist in the solitarious phase of the desert locust.     

Predator-prey interactions in a nonequilibrium context: the metapopulation approach to modeling "hide-and-seek" dynamics in a spatially explicit tri-trophic system

Predators and prey are usually heterogeneously distributed in space so that the ability of the predators to respond to the distribution of their prey may have a profound influence on the stability and persistence of a predator‐prey system. A special type of dynamics is “hide‐and‐seek” characterized by a high turnover rate of local populations of prey and predators, because once the predators have found a patch of prey they quickly overexploit it, whereupon the starving predators either should move to better places or die. Continued persistence of prey and predators thus hinges on a long‐term balance between local extinctions and founding of new subpopulations. The colonization rate depends on the rate of emigration from occupied patches and the likelihood of successfully arriving at a suitable new patch, while extinction rate depends on the local population dynamics. Since extinctions and colonizations are both discrete probabilistic events, these phenomena are most adequately modeled by means of a stochastic model. In order to demonstrate the qualitative differences between a deterministic and stochastic approach to population dynamics, a spatially explicit tritrophic predator‐prey model is developed in a deterministic and a stochastic version. The model is parameterized using data for the two‐spotted spider mite (Tetranychus urticae) and the phytoseiid mite predator Phytoseiulus persimilis inhabiting greenhouse cucumbers.
Simulations show that the deterministic and stochastic approaches yield different results. The deterministic version predicts that the populations will exhibit violent fluctuations, implying that the system is fundamentally unstable. In contrast, the stochastic version predicts that the two species will be able to coexist in spite of frequent local extinctions of both species, provided the system consists of a sufficiently large number of local populations. This finding is in agreement with experimental results. It is therefore concluded that demographic stochasticity in combination with dispersal is capable of producing and maintaining sufficient asynchrony between local populations to ensure long‐term regional (metapopulation) persistence.     

Asynchronization of local population dynamics and persistence of a metapopulation: a lesson from an endangered composite plant, Aster kantoensis

Fragmentation of a large habitat makes local populations less linked to others, and a whole population structure changes to a metapopulation. The smaller a local population is, the more strengthened extinction factors become. Then, frequent extinctions of local populations threaten persistence of the metapopulation unless recolonizations occur rapidly enough after local extinctions. Spatially structured models have been more widely used for predicting future population dynamics and for assessing the extinction risk of a metapopulation. In this article, we first review such spatially structured models that have been applied to conservation biology, focusing on effects of asynchronization among local population dynamics on persistence of the whole metapopulation. Second, we introduce our ongoing project on extinction risk assessment of an endangered composite biennial plant, Aster kantoensis, in the riverside habitat, based on a lattice model for describing its spatiotemporal population dynamics. The model predicted that the extinction risk of A. kantoensis depends on both the frequency of flood occurrence and the time to coverage of a local habitat by other competitively stronger perennials. Finally, we present a measure (Hassell and Pacala's CV ²) for quantifying the effect of asynchronization among local population dynamics on the persistence of a whole metapopulation in conservation ecology. Received: January 12, 2000 / Accepted: February 8, 2000     

Extinction, colonization, and species occupancy in tidepool fishes

Despite the increasing sophistication of ecological models with respect to the size and spatial arrangement of habitat, there is relatively little empirical documentation of how species dynamics change as a function of habitat size and the fraction of habitat occupied. In an assemblage of tidepool fishes, I used maximum-likelihood estimation to test whether models which included habitat size provided a better fit to empirical data on extinction and colonization probabilities than models that assumed constant probabilities over all habitats. I found species differences in how extinction and colonization probabilities scaled with habitat size (and hence local population size). However, there was little evidence for a relationship between extinction and colonization probabilities and the fraction of occupied tidepools, as assumed in simple metapopulation models. Instead, colonization and extinction were independent of the fraction of occupied tidepools, favoring a MacArthur-Wilson island-mainland model. When I incorporated declines in extinction probability with tidepool volume in a simple simulation model, I found that predicted occupancy could change greatly, especially when colonization was low. However, the predicted fraction of occupied patches in the simulation model changed little when I incorporated the range of values reported here for extinction and colonization and the rate at which they scale with habitat size. Quantifying extinction and colonization patterns of natural populations is fundamental to understanding how species are distributed spatially and whether metapopulation models of species occupancy provide explanatory power for field populations. Received: 14 March 1997 / Accepted: 21 September 1997     

The applicability of metapopulation theory to large mammals

Metapopulation theory has become a common framework in conservation biology and it is sometimes suggested that a metapopulation approach should be used for management of large mammals. However, it has also been suggested that metapopulation theory would not be applicable to species with long generations compared to those with short ones. In this paper, we review how and on what empirical ground metapopulation terminology has been applied to insects, small mammals and large mammals. The review showed that the metapopulation term sometimes was used for population networks which only fulfilled the broadest possible definition of a metapopulation, i.e. they were subpopulations connected by migrating individuals. We argue that the metapopulation concept should be reserved for networks that also show some kind of metapopulation dynamics. Otherwise it applies to almost all populations and loses its substance. We found much empirical support for metapopulation dynamics in both insects and small mammals, but not in large mammals. A possible reason is the methods used to confirm the existence of metapopulation dynamics. For insects and small mammals, the common approach is to study population turnover through patch occupancy data. Such data is difficult to obtain for large mammals, since longer temporal scales need to be covered to record extinctions and colonizations. Still, many populations of large mammals are exposed to habitat fragmentation and the resulting subpopulations sometimes have high risks of extinction. If there is migration between the subpopulations, the metapopulation framework could provide valuable information on their population dynamics. We suggest that a metapopulation approach can be interesting for populations of large mammals, when there are discrete breeding subpopulations and when these subpopulations have different growth rates and demographic fates. Thus, a comparison of the subpopulations’ demographic fates, rather than subpopulation turnover, can be a feasible alternative for studies of metapopulation dynamics in large mammals.     

Sperm storage reduces the strength of the mate‐finding Allee effect

Mate searching is a key component of sexual reproduction that can have important implications for population viability, especially for the mate‐finding Allee effect. Interannual sperm storage by females may be an adaptation that potentially attenuates mate limitation, but the demographic consequences of this functional trait have not been studied. Our goal is to assess the effect of female sperm storage durability on the strength of the mate‐finding Allee effect and the viability of populations subject to low population density and habitat alteration. We used an individual‐based simulation model that incorporates realistic representations of the demographic and spatial processes of our model species, the spur‐thighed tortoise (Testudo graeca). This allowed for a detailed assessment of reproductive rates, population growth rates, and extinction probabilities. We also studied the relationship between the number of reproductive males and the reproductive rates for scenarios combining different levels of sperm storage durability, initial population density, and landscape alteration. Our results showed that simulated populations parameterized with the field‐observed demographic rates collapsed for short sperm storage durability, but were viable for a durability of one year or longer. In contrast, the simulated populations with a low initial density were only viable in human‐altered landscapes for sperm storage durability of 4 years. We find that sperm storage is an effective mechanism that can reduce the strength of the mate‐finding Allee effect and contribute to the persistence of low‐density populations. Our study highlights the key role of sperm storage in the dynamics of species with limited movement ability to facilitate reproduction in patchy landscapes or during population expansion. This study represents the first quantification of the effect of sperm storage durability on population dynamics in different landscapes and population scenarios.     

Loss of dendritic connectivity in southern California's urban riverscape facilitates decline of an endemic freshwater fish

Life history adaptations and spatial configuration of metapopulation networks allow certain species to persist in extreme fluctuating environments, yet long‐term stability within these systems relies on the maintenance of linkage habitat. Degradation of such linkages in urban riverscapes can disrupt this dynamic in aquatic species, leading to increased extinction debt in local populations experiencing environment‐related demographic flux. We used microsatellites and mtDNA to examine the effects of collapsed network structure in the endemic Santa Ana sucker Catostomus santaanae of southern California, a threatened species affected by natural flood–drought cycles, “boom‐and‐bust” demography, hybridization and presumed artificial transplantation. Our results show a predominance of drift‐mediated processes in shaping population structure and that reverse mechanisms for counterbalancing the genetic effects of these phenomena have dissipated with the collapse of dendritic connectivity. We use approximate Bayesian models to support two cases of artificial transplantation and provide evidence that one of the invaded systems better represents the historic processes that maintained genetic variation within watersheds than any remaining drainages where C. santaanae is considered native. We further show that a stable dry gap in the northern range is preventing genetic dilution of pure C. santaanae persisting upstream of a hybrid assemblage involving a non‐native sucker and that local accumulation of genetic variation in the same drainage is influenced by position within the network. This work has important implications for declining species that have historically relied on dendritic metapopulation networks to maintain source–sink dynamics in phasic environments, but no longer possess this capacity in urban‐converted landscapes.     

Regional and local scale metapopulation dynamics in the interaction between Callosobruchus maculatus and Anisopteromalus calandrae

Habitat structure increases the persistence of many extinction‐prone resource–consumer interactions. Metapopulation theory is one of the leading approaches currently used to explain why local, ephemeral populations persist at a regional scale. Central to the metapopulation concept is the amount of dispersal occurring between patches, too much or too little can result in regional extinction. In this study, the role of dispersal on the metapopulation dynamics of an over‐exploitative host–parasitoid interaction is assessed. In the absence of the parasitoid the highly vagile bruchid, Callosobruchus maculatus, can maintain a similar population size regardless of the permeability of the inter‐patch matrix and exhibits strong negative density‐dependence. After the introduction of the parasitoid the size of the bruchid population decreases with a corresponding increase in the occurrence of empty patches. In this case, limiting the dispersal of both species decouples the interaction to a greater extent and results in larger regional bruchid populations. Given the disparity between the dispersal rates of the two species, it is proposed that the more dispersive host benefits from the reduction in landscape permeability by increasing the opportunity to colonise empty patches and rescue extinction prone populations. Associated with the introduction of the parasitoid is a shift in the strength of density‐dependence as the population moves from bottom–up towards top–down regulation. The importance of local and regional scale measurements is apparent when the role of individual patches on regional dynamics is considered. By only taking regional dynamics into account the importance of dispersal regime on local dynamics is overlooked. Similarly, when local dynamics were examined, patches were found to have different influences on regional dynamics depending on dispersal regime and patch location.     

Assessing the role of off‐take and source–sink dynamics in the extinction of the amphidromous New Zealand grayling (Prototroctes oxyrhynchus)

相似文献   

Large‐scale predator control increases population viability of a rare New Zealand riverine duck

The introduction of mammalian predators to oceanic islands has led to dramatic declines in the abundance of many native species. Conservation management of these species often relies on low‐cost predator control techniques that can be implemented over large scales. Assessing the effectiveness of such management techniques is difficult, but using population viability analyses (PVA), which identify the population growth rate (λ) and extinction risk of threatened species, may offer a solution. PVA provide the opportunity to compare the relative effectiveness of various management options and can identify knowledge gaps to prioritize research efforts. We used PVA to assess the population viability of whio (Hymenolaimus malacorhynchos), a rare riverine duck endemic to New Zealand. Current populations are threatened by introduced mammalian predators and are rapidly declining in both distribution and abundance. Whio conservation management is dominated by large‐scale, low‐intensity predator control, targeting introduced stoats (Mustela erminea). There is evidence that such control increases whio productivity but it is unknown if this increase is sufficient for long‐term population persistence. We undertook a stochastic PVA to assess the viability of whio populations under different management scenarios using data obtained from a 6‐year study of whio demographic responses to predator control. Populations with no predator control and low productivity will rapidly decline to extinction. Increasing productivity through predator control increased population viability but populations still showed a declining trajectory. A perturbation analysis showed that the growth rate of whio populations was largely driven by adult survival. Therefore, future research should target obtaining more robust estimates of adult survival, particularly how it is affected by predator control. Overall, our analysis indicated that large‐scale predator control increases the short‐term viability of whio populations but is insufficient for long‐term population persistence.     

植物种群生存力分析研究进展

对十多年来国外植物PVA的研究进行了综合评述;具体分析了影响植物种群生存力的各种随机性因子及确定性因子;总结了植物PVA研究的方法步骤及采用的模拟模型;探讨了植物PVA的难点,PVA对管理措施的评价效果;并提出对今后植物PVA的研究展望,认为PVA是研究濒危植物种群灭绝及评价管理或保护措施的有力工具;发展描述复杂种间关系的多种种的PVA模型以及包含多个影响因素的PVA应用模型是未来植物PVA的研究方向。     

局域种群的Allee效应和集合种群的同步性

从包含Allee效应的局域种群出发,建立了耦合映像格子模型,即集合种群模型.通过分析和计算机模拟表明:(1)当局域种群受到Allee效应强度较大时,集合种群同步灭绝;(2)而当Allee效应强度相对较弱时,通过稳定局域种群动态(减少混沌)使得集合种群发生同步波动,而这种同步波动能够增加集合种群的灭绝风险;(3)斑块间的连接程度对集合种群同步波动的发生有很大的影响,适当的破碎化有利于集合种群的续存.全局迁移和Allee效应结合起来增加了集合种群同步波动的可能,从而增加集合种群的灭绝风险.这些结果对理解同步性的机理、利用同步机理来制定物种保护策略和害虫防治都有重要的意义.     

Stochastic dynamics in exotic and native blowflies: an analysis combining laboratory experiments and a two-patch metapopulation model

In this study we explored the stochastic population dynamics of three exotic blowfly species, Chrysomya albiceps, Chrysomya megacephala and Chrysomya putoria, and two native species, Cochliomyia macellaria and Lucilia eximia, by combining a density-dependent growth model with a two-patch metapopulation model. Stochastic fecundity, survival and migration were investigated by permitting random variations between predetermined demographic boundary values based on experimental data. Lucilia eximia and Chrysomya albiceps were the species most susceptible to the risk of local extinction. Cochliomyia macellaria, C. megacephala and C. putoria exhibited lower risks of extinction when compared to the other species. The simultaneous analysis of stochastic fecundity and survival revealed an increase in the extinction risk for all species. When stochastic fecundity, survival and migration were simulated together, the coupled populations were synchronized in the five species. These results are discussed, emphasizing biological invasion and interspecific interaction dynamics.