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.

[1]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[2]  Atte Moilanen,et al.  PATCH OCCUPANCY MODELS OF METAPOPULATION DYNAMICS: EFFICIENT PARAMETER ESTIMATION USING IMPLICIT STATISTICAL INFERENCE , 1999 .

[3]  Stephen R. Baillie,et al.  Patterns of natal and breeding dispersal in birds , 1998 .

[4]  Orien M. W. Richmond,et al.  Distribution of California Black Rails in the Sierra Nevada foothills , 2008 .

[5]  T. Schoener.  Sizes of Feeding Territories among Birds , 1968 .

[6]  Ilkka Hanski,et al.  Long‐Term Dynamics in a Metapopulation of the American Pika , 1998, The American Naturalist.

[7]  J. Gamarra,et al.  Metapopulation Ecology , 2007 .

[8]  Bradley P. Carlin,et al.  Markov Chain Monte Carlo conver-gence diagnostics: a comparative review , 1996 .

[9]  J. Andrew Royle,et al.  Hierarchical Modeling and Inference in Ecology: The Analysis of Data from Populations, Metapopulations and Communities , 2008 .

[10]  P. Valpine.  Monte Carlo State-Space Likelihoods by Weighted Posterior Kernel Density Estimation , 2004 .

[11]  Benjamin B. Risk,et al.  California black rails depend on irrigation-fed wetlands in the Sierra Nevada foothills , 2010 .

[12]  I. Hanski.  A Practical Model of Metapopulation Dynamics , 1994 .

[13]  M. Legare,et al.  HOME RANGE SIZE, NEST-SITE SELECTION AND NESTING SUCCESS OF BLACK RAILS IN FLORIDA , 2001 .

[14]  Rampal S. Etienne,et al.  IMPROVED BAYESIAN ANALYSIS OF METAPOPULATION DATA WITH AN APPLICATION TO A TREE FROG METAPOPULATION , 2003 .

[15]  Marc Kéry,et al.  Hierarchical modelling and estimation of abundance and population trends in metapopulation designs. , 2010, The Journal of animal ecology.

[16]  Hannu Toivonen,et al.  BAYESIAN ANALYSIS OF METAPOPULATION DATA , 2002 .

[17]  K. Pollock.  A Capture-Recapture Design Robust to Unequal Probability of Capture , 1982 .

[18]  Atte Moilanen,et al.  Implications of empirical data quality to metapopulation model parameter estimation and application , 2002 .

[19]  E. Fleishman,et al.  An empirical evaluation of the area and isolation paradigm of metapopulation dynamics , 2007 .

[20]  R. Swihart,et al.  Modeling patch occupancy: Relative performance of ecologically scaled landscape indices , 2008, Landscape Ecology.

[21]  D. MacKenzie.  Occupancy Estimation and Modeling: Inferring Patterns and Dynamics of Species Occurrence , 2005 .

[22]  Atte Moilanen,et al.  SPOMSIM: software for stochastic patch occupancy models of metapopulation dynamics , 2004 .

[23]  Peter Green,et al.  Markov chain Monte Carlo in Practice , 1996 .

[24]  Robert J. Wilson,et al.  Minimum viable metapopulation size, extinction debt, and the conservation of a declining species. , 2007, Ecological applications : a publication of the Ecological Society of America.

[25]  D. Blumstein,et al.  Effects of patch quality and network structure on patch occupancy dynamics of a yellow-bellied marmot metapopulation. , 2006, The Journal of animal ecology.

[26]  J. L. Parra,et al.  Impact of a Century of Climate Change on Small-Mammal Communities in Yosemite National Park, USA , 2008, Science.

[27]  J. Nichols,et al.  ESTIMATING SITE OCCUPANCY, COLONIZATION, AND LOCAL EXTINCTION WHEN A SPECIES IS DETECTED IMPERFECTLY , 2003 .

[28]  B. Taylor,et al.  Rails: A Guide to the Rails, Crakes, Gallinules and Coots of the World , 1998 .

[29]  Karin Frank,et al.  RANKING METAPOPULATION EXTINCTION RISK: FROM PATTERNS IN DATA TO CONSERVATION MANAGEMENT DECISIONS , 2003 .

[30]  J. Hines,et al.  Two-species occupancy models: a new parameterization applied to co-occurrence of secretive rails. , 2010, Ecological applications : a publication of the Ecological Society of America.

[31]  I. Hanski.  Metapopulation dynamics , 1998, Nature.

[32]  K. E. Hodges,et al.  Effect of habitat area and isolation on fragmented animal populations , 2008, Proceedings of the National Academy of Sciences.

[33]  J. Andrew Royle,et al.  A Bayesian state-space formulation of dynamic occupancy models. , 2007, Ecology.

[34]  A. F. O'connell,et al.  Multi-scale occupancy estimation and modelling using multiple detection methods , 2008 .

[35]  Brian Dennis,et al.  Data cloning: easy maximum likelihood estimation for complex ecological models using Bayesian Markov chain Monte Carlo methods. , 2007, Ecology letters.