Effectiveness of Corridors Relative to Enlargement of Habitat Patches

Abstract:  The establishment of biological corridors between two otherwise isolated habitat patches is a common yet contentious strategy for conserving populations in fragmented landscapes. We compared the effectiveness of corridors with the effectiveness of an alternate conservation strategy, the enlargement of existing habitat patches. We used a spatially explicit population model that simulated population size in two kinds of patches. One patch had a corridor that connected it to a larger "source" patch and the other patch was unconnected and enlarged at the periphery by an area the same size as the corridor. Patch isolation, corridor width, patch size, and the probability that individuals would cross the border from habitat to matrix were varied independently. In general, population size was greater in enlarged patches than in connected patches when patches were relatively large and isolated. Corridor width and the probability of crossing the border from habitat to matrix did not affect the relative benefit of corridors versus patch enlargement. Although biological corridors may mitigate potential effects of inbreeding depression at long time scales, our results suggest that they are not always the best method of conserving fragmented populations.     

Influence of habitat quality, population size, patch size, and connectivity on patch-occupancy dynamics of the middle spotted woodpecker

Despite extensive research on the effects of habitat fragmentation, the ecological mechanisms underlying colonization and extinction processes are poorly known, but knowledge of these mechanisms is essential to understanding the distribution and persistence of populations in fragmented habitats. We examined these mechanisms through multiseason occupancy models that elucidated patch-occupancy dynamics of Middle Spotted Woodpeckers (Dendrocopos medius) in northwestern Spain. The number of occupied patches was relatively stable from 2000 to 2010 (15-24% of 101 patches occupied every year) because extinction was balanced by recolonization. Larger and higher quality patches (i.e., higher density of oaks >37 cm dbh [diameter at breast height]) were more likely to be occupied. Habitat quality (i.e., density of large oaks) explained more variation in patch colonization and extinction than did patch size and connectivity, which were both weakly associated with probabilities of turnover. Patches of higher quality were more likely to be colonized than patches of lower quality. Populations in high-quality patches were less likely to become extinct. In addition, extinction in a patch was strongly associated with local population size but not with patch size, which means the latter may not be a good surrogate of population size in assessments of extinction probability. Our results suggest that habitat quality may be a primary driver of patch-occupancy dynamics and may increase the accuracy of models of population survival. We encourage comparisons of competing models that assess occupancy, colonization, and extinction probabilities in a single analytical framework (e.g., dynamic occupancy models) so as to shed light on the association of habitat quality and patch geometry with colonization and extinction processes in different settings and species.     

Persistence of Different-sized Populations: An Empirical Assessment of Rapid Extinctions in Bighorn Sheep

Abstract: Theory and simulation models suggest that small populations are more susceptible to extinction than large populations, yet assessment of this idea has been hampered by lack of an empirical base. I address the problem by asking how long different-sized populations persist and present demographic and weather data spanning up to 70 years for 122 bighorn sheep ( Ovis canadensis ) populations in southwestern North America Analyses reveal that: (1) 100 percent of the populations with fewer than 50 individuals went extinct within 50 years; (2) populations with greater than 100 individuals persisted for up to 70 years; and (3) the rapid loss of populations was not likely to be caused by food shortages, severe weather, predation, or interspecific competition These data suggest that population size is a marker of persistence trajectories and they indicate that local extinction cannot be overcome because 50 individuals, even in the short term, are not a minimum viable population size for bighorn sheep.     

Populations in small, ephemeral habitat patches may drive dynamics in a Daphnia magna metapopulation

Migration is the key process to understand the dynamics and persistence of a metapopulation. Many metapopulation models assume a positive correlation between habitat patch size or stability and the number of emigrants. However, few empirical data exist, and habitat patch size and habitat stability may affect dispersal differently than they affect local persistence. Here, we studied the production of the migration stage (i.e., resting eggs called ephippia) of the cladoceran Daphnia magna in a metapopulation consisting of 530 rock pool habitat patches over 25 years. Earlier, the functioning of this metapopulation was explained with a Levins-type metapopulation model or with a mainland-island metapopulation model, based on local extinction and colonization data or time series data, respectively. We used pool volume, hydroperiod length, and number of desiccation events to calculate per-pool production of ephippia (i.e., migration stages). We estimated that populations in small and ephemeral habitat patches produced more than half of the 250 000 to 1 million ephippia that were produced in the metapopulation as a whole per year between 1982 and 2006. Furthermore, these small populations contributed approximately 90% of the ephippia exposed during desiccation events, while the contribution of the long-lived populations in large pools was minimal. We term this an "inverse mainland-island" type metapopulation and propose that populations in small, ephemeral habitat patches may also be the driving force for metapopulation dynamics in other systems.     

