A model for behavioral regulation of metapopulation dynamics

Habitat fragmentation can decrease local population persistence by reducing connectivity, which is a function of dispersal of individuals among habitat fragments. Dispersal is often treated as diffusion in population models, even though for many species it is a result of a series of behavioral decisions. We developed a metapopulation model to explore the potential importance of dispersal behaviors in driving metapopulation dynamics. We incorporated types of behavior that affect dispersal—colonization inhibiting, colonization enhancing, extinction inhibiting, extinction enhancing, rescue enhancing, rescue inhibiting—into Levins’ (1969) metapopulation model and projected occupancy rates for a variety of parameter values. Examples from the literature of behaviors associated with each of these parameters are provided. Our model simplifies into previously published metapopulation models that incorporate only a single behavior, and we present a density-dependent rescue function that leads to multiple non-zero equilibria. We found a variety of behavioral effects on metapopulations. Rescue enhancement fills patches faster than does colonization enhancement or extinction inhibition, and declines in patch occupancy are moderate with extinction enhancement, but colonization inhibition causes metapopulation extinction. We also found that with colonization and extinction inhibitions, equilibrium patch occupancy is inversely related to patch turnover rate. With density-dependent rescue, persistence depends not only on the strength of the strong rescue effect, but also on having a sufficient initial fraction of patches occupied; the stronger the rescue effect, the lower this fraction can be. This study suggests that dispersal behavior can have strong influences on metapopulation dynamics. It confirms the importance of understanding the relationship between landscape structure and dispersal behavior in understanding population persistence.     

A robust-design formulation of the incidence function model of metapopulation dynamics applied to two species of rails

The incidence function model (IFM) uses area and connectivity to predict metapopulation dynamics. However, false absences and missing data can lead to underestimates of the number of sites contributing to connectivity, resulting in overestimates of dispersal ability and turnovers (extinctions plus colonizations). We extend estimation methods for the IFM by using a hierarchical Bayesian model to account both for false absences due to imperfect detection and for missing data due to sites not surveyed in some years. We compare parameter estimates, measures of metapopulation dynamics, and forecasts using stochastic patch occupancy models (SPOMs) among three IFM models: (1) a Bayesian formulation assuming no false absences and omitting site-year combinations with missing data; (2) a hierarchical Bayesian formulation assuming no false absences but incorporating missing data; and (3) a hierarchical Bayesian formulation allowing for imperfect detection and incorporating missing data. We fit the models to multiyear data sets of occupancy for two bird species that differ in body size and presumed dispersal ability but inhabit the same network of sites: the small Black Rail (Laterallus jamaicensis) and the medium-sized Virginia Rail (Rallus limicola). Incorporating missing data affected colonization parameters and led to lower estimates of dispersal ability for the Black Rail. Detection rates were high for the Black Rail in most years but moderate for the Virginia Rail. Incorporating imperfect detection resulted in higher occupancy and lower turnover rates for both species, with largest effects for the Virginia Rail. Forecasts using SPOMs were sensitive to both missing data and false absences; persistence in models assuming no false absences was more optimistic than from robust models. Our results suggest that incorporating false absences and missing data into the IFM can improve (1) estimates of dispersal ability and the effect of connectivity on colonization, (2) the scaling of extinction risk with patch area, and (3) forecasts of occupancy and turnover rates.     

Neighborhood and habitat effects on vital rates: expansion of the Barred Owl in the Oregon coast ranges

