Introduction

This operating model is intended to be an exact replica of the age-based population dynamics of the MLE VPA stock assessment:

Swain, D.P., Savoie, L., Cox, S.P., and Aubry, E. 2015. Assessment of the southern Gulf of St. Lawrence Atlantic cod (Gadus morhua) stock of NAFO Div. 4T and 4Vn (November to April), March 2015. DFO Can. Sci. Advis. Sec. Res. Doc. 2015/080. xiv + 137 p.

Early indications are that stock depletion is entirely driven by non fishing related factors such as changes in survival, weight at age and recruitment strength. Apical fishing mortality rates (catchability, the y-axis) of 15 lead to 3.1% stock depletion (SSB in 2014 compared to unfished- x axis), apical fishing mortality rates of 0.0001 led to 3.25% stock depletion (see figure below).

This preliminary operating model assumes full stock-mixing (essentially a single area) and has no spatial structure.

Robustness tests

Future scenarios for maturity and growth

Alternative scenarios for future M (or MICE modelling with multiMSE)

Improvements (operating model specification given more time)

VPA estimates were only available from 1971 onwards despite depletion and SSB estimates starting from 1950. Where possible these inputs should match the breadth of VPA assessment estimates.

It would be desirable to obtain the bootstrapped estimates of VPA quantities for bracketting within operating model uncertainty. Right now this is achieved (probably not that badly) by uncertainty in current stock depletion (from the assessment) and a recalibration of fishing mortality rate to meet these Monte Carlo depletion estimates. Once these area available, this is trivial to implement in the build script.

The real relative abundance indices should be added to the observation model to keep these historical values constant under projections and generate future values with empirically derived statistical properties.

Obtain annual effort data to empirically derive interannual variability in catchability.

Operating Model

The OM rdata file can be downloaded from here

Download and import into R using myOM <- readRDS('OM.rdata')

Species Information

Species: Gadus morhua

Common Name: Atlantic cod

Management Agency: DFO

Region: Atlantic 4T4Vn

Latitude: 48.25

Longitude: -62

OM Parameters

OM Name: Name of the operating model: Cod_4T4Vn_DFO

nsim: The number of simulations: 192

proyears: The number of projected years: 50

interval: The assessment interval - how often would you like to update the management system? 4

pstar: The percentile of the sample of the management recommendation for each method: 0.5

maxF: Maximum instantaneous fishing mortality rate that may be simulated for any given age class: 2

reps: Number of samples of the management recommendation for each method. Note that when this is set to 1, the mean value of the data inputs is used. 1

Source: A reference to a website or article from which parameters were taken to define the operating model

Swain

Custom Parameters

The operating model makes extensive use of the cpars slot of the OM object to use the age and year specific VPA assessment estimates of:

  • length-at-age by year (survey)

  • weight-at-age by year (survey)

  • F-at-age by year (VPA)

  • maturity-at-age by year (VPA)

  • Natural mortality rate (ages 2-4 and 5+) (VPA)

Future age-selectivity was assumed to follow the average pattern estimated over the last three historical years (2012-2014).

In cases where length-at-age and weight-age-age were missing, 5 year moving averages were used to fill these (mostly older age classes in recent years). In cases where there were still NAs, a 10-year moving average was applied. Weights and lengths for fish greater than age 15 were assumed to be the same as fish age 15 (highest reported value in the VPA). Future weights and lengths at age were assumed to be the same as the final historical year (2014).

Maturity was assumed to be zero for age classes 0 and 1, and 100% for age classes > 12.

Depletion in OMx is phased in terms of current spawning stock biomass divided by the mean of the asymptotic SSB0 infered by maturity, survival, and weight at age in the first nm years where nm is the age at maturity. In this case 25550 / 1233137 which is around 2%. Arbitrarily we place a 5% CV around this specified depletion.

Stock-recruitment was assumed to be Beverton-Holt and steepness and R0 were estimated from a (quick and dirty) post-hoc MLE fit to the assessment SSB and (lagged by 2-years) estimates of recruitment. Spawning stock biomass per recruit was calculated from the maturity and mortality at age matrices. Mortality at age for years after 2010 were assumed to be the same as 2010.