Population History, Genetic Variability, and Horn Growth in Bighorn Sheep

Bighorn sheep ( Ovis canadensis ) are restricted in distribution and numbers relative to presettlement conditions. Some populations have alledgedly suffered losses of fitness resulting from small, insular populations and a breeding system that reduces effective population size. Large horns in rams, which confer breeding superiority, are absent from some populations, and this absence may result in part from loss of genetic variability. We investigated the relationship among allozyme variability, population history, and horn growth in bighorn sheep from the Rocky Mountains. Heterozygosity was higher for bighorn sheep than has been reported for Dall sheep ( O. dalli ). Heterozygosity and allelic variability were marginally related to effective population size for the proceeding 15 years. Horn growth was significantly higher in more heterozygous than in less heterozygous rams for years 6, 7, and 8 of life. By the end of year 8, more heterozygous rams had 13% higher horn volumes than less heterozygous rams. Most hunting of bighorn sheep involves selective removal of large-horned rams, which we hypothesize may reduce genetic variability of these populations and contribute to losses in fitness.     

Linking edge effects and patch size effects: Importance of matrix nest predators

The edge effect is usually considered to be the proximate cause of area sensitivity in forest birds. We tested if birds nesting in large patches are less vulnerable to the edge effect using a simple model that assumes an increase in patch size reduces the probability of a matrix predator moving to the core areas of forest and that larger perimeter/area ratios result in a higher number of matrix predators per unit of area. The probability of a nest being successful decreased asymptotically with an increase in either the patch penetration distance of predators or predator density, but those effects were reduced when patch size was increased. Large patches have a lower probability of being affected by an Allee effect and they can function as sink habitats only if penetration distance and predator density are largely increased. However, the transition from an Allee effect to a sink condition occurs with a small increase in penetration distance and predator density. Since birds nesting in large patches are less vulnerable to an increase in matrix predator populations, persistence of bird populations may be possible by increasing the size of habitat patches that can act as source populations.     

Effects of Climate Change on Population Persistence of Desert-Dwelling Mountain Sheep in California

Abstract:  Metapopulations may be very sensitive to global climate change, particularly if temperature and precipitation change rapidly. We present an analysis of the role of climate and other factors in determining metapopulation structure based on presence and absence data. We compared existing and historical population distributions of desert bighorn sheep ( Ovis canadensis ) to determine whether regional climate patterns were correlated with local extinction. To examine all mountain ranges known to hold or to have held desert bighorn populations in California and score for variables describing climate, metapopulation dynamics, human impacts, and other environmental factors, we used a geographic information system (GIS) and paper maps. We used logistic regression and hierarchical partitioning to assess the relationship among these variables and the current status of each population (extinct or extant). Parameters related to climate—elevation, precipitation, and presence of dependable springs—were strongly correlated with population persistence in the twentieth century. Populations inhabiting lower, drier mountain ranges were more likely to go extinct. The presence of domestic sheep grazing allotments was negatively correlated with population persistence. We used conditional extinction probabilities generated by the logistic-regression model to rank native, naturally recolonized, and reintroduced populations by vulnerability to extinction under several climate-change scenarios. Thus risk of extinction in metapopulations can be evaluated for global-climate-change scenarios even when few demographic data are available.     

Space use in lions, <Emphasis Type="Italic">Panthera leo</Emphasis>, in the Selous Game Reserve: social and ecological factors

