Background: Long-lived marine megavertebrates (e.g. sharks, turtles, mammals, and seabirds) are inherently vulnerable to anthropogenic mortality. Although some mathematical models have been applied successfully to manage these animals, more detailed treatments are often needed to assess potential drivers of population dynamics. In particular, factors such as age-structure, density-dependent feedbacks on reproduction, and demographic stochasticity are important for understanding population trends, but are often difficult to assess. Lemon sharks (Negaprion brevirostris) have a pelagic adult phase that makes them logistically difficult to study. However, juveniles use coastal nursery areas where their densities can be high.

Results: We use a stage-structured, Markov-chain stochastic model to describe lemon shark population dynamics from a 17-year longitudinal dataset at a coastal nursery area at Bimini, Bahamas. We found that the interaction between delayed breeding, density-dependence, and demographic stochasticity accounts for 33 to 49% of the variance in population size.

Conclusions: Demographic stochasticity contributed all random effects in this model, suggesting that the existence of unmodeled environmental factors may be driving the majority of interannual population fluctuations. In addition, we are able to use our model to estimate the natural mortality rate of older age classes of lemon sharks that are difficult to study. Further, we use our model to examine what effect the length of a time series plays on deciphering ecological patterns. We find that—even with a relatively long time series—our sampling still misses important rare events. Our approach can be used more broadly to infer population dynamics of other large vertebrates in which age structure and demographic stochasticity are important.