Estimating the complexity of bird song by using capture-recapture approaches from community ecology

Repertoire size, the number of unique song or syllable types in the repertoire, is a widely used measure of song complexity in birds, but it is difficult to calculate this exactly in species with large repertoires. A new method of repertoire size estimation applies species richness estimation procedures from community ecology, but such capture-recapture approaches have not been much tested. Here, we establish standardized sampling schemes and estimation procedures using capture-recapture models for syllable repertoires from 18 bird species, and suggest how these may be used to tackle problems of repertoire estimation. Different models, with different assumptions regarding the heterogeneity of the use of syllable types, performed best for different species with different song organizations. For most species, models assuming heterogeneous probability of occurrence of syllables (so-called detection probability) were selected due to the presence of both rare and frequent syllables. Capture-recapture estimates of syllable repertoire size from our small sample did not differ significantly from previous estimates using larger samples of count data. However, the enumeration of syllables in 15 songs yielded significantly lower estimates than previous reports. Hence, heterogeneity in detection probability of syllables should be addressed when estimating repertoire size. This is neglected using simple enumeration procedures, but is taken into account when repertoire size is estimated by appropriate capture-recapture models adjusted for species-specific song organization characteristics. We suggest that such approaches, in combination with standardized sampling, should be applied in species with potentially large repertoire size. On the other hand, in species with small repertoire size and homogenous syllable usage, enumerations may be satisfactory. Although researchers often use repertoire size as a measure of song complexity, listeners to songs are unlikely to count entire repertoires and they may rely on other cues, such as syllable detection probability.Communicated by A. Cockburn     

An incidence-based richness estimator for quadrats sampled without replacement

Most richness estimators currently in use are derived from models that consider sampling with replacement or from the assumption of infinite populations. Neither of the assumptions is suitable for sampling sessile organisms such as plants where quadrats are often sampled without replacement and the area of study is always limited. In this paper, we propose an incidence-based parametric richness estimator that considers quadrat sampling without replacement in a fixed area. The estimator is derived from a zero-truncated binomial distribution for the number of quadrats containing a given species (e.g., species i) and a modified beta distribution for the probability of presence-absence of a species in a quadrat. The maximum likelihood estimate of richness is explicitly given and can be easily solved. The variance of the estimate is also obtained. The performance of the estimator is tested against nine other existing incidence-based estimators using two tree data sets where the true numbers of species are known. Results show that the new estimator is insensitive to sample size and outperforms the other methods as judged by the root mean squared errors. The superiority of the new method is particularly noticeable when large quadrat size is used, suggesting that a few large quadrats are preferred over many small ones when sampling diversity.     

Multiple data sources improve DNA-based mark-recapture population estimates of grizzly bears.

A fundamental challenge to estimating population size with mark-recapture methods is heterogeneous capture probabilities and subsequent bias of population estimates. Confronting this problem usually requires substantial sampling effort that can be difficult to achieve for some species, such as carnivores. We developed a methodology that uses two data sources to deal with heterogeneity and applied this to DNA mark-recapture data from grizzly bears (Ursus arctos). We improved population estimates by incorporating additional DNA "captures" of grizzly bears obtained by collecting hair from unbaited bear rub trees concurrently with baited, grid-based, hair snag sampling. We consider a Lincoln-Petersen estimator with hair snag captures as the initial session and rub tree captures as the recapture session and develop an estimator in program MARK that treats hair snag and rub tree samples as successive sessions. Using empirical data from a large-scale project in the greater Glacier National Park, Montana, USA, area and simulation modeling we evaluate these methods and compare the results to hair-snag-only estimates. Empirical results indicate that, compared with hair-snag-only data, the joint hair-snag-rub-tree methods produce similar but more precise estimates if capture and recapture rates are reasonably high for both methods. Simulation results suggest that estimators are potentially affected by correlation of capture probabilities between sample types in the presence of heterogeneity. Overall, closed population Huggins-Pledger estimators showed the highest precision and were most robust to sparse data, heterogeneity, and capture probability correlation among sampling types. Results also indicate that these estimators can be used when a segment of the population has zero capture probability for one of the methods. We propose that this general methodology may be useful for other species in which mark-recapture data are available from multiple sources.     

Optimization of capture–recapture monitoring of elusive species illustrated with a threatened grasshopper