To calculate recruitment (age 0 N), M estimates for age classes 2-4 were assumed for age classes 0 and 1. These were only available from 1971. It follows that corresponding M age 0 (1969) and M age 1 (1970) to calculate N age 0 from N age 2 in 1971 are not available. In this case the most recent available M estimates (1971) were used.

Modelling fish from age 0 is necessary in this operating model as it may be desirable to estimate the effect of herring on age-0 cod survival.

Stock Parameters

Mortality and age: maxage, R0, M, M2, Mexp, Msd

maxage: The maximum age of individuals that is simulated (there is no plus group ). Single value. Positive integer

Specified Value(s): 30

Following catch-at-age data this is a maximum of 30 to ensure calculations include all older individuals (there is no plus-group in the OMx operating model of DLMtool and MSEtool).

R0: The magnitude of unfished recruitment. Single value. Positive real number

Specified Value(s): 1945012.84

Defined in cpars. This is crudely calculated by MLE estimation in the build script. Note that the correct R0 is only require in MSE if absolute quantities of TAC advice are to be used. Ie if you want to know how well a particular TAC (e.g. 10000 tonnes) may work in the future, as opposed to a scale-less MP which generally provide advice scaled to the catches of any stock.

M: Natural mortality rate. Uniform distribution lower and upper bounds. Positive real number

Specified Value(s): 0.21, 0.21

Defined in cpars. Time varying age-specific M was specified based on the values presented in the stock assessment.

M2: (Optional) Natural mortality rate at age. Vector of length maxage . Positive real number

Slot not used.

Mexp: Exponent of the Lorenzen function assuming an inverse relationship between M and weight. Uniform distribution lower and upper bounds. Real numbers <= 0.

Specified Value(s): 0, 0

NA

Msd: Inter-annual variability in natural mortality rate expressed as a coefficient of variation. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0, 0

Defined in cpars. Time varying age-specific M was specified based on the values presented in the stock assessment.

Natural Mortality Parameters

Sampled Parameters

Histograms of 48 simulations of M, Mexp, and Msd parameters, with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

The average natural mortality rate by year for adult fish for 3 simulations. The vertical dashed line indicates the end of the historical period:

M-at-Age

Natural mortality-at-age for 3 simulations in the first historical year, the last historical year (i.e., current year), and the last projected year:

M-at-Length

Natural mortality-at-length for 3 simulations in the first historical year, the last historical year (i.e., current year), and the last projected year:

Recruitment: h, SRrel, Perr, AC

h: Steepness of the stock recruit relationship. Uniform distribution lower and upper bounds. Values from 1/5 to 1

Specified Value(s): 1, 1

We model an SSB-invariant recruitment relationship with steepness very close to 1.

SRrel: Type of stock-recruit relationship. Single value, switch (1) Beverton-Holt (2) Ricker. Integer

Specified Value(s): 1

A value of 1 represents the Beverton-Holt stock recruitment curve but this is essentially ignored and a mean recruitment with deviations model is uses instead.

Perr: Process error, the CV of lognormal recruitment deviations. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0.85, 0.85

Specified in cpars.

AC: Autocorrelation in recruitment deviations rec(t)=ACrec(t-1)+(1-AC)sigma(t). Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0.84, 0.84

Specified in cpars.

Recruitment Parameters

Sampled Parameters

Histograms of 48 simulations of steepness (h), recruitment process error (Perr) and auto-correlation (AC) for the Beverton-Holt stock-recruitment relationship, with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Non-stationarity in stock productivity: Period, Amplitude

Period: (Optional) Period for cyclical recruitment pattern in years. Uniform distribution lower and upper bounds. Non-negative real numbers

Slot not used.

Amplitude: (Optional) Amplitude in deviation from long-term average recruitment during recruitment cycle (eg a range from 0 to 1 means recruitment decreases or increases by up to 100% each cycle). Uniform distribution lower and upper bounds. 0 < Amplitude < 1

Slot not used.

Growth: Linf, K, t0, LenCV, Ksd, Linfsd

Linf: Maximum length. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 105, 105

Back-calculated from the length at age matrices specified in cpars.

K: von Bertalanffy growth parameter k. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.2, 0.2