In this paper, we modify dynamic occupancy models developed for detection-nondetection data to allow for the dependence of local vital rates on neighborhood occupancy, where neighborhood is defined very flexibly. Such dependence of occupancy dynamics on the status of a relevant neighborhood is pervasive, yet frequently ignored. Our framework permits joint inference about the importance of neighborhood effects and habitat covariates in determining colonization and extinction rates. Our specific motivation is the recent expansion of the Barred Owl (Strix varia) in western Oregon, USA, over the period 1990-2010. Because the focal period was one of dramatic range expansion and local population increase, the use of models that incorporate regional occupancy (sources of colonists) as determinants of dynamic rate parameters is especially appropriate. We began our analysis of 21 years of Barred Owl presence/nondetection data in the Tyee Density Study Area (TDSA) by testing a suite of six models that varied only in the covariates included in the modeling of detection probability. We then tested whether models that used regional occupancy as a covariate for colonization and extinction outperformed models with constant or year-specific colonization or extinction rates. Finally we tested whether habitat covariates improved the AIC of our models, focusing on which habitat covariates performed best, and whether the signs of habitat effects are consistent with a priori hypotheses. We conclude that all covariates used to model detection probability lead to improved AIC, that regional occupancy influences colonization and extinction rates, and that habitat plays an important role in determining extinction and colonization rates. As occupancy increases from low levels toward equilibrium, colonization increases and extinction decreases, presumably because there are more and more dispersing juveniles. While both rates are affected, colonization increases more than extinction decreases. Colonization is higher and extinction is lower in survey polygons with more riparian forest. The effects of riparian forest on extinction rates are greater than on colonization rates. Model results have implications for management of the invading Barred Owl, both through habitat alteration and removal.     

Metapopulation dynamics in the butterfly Hipparchia semele changed decades before occupancy declined in The Netherlands

The survival of many species in human-dominated, fragmented landscapes depends on metapopulation dynamics, i.e., on a dynamic equilibrium of extinctions and colonizations in patches of suitable habitat. To understand and predict distributional changes, knowledge of these dynamics can be essential, and for this, metapopulation studies are preferably based on long-time-series data from many sites. Alas, such data are very scarce. An alternative is to use opportunistic data (i.e., collected without applying standardized field methods), but these data suffer from large variations in field methods and search intensity between sites and years. Dynamic site-occupancy models offer a general approach to adjust for variable survey effort. These models extend classical metapopulation models to account for imperfect detection of species and yield estimates of the probabilities of occupancy, colonization, and survival of species at sites. By accounting for detection, they fully correct for among-year variability in search effort. As an illustration, we fitted a dynamic site-occupancy model to 60 years of presence-absence data (more precisely, detection-nondetection) of the heathland butterfly Hipparchia semele in The Netherlands. Detection records were obtained from a database containing volunteer-based data from 1950-2009, and nondetection records were deduced from database records of other butterfly species. Our model revealed that metapopulation dynamics of Hipparchia had changed decades before the species' distribution began to contract. Colonization probability had already started to decline from 1950 onward, but this was counterbalanced by an increase in the survival of existing populations, the result of which was a stable distribution. Only from 1990 onward was survival not sufficient to compensate for the further decrease of colonization, and occupancy started to decline. Hence, it appears that factors acting many decades ago triggered a change in the metapopulation dynamics of this species, which ultimately led to a severe decline in occupancy that only became apparent much later. Our study emphasizes the importance of knowledge of changes in survival and colonization of species in modern landscapes over a very long time scale. It also demonstrates the power of site-occupancy modeling to obtain important population dynamics information from databases containing opportunistic sighting records.     

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.     

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.     

Minimum viable metapopulation size, extinction debt, and the conservation of a declining species.