Information on population sizes and trends of threatened species is essential for their conservation, but obtaining reliable estimates can be challenging. We devised a method to improve the precision of estimates of population size obtained from capture–recapture studies for species with low capture and recapture probabilities and short seasonal activity, illustrated with population data of an elusive grasshopper (Prionotropis rhodanica). We used data from 5 capture–recapture studies to identify methodological and environmental factors affecting capture and recapture probabilities and estimates of population size. In a simulation, we used the population size and capture and recapture probability estimates obtained from the field studies to identify the minimum number of sampling occasions needed to obtain unbiased and robust estimates of population size. Based on these results we optimized the capture–recapture design, implemented it in 2 additional studies, and compared their precision with those of the nonoptimized studies. Additionally, we simulated scenarios based on thresholds of population size in criteria C and D of the International Union for Conservation of Nature (IUCN) Red List to investigate whether estimates of population size for elusive species can reliably inform red-list assessments. Identifying parameters that affect capture and recapture probabilities (for the grasshopper time since emergence of first adults) and optimizing field protocols based on this information reduced study effort (−6% to −27% sampling occasions) and provided more precise estimates of population size (reduced coefficient of variation) compared with nonoptimized studies. Estimates of population size from the scenarios based on the IUCN thresholds were mostly unbiased and robust (only the combination of very small populations and little study effort produced unreliable estimates), suggesting capture–recapture can be considered reliable for informing red-list assessments. Although capture–recapture remains difficult and costly for elusive species, our optimization procedure can help determine efficient protocols to increase data quality and minimize monitoring effort.     

Detecting Rare Species with Random or Subjective Sampling: a Case Study of Red-Listed Saproxylic Beetles in Boreal Sweden

Abstract:  Efficient sampling design in field studies is important for economical and statistical reasons. We compared two ways to distribute sampling effort over an area, either randomly or subjectively. We searched for red-listed saproxylic (wood-living) beetles in 30 spruce stands in boreal Sweden by sifting wood from dead trees. We randomly selected positions within each stand with a geographic positioning system and sampled the nearest dead tree (random sample). In the same stand we also sampled dead trees that, based on literature, were likely to host such species (subjective sampling). The subjective sampling (two to five samples per stand, depending on stand size) was compared with the higher, random sampling effort (fixed level of 12 samples/stand). Subjective sampling was significantly more efficient. Red-listed species were found in 36% of the subjective samples and in 16% of the random samples. Nevertheless, the larger random effort resulted in a comparable number of red-listed species per stand and in 13 detected species in total (vs. 12 species with subjective sampling). Random sampling was less efficient, but provided an unbiased alternative more suitable for statistical purposes, as needed in, for example, monitoring programs. Moreover, new species-specific knowledge can be gained through random searches.     

Effects of sampling completeness on the structure of plant-pollinator networks

Plant-animal interaction networks provide important information on community organization. One of the most critical assumptions of network analysis is that the observed interaction patterns constitute an adequate sample of the set of interactions present in plant-animal communities. In spite of its importance, few studies have evaluated this assumption, and in consequence, there is no consensus on the sensitivity of network metrics to sampling methodological shortcomings. In this study we examined how variation in sampling completeness influences the estimation of six network metrics frequently used in the literature (connectance, nestedness, modularity, robustness to species loss, path length, and centralization). We analyzed data of 186 flowering plants and 336 pollinator species in 10 networks from a forest-fragmented system in central Chile. Using species-based accumulation curves, we estimated the deviation of network metrics in undersampled communities with respect to exhaustively sampled communities and the effect of network size and sampling evenness on network metrics. Our results indicate that: (1) most metrics were affected by sampling completeness but differed in their sensitivity to sampling effort; (2) nestedness, modularity, and robustness to species loss were less influenced by insufficient sampling than connectance, path length, and centralization; (3) robustness was mildly influenced by sampling evenness. These results caution studies that summarize information from databases with high, or unknown, heterogeneity in sampling effort per species and should stimulate researchers to report sampling intensity to standardize its effects in the search for broad patterns in plant-pollinator networks.     

Nonparametric estimation of Shannon’s index of diversity when there are unseen species in sample