Back-calculated from the length at age matrices specified in cpars.

t0: von Bertalanffy theoretical age at length zero. Uniform distribution lower and upper bounds. Non-positive real numbers

Specified Value(s): 0, 0

Back-calculated from the length at age matrices specified in cpars.

LenCV: Coefficient of variation of length-at-age (assumed constant for all age classes). Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.1, 0.1

Back-calculated from the length at age matrices specified in cpars.

Ksd: Inter-annual variability in growth parameter k expressed as coefficient of variation. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0, 0

No variability in growth is assumed here.

Linfsd: Inter-annual variability in maximum length expressed as a coefficient of variation. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0, 0

No variability in growth is assumed here.

Growth Parameters

Sampled Parameters

Histograms of 48 simulations of von Bertalanffy growth parameters Linf, K, and t0, and inter-annual variability in Linf and K (Linfsd and Ksd), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

The Linf and K parameters in each year for 3 simulations. The vertical dashed line indicates the end of the historical period:

Growth Curves

Sampled length-at-age curves for 3 simulations in the first historical year, the last historical year, and the last projection year.

Maturity: L50, L50_95

L50: Length at 50 percent maturity. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 39, 49

Back-calculated from the maturity at age matrices specified in cpars.

L50_95: Length increment from 50 percent to 95 percent maturity. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 2.6, 3.4

Back-calculated from the maturity at age matrices specified in cpars.

Maturity Parameters

Sampled Parameters

Histograms of 48 simulations of L50 (length at 50% maturity), L95 (length at 95% maturity), and corresponding derived age at maturity parameters (A50 and A95), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Maturity at Age and Length

Maturity-at-age and -length for 3 simulations in the first historical year, the last historical year (i.e., current year), and the last projected year:

Stock depletion and Discard Mortality: D, Fdisc

D: Current level of stock depletion SSB(current)/SSB(unfished). Uniform distribution lower and upper bounds. Fraction

Specified Value(s): 0.03, 0.03

In incredibly narrow range of depletions can be obtained by varying fishing from almost taking all vulnerable fish to zero fish. In this case SSB(2014) / SSB0 ranged from 0.0315 to 0.0322 (3%). Stock depletion in 2014 is entirely driven by non-fishing related factors such as changes in survival rate, maturity, weight and recruitment.

2014 SSB trend VPA

Fdisc: Fraction of discarded fish that die. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 1, 1

Discard mortality rate is assumed to be 100%.

Depletion and Discard Mortality

Sampled Parameters

Histograms of 48 simulations of depletion (spawning biomass in the last historical year over average unfished spawning biomass; D) and the fraction of discarded fish that are killed by fishing mortality (Fdisc), with vertical colored lines indicating 3 randomly drawn values.

Length-weight conversion parameters: a, b

a: Length-weight parameter alpha. Single value. Positive real number

Specified Value(s): 0

Back-calculated from the length and weight at age matrices specified in cpars.

b: Length-weight parameter beta. Single value. Positive real number

Specified Value(s): 3.08

Back-calculated from the length and weight at age matrices specified in cpars.

Spatial distribution and movement: Size_area_1, Frac_area_1, Prob_staying

Size_area_1: The size of area 1 relative to area 2. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.5, 0.5

No justification provided.

Frac_area_1: The fraction of the unfished biomass in stock 1. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.5, 0.5

A fully mixed stock is simulated.

Prob_staying: The probability of inviduals in area 1 remaining in area 1 over the course of one year. Uniform distribution lower and upper bounds. Positive fraction.

Specified Value(s): 0.5, 0.5

A fully mixed stock is simulated.

Spatial & Movement

Sampled Parameters

Histograms of 48 simulations of size of area 1 (Size_area_1), fraction of unfished biomass in area 1 (Frac_area_1), and the probability of staying in area 1 in a year (Frac_area_1), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Fleet Parameters

Historical years of fishing, spatial targeting: nyears, Spat_targ

nyears: The number of years for the historical spool-up simulation. Single value. Positive integer

Specified Value(s): 44

1971-2014, a total of 44 years.

Spat_targ: Distribution of fishing in relation to spatial biomass: fishing distribution is proportional to B^Spat_targ. Uniform distribution lower and upper bounds. Real numbers

