How do changing somatic growth rates affect estimation of management quantities in fisheries stock assessments, and can such growth changes be detected?

Christine Stawitz1
1SAFS
November 08, 2016 9:00 (PST): FSH 203

How do changing somatic growth rates affect estimation of management quantities in fisheries stock assessments, and can such growth changes be detected?

A substantial body of work suggests that somatic growth in marine fish is plastic, and thus growth rates may change substantially in response to both extrinsic (i.e. climate) and intrinsic (i.e. competition) factors. Additionally, such changes in somatic growth rate may substantially affect population dynamics of marine fish. However, these changes in growth rate are not incorporated into many fisheries stock assessment models and, consequently, management advice. Growth rate changes are difficult to incorporate in management models, since such models are already highly parameterized and often contain temporal variation in observational parameters which are confounded with somatic growth parameters, such as fisheries selectivity. In this analysis, we evaluate the detectability of temporal changes in somatic growth and how such changes may impact fisheries management quantities. We examine these questions using a simulation framework, employing an operating and estimation model closely modeled on stock assessment models for North Pacific groundfish and the Stock Synthesis 3 (SS3) program. Simulated variation in population processes, including growth, recruitment, and fishing, is derived directly from empirical estimates. We test, first, if regime-like growth shifts can be detected using realistic quantities of time series data. We find regime-like patterns in growth rates are not detectable in population models without many decades of high-quality composition data, even when process error is low. Secondly, we find such changes, while difficult to detect, can introduce substantial bias into management reference points, particularly current depletion levels. This suggests more work is needed to increase estimability of temporal growth changes and incorporate them into fisheries stock assessment models.

Resolving issues of hook saturation, hook competition, and fixed-site design in the Southern California hook-and-line survey

Peter Kuriyama
1SAFS
October 11, 2016 9:00 (PST): FSH 203

Resolving issues of hook saturation, hook competition, and fixed-site design in the Southern California hook-and-line survey

The Southern California hook-and-line survey has been conducted by the Northwest Fisheries Science Center since 2004 to monitor the untrawlable habitat of the Southern California Bight. Data from the survey have been used in stock assessments and supporting research for a number of shelf rockfish species, such as bocaccio (Sebastes pauicispinis) and vermilion rockfish (S. miniatus). However, an index of abundance estimated from hook-and-line data may be biased due to the fixed-site design of the survey and issues with hook saturation and hook competition. Here, I will present empirical results from the hook-and-line data and results of a simulation study exploring the biases associated with aspects of the survey. Bocaccio are the most sampled species in the survey, and sites with low catch rates of bocaccio have high catch rates of vermilion rockfish. Preliminary results from the simulations indicate that hook saturation causes estimates of abundance to be negatively biased at large population sizes. Additionally, hook competition leads to positively biased indices of abundance, and weighting catch rates by site leads to the least biased index of abundance. These results identify methods of incorporating hook-and-line data from untrawlable habitat into stock assessments and identify biases that are applicable to general hook-and-line survey methods.

Progressing scientific advice and management of mixed demersal fisheries in Europe

Paul Dolder1
1 MARES Doctoral Programme on Marine Ecosystem Health and Conservation

October 04, 2016 9:00 (PST): FSH 203

Progressing scientific advice and management of mixed demersal fisheries in
Europe

In Europe, scientific advice for managing multi-gear, multi-country mixed
fisheries is rapidly developing. Reforms to the EU’s Common Fisheries Policy
have given priority to the development and implementation of fishery-based
management plans and emphasize the importance of understanding how technical
interactions impede effective management of all stocks caught in mixed
fisheries. Modelling approaches which i) characterize these interactions and,
ii) can simulate a range of management options to contrast their bio-economic
impact, are gaining prominence. In this talk I’ll give an overview of the
developments happening within the EU and the ICES frameworks and highlight how
my PhD research aims to improve our ability to understand fleet behavioural
dynamics in mixed fisheries so that managers can better evaluate different
management options.

Optimizing otoliths sampling design in fishery-independent surveys for stock assessments

Dr. Gwladys Lambert1
1AFSC
May 31, 2016 9:00 (PST): FSH 105

Optimizing otoliths sampling design in fishery-independent surveys for stock assessments

