Extinction and colonization processes in subpopulations of five neotropical small mammal species

We report the extinction and colonization rates of five sympatric small mammal species at a semiarid locality in north central Chile. We provide information based on 6 years of monitoring on how colonization and extinction rates change according to landscape features (slope aspect) and on their relationship to populations size, population variability, and body size. We found that: (1) for all species in the assemblage, extinction rates of subpopulations from equatorial-facing slopes were significantly lower than those in polar-facing slopes, (2) population size was the most important factor determining extinction rates, (3) colonization rates did not vary between slopes, and were affected by population size only in equatorial-facing slopes, and (4) most species had higher extinction than colonization rates. Persistence of the metapopulation system for all five small mammal species appears to be fueled by repeated colonization events.     

Comparison of reintroduction and enhancement effects on metapopulation viability

Metapopulation viability depends upon a balance of extinction and colonization of local habitats by a species. Mechanisms that can affect this balance include physical characteristics related to natural processes (e.g. succession) as well as anthropogenic actions. Plant restorations can help to produce favorable metapopulation dynamics and consequently increase viability; however, to date no studies confirm this is true. Population viability analysis (PVA) allows for the use of empirical data to generate theoretical future projections in the form of median time to extinction and probability of extinction. In turn, PVAs can inform and aid the development of conservation, recovery, and management plans. Pitcher's thistle (Cirsium pitcheri) is a dune endemic that exhibited metapopulation dynamics. We projected viability of three natural and two restored populations with demographic data spanning 15–23 years to determine the degree the addition of reintroduced population affects metapopulation viability. The models were validated by comparing observed and projected abundances and adjusting parameters associated with demographic and environmental stochasticity to improve model performance. Our chosen model correctly predicted yearly population abundance for 60% of the population‐years. Using that model, 50‐year projections showed that the addition of reintroductions increases metapopulation viability. The reintroduction that simulated population performance in early‐successional habitats had the maximum benefit. In situ enhancements of existing populations proved to be equally effective. This study shows that restorations can facilitate and improve metapopulation viability of species dependent on metapopulation dynamics for survival with long‐term persistence of C. pitcheri in Indiana likely to depend on continued active management.     

Long-term demographic surveys reveal a consistent relationship between average occupancy and abundance within local populations of a butterfly metapopulation

Species distribution models are the tool of choice for large-scale population monitoring, environmental association studies and predictions of range shifts under future environmental conditions. Available data and familiarity of the tools rather than the underlying population dynamics often dictate the choice of specific method – especially for the case of presence–absence data. Yet, for predictive purposes, the relationship between occupancy and abundance embodied in the models should reflect the actual population dynamics of the modelled species. To understand the relationship of occupancy and abundance in a heterogeneous landscape at the scale of local populations, we built a spatio-temporal regression model of populations of the Glanville fritillary butterfly Melitaea cinxia in a Baltic Sea archipelago. Our data comprised nineteen years of habitat surveys and snapshot data of land use in the region. We used variance partitioning to quantify relative contributions of land use, habitat quality and metapopulation covariates. The model revealed a consistent and positive, but noisy relationship between average occupancy and mean abundance in local populations. Patterns of abundance were highly variable across years, with large uncorrelated random variation and strong local population stochasticity. In contrast, the spatio-temporal random effect, habitat quality, population connectivity and patch size explained variation in occupancy, vindicating metapopulation theory as the basis for modelling occupancy patterns in fragmented landscapes. Previous abundance was an important predictor in the occupancy model, which points to a spillover of abundance into occupancy dynamics. While occupancy models can successfully model large-scale population structure and average occupancy, extinction probability estimates for local populations derived from occupancy-only models are overconfident, as extinction risk is dependent on actual, not average, abundance.     

Daily age determination and growth rates of freshwater fish throughout a regulated lotic system of the Murray‐Darling Basin Australia