A biological community usually has a large number of species with relatively small abundances. When a random sample of individuals is selected and each individual is classified according to species identity, some rare species may not be discovered. This paper is concerned with the estimation of Shannons index of diversity when the number of species and the species abundances are unknown. The traditional estimator that ignores the missing species underestimates when there is a non-negligible number of unseen species. We provide a different approach based on unequal probability sampling theory because species have different probabilities of being discovered in the sample. No parametric forms are assumed for the species abundances. The proposed estimation procedure combines the Horvitz–Thompson (1952) adjustment for missing species and the concept of sample coverage, which is used to properly estimate the relative abundances of species discovered in the sample. Simulation results show that the proposed estimator works well under various abundance models even when a relatively large fraction of the species is missing. Three real data sets, two from biology and the other one from numismatics, are given for illustration.     

Complex national sampling design for long-term monitoring of protected dry grasslands in Switzerland

We describe a probabilistic sampling design of circular permanent plots for the long-term monitoring of protected dry grasslands in Switzerland. The population under study is defined by the perimeter of a national inventory. The monitoring focus is on the species composition of the protected grassland vegetation and derived conservation values. Efficient trend estimations are required for the whole country and for some predefined target groups (six biogeographical regions and eleven vegetation types). The target groups are equally important regardless of their size. Consequently, intensified sampling of the less frequent groups is essential for sample efficiency. The prior information needed to draw a targeted sample is obtained from the sampling frame and external databases. The logistics and generalized delineation of the target population may pose further problems. Thus, investments in fieldwork and travel time should be well balanced by selecting a cluster sample. Second, any access problems in the field and non-target units in the sample should be compensated for by selecting reserve plots as they otherwise may considerably reduce the effective sample size. Finally, the design has to be flexible as the sampling frame may change over time and sampling intensity might have to be adjusted to redefined budgets or requirements. Likewise, the variables and biological items of interest may change. To fulfil all these constraints and to optimally use the available prior information, we propose a multi-stage self-weighted unequal probability sampling design. The design uses modern techniques such as: balanced sampling, spreading, stratified balancing, calibration, unequal probability sampling and power allocation. This sampling design meets the numerous requirements of this study and provides a very efficient estimator.     

Species abundance distributions result from body size-energetics relationships

Abundance distributions are a central characteristic of ecosystems. Certain distributions have been derived from theoretical models of community organization, and therefore the fit of data to these distributions has been proposed as a test of these theories. However, it is shown here that the geometric sequence distribution can be derived directly from the empirical relationship between population density and body size, with the assumption of random or uniform body size distributions on a log scale (as holds at local scales). The geometric sequence model provides a good to excellent fit to empirical data. The presence of noise in the relationship between population density and body size creates a curve that begins to approximate a lognormal species abundance distribution as the noise term increases. For continental-scale data in which the body size distribution is not flat, the result of sampling tends again toward the lognormal. Repeat sampling over time smooths out species population fluctuations and damps out the noise, giving a more precise geometric sequence abundance distribution. It is argued that the direct derivation of this distribution from empirical relationships gives it priority over distributions derived from complex theoretical community models.     

A combination line transect and capture-recapture sampling model for multiple observers in aerial surveys

We present a robust sampling methodology to estimate population size using line transect and capture-recapture procedures for aerial surveys. Aerial surveys usually underestimate population density due to animals being missed. A combination of capture-recapture and line transect sampling methods with multiple observers allows violation of the assumption that all animals on the centreline are sighted from the air. We illustrate our method with an example of inanimate objects which shows evidence of failure of the assumption that all objects on the centreline have probability 1 of being detected. A simulation study is implemented to evaluate the performance of three variations of the Lincoln-Petersen estimator: the overall estimator, the stratified estimator, and the general stratified estimator based on the combined likelihood proposed in this paper. The stratified Lincoln-Petersen estimator based on the combined likelihood is found to be generally superior to the other estimators.     

Chilean Blue Whales as a Case Study to Illustrate Methods to Estimate Abundance and Evaluate Conservation Status of Rare Species