A key question facing conservation biologists is whether declines in species' distributions are keeping pace with landscape change, or whether current distributions overestimate probabilities of future persistence. We use metapopulations of the marsh fritillary butterfly Euphydryas aurinia in the United Kingdom as a model system to test for extinction debt in a declining species. We derive parameters for a metapopulation model (incidence function model, IFM) using information from a 625-km2 landscape where habitat patch occupancy, colonization, and extinction rates for E. aurinia depend on patch connectivity, area, and quality. We then show that habitat networks in six extant metapopulations in 16-km2 squares were larger, had longer modeled persistence times (using IFM), and higher metapopulation capacity (lambdaM) than six extinct metapopulations. However, there was a > 99% chance that one or more of the six extant metapopulations would go extinct in 100 years in the absence of further habitat loss. For 11 out of 12 networks, minimum areas of habitat needed for 95% persistence of metapopulation simulations after 100 years ranged from 80 to 142 ha (approximately 5-9% of land area), depending on the spatial location of habitat. The area of habitat exceeded the estimated minimum viable metapopulation size (MVM) in only two of the six extant metapopulations, and even then by only 20%. The remaining four extant networks were expected to suffer extinction in 15-126 years. MVM was consistently estimated as approximately 5% of land area based on a sensitivity analysis of IFM parameters and was reduced only marginally (to approximately 4%) by modeling the potential impact of long-distance colonization over wider landscapes. The results suggest a widespread extinction debt among extant metapopulations of a declining species, necessitating conservation management or reserve designation even in apparent strongholds. For threatened species, metapopulation modeling is a potential means to identify landscapes near to extinction thresholds, to which conservation measures can be targeted for the best chance of success.     

Transient dynamics of invasive competition: barred owls, spotted owls, habitat, and the demons of competition present

The recent range expansion of Barred Owls (Strix varia) into the Pacific Northwest, where the species now co-occurs with the endemic Northern Spotted Owl (Strix occidentalis caurina), resulted in a unique opportunity to investigate potential competition between two congeneric, previously allopatric species. The primary criticism of early competition research was the use of current species' distribution patterns to infer past processes; however, the recent expansion of the Barred Owl and the ability to model the processes that result in site occupancy (i.e., colonization and extinction) allowed us to address the competitive process directly rather than inferring past processes through current patterns. The purpose of our study was to determine whether Barred Owls had any negative effects on occupancy dynamics of nesting territories by Northern Spotted Owls and how these effects were influenced by habitat characteristics of Spotted Owl territories. We used single-species, multi-season occupancy models and covariates quantifying Barred Owl detections and habitat characteristics to model extinction and colonization rates of Spotted Owl pairs in southern Oregon, USA. We observed a strong, negative association between Barred Owl detections and colonization rates and a strong positive effect of Barred Owl detections on extinction rates of Spotted Owls. We observed increased extinction rates in response to decreased amounts of old forest at the territory core and higher colonization rates when old-forest habitat was less fragmented. Annual site occupancy for pairs reflected the strong effects of Barred Owls on occupancy dynamics with much lower occupancy rates predicted for territories where Barred Owls were detected. The strong Barred Owl and habitat effects on occupancy dynamics of Spotted Owls provided evidence of interference competition between the species. These effects increase the importance of conserving large amounts of contiguous, old-forest habitat to maintain Northern Spotted Owls in the landscape.     

Colonization and extinction in dynamic habitats: an occupancy approach for a Great Plains stream fish assemblage

Despite the importance of habitat in determining species distribution and persistence, habitat dynamics are rarely modeled in studies of metapopulations. We used an integrated habitat-occupancy model to simultaneously quantify habitat change, site fidelity, and local colonization and extinction rates for larvae of a suite of Great Plains stream fishes in the Arikaree River, eastern Colorado, USA, across three years. Sites were located along a gradient of flow intermittency and groundwater connectivity. Hydrology varied across years: the first and third being relatively wet and the second dry. Despite hydrologic variation, our results indicated that site suitability was random from one year to the next. Occupancy probabilities were also independent of previous habitat and occupancy state for most species, indicating little site fidelity. Climate and groundwater connectivity were important drivers of local extinction and colonization, but the importance of groundwater differed between periods. Across species, site extinction probabilities were highest during the transition from wet to dry conditions (range: 0.52-0.98), and the effect of groundwater was apparent with higher extinction probabilities for sites not fed by groundwater. Colonization probabilities during this period were relatively low for both previously dry sites (range: 0.02-0.38) and previously wet sites (range: 0.02-0.43). In contrast, no sites dried or remained dry during the transition from dry to wet conditions, yielding lower but still substantial extinction probabilities (range: 0.16-0.63) and higher colonization probabilities (range: 0.06-0.86), with little difference among sites with and without groundwater. This approach of jointly modeling both habitat change and species occupancy will likely be useful to incorporate effects of dynamic habitat on metapopulation processes and to better inform appropriate conservation actions.     