The sampling design of composition data can have a major impact on our perception of the stock status and the estimated reference points for management. Composition data as sampled from bottom trawl surveys are famously clustered which creates major challenges when it comes to estimating length and age composition of the population of interest and in particular with respect to effective sample size/data weighting in stock assessment models. Here I present the approach I have been exploring in order to prepare a simulation framework to compare otoliths sampling strategies and their impact on stock assessment outputs. First I will talk about the questionnaires I have designed as part of the Otolith Sampling Size Working Group that aimed at informing the development of an operating model. Then I will discuss the implications of removing age data in data rich stock assessments based on a couple of examples and I will show some early simulation work using Stock Synthesis and developments from the R package ss3sim. Ss3sim developments include the addition of an option to test for the effect of ageing error on stock assessment outputs. The main focus will be on another option I have been experimenting with which aimed at simulating the sampling of clustered length data (that resemble real survey data) from random or length-stratified sampling in order to investigate the effect of different sampling strategies on the effective sample size estimated in stock assessments. I will also briefly talk about the Stock Synthesis model for GOA Pacific Ocean Perch I have been working on in order to use in the simulations described above.

Coupled biophysical/individual-based models: can they inform recruitment, stock assessment and/or fishery management? DisMELS lessons from the GOAIERP

Dr. William T. Stockhausen1
1AFSC
May 24, 2016 9:00 (PST): FSH 105

Coupled biophysical/individual-based models: can they inform recruitment, stock assessment and/or fishery management? DisMELS lessons from the GOAIERP

One of the key hypotheses underlying NPRB’s recent Gulf of Alaska (GOA) Integrated Ecosystem Research Project (IERP) was that year-class strength for the five commercially- and ecologically-important groundfish stocks at the focus of the project was determined by success or failure in running the early life, pelagic “gauntlet” from adult spawning areas to young-of-the-year nursery areas. As part of the multi-institution, multidisciplinary effort to address this hypothesis for pollock, Pacific cod, arrowtooth flounder, sablefish, and Pacific ocean perch, a group of modelers (including the author) constructed early life, coupled biophysical/individual-based models (CBP-IBMs) for each stock. For marine species with pelagic early life stages, CBP-IBMs provide one of the best approaches to synthesizing disparate information on individual life stages (e.g., eggs, larvae) and environmental variability into an integrated picture of early life processes. Among our goals for these models was to be able to improve current stock assessments and enhance fishery management by identifying the critical drivers of recruitment variability and hindcasting (or potentially forecasting) recent trends in recruitment for the five focal species. In regard to identifying critical drivers of recruitment variability, I think we can say that variability in physical transport, i.e. where young-of-the-year end up, is less important to recruitment variability than what happens along the way. However, in regard to predicting recruitment variability itself, we were certainly not terribly successful—except perhaps for pollock, the most data-informed of all the IBMs. As an extended introduction, I’ll briefly 1) introduce the concepts behind CBP-IBMs; 2) discuss DisMELS (the Dispersal Model for Early Life Stages), a framework for creating and running CBP-IBMs; and 3) review the five species-specific IBMs created for the GOAIERP. Then I’ll discuss results from the models and what they tell us about early life processes in the GOA. Finally, I’ll present a couple of practical suggestions for some “best practices” using these types of models.

Simulation Testing Hierarchical Multi-Stock Assessment Models for a BC Flatfish Complex

Samuel Johnson1
Simon Fraser University
Novermbet 22, 2016 9:00 (PST): FSH 203

Simulation Testing Hierarchical Multi-Stock Assessment Models for a BC Flatfish Complex

I am interested in finding the conditions that favour joint stock assessments for multiple flatfish species using hierarchical models. Hierarchical multi-level modeling approaches are being increasingly applied in stock asssessment to share information across stocks and species. These so-called “Robin Hood” approaches borrow information from data-rich to potentially improve assessments for data poor species; however, little research exists to demonstrate the boundary conditions where hierarchical approaches outperform single species models. In this talk, I present a hierarchical surplus production modeling approach for a multi-species flatfish complex in British Columbia, Canada. If applicable, the hierarchical modeling approach could help provide more up-to-date assessments for species last assessed in the 1990s!

The genetic basis of pathogen resistance across outbred salmon populations: a comparative analysis using Random Forest

Dr. Kerry Naish1, Dr. Marine Brieuc1, Charles Waters1, Dr. Maureen Purcell2
1SAFS
2Western Fisheries Research Center, USGS
May 17, 2016 9:00 (PST): FSH 105

The genetic basis of pathogen resistance across outbred salmon populations: a comparative analysis using Random Forest

