Joint estimation of habitat dynamics and species interactions: disturbance reduces co-occurrence of non-native predators with an endangered toad.

1. Ecologists have long been interested in the processes that determine patterns of species occurrence and co-occurrence. Potential short-comings of many existing empirical approaches that address these questions include a reliance on patterns of occurrence at a single time point, failure to account properly for imperfect detection and treating the environment as a static variable. 2. We fit detection and non-detection data collected from repeat visits using a dynamic site occupancy model that simultaneously accounts for the temporal dynamics of a focal prey species, its predators and its habitat. Our objective was to determine how disturbance and species interactions affect the co-occurrence probabilities of an endangered toad and recently introduced non-native predators in stream breeding habitats. For this, we determined statistical support for alternative processes that could affect co-occurrence frequency in the system. 3. We collected occurrence data at stream segments in two watersheds where streams were largely ephemeral and one watershed dominated by perennial streams. Co-occurrence probabilities of toads with non-native predators were related to disturbance frequency, with low co-occurrence in the ephemeral watershed and high co-occurrence in the perennial watershed. This occurred because once predators were established at a site, they were rarely lost from the site except in cases when the site dried out. Once dry sites became suitable again, toads colonized them much more rapidly than predators, creating a period of predator-free space. 4. We attribute the dynamics to a storage effect, where toads persisting outside the stream environment during periods of drought rapidly colonized sites when they become suitable again. Our results support that even in highly connected stream networks, temporal disturbance can structure frequencies with which breeding amphibians encounter non-native predators. 5. Dynamic multi-state occupancy models are a powerful tool for rigorously examining hypotheses about inter-species and species-habitat interactions. In contrast to previous methods that infer dynamic processes based on static patterns in occupancy, the approach we took allows the dynamic processes that determine species-species and species-habitat interactions to be directly estimated.

[1]  Larissa L. Bailey,et al.  Modeling co-occurrence of northern spotted and barred owls: Accounting for detection probability differences , 2009 .

[2]  J. Lawton,et al.  Enemy free space and the structure of ecological communities , 1984 .

[3]  Joseph H. Connell,et al.  Species coexistence, keystone species, and succession: a sensitivity analysis , 1994 .

[4]  J. Wootton,et al.  PREDICTION IN COMPLEX COMMUNITIES: ANALYSIS OF EMPIRICALLY DERIVED MARKOV MODELS , 2001 .

[5]  Joel E. Cohen,et al.  Communities in Patchy Environments: A Model of Disturbance, Competition, and Heterogeneity , 1991 .

[6]  G. Fellers,et al.  Collapse of a Regional Frog Fauna in the Yosemite Area of the California Sierra Nevada, USA , 1996 .

[7]  N. Gotelli Null model analysis of species co-occurrence patterns , 2000 .

[8]  J. Connell Diversity in tropical rain forests and coral reefs. , 1978, Science.

[9]  J. C. Grubb Differential Predation by Gambusia affinis on the Eggs of Seven Species of Anuran Amphibians , 1972 .

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

[11]  L. Kats,et al.  Alien predators and amphibian declines: review of two decades of science and the transition to conservation , 2003 .

[12]  David K. Skelly,et al.  Pond Drying, Predators, and the Distribution of Pseudacris Tadpoles , 1996 .

[13]  H. Wilbur,et al.  Regulation of Structure in Complex Systems: Experimental Temporary Pond Communities , 1987 .

[14]  Colin R. Townsend,et al.  The Patch Dynamics Concept of Stream Community Ecology , 1989, Journal of the North American Benthological Society.

[15]  J. Burger,et al.  Relationships among Isolated Wetland Size, Hydroperiod, and Amphibian Species Richness: Implications for Wetland Regulations , 2000 .

[16]  M. V. Jiménez-Franco,et al.  Territorial occupancy dynamics in a forest raptor community , 2011, Oecologia.

[17]  Stephen H. Roxburgh,et al.  THE INTERMEDIATE DISTURBANCE HYPOTHESIS: PATCH DYNAMICS AND MECHANISMS OF SPECIES COEXISTENCE , 2004 .

[18]  Julien Martin,et al.  Perturbation analysis for patch occupancy dynamics. , 2009, Ecology.

[19]  David A W Miller,et al.  General methods for sensitivity analysis of equilibrium dynamics in patch occupancy models. , 2012, Ecology.

[20]  P. Waggoner,et al.  Transition Probabilities for a Forest , 1970, Nature.

[21]  H. Caswell,et al.  Markov Chain Analysis of Succession in a Rocky Subtidal Community , 2004, The American Naturalist.

[22]  T. Beebee,et al.  The amphibian decline crisis: A watershed for conservation biology? , 2005 .

[23]  M. Huston A General Hypothesis of Species Diversity , 1979, The American Naturalist.

[24]  Andrew M Hein,et al.  Predators, prey, and transient states in the assembly of spatially structured communities. , 2011, Ecology.

[25]  James D. Nichols,et al.  Exploring sensitivity of a multistate occupancy model to inform management decisions , 2011 .

[26]  Terrestrial habitat preferences of adult arroyo southwestern toads , 2001 .

[27]  P. Moyle Effects of introduced bullfrogs, Rana catesbeiana, on the native frogs of the San Joaquin Valley, California , 1973 .

[28]  S. Walde Effect of an abiotic disturbance on a lotic predator-prey interaction , 1986, Oecologia.

[29]  R. Nisbet,et al.  Bullfrogs, disturbance regimes, and the persistence of California red-legged frogs , 2003 .

[30]  L. Kats,et al.  Aggression by Non‐Native Crayfish Deters Breeding in California Newts , 1997 .

[31]  Simon Chamaillé-Jammes,et al.  Simultaneous modeling of habitat suitability, occupancy, and relative abundance: African elephants in Zimbabwe. , 2010, Ecological applications : a publication of the Ecological Society of America.

[32]  David R. Anderson,et al.  Model selection and multimodel inference : a practical information-theoretic approach , 2003 .

[33]  James B. Johnson,et al.  Accelerated Hatching of Southern Leopard Frog (Rana sphenocephala) Eggs in Response to the Presence of a Crayfish (Procambarus nigrocinctus) Predator , 2003, Copeia.

[34]  E. Fleishman,et al.  The relationship between wetland hydroperiod and nestedness patterns in assemblages of larval amphibians and predatory macroinvertebrates , 2004 .

[35]  Robert N. Fisher,et al.  The Decline of Amphibians in California’s Great Central Valley , 1996 .

[36]  J. Nichols,et al.  An integrated model of habitat and species occurrence dynamics , 2011 .

[37]  L. Kats,et al.  Effects of Urbanization on the Distribution and Abundance of Amphibians and Invasive Species in Southern California Streams , 2005 .

[38]  M. J. Adams,et al.  POND PERMANENCE AND THE EFFECTS OF EXOTIC VERTEBRATES ON ANURANS , 2000 .

[39]  J. Nichols,et al.  Investigating species co-occurrence patterns when species are detected imperfectly , 2004 .

[40]  Viviana Ruiz-Gutierrez,et al.  Occupancy dynamics in a tropical bird community: unexpectedly high forest use by birds classified as non‐forest species , 2010 .

[41]  E. Grant Structural complexity, movement bias, and metapopulation extinction risk in dendritic ecological networks , 2011 .

[42]  Darryl I MacKenzie,et al.  Modeling species occurrence dynamics with multiple states and imperfect detection. , 2009, Ecology.