Habitat qualities, such as food supply or access to refuges, often influence home-range size. Furthermore, such qualitative differences usually lead to conspecific competition over space, which can be an important factor in determining the distribution of individuals within populations. In carnivores, patterns of resource dispersion are hypothesized to determine home range-size and group size. But in contests over space (or other resources), larger groups usually dominate smaller ones, and group size should therefore also affect home-range size. Here I describe the space use of lions, Panthera leo, in the Selous Game Reserve, Tanzania, and ask whether space use is related to pride size, habitat, or relatedness. Home ranges varied in size, but size showed no correlation to number of adult females in the pride or to habitat type. Lions exhibited a significant preference for riverine and short-grass habitat, and a significant avoidance of acacia woodland. Habitat preference ratios largely reflected prey availability in each habitat. Outer areas, as well as core areas of home ranges, were often used by two or more prides. Overlaps showed no correlation to relatedness among prides or habitat type. Thus, whereas home-range sizes and overlaps were determined by factors that could not be revealed by demographic factors or analyses of habitat composition or genetic structure, lion space use within each home range seemed driven mostly by prey availability, which mainly varies with habitat type.     

Inferring Metapopulation Dynamics from Patch-Level Incidence of Florida Scrub Plants

Spatial structure and dynamics of multiple populations may explain species distribution patterns in patchy communities with heterogeneous disturbance regimes, especially when species have poor dispersal. The endemic-rich Florida (U.S.A.) rosemary scrub occupies about 4% of the west portion of Archbold Biological Station and occurs scattered within a matrix of less xeric vegetation. Longer fire-return times and higher frequency of open patches in rosemary scrub provide favorable habitat for many plant species. Occupancy of 123 species of vascular plants and ground lichens in 89 patches was determined by repeated site surveys. About two-thirds of the species occurring at more than 14 patches had a significant logistic regression of presence on time-since-fire, patch size, patch isolation, or their interactions. Species with presence related to the interaction between patch isolation and patch size were primarily herbs and small shrubs specializing in rosemary scrub. These results suggest the importance of spatial characteristics of the landscape for population turnover of these species. An incidence-based metapopulation model was used to predict extinction and colonization probabilities of those species with presence in rosemary scrub patches related to the studied spatial variables. This is the first attempt to apply incidence-based metapopulation models to plants. The results showed stronger effects of patch size and patch isolation on extinction probabilities of herbs than on those of woody species. Because of their effect on spatial heterogeneity and habitat availability, fire suppression and habitat destruction may decrease persistence probabilities for these rosemary scrub specialists, many of which are endangered species.     

Population Demographics and Genetic Diversity in Remnant and Translocated Populations of Sea Otters

Abstract: The effects of small population size on genetic diversity and subsequent population recovery are theoretically predicted, but few empirical data are available to describe those relations. We use data from four remnant and three translocated sea otter ( Enhydra lutris ) populations to examine relations among magnitude and duration of minimum population size, population growth rates, and genetic variation. Mitochondrial (mt)DNA haplotype diversity was correlated with the number of years at minimum population size ( r _s = −0.741, p = 0.038) and minimum population size ( r _s = 0.709, p = 0.054). We found no relation between population growth and haplotype diversity, although growth was significantly greater in translocated than in remnant populations. Haplotype diversity in populations established from two sources was higher than in a population established from a single source and was higher than in the respective source populations. Haplotype frequencies in translocated populations of founding sizes of 4 and 28 differed from expected, indicating genetic drift and differential reproduction between source populations, whereas haplotype frequencies in a translocated population with a founding size of 150 did not. Relations between population demographics and genetic characteristics suggest that genetic sampling of source and translocated populations can provide valuable inferences about translocations.     

Captive Breeding and Reintroduction Evaluation Criteria: a Case Study of Peninsular Bighorn Sheep

Abstract: Captive breeding and reintroduction programs are rarely evaluated, and assessment criteria vary widely. We used the following criteria to evaluate a bighorn sheep ( Ovis canadensis ) augmentation program: (1) survival and recruitment rates in the captive population, (2) survival of released animals, (3) recruitment of released animals, (4) growth rate of the reintroduced or augmented population, and (5) establishment of a viable wild population. Captive bighorn survival and recruitment was high, averaging 0.98 (SD = 0.05) and 71.0% (SD = 19.4), respectively. Annual survival of free-ranging captive-reared bighorn ( n = 73, x = 0.80, SD = 0.11) did not differ (   Z = −0.85, p = 0.40; n = 14) from survival of wild-reared bighorn ( n = 43, x = 0.81, SD = 0.12). Recruitment was unusually low for both captive-reared (  x = 13.7%, SD = 0.24) and wild-reared ewes (  x = 13.7%, SD = 0.20). Although reintroduction did not result in population growth or establishment of a viable population, it helped prevent extirpation of the reinforced deme, preserved metapopulation linkage, and aided habitat preservation. Chronic low recruitment and low adult survivorship precluded achievement of criteria 3–5. Environmental conditions in the release area also appeared to hinder program success. Standard evaluation criteria for ongoing reintroductions allow for informative assessments and facilitate comparisons needed to refine reintroduction science as a recovery tool for threatened or endangered populations.     