Metapopulation dynamics of a perennial plant, Succisa pratensis, in an agricultural landscape

Most metapopulation models neglect the local dynamics, and systems characterized by slow population turnover, time lags and non-equilibrium, are only rarely examined within a metapopulation context. In this study we used a realistic, spatially explicit, dynamic metapopulation model of a long-lived grassland plant, Succisa pratensis, to examine the relative importance of local population dynamics, and short and long-distance dispersal of seeds.     

General methods for sensitivity analysis of equilibrium dynamics in patch occupancy models

Sensitivity analysis is a useful tool for the study of ecological models that has many potential applications for patch occupancy modeling. Drawing from the rich foundation of existing methods for Markov chain models, I demonstrate new methods for sensitivity analysis of the equilibrium state dynamics of occupancy models. Estimates from three previous studies are used to illustrate the utility of the sensitivity calculations: a joint occupancy model for a prey species, its predators, and habitat used by both; occurrence dynamics from a well-known metapopulation study of three butterfly species; and Golden Eagle occupancy and reproductive dynamics. I show how to deal efficiently with multistate models and how to calculate sensitivities involving derived state variables and lower-level parameters. In addition, I extend methods to incorporate environmental variation by allowing for spatial and temporal variability in transition probabilities. The approach used here is concise and general and can fully account for environmental variability in transition parameters. The methods can be used to improve inferences in occupancy studies by quantifying the effects of underlying parameters, aiding prediction of future system states, and identifying priorities for sampling effort.     

Extinction debt of forest plants persists for more than a century following habitat fragmentation

Following habitat fragmentation individual habitat patches may lose species over time as they pay off their "extinction debt." Species with relatively low rates of population extinction and colonization ("slow" species) may maintain extinction debts for particularly prolonged periods, but few data are available to test this prediction. We analyzed two unusually detailed data sets on forest plant distributions and land-use history from Lincolnshire, United Kingdom, and Vlaams-Brabant, Belgium, to test for an extinction debt in relation to species-specific extinction and colonization rates. Logistic regression models predicting the presence-absence of 36 plant species were first parameterized using data from Lincolnshire, where forest cover has been relatively low (approximately 5-8%) for the past 1000 years. Consistent with extinction debt theory, for relatively slow species (but not fast species) these models systematically underpredicted levels of patch occupancy in Vlaams-Brabant, where forest cover was reduced from approximately 25% to <10% between 1775 and 1900 (it is presently 6.5%). As a consequence, the ability of the Lincolnshire models to predict patch occupancy in Vlaams-Brabant was worse for slow than for fast species. Thus, more than a century after forest fragmentation reached its current level an extinction debt persists for species with low rates of population turnover.     

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.     

The role of habitat disturbance and recovery in metapopulation persistence

Classical metapopulation theory assumes a static landscape. However, empirical evidence indicates many metapopulations are driven by habitat succession and disturbance. We develop a stochastic metapopulation model, incorporating habitat disturbance and recovery, coupled with patch colonization and extinction, to investigate the effect of habitat dynamics on persistence. We discover that habitat dynamics play a fundamental role in metapopulation dynamics. The mean number of suitable habitat patches is not adequate for characterizing the dynamics of the metapopulation. For a fixed mean number of suitable patches, we discover that the details of how disturbance affects patches and how patches recover influences metapopulation dynamics in a fundamental way. Moreover, metapopulation persistence is dependent not only on the average lifetime of a patch, but also on the variance in patch lifetime and the synchrony in patch dynamics that results from disturbance. Finally, there is an interaction between the habitat and metapopulation dynamics, for instance declining metapopulations react differently to habitat dynamics than expanding metapopulations. We close, emphasizing the importance of using performance measures appropriate to stochastic systems when evaluating their behavior, such as the probability distribution of the state of the metapopulation, conditional on it being extant (i.e., the quasistationary distribution).     