To facilitate future research in freshwater fish recruitment response to environmental flow delivery, size‐at‐age and growth models are presented for eight fish species occurring in south‐eastern Australia; three small‐bodied species and the juvenile 0+ age classes of five large‐bodied species. Otolith increments were used to estimate the daily age of golden perch Macquaria ambigua, bony bream Nematalosa erebi, common carp Cyprinus carpio; Murray cod Maccullochella peelii, freshwater eel‐tailed catfish Tandanus tandanus, Australian smelt Retropinna semoni, un‐specked hardyhead Craterocephalus stercusmuscarum fulvus and Murray–Darling rainbowfish Melanotaenia fluviatilis. Linear growth models provided the best fit for length‐at‐age data of juvenile 0+ age large‐bodied species; whereas von Bertalanffy growth functions provided the best fit to length‐at‐age data of small‐bodied species. The results provide novel baseline data for future research in this area.     

Metapopulation dynamics: Does it help to have more of the same?

With the interest in conservation biology shifting towards processes from patterns, and to populations from communities, the theory of metapopulation dynamics is replacing the equilibrium theory of island biogeography as the population ecology paradigm in conservation biology. The simplest models of metapopulation dynamics make predictions about the effects of habitat fragmentation - size and isolation of habitat patches - on metapopulation persistence. The simple models may be enriched by considerations of the effects of demographic and environmental stochasticity on the size and extinction probability of local populations. Environmental stochasticity affects populations at two levels: it makes local extinctions more probable, and it also decreases metapopulation persistence time by increasing the correlation of extinction events across populations. Some controversy has arisen over the significance of correlated extinctions, and how they may affect the optimal subdivision of metapopulations to maximize their persistence time.     

Trophic interactions affect the population dynamics and risk of extinction of basal species in food webs

This paper addresses effects of trophic complexity on basal species, in a Lotka–Volterra model with stochasticity. We use simple food web modules, with three trophic levels, and expose every species to random environmental stochasticity and analyze (1) the effect of the position of strong trophic interactions on temporal fluctuations in basal species’ abundances and (2) the relationship between fluctuation patterns and extinction risk. First, the numerical simulations showed that basal species do not simply track the environment, i.e. species dynamics do not simply mirror the characteristics of the applied environmental stochasticity. Second, the extinction risk of species was related to the fluctuation patterns of the species.More specifically, we show (i) that despite being forced by random stochasticity without temporal autocorrelation (i.e. white noise), there is significant temporal autocorrelation in the time series of all basal species’ abundances (i.e. the spectra of basal species are red-shifted), (ii) the degree of temporal autocorrelation in basal species time series is affected by food web structure and (iii) the degree of temporal autocorrelation tend to be correlated to the extinction risks of basal species.Our results emphasize the role of food web structure and species interactions in modifying the response of species to environmental variability. To shed some light on the mechanisms we compare the observed pattern in abundances of basal species with analytically predicted patterns and show that the change in the predicted pattern due to the addition of strong trophic interactions is correlated to the extinction risk of the basal species. We conclude that much remain to be understood about the mechanisms behind the interaction among environmental variability, species interactions, population dynamics and vulnerability before we quantitatively can predict, for example, effects of climate change on species and ecological communities. Here, however, we point out a new possible approach for identifying species that are vulnerable to environmental stochasticity by checking the degree of temporal autocorrelation in the time series of species. Increased autocorrelation in population fluctuations can be an indication of increased extinction risk.     

The equilibrium assumption in estimating the parameters of metapopulation models