Specified Value(s): 1, 1

Spatial targetting. We’re going to stick to the default level of 1 (effort distributed in proportion to density)

Trend in historical fishing effort (exploitation rate), interannual variability in fishing effort: EffYears, EffLower, EffUpper, Esd

EffYears: Years representing join-points (vertices) of time-varying effort. Vector. Non-negative real numbers

Specified by Find in custom parameters

EffLower: Lower bound on relative effort corresponding to EffYears. Vector. Non-negative real numbers

Specified by Find in custom parameters

EffUpper: Upper bound on relative effort corresponding to EffYears. Vector. Non-negative real numbers

Specified by Find in custom parameters

EffYears EffLower EffUpper
1971 0.829 0.829
1972 0.642 0.642
1973 0.593 0.593
1974 0.752 0.752
1975 1.050 1.050
1976 0.614 0.614
1977 0.434 0.434
1978 0.546 0.546
1979 0.754 0.754
1980 0.500 0.500
1981 0.603 0.603
1982 0.355 0.355
1983 0.500 0.500
1984 0.422 0.422
1985 0.396 0.396
1986 0.516 0.516
1987 0.379 0.379
1988 0.528 0.528
1989 0.630 0.630
1990 0.860 0.860
1991 0.725 0.725
1992 0.822 0.822
1993 0.227 0.227
1994 0.046 0.046
1995 0.026 0.026
1996 0.025 0.025
1997 0.039 0.039
1998 0.051 0.051
1999 0.133 0.133
2000 0.169 0.169
2001 0.182 0.182
2002 0.192 0.192
2003 0.020 0.020
2004 0.096 0.096
2005 0.098 0.098
2006 0.153 0.153
2007 0.082 0.082
2008 0.084 0.084
2009 0.015 0.015
2010 0.008 0.008
2011 0.011 0.011
2012 0.019 0.019
2013 0.020 0.020
2014 0.021 0.021

Esd: Additional inter-annual variability in fishing mortality rate. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0, 0

Specified by Find in custom parameters

Historical Effort

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in historical fishing mortality (Esd), with vertical colored lines indicating 3 randomly drawn values used in the time-series plot:

Time-Series

Time-series plot showing 3 trends in historical fishing mortality (OM@EffUpper and OM@EffLower or OM@cpars$Find):

Annual increase in catchability, interannual variability in catchability: qinc, qcv

qinc: Average percentage change in fishing efficiency (applicable only to forward projection and input controls). Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): -0.1, 0.1

Trend in fishing efficiency (projections). In the future, is a unit of effort going to lead to higher fishing mortality (positive qinc) or lower fishing mortality (negative qinc). Since there isn’t a compelling reason to expect fishing to become more or less efficient, we set the % annual increase to be very close to zero, -0.1 to 0.1.

qcv: Inter-annual variability in fishing efficiency (applicable only to forward projection and input controls). Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0.1, 0.1

This should be calculated from assessment fishing mortality rate divided by observed fishing effort ie sd(F(y)/E(y)). In the absence of effort data (E) I’m guessing this is moderate at 10% interannual variability.

Future Catchability

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in fishing efficiency (qcv) and average annual change in fishing efficiency (qinc), with vertical colored lines indicating 3 randomly drawn values used in the time-series plot:

Time-Series

Time-series plot showing 3 trends in future fishing efficiency (catchability):

Fishery gear length selectivity: L5, LFS, Vmaxlen, isRel

L5: Shortest length corresponding to 5 percent vulnerability. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 1, 1

From the VPA

LFS: Shortest length that is fully vulnerable to fishing. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 1, 1

From the VPA

Vmaxlen: The vulnerability of fish at . Uniform distribution lower and upper bounds. Fraction

Specified Value(s): 1, 1

From the VPA

isRel: Selectivity parameters in units of size-of-maturity (or absolute eg cm). Single value. Boolean.

Specified Value(s): FALSE

Selectivity is in terms of absolute length, not relative to length at maturity.

Fishery length retention: LR5, LFR, Rmaxlen, DR

LR5: Shortest length corresponding ot 5 percent retention. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0, 0

Retention follows selectivity.

LFR: Shortest length that is fully retained. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0, 0

