JABBA-Select: Incorporating life history and fisheries’ selectivity into surplus production models

Abstract Age-structured production models (ASPMs) are often preferred over biomass aggregated surplus production models (SPMs) as the former can track the propagation of cohorts and explicitly account for the effects of selective fishing, even in the absence of reliable size- or age data. Here, we introduce ‘JABBA-Select’, an extension of the JABBA software (Just Another Bayesian Biomass Assessment; Winker et al., 2018), that is able to overcome some of the shortcomings of conventional SPMs and allows a direct comparison to ASPMs. JABBA-Select incorporates life history parameters and fishing selectivity and distinguishes between exploitable biomass (used to fit indices given fishery selectivity) and spawning biomass (used to predict surplus production). Applying JABBA-Select involves using an age-structured equilibrium model to convert the input parameters into multivariate normal priors for surplus-production productivity parameters. We illustrate the main elements of JABBA-Select using the stock parameters of South African silver kob (Argyrosomus inodorus, Scienidae) as a case study. This species is exploited by multiple fisheries and was selected as an example of a data moderate fishery that features strong contrast in selectivity over time and across fleets. For proof-of-concept, we use an age-structured simulation framework to compare the performance of JABBA-Select to: 1) a conventional Bayesian state-space model using a Pella-Tomlinson (PT) production function, 2) an ASPM with deterministic recruitment; and 3) an ASPM with stochastic recruitment. The PT model produced biased estimates of relative and absolute spawning biomass trajectories and associated reference points, which by contrast could be fairly accurately estimated by JABBA-Select. JABBA-Select performed at least as well as the ASPMs in accuracy for most of the performance metrics and best characterized the stock status uncertainty. The results indicate that JABBA-Select is able to accurately account for moderate changes in selectivity and fleet dynamics over time and to provide a robust tool for data-moderate stock assessments.

[1]  M. B. Schaefer,et al.  A study of the dynamics of the fishery for yellowfin tuna in the Eastern Tropical Pacific Ocean , 1957 .

[2]  Philip Heidelberger,et al.  Simulation Run Length Control in the Presence of an Initial Transient , 1983, Oper. Res..

[3]  André E. Punt,et al.  Model performance analysis for Bayesian biomass dynamics models using bias, precision and reliability metrics , 2012 .

[4]  K. F. Johnson,et al.  Can autocorrelated recruitment be estimated using integrated assessment models and how does it affect population forecasts , 2016 .

[5]  Patrick J. Farrell,et al.  Local influence diagnostics for the retrospective problem in sequential population analysis , 2005 .

[6]  Richard D. Methot,et al.  Stock synthesis: A biological and statistical framework for fish stock assessment and fishery management , 2013 .

[7]  J. Thorson,et al.  Hierarchical analysis of taxonomic variation in intraspecific competition across fish species. , 2016, Ecology.

[8]  Adrian R. Hordyk,et al.  A simple length-structured model based on life history ratios and incorporating size-dependent selectivity: application to spawning potential ratios for data-poor stocks , 2016 .

[9]  David C. Smith,et al.  Selecting relative abundance proxies for BMSY and BMEY , 2014 .

[10]  Kevin R. Piner,et al.  Can diagnostic tests help identify model misspecification in integrated stock assessments , 2017 .

[11]  Toshihide Kitakado,et al.  Evaluation of the prediction skill of stock assessment using hindcasting , 2016 .

[12]  John Shepherd,et al.  An Alternative Perspective on Recruitment Overfishing and Biological Reference Points , 1987 .

[13]  Felipe Carvalho,et al.  JABBA: Just Another Bayesian Biomass Assessment , 2018, Fisheries Research.

[14]  R. Deriso.  Harvesting Strategies and Parameter Estimation for an Age-Structured Model , 1980 .

[15]  C. Legault,et al.  A Stochastic Implementation of an Age-Structured Production Model , 1998 .

[16]  F. Carvalho,et al.  Incorporating specific change points in catchability in fisheries stock assessment models: An alternative approach applied to the blue shark (Prionace glauca) stock in the south Atlantic Ocean , 2014 .

[17]  James N. Ianelli,et al.  A Statistical Model for Stock Assessment of Southern Bluefin Tuna with Temporal Changes in Selectivity , 2003 .

[18]  Richard D. Methot,et al.  A review of stock assessment packages in the United States , 2016 .

[19]  J. Schnute.  A General Theory for Analysis of Catch and Effort Data , 1985 .

[20]  Stephan B Munch,et al.  Predicting life history parameters for all fishes worldwide. , 2017, Ecological applications : a publication of the Ecological Society of America.

[21]  Athanassios C. Tsikliras,et al.  On the pile‐up effect and priors for Linf and M/K: response to a comment by Hordyk et al. on “A new approach for estimating stock status from length frequency data” , 2018, ICES Journal of Marine Science.

[22]  Kotaro Ono,et al.  A Bayesian approach to identifying and compensating for model misspecification in population models. , 2014, Ecology.

[23]  A. Punt,et al.  Extending integrated stock assessment models to use non-depensatory three-parameter stock-recruitment relationships , 2017, Fisheries Research.

[24]  Cody S. Szuwalski,et al.  How well can FMSY and BMSY be estimated using empirical measures of surplus production , 2012 .

[25]  James T. Thorson,et al.  Catch curve stock-reduction analysis: An alternative solution to the catch equations , 2015 .