Abstract: Often abundance of rare species cannot be estimated with conventional design‐based methods, so we illustrate with a population of blue whales (Balaenoptera musculus) a spatial model‐based method to estimate abundance. We analyzed data from line‐transect surveys of blue whales off the coast of Chile, where the population was hunted to low levels. Field protocols allowed deviation from planned track lines to collect identification photographs and tissue samples for genetic analyses, which resulted in an ad hoc sampling design with increased effort in areas of higher densities. Thus, we used spatial modeling methods to estimate abundance. Spatial models are increasingly being used to analyze data from surveys of marine, aquatic, and terrestrial species, but estimation of uncertainty from such models is often problematic. We developed a new, broadly applicable variance estimator that showed there were likely 303 whales (95% CI 176–625) in the study area. The survey did not span the whales' entire range, so this is a minimum estimate. We estimated current minimum abundance relative to pre‐exploitation abundance (i.e., status) with a population dynamics model that incorporated our minimum abundance estimate, likely population growth rates from a meta‐analysis of rates of increase in large baleen whales, and two alternative assumptions about historic catches. From this model, we estimated that the population was at a minimum of 9.5% (95% CI 4.9–18.0%) of pre‐exploitation levels in 1998 under one catch assumption and 7.2% (CI 3.7–13.7%) of pre‐exploitation levels under the other. Thus, although Chilean blue whales are probably still at a small fraction of pre‐exploitation abundance, even these minimum abundance estimates demonstrate that their status is better than that of Antarctic blue whales, which are still <1% of pre‐exploitation population size. We anticipate our methods will be broadly applicable in aquatic and terrestrial surveys for rarely encountered species, especially when the surveys are intended to maximize encounter rates and estimate abundance.     

Spatially Explicit Power Analyses for Occupancy‐Based Monitoring of Wolverine in the U.S. Rocky Mountains

Conservation scientists and resource managers often have to design monitoring programs for species that are rare or patchily distributed across large landscapes. Such programs are frequently expensive and seldom can be conducted by one entity. It is essential that a prospective power analysis be undertaken to ensure stated monitoring goals are feasible. We developed a spatially based simulation program that accounts for natural history, habitat use, and sampling scheme to investigate the power of monitoring protocols to detect trends in population abundance over time with occupancy‐based methods. We analyzed monitoring schemes with different sampling efforts for wolverine (Gulo gulo) populations in 2 areas of the U.S. Rocky Mountains. The relation between occupancy and abundance was nonlinear and depended on landscape, population size, and movement parameters. With current estimates for population size and detection probability in the northern U.S. Rockies, most sampling schemes were only able to detect large declines in abundance in the simulations (i.e., 50% decline over 10 years). For small populations reestablishing in the Southern Rockies, occupancy‐based methods had enough power to detect population trends only when populations were increasing dramatically (e.g., doubling or tripling in 10 years), regardless of sampling effort. In general, increasing the number of cells sampled or the per‐visit detection probability had a much greater effect on power than the number of visits conducted during a survey. Although our results are specific to wolverines, this approach could easily be adapted to other territorial species. Poder de Análisis Espacialmente Explícito para el Monitoreo Basado en Ocupación del Glotón (Gulo gulo) en las Montañas Rocallosas de Estados Unidos     

Effects of sampling error and temporal correlations in population growth on process variance estimators

Estimates of a population’s growth rate and process variance from time-series data are often used to calculate risk metrics such as the probability of quasi-extinction, but temporal correlations in the data from sampling error, intrinsic population factors, or environmental conditions can bias process variance estimators and detrimentally affect risk predictions. It has been claimed (McNamara and Harding, Ecol Lett 7:16–20, 2004) that estimates of the long-term variance that incorporate observed temporal correlations in population growth are unaffected by sampling error; however, no estimation procedures were proposed for time-series data. We develop a suite of such long-term variance estimators, and use simulated data with temporally autocorrelated population growth and sampling error to evaluate their performance. In some cases, we get nearly unbiased long-term variance estimates despite ignoring sampling error, but the utility of these estimators is questionable because of large estimation uncertainty and difficulties in estimating correlation structure in practice. Process variance estimators that ignored temporal correlations generally gave more precise estimates of the variability in population growth and of the probability of quasi-extinction. We also found that the estimation of probability of quasi-extinction was greatly improved when quasi-extinction thresholds were set relatively close to population levels. Because of precision concerns, we recommend using simple models for risk estimates despite potential biases, and limiting inference to quantifying relative risk; e.g., changes in risk over time for a single population or comparative risk among populations.     

Adaptive cluster sampling with clusters selected without replacement and stopping rule