Retention follows selectivity.

Rmaxlen: The retention of fish at . Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 1, 1

Retention follows selectivity.

DR: Discard rate - the fraction of caught fish that are discarded. Uniform distribution lower and upper bounds. Fraction

Specified Value(s): 0, 0

No general discarding rate.

Time-varying selectivity: SelYears, AbsSelYears, L5Lower, L5Upper, LFSLower, LFSUpper, VmaxLower, VmaxUpper

SelYears: (Optional) Years representing join-points (vertices) at which historical selectivity pattern changes. Vector. Positive real numbers

Slot not used.

AbsSelYears: (Optional) Calendar years corresponding with SelYears (eg 1951, rather than 1), used for plotting only. Vector (of same length as SelYears). Positive real numbers

Slot not used.

L5Lower: (Optional) Lower bound of L5 (use ChooseSelect function to set these). Vector. Non-negative real numbers

Slot not used.

L5Upper: (Optional) Upper bound of L5 (use ChooseSelect function to set these). Vector. Non-negative real numbers

Slot not used.

LFSLower: (Optional) Lower bound of LFS (use ChooseSelect function to set these). Vector. Non-negative real numbers

Slot not used.

LFSUpper: (Optional) Upper bound of LFS (use ChooseSelect function to set these). Vector. Non-negative real numbers

Slot not used.

VmaxLower: (Optional) Lower bound of Vmaxlen (use ChooseSelect function to set these). Vector. Fraction

Slot not used.

VmaxUpper: (Optional) Upper bound of Vmaxlen (use ChooseSelect function to set these). Vector. Fraction

Slot not used.

Current Year: CurrentYr

CurrentYr: The current calendar year (final year) of the historical simulations (eg 2011). Single value. Positive integer.

Specified Value(s): 2014

The last year of F estimates from the update assessment

Existing Spatial Closures: MPA

MPA: (Optional) Matrix specifying spatial closures for historical years.

Slot not used.

Obs Parameters

The observation model parameter are taken from the Generic_Obs model subject to a few addtional changes which are documented here.

Catch statistics: Cobs, Cbiascv, CAA_nsamp, CAA_ESS, CAL_nsamp, CAL_ESS

Cobs: Log-normal catch observation error expressed as a coefficient of variation. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0.05, 0.1

Catches are observed more precisely than the Precise_Unbiased object with a CV of between 5 and 10 per cent.

Cbiascv: Log-normal coefficient of variation controlling the sampling of bias in catch observations for each simulation. Uniform distribution lower and upper bounds. Non-negative real numbers

Specified Value(s): 0.02

Mean bias (under / over reporting) in catches is assumed to be small with a CV of 0.025 - 95% of simulations are between 95% and 105% of true simulated catches.

CAA_nsamp: Number of catch-at-age observation per time step. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 100, 200

As Precise_Unbiased

CAA_ESS: Effective sample size (independent age draws) of the multinomial catch-at-age observation error model. Uniform distribution lower and upper bounds. Positive integers

Specified Value(s): 25, 50

As Precise_Unbiased

CAL_nsamp: Number of catch-at-length observation per time step. Uniform distribution lower and upper bounds. Positive integers

Specified Value(s): 100, 200

As Precise_Unbiased

CAL_ESS: Effective sample size (independent length draws) of the multinomial catch-at-length observation error model. Uniform distribution lower and upper bounds. Positive integers

Specified Value(s): 25, 50

As Precise_Unbiased

Index imprecision, bias and hyperstability: Iobs, Ibiascv, Btobs, Btbiascv, beta

Iobs: Observation error in the relative abundance indices expressed as a coefficient of variation. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.1, 0.15

Relative abundance indices are assumed to be observed somewhat more precisely than the Precise_Unbiased model with a CV of between 10 and 15 per cent.

Ibiascv: Not Used. Log-normal coefficient of variation controlling error in observations of relative abundance index. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.2

This parameter is not used in this version of DLMtool.

Btobs: Log-normal coefficient of variation controlling error in observations of current stock biomass among years. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.2, 0.5

As Precise_Unbiased

Btbiascv: Uniform-log bounds for sampling persistent bias in current stock biomass. Uniform-log distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.5, 2