[26]  FrancisR.I.C. Chris,et al.  Data weighting in statistical fisheries stock assessment models , 2011 .

[27]  André E. Punt,et al.  Extending production models to include process error in the population dynamics , 2003 .

[28]  J. Hoenig,et al.  Comparative Performance of Three Length-Based Mortality Estimators , 2018, Marine and Coastal Fisheries.

[29]  Michael H. Prager,et al.  A suite of extensions to a nonequilibrium surplus-production model , 1994 .

[30]  D. Butterworth,et al.  Does selectivity matter? A fisheries management perspective , 2014 .

[31]  M. Maunder,et al.  Selectivity's distortion of the production function and its influence on management advice from surplus production models , 2014 .

[32]  H. Winker,et al.  Comparison of two approaches to standardize catch-per-unit-effort for targeting behaviour in a multispecies hand-line fishery , 2013 .

[33]  Russell B. Millar,et al.  BUGS in Bayesian stock assessments , 1999 .

[34]  André E. Punt,et al.  Evaluating methods for setting catch limits in data-limited fisheries , 2014 .

[35]  A. AaltoEmilius,et al.  Separating recruitment and mortality time lags for a delay-difference production model , 2015 .

[36]  Murdoch K. McAllister,et al.  A Bayesian hierarchical analysis of stock-recruit data: quantifying structural and parameter uncertainties , 2004 .

[37]  Martyn Plummer,et al.  JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling , 2003 .

[38]  André E. Punt,et al.  Evaluating the performance of data-moderate and catch-only assessment methods for U.S. west coast groundfish , 2015 .

[39]  J. Pella,et al.  A generalized stock production model , 1969 .

[40]  William W. Fox,et al.  An Exponential Surplus-Yield Model for Optimizing Exploited Fish Populations , 1970 .

[41]  Merrill B. Rudd,et al.  Accounting for variable recruitment and fishing mortality in length-based stock assessments for data-limited fisheries , 2018, Canadian Journal of Fisheries and Aquatic Sciences.

[42]  M. Griffiths,et al.  Long-term trends in catch and effort of commercial linefish off South Africa's Cape Province: snapshots of the 20th century , 2000 .

[43]  G. Ishimura,et al.  Development of Bayesian production models for assessing the North Pacific swordfish population , 2010, Fisheries Science.

[44]  Mark N. Maunder,et al.  Evaluating fishery impacts: Application to bigeye tuna (Thunnus obesus) in the eastern Pacific Ocean , 2009 .

[45]  Ray Hilborn,et al.  Comparison of Fisheries Control Systems That Utilize Catch and Effort Data , 1979 .

[46]  A. Punt,et al.  Evidence for rebuilding in the panga stock on the Agulhas Bank, South Africa , 1998 .

[47]  Mark N. Maunder,et al.  Is it time to discard the Schaefer model from the stock assessment scientist's toolbox? , 2003 .

[48]  J. Thorson.  Guidance for decisions using the Vector Autoregressive Spatio-Temporal (VAST) package in stock, ecosystem, habitat and climate assessments , 2019, Fisheries Research.

[49]  André E. Punt,et al.  A Bayesian Approach to Stock Assessment and Harvest Decisions Using the Sampling/Importance Resampling Algorithm , 1994 .

[50]  André E. Punt,et al.  Model selection for selectivity in fisheries stock assessments , 2014 .

[51]  C. Walters,et al.  Quantitative fisheries stock assessment: Choice, dynamics and uncertainty , 2004, Reviews in Fish Biology and Fisheries.

[52]  M. Mangel,et al.  A perspective on steepness, reference points, and stock assessment , 2013 .

[53]  Kasper Kristensen,et al.  Implementing a generic method for bias correction in statistical models using random effects, with spatial and population dynamics examples , 2016 .

[54]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[55]  Richard D. Methot,et al.  Can steepness of the stock–recruitment relationship be estimated in fishery stock assessment models? , 2012 .

[56]  John Geweke,et al.  Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments , 1991 .

[57]  Keisuke Satoh,et al.  Get the biology right, or use size-composition data at your own risk , 2017 .

[58]  Philippe Grosjean,et al.  FLR: an open-source framework for the evaluation and development of management strategies , 2007 .

[59]  J. Thorson,et al.  A new role for effort dynamics in the theory of harvested populations and data-poor stock assessment , 2013 .

[60]  Murdoch K. McAllister,et al.  Using demographic methods to construct Bayesian priors for the intrinsic rate of increase in the Schaefer model and implications for stock rebuilding , 2001 .

[61]  Merrill B. Rudd,et al.  The case for estimating recruitment variation in data-moderate and data-poor age-structured models , 2019, Fisheries Research.

[62]  J. Thorson,et al.  Spawning biomass reference points for exploited marine fishes, incorporating taxonomic and body size information , 2012 .

[63]  André E. Punt,et al.  Assessing billfish stocks: A review of current methods and some future directions , 2015 .

[64]  A. Booth,et al.  Maximum likelihood and Bayesian approaches to stock assessment when data are questionable , 2006 .

[65]  Laurence T. Kell,et al.  Accounting for uncertainty due to data processing in virtual population analysis using Bayesian multiple imputation , 2018, Canadian Journal of Fisheries and Aquatic Sciences.

[66]  Martin Wæver Pedersen,et al.  A stochastic surplus production model in continuous time , 2017 .

[67]  H. Winker,et al.  Proof of concept for a novel procedure to standardize multispecies catch and effort data , 2014 .