Potential Effects of the United States-Mexico Border Fence on Wildlife

Abstract:  Security infrastructure along international boundaries threatens to degrade connectivity for wildlife. To explore potential effects of a fence under construction along the U.S.–Mexico border on wildlife, we assessed movement behavior of two species with different life histories whose regional persistence may depend on transboundary movements. We used radiotelemetry to assess how vegetation and landscape structure affect flight and natal dispersal behaviors of Ferruginous Pygmy-Owls ( Glaucidium brasilianum ), and satellite telemetry, gene-flow estimates, and least-cost path models to assess movement behavior and interpopulation connectivity of desert bighorn sheep ( Ovis canadensis mexicana ). Flight height of Pygmy-Owls averaged only 1.4 m (SE 0.1) above ground, and only 23% of flights exceeded 4 m. Juvenile Pygmy-Owls dispersed at slower speeds, changed direction more, and had lower colonization success in landscapes with larger vegetation openings or higher levels of disturbance ( p ≤ 0.047), which suggests large vegetation gaps coupled with tall fences may limit transboundary movements. Female bighorn sheep crossed valleys up to 4.9 km wide, and microsatellite analyses indicated relatively high levels of gene flow and migration (95% CI for F_ST= 0.010–0.115, Nm = 1.9–24.8, M = 10.4–15.4) between populations divided by an 11-km valley. Models of gene flow based on regional topography and movement barriers suggested that nine populations of bighorn sheep in northwestern Sonora are linked by dispersal with those in neighboring Arizona. Disruption of transboundary movement corridors by impermeable fencing would isolate some populations on the Arizona side. Connectivity for other species with similar movement abilities and spatial distributions may be affected by border development, yet mitigation strategies could address needs of wildlife and humans.     

Local Population Dynamics in Metapopulation Models: Implications for Conservation

Abstract: Habitat fragmentation and the division of populations into spatially separated units have led to the increasing use of metapopulation models to characterize these populations. One prominent model that has served as a heuristic tool was introduced by Levins and is based on a collection of simplifying assumptions that exclude information on the dynamics and spatial distribution of local populations. Levins's and similar models predict the proportion of occupied habitat patches at equilibrium and the conditions needed to avoid total extinction. There are many obvious concerns about using such models, including how realistic alterations might change the predictions and whether occupancy has any relationship to population-level processes. Although many of the assumptions of these simple models are known to be unrealistic, we do not know how the assumptions affect model predictions. We simulated a metapopulation, and our results show that assumptions such as homogeneity of habitat patches, random migration among patches, equivalent extinction probabilities in all patches, and a large number of patches can lead to large overestimations of habitat occupancy. But when we explicitly modeled the underlying population dynamics within each patch, we found (1) that there was a strong correlation between proportion of occupied patches and total metapopulation size and (2) that the distribution of individuals among patches was relatively insensitive to model assumptions. Thus, our results show that although realistic modifications will change model predictions for occupancy, occupancy and population trends will be correlated. These correlations between occupancy and population size suggest that occupancy models may have some utility in conservation applications.     

Spatial characteristics of habitat patches and population survival

The effect of the spatial arrangement of habitat patches on the survival of resident populations were considered in a stochastic model using population parameters appropriate to Peromyscus leucopus. The 34 possible arrangements of connections among five otherwise identical patches were simulated in order to determine the survival probabilities and population sizes.The main findings are that populations in completely isolated patches have lower survival probabilities than those in patches that are connected to other patches, are connected to the survival probabilities of populations in connected patches increases with the size of the largest geometric figure of which the patch is a part. The results are discussed in the context of resource management.     

Translating effects of inbreeding depression on component vital rates to overall population growth in endangered bighorn sheep