Epiphyte metapopulation dynamics are explained by species traits, connectivity, and patch dynamics

The colonization-extinction dynamics of many species are affected by the dynamics of their patches. For increasing our understanding of the metapopulation dynamics of sessile species confined to dynamic patches, we fitted a Bayesian incidence function model extended for dynamic landscapes to snapshot data on five epiphytic lichens among 2083 mapped oaks (dynamic patches). We estimate the age at which trees become suitable patches for different species, which defines their niche breadth (number of suitable trees). We show that the colonization rates were generally low, but increased with increasing connectivity in accordance with metapopulation theory. The rates were related to species traits, and we show, for the first time, that they are higher for species with wide niches and small dispersal propagules than for species with narrow niches or large propagules. We also show frequent long-distance dispersal in epiphytes by quantifying the relative importance of local dispersal and background deposition of dispersal propagules. Local stochastic extinctions from intact trees were negligible in all study species, and thus, the extinction rate is set by the rate of patch destruction (tree fall). These findings mean that epiphyte metapopulations may have slow colonization-extinction dynamics that are explained by connectivity, species traits, and patch dynamics.     

Geometrical Factors and Metapopulation Dynamics of the Bush Cricket, Metrioptera bicolor Philippi (Orthoptera: Tettigoniidae)

This paper presents a metapopulation study of the bush cricket, Metrioptera bicolor , living in a recently fragmented landscape. The species inhabits grass and heathland patches of varying area and isolation. Analyses are made of how these geometrical factors affect local population size and density, distribution pattern, and the probability of local extinction and colonization. The proportion of available patches occupied varied between 72 and 79% during 1985–1990. Unoccupied patches were smaller and more isolated than those that were occupied. Patches where populations became extinct during this period were smaller than those with persisting populations. Since local population size was well correlated with patch area, it was clear that stochastic extinctions only occurred in small populations. Critical patch size for population extinction was approximately half a hectare. Colonized patches were less isolated than those that had not been colonized. Critical inter-patch distance for colonization was about 100 meters. The turnover was restricted to an identifiable share of the available patches. Only 33% of the patches were so small that extinction due to stochastic causes could be considered highly probable. This metapopulation will therefore most likely persist over a considerable period in its present spatial structure. There are apparent threats of further fragmentation, however, and nothing is known about the likelihood of large-scale extinctions resulting from extremely unfavorable weather conditions. Nevertheless, our results show that it is appropriate to include geometrical factors in metapopulation models.     

Evaluating the long-term metacommunity dynamics of tree hole mosquitoes

Four different conceptual models of metacommunities have been proposed, termed "patch dynamics," "species sorting", "mass effect", and "neutral". These models simplify thinking about metacommunities and improve our understanding of the role of spatial dynamics both in structuring communities and in determining local and regional diversity. We tested whether mosquito communities inhabiting water-filled tree holes in southeastern Florida, U.S.A., displayed any of the characteristics and dynamics predicted by the four models. The densities of the five most common species in 3-8 tree holes were monitored every two weeks during 1978-2003. We tested relationships between habitat variables and species densities, spatial synchrony, the presence of life history trade-offs, and species turnover. Dynamics showed strong elements of species sorting, but with considerable turnover, as predicted by the patch dynamics model. Consistent with patch dynamics, there was substantial asynchrony in dynamics for different tree holes, substantial species turnover in space and time, and an occupancy/colonization trade-off. Substantial correlations of density and occupancy with tree hole volume were consistent with the species-sorting model, but unlike this model, species did not have permanent refuges. No evidence of mass effects was found, and correlations between habitat variables and dynamics were inconsistent with neutral models. Our results did not match a single model and therefore caution against overly simplifying metacommunity dynamics by using one dynamical characteristic to select a particular metacommunity perspective.     

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.     