Adaptive cluster sampling (ACS) has received much attention in recent years since it yields more precise estimates than conventional sampling designs when applied to rare and clustered populations. These results, however, are impacted by the availability of some prior knowledge about the spatial distribution and the absolute abundance of the population under study. This prior information helps the researcher to select a suitable critical value that triggers the adaptive search, the neighborhood definition and the initial sample size. A bad setting of the ACS design would worsen the performance of the adaptive estimators. In particular, one of the greatest weaknesses in ACS is the inability to control the final sampling effort if, for example, the critical value is set too low. To overcome this drawback one can introduce ACS with clusters selected without replacement where one can fix in advance the number of distinct clusters to be selected or ACS with a stopping rule which stops the adaptive sampling when a predetermined sample size limit is reached or when a given stopping rule is verified. However, the stopping rule breaks down the theoretical basis for the unbiasedness of the ACS estimators introducing an unknown amount of bias in the estimates. The current study improves the performance of ACS when applied to patchy and clustered but not rare populations and/or less clustered populations. This is done by combining the stopping rule with ACS without replacement of clusters so as to further limit the sampling effort in form of traveling expenses by avoiding repeat observations and by reducing the final sample size. The performance of the proposed design is investigated using simulated and real data.     

Population Viability Analysis with Species Occurrence Data from Museum Collections

Abstract: The most comprehensive data on many species come from scientific collections. Thus, we developed a method of population viability analysis (PVA) in which this type of occurrence data can be used. In contrast to classical PVA, our approach accounts for the inherent observation error in occurrence data and allows the estimation of the population parameters needed for viability analysis. We tested the sensitivity of the approach to spatial resolution of the data, length of the time series, sampling effort, and detection probability with simulated data and conducted PVAs for common, rare, and threatened species. We compared the results of these PVAs with results of standard method PVAs in which observation error is ignored. Our method provided realistic estimates of population growth terms and quasi‐extinction risk in cases in which the standard method without observation error could not. For low values of any of the sampling variables we tested, precision decreased, and in some cases biased estimates resulted. The results of our PVAs with the example species were consistent with information in the literature on these species. Our approach may facilitate PVA for a wide range of species of conservation concern for which demographic data are lacking but occurrence data are readily available.     

Combining multistate capture-recapture data with tag recoveries to estimate demographic parameters

Matrix population models that allow an animal to occupy more than one state over time are important tools for population and evolutionary ecologists. Definition of state can vary, including location for metapopulation models and breeding state for life history models. For populations whose members can be marked and subsequently reencountered, multistate mark-recapture models are available to estimate the survival and transition probabilities needed to construct population models. Multistate models have proved extremely useful in this context, but they often require a substantial amount of data and restrict estimation of transition probabilities to those areas or states subjected to formal sampling effort. At the same time, for many species, there are considerable tag recovery data provided by the public that could be modeled in order to increase precision and to extend inference to a greater number of areas or states. Here we present a statistical model for combining multistate capture-recapture data (e.g., from a breeding ground study) with multistate tag recovery data (e.g., from wintering grounds). We use this method to analyze data from a study of Canada Geese (Branta canadensis) in the Atlantic Flyway of North America. Our analysis produced marginal improvement in precision, due to relatively few recoveries, but we demonstrate how precision could be further improved with increases in the probability that a retrieved tag is reported.     

Toward estimation of map accuracy without a probability test sample

The time and effort required of probability sampling for accuracy assessment of large-scale land cover maps often means that probability test samples are not collected. Yet, map usefulness is substantially reduced without reliable accuracy estimates. In this article, we introduce a method of estimating the accuracy of a classified map that does not utilize a test sample in the usual sense, but instead estimates the probability of correct classification for each map unit using only the classification rule and the map unit covariates. We argue that the method is an improvement over conventional estimators, though it does not eliminate the need for probability sampling. The method also provides a new and simple method of constructing accuracy maps. We illustrate some of problems associated with accuracy assessment of broad-scale land cover maps, and our method, with a set of nine Landsat Thematic Mapper satellite image-based land cover maps from Montana and Wyoming, USA.     

Repeated geographic divergence in behavior: a case study employing phenotypic trajectory analyses