Key questions in both aquaculture and conservation of wild populations are whether genomic regions associated with specific traits are conserved across lines or populations, and to what extent trait-linked markers might be used to predict population responses to selection. In aquaculture, the ability to improve unrelated lines using common trait-linked markers provides the industry with tools necessary to maintain several diverse breeding populations. In conservation, predicting differential population responses to anthropogenic change provides a way to evaluate different remedial actions. Infectious hematopoietic necrosis virus (IHNV) is a rhabdovirus endemic to the northeast Pacific that can cause morbidity and mortality in the commercial rainbow trout and Atlantic salmon aquaculture industries, as well as in supplementation hatcheries. We have investigated the genetic basis of IHNV resistance across different populations of steelhead salmon Oncorhynchus mykiss, in a broader effort to characterize the role of differential host response in the epidemiology of the pathogen. Using the rainbow trout 50K SNP chip, we performed genome wide association analyses across five outbred populations challenged with the pathogen. Traits within these populations may be encoded by many segregating loci of small effect, whose interactions may be influenced by dominance and epistasis. Therefore, we used Random Forest because this non-parametric approach has the advantage of identifying suites of interacting loci that explain phenotypic variation collectively, but may not have large effect sizes on their own. Here, we report several genomic regions associated with resistance and evaluate their conservation across populations. We also review the key steps, analytical approaches, and lessons we have learned in the application of Random Forest to genome wide association studies in outbred salmon populations, with a view to improving genomic resources available to the community as a whole.

Incorporating climate in to a multispecies size spectrum model of the Eastern Bering Sea

Dr. Jon Reum1
1AFSC
May 03, 2016 9:00 (PST): FSH 105

Incorporating climate in to a multispecies size spectrum model of the Eastern Bering Sea

I’ve been working with Kirstin Holsman, Kerim Aydin, and Anne Hollowed at the AFSC on developing a multispecies size-based food web model of the Eastern Bering Sea. I’ll give an overview of the model and how I’m trying to parameterize it. The model is part of a larger effort to evaluate food web and fisheries responses (and associated uncertainty) to climate change by using a variety of ecological models that vary in complexity and structural assumptions.

Evaluating alternative harvest strategies for Alaska groundfish fisheries based on estimates of the probability of overfishing

Dr. Carey R. McGilliard1
1AKFSC
April 19, 2016 9:00 (PST): FSH 105

Evaluating alternative harvest strategies for Alaska groundfish fisheries based on estimates of the probability of overfishing

Increasingly, management systems require that catch limits account for scientific uncertainty by specifying a buffer whereby the catch limit is lower than that determined based on an assumption of perfect information. A “P* approach” is a method whereby catch limits are specified such that the probability that a future fishing mortality (F) exceeds the F corresponding to maximum sustainable yield (FMSY) remains at or below a pre-specified value (P*). Some sources of scientific uncertainty are typically ignored when calculating catch limits, leading to underestimates of uncertainty about population biomass. A management strategy evaluation was developed to examine the performance of a P* approach that included accounting for typically ignored uncertainties. The P* approach was considered because it seems mandated by U.S. law; however, our results demonstrate that the approach has concerning properties. First, the amount of uncertainty that is taken into account is arbitrary and influences resulting catch limits. Secondly, small changes in the acceptable probability of overfishing can result in large changes to resulting catch limits. Third, the estimated and true probability of overfishing diverge as additional uncertainty is taken into account. Simulated cases will be shown to demonstrate each of these properties. Alternative harvest policies will be discussed.

The Alaska arrowtooth flounder stock assessment: how can we make it better?

Dr. Ingrid Spies1
1AKFSC
February 16, 2016 9:00 (PST): FSH 203

The Alaska arrowtooth flounder stock assessment: how can we make it better?

Coauthors: Steve Barbeaux and Miriam Doyle

A single stock assessment model written in AD Model Builder is used for the Gulf of Alaska and Bering Sea/Aleutian Island stocks of arrowtooth flounder. The age-structured model begins in 1961 in the GOA and 1976 in the BSAI. The modeled population includes individuals from age 1 to age 21, with age 21 defined as a “plus” group (all individuals age 21 and older). The age composition of the species shows fewer males relative to females as fish increase in age, which suggests higher natural mortality (M) for males. To account for this process, natural mortality is fixed at 0.2 for females and 0.35 for males in the model. The model is fit to time series of catch, survey indices of abundance (3 surveys in the BSAI: Bering Sea slope, Bering Sea shelf, and Aleutian Islands), and estimates of age and length composition from fishery and survey. Model parameters are estimated by maximizing the log likelihood of the data, and there are likelihood components for: survey biomass, lengths from the fishery and the survey, survey ages, catch, recruitment, and female and male fishery selectivity.

There are several aspects of arrowtooth distribution and life history that have not been considered in the stock assessment, and this presentation is designed to stimulate a discussion on how to improve the assessment. Models suggest that the population sizes have increased 8-fold in the BSAI and four-fold in the GOA over the past few decades. The model does not incorporate carrying capacity, but there are some signs that population growth is decreasing. In addition, the maximum age and sex ratio of arrowtooth flounder appears different by region, which may contradict the assumption of different natural mortality rates for males and females.