The bias in the absolute abundance index is assumed to be reasonably high and could be 1/2 to 2 times the true value.

beta: A parameter controlling hyperstability/hyperdepletion where values below 1 lead to hyperstability (an index that decreases slower than true abundance) and values above 1 lead to hyperdepletion (an index that decreases more rapidly than true abundance). Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.66, 1.5

Since the survey is carried out according to a systematic design we assume that it varies roughly proportionally to real abundance and specify a beta parameter between 2/3 and 3/2.

Bias in maturity, natural mortality rate and growth parameters: LenMbiascv, Mbiascv, Kbiascv,t0biascv, Linfbiascv

LenMbiascv: Log-normal coefficient of variation for sampling persistent bias in length at 50 percent maturity. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.1

As Precise_Unbiased

Mbiascv: Log-normal coefficient of variation for sampling persistent bias in observed natural mortality rate. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.2

As Precise_Unbiased

Kbiascv: Log-normal coefficient of variation for sampling persistent bias in observed growth parameter K. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.05

As Precise_Unbiased

t0biascv: Log-normal coefficient of variation for sampling persistent bias in observed t0. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0

Unbiased

Linfbiascv: Log-normal coefficient of variation for sampling persistent bias in observed maximum length. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.02

Twice as accurate as Generic_Obs with L-infinity values within plus or minus 5% of true value.

Bias in length at first capture, length at full selection: LFCbiascv, LFSbiascv

LFCbiascv: Log-normal coefficient of variation for sampling persistent bias in observed length at first capture. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.05

As Precise_Unbiased

LFSbiascv: Log-normal coefficient of variation for sampling persistent bias in length-at-full selection. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.05

As Generic_Obs

Bias in fishery reference points, unfished biomass, FMSY, FMSY/M ratio, biomass at MSY relative to unfished: FMSYbiascv, FMSY_Mbiascv, BMSY_B0biascv

FMSYbiascv: Not used. Log-normal coefficient of variation for sampling persistent bias in FMSY. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.4

Highly dependent on M for which future levels are highly uncertain.

FMSY_Mbiascv: Log-normal coefficient of variation for sampling persistent bias in FMSY/M. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.15

More stable and better known than FMSY.

BMSY_B0biascv: Log-normal coefficient of variation for sampling persistent bias in BMSY relative to unfished. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.2

Again, highly dependent on M which is very uncertain in the future.

Management targets in terms of the index (i.e., model free), the total annual catches and absolute biomass levels: Irefbiascv, Crefbiascv, Brefbiascv

Irefbiascv: Log-normal coefficient of variation for sampling persistent bias in relative abundance index at BMSY. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.5

Uncertain due to uncertainty in BMSY/B0

Crefbiascv: Log-normal coefficient of variation for sampling persistent bias in MSY. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.5

Less likely to be biased than FMSY as high estimated Ms lead to low estimated B0 (MSY is the proportional to the product more or less)

Brefbiascv: Log-normal coefficient of variation for sampling persistent bias in BMSY. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.5

Arbitrarily we make this twice as potentially biased as Iref and Cref.

Depletion bias and imprecision: Dbiascv, Dobs

Dbiascv: Log-normal coefficient of variation for sampling persistent bias in stock depletion. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.5

Could be very seriously biased given changing M scenarios

Dobs: Log-normal coefficient of variation controlling error in observations of stock depletion among years. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.05, 0.1

In a data-limited situation it is unlikely that radically new data would become available regarding depletion meaning that while estimates may be biased, they are likely to be relatively precise. We assign a level of imprecision consistent with observations of catch rate data among years at between 0.05- 0.1.

Recruitment compensation and trend: hbiascv, Recbiascv

hbiascv: Log-normal coefficient of variation for sampling persistent bias in steepness. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.1

Steepness appears consistently low according to S-R estimates from the 2008 VPA

Recbiascv: Log-normal coefficient of variation for sampling persistent bias in recent recruitment strength. Uniform distribution lower and upper bounds. Positive real numbers

Specified Value(s): 0.05, 0.1

As Precise_Unbiased

Obs Plots

Observation Parameters