Use of fecundity measured directly throughout the breeding season to test a source-sink demographic model

Populations of landbirds (bird species that occupy terrestrial habitats for most of their life cycle) are declining throughout North America (north of Mexico) and Europe, yet little is known about how demography is driving this trend. A recent model of 5 geographically separated populations of Cerulean Warblers (Dendroica cerulea) that was based on within-season sampling of nest survival and fledgling success shows that all populations are sinks (annual reproduction is consistently less than annual adult mortality). I tested this indirect model by directly measuring fecundity (number of female fledglings/female) during the breeding season for 2 years in a Cerulean Warbler population occupying a mature forest in southwestern Michigan (U.S.A.) I determined territories of male birds on the basis of male plumage characters and phases of the nesting cycle (2007) and on uniquely color-banded males (2008). I transferred locations of identified males to topographic maps. I counted all fledglings in territories from May to July each year. The model I tested may apply only to single-brooded species; therefore, I searched the literature to estimate the percentage of single-brooded species in North America. The breeding season of Cerulean Warblers was short- nearly all nests were initiated from mid-May to late June. Nest predation and brood parasitism were primary and rare causes of nest failure, respectively. Significantly fewer Cerulean Warblers fledged from parasitized than from nonparasitized nests. Fledgling survival required to maintain the population size was well above previously published values for Neotropical migrants. Single-brooded species comprise 62% of North American breeding bird species for which the number of broods per year is known; I believe my results may apply to these species. The consistency between identification of populations as sources or sinks on the basis of either model estimates or direct measurements suggests that a demographic model relying on within-season sampling of fecundity is adequate to determine population status of single-brooded avian populations. In addition, on the basis of results of previous studies, annual adult survival rate of the Cerulean Warbler is typical of parulid warblers that are not declining. Thus, low fecundity, here determined with different quantitative methods, can drive status of landbird species with high-observed survival.     

Assessing hypotheses about nesting site occupancy dynamics

Hypotheses about habitat selection developed in the evolutionary ecology framework assume that individuals, under some conditions, select breeding habitat based on expected fitness in different habitat. The relationship between habitat quality and fitness may be reflected by breeding success of individuals, which may in turn be used to assess habitat quality. Habitat quality may also be assessed via local density: if high-quality sites are preferentially used, high density may reflect high-quality habitat. Here we assessed whether site occupancy dynamics vary with site surrogates for habitat quality. We modeled nest site use probability in a seabird subcolony (the Black-legged Kittiwake, Rissa tridactyla) over a 20-year period. We estimated site persistence (an occupied site remains occupied from time t to t+1) and colonization through two subprocesses: first colonization (site creation at the timescale of the study) and recolonization (a site is colonized again after being deserted). Our model explicitly incorporated site-specific and neighboring breeding success and conspecific density in the neighborhood. Our results provided evidence that reproductively "successful" sites have a higher persistence probability than "unsuccessful" ones. Analyses of site fidelity in marked birds and of survival probability showed that high site persistence predominantly reflects site fidelity, not immediate colonization by new owners after emigration or death of previous owners. There is a negative quadratic relationship between local density and persistence probability. First colonization probability decreases with density, whereas recolonization probability is constant. This highlights the importance of distinguishing initial colonization and recolonization to understand site occupancy. All dynamics varied positively with neighboring breeding success. We found evidence of a positive interaction between site-specific and neighboring breeding success. We addressed local population dynamics using a site occupancy approach integrating hypotheses developed in behavioral ecology to account for individual decisions. This allows development of models of population and metapopulation dynamics that explicitly incorporate ecological and evolutionary processes.