1.  The construction of a predictive metapopulation model includes three steps: the choice of factors affecting metapopulation dynamics, the choice of model structure, and finally parameter estimation and model testing.
2.  Unless the assumption is made that the metapopulation is at stochastic quasi-equilibrium and unless the method of parameter estimation of model parameters uses that assumption, estimates from a limited amount of data will usually predict a trend in metapopulation size.
3.  This implicit estimation of a trend occurs because extinction-colonization stochasticity, possibly amplified by regional stochasticity, leads to unequal numbers of observed extinction and colonization events during a short study period.
4.  Metapopulation models, such as those based on the logistic regression model, that rely on observed population turnover events in parameter estimation are sensitive to the implicit estimation of a trend.
5.  A new parameter estimation method, based on Monte Carlo inference for statistically implicit models, allows an explicit decision about whether metapopulation quasi-stability is assumed or not.
6. Our confidence in metapopulation model parameter estimates that have been produced from only a few years of data is decreased by the need to know before parameter estimation whether the metapopulation is in quasi-stable state or not.
7. The choice of whether metapopulation stability is assumed or not in parameter estimation should be done consciously. Typical data sets cover only a few years and rarely allow a statistical test of a possible trend. While making the decision about stability one should consider any information about the landscape history and species and metapopulation characteristics.     

甘南高寒草甸坡向梯度对植物群落功能多样性的影响

植物群落功能多样性对于维持生态系统功能具有十分重要的意义。该研究以青藏高原高寒草甸植物群落为研究对象,运用R软件程序包(FD、Vegan)与单因素方差分析法,分析高寒草甸不同坡向植物群落功能多样性的变化及其与环境因子的关系,以揭示微生境梯度上植物对环境的适应策略以及群落种间功能特征及生态系统内部结构和功能的变化机制。结果显示：(1)青藏高原高寒草甸的北坡和西北坡的物种丰富度、植物株高、比叶面积、叶片有机碳、叶片氮以及叶片磷显著高于其他坡向,且坡向间差异显著(P<0.05)。(2)在北坡到南坡的变化中,功能丰富度差异比较明显(P<0.05),株高、比叶面积、叶片有机碳、叶片氮及叶片磷等功能丰富度均呈递减趋势。(3)北坡的株高功能均匀度、西北坡比叶面积、叶片有机碳、叶片氮及叶片磷等功能均匀度均显著高于南坡,且多元性状功能均匀度在坡向间差异显著(P<0.05)。(4)北坡的比叶面积、叶片氮功能离散度及西北坡的株高、叶片有机碳、叶片磷等功能离散度均高于南坡,且坡向间差异显著(P<0.05)。(5)植物群落功能多样性与土壤含水量、土壤有机碳、土壤全氮、土壤全磷呈显著正相关关系(P<0.05),与坡度、土壤pH、土温、照度呈显著负相关关系(P<0.05)。     

Life-history models of extinction: a test with island spiders

This study analyzes extinction patterns for two species of orb spiders monitored annually on 77 islands over a continuous 20-yr period. One species, Argiope argentata, has large populations sometimes crashing quickly to extinction and a much weaker relation of extinction likelihood to population size than does the other species, Metepeira datona. Demographic models were built for both species and matched against observations. Differences between the species in life-history traits-estimated with measurements from the field-together with incorporation of demographic stochasticity, a population ceiling, and environmental stochasticity, were necessary to fit the observed extinction curves. As predicted from life-history patterns, long-term population growth rates (and hence predicted extinction probabilities) are relatively very sensitive to values of juvenile survivorship. Models are also sensitive to variation in the population ceiling and environmental noise, which tend to act in a complementary manner. A simple model with no age structure was able to fit the data on large initial population sizes but not on small initial population sizes, showing that life cycle characteristics interact with the various sources of stochasticity and hence have to be taken into account to produce a precise model of the extinction process.     

Species–area relationships in Mediterranean-climate plant communities