Environmental effects on behavior have long been a focus of behavioral ecologists. Among the important drivers of behavior is predation environment, which can include the presence/absence of predators, differences in resource availability, and variation in individual density. Environments with predators are often more ecologically complex and “risky” than those without predators. Populations from these environments are sometimes more active and explorative than populations from low-risk, less complex environments. To date, most comparative studies of behavior are limited to within-species comparisons of populations from divergent environments, but neglect comparisons between species following speciation, thus limiting our understanding of post-speciation behavioral evolution. Brachyrhaphis fishes provide an ideal system for studying correlations between divergent environments and behavior within and between species. Here, we test for differences in two behavioral traits—activity and exploration —between sister species Brachyrhaphis roseni and Brachyrhaphis terrabensis that occur in divergent predation environments. Species differed in activity and exploration, with higher activity and exploration levels in populations that co-occur with predators. Furthermore, we found drainage-by-species interactions, indicating that the nature of divergence varied geographically. Using the recently developed phenotypic trajectory analysis (PTA), we quantified this difference and found that, while the geographically isolated populations of sister species tended to evolve in parallel, the magnitude of divergence between species differed between drainages. Our results highlight the utility of PTA for multivariate behavioral data and corroborate past predictions that complex and risky environments are correlated with increased activity and exploration levels and that divergence continues post-speciation.     

Development of the gap model ZELIG-CFS to predict the dynamics of North American mixed forest types with complex structures

When the development of gap models began about three decades ago, they became a new category of forest productivity models. Compared with traditional growth and yield models, which aim at deriving empirical relationships that best fit data, gap models use semi-theoretical relationships to simulate biotic and abiotic processes in forest stands, including the effects of photosynthetic active radiation interception, site fertility, temperature and soil moisture on tree growth and seedling establishment. While growth and yield models are appropriate to predict short-term stemwood production, gap models may be used to predict the natural course of species replacement for several generations. Because of the poor availability of historical data and knowledge on species-specific allometric relationships, species replacement and death rate, it has seldom been possible to develop and evaluate the most representative algorithms to predict growth and mortality with a high degree of accuracy. For this reason, the developers of gap models focused more on developing simulation tools to improve the understanding of forest succession than predicting growth and yield accurately.In a previous study, the predictions of simulations in two southeastern Canadian mixed ecosystem types using the ZELIG gap model were compared with long-term historical data. This exercise highlighted model components that needed modifications to improve the predictive capacity of ZELIG. The updated version of the model, ZELIG-CFS, includes modifications in the modelling of crown interaction effects, survival rate and regeneration. Different algorithms representing crown interactive effects between crowns were evaluated and species-specific model components that compute individual-tree mortality probability rate were derived. The results of the simulations were compared using long-term remeasurement data obtained from sample plots located in La Mauricie National Park of Canada in Quebec. In the present study, three forest types were studied: (1) red spruce-balsam fir-yellow birch, (2) yellow birch-sugar maple-balsam fir, and (3) red spruce-balsam fir-white birch mixed ecosystems. Among the seven algorithms that represented individual crown interactions, two better predicted the changes in basal area and individual-tree growth: (1) the mean available light growing factor (ALGF), which is computed from the proportion of light intercepted at different levels of individual crowns adjusted by the species-specific shade tolerance index, and (2) the ratio of mean ALGF to crown width. The long-term predicted patterns of change in basal area were consistent with the life history of the different species.     

Inbreeding in a lek-mating ant species, Pogonomyrmex occidentalis

In this paper we have two goals. First, we examine the effects of sample size on the statistical power to detect a given amount of inbreeding in social insect populations. The statistical power to detect a given level of inbreeding is largely a function of the number of colonies sampled. We explore two sampling schemes, one in which a single individual per colony is sampled for different sample sizes and a second sampling scheme in which constant sampling effort is maintained (the product of the number of colonies and the number of workers per colony is constant). We find that adding additional workers to a sample from a colony makes it easier to detect inbreeding in samples from given number of colonies; however, adding more colonies rather than more workers per colony always gives greater power to detect inbreeding. Because even relatively large amounts of sib-mating generate relatively small inbreeding coefficients, detection of even substantial deviations from random mating will require very large samples. Second, we look at the amount of inbreeding in a large population of the western harvest ant, Pogonomyrmex occidentalis. We find deviations from Hardy-Weinberg equilibrium equivalent to approximately 27% sib-mating in our population ( f = 0.09). Review of past studies on the population structure of other Pogonomyrmex species suggests that inbreeding may be a regular feature of the mating system of these ants. Although P. occidentalisis a swarm-mating species, there are a number of features of its population biology which suggest that the effective population size may be small. These include topographical variation that potentially breaks the population into demes, variation in the reproductive output of colonies, and variation in the size of reproductives produced by colonies. Received: 6 May 1996 / Accepted after revision: 6 October 1996