Breeding Population Dynamics of Threatened Crawfish Frogs Inform Targets for Habitat Management
We used data from a five-year study at two focal breeding wetlands of a single Crawfish Frog (Rana areolata) population to demonstrate the relative importance and potential of tadpole survival as a target for population management. We concurrently estimated 11 vital rates in the field and used a females-only matrix projection model to estimate elasticities of demographic rates. We then simulated stochastic population dynamics at each wetland with and without immigration to estimate the intrinsic capacity for each focal breeding wetland to sustain a Crawfish Frog breeding population and the likelihood of breeding population persistence at each wetland. Elasticity of tadpole survival was second only to juvenile survival elasticity and 1.34–2.04 times greater than adult survival elasticity. Projections indicated that the Crawfish Frog population was not at risk of extinction but only one breeding site was capable of self-sustaining a breeding population. Because of low tadpole survival, the other breeding site was completely dependent on immigration to persist and was functioning as a population sink. Despite higher variability compared to terrestrial vital rates, larval survival did have a strong effect on population growth. Tadpole survival at the more productive breeding site was density dependent and likely related to wetland vegetation and predator and competitor abundance. Two additional findings were that annual survival of frogs following their first known breeding event was 48% lower compared to survival of frogs that had bred in two or more prior years, and adult temporary emigration from the breeding population was moderately high. Our study demonstrates the benefits of using population models that integrate density-dependent processes, temporary emigration from the breeding population, and state-specific adult survival, to identify larval habitats that function as population sinks and limit current population size and persistence probability. We contend that tadpole survival is an important and feasible habitat management target within broader conservation strategies for Crawfish Frogs and other amphibian species.

(A) Post-birth, females-only, stage matrix model of Crawfish Frogs that accommodates age and state-specific transition probabilities and the “temporary emigration” of individuals that skip breeding or breed at other sites. (B) Projection matrix of the model used to estimate elasticities and sensitivities of vital rates. Parameters in the matrix model: embryonic survival (Se), tadpole survival (St), metamorph and juvenile survival (Sj), adult survival following the first breeding event (Sa1), adult survival following two or more breeding events (Sa2), probability that a one-year-old juvenile matures at age 2 (p), probability that a two-year-old juvenile matures at age 3 (q), the probability that a mature individual emigrates to become a non-breeding individual (γ˝), probability that a non-breeding individual remains a non-breeding individual (γ'), and clutch size multiplied by 0.5 to represent half of all eggs developing into females (C). We assumed that all individuals mature by age 4. Individuals that emigrate to become a non-breeding individual may either skip reproduction or breed at another site.

Percent change in the instantaneous population growth rate (λ) in response to a percent change in female Crawfish Frog tadpole, juvenile, or adult survival rates for a population located in southwest Indiana, USA, between 2009–2013. Percent change in each vital rate represents the general range of that vital rate observed in this study. Change was estimated by varying the focal rate while holding all other rates constant.

Correlation between tadpole survival and the instantaneous population growth rate (λ) for 1,000 points drawn randomly from 1,000 simulations of stochastic population simulations of Nate's Pond with (left panel) and without (right panel) immigration included in the simulation model. The white line is a third order polynomial function fitted to the data, and the gray band represents the 95% confidence interval. The horizontal dashed line indicates a λ of 1. The distributions of tadpole survival rates and λ values are provided as marginal plots.
Contributor Notes
Associate Editor: J. Lamb.