Evidence of inbreeding depression is commonly detected from the fitness traits of animals, yet its effects on population growth rates of endangered species are rarely assessed. We examined whether inbreeding depression was affecting Sierra Nevada bighorn sheep (Ovis canadensis sierrae), a subspecies listed as endangered under the U.S. Endangered Species Act. Our objectives were to characterize genetic variation in this subspecies; test whether inbreeding depression affects bighorn sheep vital rates (adult survival and female fecundity); evaluate whether inbreeding depression may limit subspecies recovery; and examine the potential for genetic management to increase population growth rates. Genetic variation in 4 populations of Sierra Nevada bighorn sheep was among the lowest reported for any wild bighorn sheep population, and our results suggest that inbreeding depression has reduced adult female fecundity. Despite this population sizes and growth rates predicted from matrix-based projection models demonstrated that inbreeding depression would not substantially inhibit the recovery of Sierra Nevada bighorn sheep populations in the next approximately 8 bighorn sheep generations (48 years). Furthermore, simulations of genetic rescue within the subspecies did not suggest that such activities would appreciably increase population sizes or growth rates during the period we modeled (10 bighorn sheep generations, 60 years). Only simulations that augmented the Mono Basin population with genetic variation from other subspecies, which is not currently a management option, predicted significant increases in population size. Although we recommend that recovery activities should minimize future losses of genetic variation, genetic effects within these endangered populations-either negative (inbreeding depression) or positive (within subspecies genetic rescue)-appear unlikely to dramatically compromise or stimulate short-term conservation efforts. The distinction between detecting the effects of inbreeding depression on a component vital rate (e.g., fecundity) and the effects of inbreeding depression on population growth underscores the importance of quantifying inbreeding costs relative to population dynamics to effectively manage endangered populations.     

Local extinction of grassland plants: the landscape matrix is more important than patch attributes

Past studies of local extinctions in fragmented habitats most often tested the influence of fragment size and isolation while ignoring how differences in the surrounding landscape matrix may govern extinction. We assessed how both the spatial attributes of remnant patches (area and isolation) and landscape factors (extent of urbanization and maximum inter-fire interval) influence the persistence of native plant species in grasslands in western Victoria, Australia. Persistence was determined in 2001 by resurveying 30 remnants first surveyed in the 1980s, and correlates of extinction were assessed using Bayesian logistic regression models. On average, 26% of populations of native species became locally extinct over two decades. Area and isolation had little effect on the probability of local extinction, but urbanization and longer maximum inter-fire intervals increased extinction risk. These findings suggest that the native grasslands studied are relatively insensitive to area- and isolation-based fragmentation effects and that short-term persistence of plant populations requires the maintenance of habitat quality. The latter is strongly influenced by the landscape matrix surrounding remnant patches through changes in fire regimes and increased exogenous disturbance.     

Dispersal in Spatially Explicit Population Models

Abstract: Ruckelshaus et al. (1997) outlined a simulation model of dispersal between patches in a fragmented landscape. They showed that dispersal success—the proportion of dispersers successfully locating a patch—was particularly sensitive to errors in dispersal mortality and concluded that this limits the utility of spatially explicit population models in conservation biology. I contend that, although they explored error propagation in a simple dispersal model, they did not explore how errors are propagated in spatially explicit population models, as no consideration of population processes was included. I developed a simple simulation model to investigate the effect of varying dispersal success on predictions of patch occupancy and population viability, the conventional outputs of spatially explicit population models. The model simulates births and deaths within habitat patches and dispersal as the transfer of individuals between them. Model predictions were sensitive to changes in dispersal success across a restricted range of within-patch growth rates, which depended on the dispersal initiation mechanism, patch carrying capacities, and number of generations simulated. Predictions of persistence and patch occupancy were generally more sensitive to changes in dispersal success (1) under presaturation rather than saturation dispersal; (2) at lower patch carrying capacities; and (3) over longer time periods. The framework I present provides a means of assessing, quantitatively, the regions of parameter space for which differences in dispersal success are likely to have a large effect on population model outputs. Investigating the effect of the representation of dispersal behavior within the demographic and landscape context provides a more useful assessment of whether our lack of knowledge is likely to cause unacceptable uncertainty in the predictions of spatially explicit population models.     

Patch Occupancy and Potential Metapopulation Dynamics of Three Forest Mammals in Fragmented Afromontane Forest in South Africa