Aim To determine the best‐fit model of species–area relationships for Mediterranean‐type plant communities and evaluate how community structure affects these species–area models. Location Data were collected from California shrublands and woodlands and compared with literature reports for other Mediterranean‐climate regions. Methods The number of species was recorded from 1, 100 and 1000 m² nested plots. Best fit to the power model or exponential model was determined by comparing adjusted r² values from the least squares regression, pattern of residuals, homoscedasticity across scales, and semi‐log slopes at 1–100 m² and 100–1000 m². Dominance–diversity curves were tested for fit to the lognormal model, MacArthur's broken stick model, and the geometric and harmonic series. Results Early successional Western Australia and California shrublands represented the extremes and provide an interesting contrast as the exponential model was the best fit for the former, and the power model for the latter, despite similar total species richness. We hypothesize that structural differences in these communities account for the different species–area curves and are tied to patterns of dominance, equitability and life form distribution. Dominance–diversity relationships for Western Australian heathlands exhibited a close fit to MacArthur's broken stick model, indicating more equitable distribution of species. In contrast, Californian shrublands, both postfire and mature stands, were best fit by the geometric model indicating strong dominance and many minor subordinate species. These regions differ in life form distribution, with annuals being a major component of diversity in early successional Californian shrublands although they are largely lacking in mature stands. Both young and old Australian heathlands are dominated by perennials, and annuals are largely absent. Inherent in all of these ecosystems is cyclical disequilibrium caused by periodic fires. The potential for community reassembly is greater in Californian shrublands where only a quarter of the flora resprout, whereas three quarters resprout in Australian heathlands. Other Californian vegetation types sampled include coniferous forests, oak savannas and desert scrub, and demonstrate that different community structures may lead to a similar species–area relationship. Dominance–diversity relationships for coniferous forests closely follow a geometric model whereas associated oak savannas show a close fit to the lognormal model. However, for both communities, species–area curves fit a power model. The primary driver appears to be the presence of annuals. Desert scrub communities illustrate dramatic changes in both species diversity and dominance–diversity relationships in high and low rainfall years, because of the disappearance of annuals in drought years. Main conclusions Species–area curves for immature shrublands in California and the majority of Mediterranean plant communities fit a power function model. Exceptions that fit the exponential model are not because of sampling error or scaling effects, rather structural differences in these communities provide plausible explanations. The exponential species–area model may arise in more than one way. In the highly diverse Australian heathlands it results from a rapid increase in species richness at small scales. In mature California shrublands it results from very depauperate richness at the community scale. In both instances the exponential model is tied to a preponderance of perennials and paucity of annuals. For communities fit by a power model, coefficients z and log c exhibit a number of significant correlations with other diversity parameters, suggesting that they have some predictive value in ecological communities.     

Biodiversity and persistence of ecological communities in variable environments

Recent analyses of climate data indicate that the intensity and frequency of different weather extremes have increased. Such increased environmental variability may lead to increased species extinction rates and hence have important consequences for the long-term persistence of ecological communities. Here we use model communities in order to investigate the relationship between species richness and community persistence in a fluctuating environment. We model two scenarios: (1) correlated species responses to environmental fluctuations and (2) uncorrelated (independent) species responses. We quantify the risk and extent of species extinctions using the so-called community viability analysis. It is shown that species-rich communities are more sensitive to environmental stochasticity than species-poor communities. Specifically, per species risk of extinction is higher in species-rich communities than in species-poor ones. Moreover, for a given species richness, communities with uncorrelated species responses to environmental variation run a considerable higher risk of losing a fixed proportion of species compared with communities with correlated species responses. We discuss the compatibility of these results with the ecological insurance hypothesis.     

Relative importance of host plant patch geometry and habitat quality on the patterns of occupancy,extinction and density of the monophagous butterfly <Emphasis Type="Italic">Iolana iolas</Emphasis>

Habitat fragmentation is a major cause of species rarity and decline because it increases local population extinctions and reduces recolonisation rates of remnant patches. Although two major patch characteristics (area and connectivity) have been used to predict distribution patterns in fragmented landscapes, other factors can affect the occurrence of a species as well as the probability of it becoming extinct. In this paper, we study the spatial structure and dynamics of the butterfly Iolana iolas in a 75-patch network of its host plant (Colutea hispanica) to determine the relative importance of patch area, connectivity and habitat quality characteristics on occupancy, extinction and density over the period 2003–2006. Occupancy in 2003, incidence (proportion of years occupied) and probability of extinction were mostly affected by patch area. Smaller patches were less likely to be occupied because they had a higher probability of extinction, partly due to environmental stochasticity. The density of I. iolas was negatively related to patch area in all study years. Only in 2004 was the density of I. iolas positively influenced by fruit production per plant. Our results suggest that for I. iolas, and probably for other specialist butterflies with clearly delimited resource requirements, metapopulation dynamics can be satisfactorily predicted using only geometric variables because most habitat characteristics are subsumed in patch area. However, this hypothesis should be subject to further testing under diverse environmental conditions to evaluate the extent of its generalisation.     