Catch Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in catch observations (Csd) and persistent bias in observed catch (Cbias), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Depletion Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in depletion observations (Dobs) and persistent bias in observed depletion (Dbias), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Abundance Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in abundance observations (Btobs) and persistent bias in observed abundance (Btbias), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Index Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in index observations (Iobs) and hyper-stability/depletion in observed index (beta), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Time-series plot of 3 samples of index observation error:

Plot showing an example true abundance index (blue) with 3 samples of index observation error and the hyper-stability/depletion parameter (beta):

Recruitment Observations

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in index observations (Recsd) , with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Composition Observations

Sampled Parameters

Histograms of 48 simulations of catch-at-age effective sample size (CAA_ESS) and sample size (CAA_nsamp) and catch-at-length effective (CAL_ESS) and actual sample size (CAL_nsamp) with vertical colored lines indicating 3 randomly drawn values:

Parameter Observations

Sampled Parameters

Histograms of 48 simulations of bias in observed natural mortality (Mbias), von Bertalanffy growth function parameters (Linfbias, Kbias, and t0bias), length-at-maturity (lenMbias), and bias in observed length at first capture (LFCbias) and first length at full capture (LFSbias) with vertical colored lines indicating 3 randomly drawn values:

Reference Point Observations

Sampled Parameters

Histograms of 48 simulations of bias in observed FMSY/M (FMSY_Mbias), BMSY/B0 (BMSY_B0bias), reference index (Irefbias), reference abundance (Brefbias) and reference catch (Crefbias), with vertical colored lines indicating 3 randomly drawn values:

Imp Parameters

Output Control Implementation Error: TACFrac, TACSD

TACFrac: Mean fraction of TAC taken. Uniform distribution lower and upper bounds. Positive real number.

Specified Value(s): 1, 1.05

Here we assume that the actual catches can be up to 5% higher than the recommended TAC.

TACSD: Log-normal coefficient of variation in the fraction of Total Allowable Catch (TAC) taken. Uniform distribution lower and upper bounds. Non-negative real numbers.

Specified Value(s): 0.02, 0.05

We assume that the bias in the actual catch is relatively consistent between years and set the range for this parameter to a low value.

Effort Control Implementation Error: TAEFrac, TAESD

TAEFrac: Mean fraction of TAE taken. Uniform distribution lower and upper bounds. Positive real number.

Specified Value(s): 1, 1.05

We have little information to inform this parameter, and set the implementation error in effort equal to the TAC implementation error.

TAESD: Log-normal coefficient of variation in the fraction of Total Allowable Effort (TAE) taken. Uniform distribution lower and upper bounds. Non-negative real numbers.

Specified Value(s): 0.02, 0.05

We assume that the bias in the effort is relatively consistent between years and set the range for this parameter to a low value.

Size Limit Control Implementation Error: SizeLimFrac, SizeLimSD

SizeLimFrac: The real minimum size that is retained expressed as a fraction of the size. Uniform distribution lower and upper bounds. Positive real number.

Specified Value(s): 1, 1.05

We assume that, on average, a size limit would be well-implemented.

SizeLimSD: Log-normal coefficient of variation controlling mismatch between a minimum size limit and the real minimum size retained. Uniform distribution lower and upper bounds. Non-negative real numbers.

Specified Value(s): 0.02, 0.05

We assume that the implementation of the size limit is relatively consistent between years.

Imp Plots

Implementation Parameters

TAC Implementation

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in TAC implementation error (TACSD) and persistent bias in TAC implementation (TACFrac), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

TAE Implementation

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in TAE implementation error (TAESD) and persistent bias in TAC implementation (TAEFrac), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Size Limit Implementation

Sampled Parameters

Histograms of 48 simulations of inter-annual variability in size limit implementation error (SizeLimSD) and persistent bias in size limit implementation (SizeLimFrac), with vertical colored lines indicating 3 randomly drawn values used in other plots:

Time-Series

Historical Simulation Plots

Historical Time-Series

Spawning Biomass

Depletion

Absolute

Vulnerable Biomass

Depletion

Absolute

Total Biomass

Depletion

Absolute

Recruitment

Relative

Absolute

Catch

Relative

Absolute

Historical Fishing Mortality

Historical Time-Series

References