Abstract: We investigated the persistence of three medium-sized (2–9 kg), rare forest mammals in the fragmented mist-belt Podocarpus forests of the midlands of KwaZulu-Natal Province, South Africa. We recorded patch occupancy of blue duiker (   Philantomba monticola ), tree hyrax (   Dendrohyrax arboreus ), and samango monkey ( Cercopithecus mitis labiatus ) in 199 forest patches. Their rarity is ascribed to the fragmentation and destruction of their forest habitat. Incidence functions, derived from presence and absence data, were formulated as generalized linear models, and environmental effects were included in the fitted logistic models. The small and mostly solitary hyrax and duiker persisted in smaller patches than the large and social monkey. Although this result follows expectations based on relative home-range sizes of each species, the incidence probability of the samango monkey was invariant with increasing isolation, whereas a gradual decrease with increasing isolation was observed for the hyrax and duiker. Group dynamics may inhibit dispersal and increase the isolation effect in social species such as samango monkeys. A mainland-island metapopulation model adequately describes patterns of patch occupancy by the hyrax and duiker, but the monkeys' poor dispersal ability and obvious area-dependent extirpation suggest that they exist in transient, nonequilibrium (declining) metapopulations. Through identification of large forest patches for careful protection and management, the survival of all three species—especially the monkey—could be prolonged. Because no functional metapopulation may exist for the monkey, however, this is an emergency measure. For the duiker and hyrax, larger patches should form part of a network of smaller and closer patches in a natural matrix.     

Long-term drivers of persistence and colonization dynamics in spatially structured amphibian populations

Many organisms live in networks of local populations connected by dispersing individuals, called spatially structured populations (SSPs), where the long-term persistence of the entire network is determined by the balance between 2 processes acting at the scale of local populations: extinction and colonization. When multiple threats act on an SSP, a comparison of the different factors determining local extinctions and colonizations is essential to plan sound conservation actions. We assessed the drivers of long-term population dynamics of multiple amphibian species at the regional scale. We used dynamic occupancy models within a Bayesian framework to identify the factors determining persistence and colonization of local populations. Because connectivity among patches is fundamental to SSPs dynamics, we considered 2 measures of connectivity acting on each focal patch: incidence of the focal species and incidence of invasive crayfish. We used meta-analysis to summarize the effect of different drivers at the community level. Persistence and colonization of local populations were jointly determined by factors acting at different scales. Persistence probability was positively related to the area and the permanence of wetlands, whereas it was negatively related to occurrence of fish. Colonization probability was highest in semipermanent wetlands and in sites with a high incidence of the focal species in nearby sites, whereas it showed a negative relationship with the incidence of invasive crayfish in the landscape. By analyzing long-term data on amphibian population dynamics, we found a strong effect of some classic features commonly used in SSP studies, such as patch area and focal species incidence. The presence of an invasive non-native species at the landscape scale emerged as one of the strongest drivers of colonization dynamics, suggesting that studies on SSPs should consider different connectivity measures more frequently, such as the incidence of predators, especially when dealing with biological invasions.     

Local population size in a flightless insect: importance of patch structure-dependent mortality

In spatially heterogeneous systems, utilizing population models to integrate the effects of multiple population rates can yield powerful insights into the relative importance of the component rates. The relative importance of demographic rates and dispersal in shaping the distribution of the western tussock moth (Orgyia vetusta) among patches of its host plant was explored using stage-structured population models. Tussock moth dispersal occurs passively in first-instar larvae and is poor or absent in all other life stages. Spatial surveys suggested, however, that moth distribution is not well explained by passive dispersal; moth populations were greater on small patches and on isolated ones. Further analysis showed that several local demographic rates varied significantly with patch characteristics. Two mortality factors in particular may explain the observed patterns. First, crawler mortality both increased with patch size and was density-dependent. A single-patch difference equation model showed mortality related to patch size is strong enough to overcome the homogenizing effect of density dependence; greater equilibrium densities were predicted for smaller patches. Second, although three rates were found to vary with local patch density, only pupal parasitism by a chalcid wasp could potentially account for higher moth abundances on isolated patches. A spatially explicit simulation model of the multiple-patch system showed that spatial variation in pupal parasitism is indeed strong enough to generate such a pattern. These results demonstrate that habitat spatial structure can affect multiple population processes simultaneously, and even relatively low attack rates imposed on a reproductively valuable life stage of the host can have a dominant effect on population distribution among habitat patches.