Small island biogeography in the Gulf of California: lizards, the subsidized island biogeography hypothesis, and the small island effect

Aim We used insular lizard communities to test the predictions of two hypotheses that attempt to explain patterns of species richness on small islands. We first address the subsidized island biogeography (SIB) hypothesis, which predicts that spatial subsidies may cause insular species richness to deviate from species–area predictions, especially on small islands. Next, we examine the small island effect (SIE), which suggests small islands may not fit the traditional log‐linear species–area curve. Location Islands with arthropodivorous lizard communities throughout the Gulf of California. Methods To evaluate the SIB hypothesis, we first identified subsidized and unsubsidized islands based on surrogate measures of allochthonous productivity (i.e. island size and bird presence). Subsequently, we created species–area curves from previously published lizard species richness and island area data. We used the residuals and slopes from these analyses to compare species richness on subsidized and unsubsidized islands. To test for an SIE, we used breakpoint regression to model the relationship between lizard species richness and island area. We compared results from this model to results from the log‐linear regression model. Results Subsidized islands had a lower slope than unsubsidized islands, and the difference between these groups was significant when small islands were defined as < 1 km². In addition to comparing slopes, we tested for differences in the magnitude of the residuals (from the species–area regression of all islands) for subsidized vs. unsubsidized islands. We found no significant patterns in the residual values for small vs. large islands, or between islands with and without seabirds. The SIE was found to be a slightly better predictor of lizard species richness than the traditional log‐linear model. Main conclusions Predictions of the SIB hypothesis were partially supported by the data. The absence of a significant SIE may be a result of spatial subsidies as explained by the SIB hypothesis and data presented here. We conclude by suggesting potential scenarios to test for interactions between these two small island hypotheses. Future studies considering factors affecting species richness should examine the possible role of spatial subsidies, an SIE, or a synergistic effect of the two in data sets with small islands.     

Dances with fire: Tracking metapopulation dynamics of<Emphasis Type="Italic">polygonella basiramia</Emphasis> in Florida scrub (USA)

The regional persistence of species subject to local population colonization and extinction necessarily depends on how landscape features and disturbance affect metapopulation dynamics. Here, we characterize the metapopulation structure and short-term dynamics ofPolygonella basiramia. This rare, short-lived perennial herb is endemic to Florida scrublands and lacks a seed bank. Fires create the open sand gaps within a shrub matrix that support this species but also kill established plants. Thus, persistence depends on frequent colonization of unoccupied gaps. We are monitoring population dynamics within and among 1204 gaps distributed among 19 shrub patches. Considerable subpopulation turnover is evident at the gap level with rates of gap extinction exceeding rates of colonization in the first year. Whether declines in overall abundance continue is likely to depend on patterns of disturbance and regional stochasticity in this dynamic landscape.Polygonella is more likely to occupy larger and less isolated gaps, demonstrating that landscape features and disturbance strongly affect metapopulation dynamics. BecausePolygonella basiramia displays characteristics, occupancy patterns, and turnover dynamics consistent with metapopulation theory, it represents a model system for studying plant metapopulations.     

Single-species metapopulation dynamics: concepts, models and observations

