Estimating reference fishing mortality rates from noisy spawnerrecruit data
We review and evaluate methods of estimating reference fishing mortality rates from spawnerrecruit (SR) data to obtain maximum sustainable yield. Using Monte Carlo simulations, we found that a reference fishing mortality rate derived from the maximum likelihood estimates of the SR parameters was less biased than reference fishing mortality rates obtained using the mode of the marginal probability distribution for the maximum rate that spawners produce recruits or by finding the fishing mortality rate that maximizes the expected yield. However, the maximum likelihood method produced the most variable estimates, at times leading to substantial under- or over-exploitation of the population. In contrast, the decision theoretic method of maximizing the expected yield exhibited less variability, produced higher yields, and substantially reduced the risk of overexploiting the population. We show how these methods can be extended to include information from other populations. Bayesian priors for the SR parameters, obtained through meta-analyses of population dynamics at some higher organizational level (e.g., the species), may be used to assess the plausibility of parameter estimates obtained for a single population or combined with the data for the population of interest. Reference fishing mortality rates are then estimated from the resulting joint posterior distribution.