This paper outlines a conceptual and theoretical framework for single-species metapopulation dynamics based on the Levins model and its variants. The significance of the following factors to metapopulation dynamics are explored: evolutionary changes in colonization ability; habitat patch size and isolation; compensatory effects between colonization and extinction rates; the effect of immigration on local dynamics (the rescue effect); and heterogeneity among habitat patches. The rescue effect may lead to alternative stable equilibria in metapopulation dynamics. Heterogeneity among habitat patches may give rise to a bimodal equilibrium distribution of the fraction of patches occupied in an assemblage of species (the core-satellite distribution). A new model of incidence functions is described, which allows one to estimate species' colonization and extinction rates on islands colonized from mainland. Four distinct kinds of stochasticity affecting metapopulation dynamics are discussed with examples. The concluding section describes four possible scenarios of metapopulation extinction.     

Correlated extinctions, colonizations and population fluctuations in a highly connected ringlet butterfly metapopulation

 The persistence of metapopulations is likely to be highly dependent on whether population dynamics are correlated among habitat patches as a result of migration between patches and spatially-correlated environmental stochasticity (weather effects). We examined whether population dynamics of the ringlet butterfly, Aphantopus hyperantus, were synchronous in an area of approximately 0.5 km², with respect to extinction, colonization and population fluctuations. Monks Wood Butterfly Monitoring Scheme transect count data from 1973 to 1995, revealed (A) a major environmental perturbation, the drought of 1976, which caused synchronized extinctions of A. hyperantus in subsequent years, (B) synchronized recolonization in years following the large number of apparent extinctions, and (C) population changes by A. hyperantus were highly correlated in many of the 14 sections of the transect, presumably reflecting similar responses to environmental stochasticity, and the exchange of individuals among sections. However, extinction and population synchrony depended on habitat type. Following the 1976 drought, A. hyperantus apparently became extinct from the most open and most shady habitats it occupied, with some persistence in habitats of intermediate shading, thus showing retraction to core populations in central parts of an environmental gradient, albeit with an average shift to relatively open habitat. Populations at extreme ends of the environmental gradient occupied by A. hyperantus fluctuated least synchronously, suggesting a potential buffering effect of habitat heterogeneity, but this was not crucial to survival after the 1976 drought. Thus, not all habitats are equally important to persistence. Correlated temporal dynamics, variation in habitat quality and the interaction between habitat quality and temporal environmental stochasticity are important determinants of metapopulation persistence and should be incorporated in metapopulation models. Received: 26 April 1996 / Accepted: 17 July 1996     

Extinctions in competitive communities forced by coloured environmental variation

Understanding the relationships between environmental fluctuations, population dynamics and species interactions in natural communities is of vital theoretical and practical importance. This knowledge is essential in assessing extinction risks in communities that are, for example, pressed by changing environmental conditions and increasing exploitation. We developed a model of density dependent population renewal, in a Lotka–Volterra competitive community context, to explore the significance of interspecific interactions, demographic stochasticity, population growth rate and species abundance on extinction risk in populations under various autocorrelation (colour) regimes of environmental forcing. These factors were evaluated in two cases, where either a single species or the whole community was affected by the external forcing. Species' susceptibility to environmental noise with different autocorrelation structure depended markedly on population dynamics, species' position in the abundance hierarchy and how similarly community members responded to external forcing. We also found interactions between demographic stochasticity and environmental noise leading to a reversal in extinction probabilities from under- to overcompensatory dynamics. We compare our results with studies of single species populations and contrast possible mechanisms leading to extinctions. Our findings indicate that abundance rank, the form of population dynamics, and the colour of environmental variation interact in affecting species extinction risk. These interactions are further modified by interspecific interactions within competitive communities as the interactions filter and modulate the environmental noise.     

Metapopulation shift and survival of woodland birds under climate change: will species be able to track?

Climate change has been widely recognized as a key factor driving changes in species distributions. In this study we use a metapopulation model, with a window of suitable climate moving polewards, to explore population shifts and survival of woodland birds under different climate change scenarios and landscape configurations. Extinction vulnerability and expansion ability are predicted for the middle spotted woodpecker Dendrocopus medius and two alternative r‐K strategies under west European climate change scenarios of 1, 2 and 4°C temperature increase per century, corresponding to isotemperature velocities of ca 2, 4 and 8 km yr^?1. The simulated northward expansion of the bird's distribution is typically in the range of only 0–3 km yr^?1, in spite of 10–20 times larger maximum dispersal distances. This is too slow to track the climate change‐driven range contraction of 4 or 8 km yr^?1 in the south resulting in metapopulation extinction. Especially K‐selected (large‐bodied) species are vulnerable in the simulations. With a temperature increase of 4°C per century bird species go extinct within 104–178 yr. We present a simple approximation formula to predict the mean time to metapopulation extinction using 1) the rate of climate change, which determines the speed of range contraction in the south, 2) the size of the distribution range, which serves as a buffer against extinction, and 3) the northward expansion velocity, determined by species traits and landscape properties. Finally, our results indicate that the northward expansion rate is not constant. It will be initially lagged suggesting that recently observed expansion rates might be underestimations of future northward expansion rates.     

Metapopulation extinction risk: Dispersal’s duplicity

Metapopulation extinction risk is the probability that all local populations are simultaneously extinct during a fixed time frame. Dispersal may reduce a metapopulation’s extinction risk by raising its average per-capita growth rate. By contrast, dispersal may raise a metapopulation’s extinction risk by reducing its average population density. Which effect prevails is controlled by habitat fragmentation. Dispersal in mildly fragmented habitat reduces a metapopulation’s extinction risk by raising its average per-capita growth rate without causing any appreciable drop in its average population density. By contrast, dispersal in severely fragmented habitat raises a metapopulation’s extinction risk because the rise in its average per-capita growth rate is more than offset by the decline in its average population density. The metapopulation model used here shows several other interesting phenomena. Dispersal in sufficiently fragmented habitat reduces a metapopulation’s extinction risk to that of a constant environment. Dispersal between habitat fragments reduces a metapopulation’s extinction risk insofar as local environments are asynchronous. Grouped dispersal raises the effective habitat fragmentation level. Dispersal search barriers raise metapopulation extinction risk. Nonuniform dispersal may reduce the effective fraction of suitable habitat fragments below the extinction threshold. Nonuniform dispersal may make demographic stochasticity a more potent metapopulation extinction force than environmental stochasticity.     

Simulating viability of a spadefoot toad Pelobates fuscus metapopulation in a landscape fragmented by a road

A stochastic simulation model allowing for demographic and environmental stochasticity was developed in order to predict the dynamics of a Pelobates fuscus metapopulation. The metapopulation (consisting of ca 1000 adult individuals) was intensively studied in the field for a period of four years and the simulation model was parameterised (sex ratio, age‐specific survival rates, fecundity and dispersal between subpopulations) using field data. A sensitivity analysis revealed that a change in the juvenile yearly survival rate has a relatively larger effect on subpopulation persistence than adult survival rate and fecundity rate have. The probability of subpopulation persistence for one hundred years increased from 0 to 0.6 with a change in juvenile yearly survival rate from 0.35 to 0.40. Varying dispersal rate (from 0 to 1% of the individuals in a subpopulation moving to another subpopulation within a year) showed that four of the five subpopulations are dependent on the last one for persistence, indicating a source‐sink structure of the metapopulation. The subpopulation with the highest estimated juvenile survival has a far higher persistence probability than the others; they in turn would go extinct were it not for the occasional input of individuals from the source subpopulation. The source‐sink structure was also apparent when simulating the isolation effect of the road: persistence of the subpopulation isolated by the road decreases markedly with only a 20% decrease in the number of individuals dispersing to this pond. Environmental stochasticity decreases persistence time of the source subpopulation, but increases persistence time of the unstable ones. This is probably due to the presence of overcompensating density‐dependent factors affecting the subpopulations and to the more general effect of stochasticity: it may temporarily change reproduction and mortality rates, increase time to extinction for unstable subpopulations through a rescue effect; and decrease time to extinction for the more stable subpopulations.