Typically, fish that are landed by fishermen are not measured for length. However samples of the landed fish are taken through time by port samplers at the request of the assessment scientist. The size of the samples (the number of fish sampled) and the frequency of the samples change through time for several reasons. As a result some species are more heavily sampled than others. In addition some species are more heavily sampled during certain times of the year. In many cases the lengthed fish are assigned a market category, for example, Large, Small, Jumbo. Other times they are not and are simple considered Unclassified
For many stock assessment models and multispecies models the catch is often required as a composition of lengths. Due to the nature of the sampling program, length samples are not always available at the level the model requires. In this case there are several options. One of them being to “borrow” length samples from adjacent time periods. This assumes that fish caught from one time period to the next have similar size (length) distributions.
Expanding catch to a composition of lengths
Hydra, the multispecies length based model, is the only model that requires this type of data. Currently Hydra incorporates three distinct fleets, a demersal fleet, a pelagic fleet, and a fixed gear fleet. This data product categorizes catch by these three fleet types.
The annual catch (measured in metric tons) for each of the focal species is expanded to a composition of length. The result being that the sum of the catch over all lengths in a given year (or any other time interval) will equal the total catch for that year (or equivalent time interval). Fish are measured in 1 cm increments.
Commercial fisheries landings were pulled using R package
comlandr
The length data were pulled from the commercial fisheries database, specifically the
stockeff.mv_cf_len
tableThe R package
mscatch
is used to expand the annual catch data. Detailed methods are explained there. Species specific rules are used in the expansion to account for life history traits. See Data sources for more information-
Two data sets of expanded catch are produced:
- One where length samples are borrowed to fill in years/season/market
code with missing length samples,
realFisheryLencomp
- Another where length samples remain missing. In this cases the
landings can not be partitioned into landings at length. For length
based models this will represent missing data.
realFisheryLencompRaw
- One where length samples are borrowed to fill in years/season/market
code with missing length samples,
Species rules
For each species, a set of rules are applied (similar to the approach taken in a stock assessment), to determine how the landings data and discard data (not currently incorporated) are aggregated based on the availability of length samples. Biomass (metric tons) is then estimated at each LENGTH (1 cm bins) category for a specified time interval, market category, and gear type. The expanded catch is then aggregated to annual time steps and fleets.
To create the data set WITH borrowed length samples:
create_real_fishery_lencomp()
A sample of the output for Atlantic cod is shown below.
ModSim |
year |
Code |
Name |
fishery |
lenbin |
variable |
value |
units |
---|---|---|---|---|---|---|---|---|
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
36 |
biomass |
0.3934795 |
metric tons |
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
36 |
abundance |
881,464.4778893 |
numbers |
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
37 |
biomass |
1.1127373 |
metric tons |
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
37 |
abundance |
2,323,974.7067534 |
numbers |
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
38 |
biomass |
3.9541105 |
metric tons |
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
38 |
abundance |
7,403,875.3661190 |
numbers |
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
39 |
biomass |
7.2900553 |
metric tons |
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
39 |
abundance |
12,669,141.9442785 |
numbers |
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
40 |
biomass |
7.3784587 |
metric tons |
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
40 |
abundance |
11,895,843.5810585 |
numbers |
To create the data set WITHOUT borrowed length samples:
create_real_fishery_lencomp_raw()
A sample of the output for Atlantic cod is shown below.
ModSim |
year |
Code |
Name |
fishery |
lenbin |
variable |
value |
units |
---|---|---|---|---|---|---|---|---|
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
biomass |
9,363.9486 |
metric tons |
|
Actual |
1964 |
164712 |
Atlantic cod |
demersal |
abundance |
numbers |
||
Actual |
1964 |
164712 |
Atlantic cod |
fixedGear |
biomass |
1,003.1213 |
metric tons |
|
Actual |
1964 |
164712 |
Atlantic cod |
fixedGear |
abundance |
numbers |
||
Actual |
1965 |
164712 |
Atlantic cod |
demersal |
biomass |
8,674.7454 |
metric tons |
|
Actual |
1965 |
164712 |
Atlantic cod |
demersal |
abundance |
numbers |
||
Actual |
1965 |
164712 |
Atlantic cod |
fixedGear |
biomass |
930.8174 |
metric tons |
|
Actual |
1965 |
164712 |
Atlantic cod |
fixedGear |
abundance |
numbers |
||
Actual |
1966 |
164712 |
Atlantic cod |
demersal |
biomass |
8,517.8823 |
metric tons |
|
Actual |
1966 |
164712 |
Atlantic cod |
demersal |
abundance |
numbers |
Fleet/Gear Types
NEGEAR2 codes observed in the catch data are attributed to one of three fleets. The designation “all” implies all other gear codes. Note that some species are not caught by all three fleets.
species |
Fleet |
Negear2 |
---|---|---|
ATLANTIC COD |
demersal |
05,13,16,36 |
fixedGear |
all |
|
ATLANTIC HERRING |
pelagic |
all |
ATLANTIC MACKEREL |
demersal |
05 |
fixedGear |
all |
|
pelagic |
12,17,37 |
|
GOOSEFISH |
demersal |
05,13 |
fixedGear |
all |
|
HADDOCK |
demersal |
05,13 |
fixedGear |
all |
|
SILVER HAKE |
demersal |
all |
SPINY DOGFISH |
demersal |
all |
fixedGear |
01,02,10,50,52 |
|
WINTER FLOUNDER |
demersal |
all |
WINTER SKATE |
demersal |
all |
YELLOWTAIL FLOUNDER |
demersal |
all |
Market Codes
In a similar fashion to fleets, the market codes observed in the catch data are relabelled according to assessment scientist recommendations. The designation “all” implies all other market codes
species |
Relablled |
Market Codes |
---|---|---|
ATLANTIC COD |
LG |
LG,ST,WH |
MK |
MK |
|
SK |
SK,ST |
|
UN |
all |
|
ATLANTIC HERRING |
UN |
all |
ATLANTIC MACKEREL |
LG |
JB,LG,XG |
MD |
MD |
|
SM |
ES,SQ,SV |
|
UN |
all |
|
GOOSEFISH |
LG |
LG |
SQ |
PW,SQ |
|
UN |
all |
|
HADDOCK |
LG |
LG,XG |
SK |
SK |
|
SR |
SR |
|
UN |
all |
|
SILVER HAKE |
UN |
all |
SPINY DOGFISH |
UN |
all |
WINTER FLOUNDER |
LG |
LG,LM |
LS |
LS,XG |
|
SM |
MD,PW,SQ |
|
UN |
all |
|
WINTER SKATE |
UN |
all |
YELLOWTAIL FLOUNDER |
LG |
LG |
SM |
MD,PW,SQ |
|
UN |
all |
Additional preferences
Additional decision points were considered for each species in accordance with stock scientist recommendations. The following table outlines these decision points:
-
temporalAggregation
= The time scale that catch is expanded to lengths. Quarterly, by semester or annual time scales. -
howAggregate
= The method for dealing with missing length samples (at above time scales). Note: not all stock scientists borrow length samples from other sampled times. This will be addressed at a later date -
LengthWeightData
= The data source for obtaining length and weight data. This data is used to fit length-weight relationships at the level defined byLengthWeightRelationship
. -
LengthWeightRelationship
= The time scale used to fit length-weight relationships. Note: some stock scientists fit length-weight relationships for each semester in each year. Some fit over a specified time range of survey data. These modifications will be addressed at a later date. -
AgeData
= The data source for obtaining age data. Note that some species do not have age data available and subsequently age based stock assessments are not carried out. -
AgeLengthKey
= The time scale used to calculate the Age length key. Note: some stock scientists use year and quarter, some also incorporate gear type. These modifications will be addressed at a later date. -
maxAge
= The approximate observed maximum age of the species. This is used in the construction of the age length key.
## All species catch data are pulled from the statistical areas, 521,522,523,524,525,526,538,551,552,561,562,537
species |
temporalAggregation |
howAggregate |
LengthWeightData |
LengthWeightRelationships |
AgeData |
AgeLengthKey |
maxAge |
---|---|---|---|---|---|---|---|
GOOSEFISH |
semester |
borrow |
survey |
semester |
survey,commerical |
year,semester |
10 |
ATLANTIC COD |
quarter |
borrow |
survey |
semester |
survey,commerical |
year,semester |
15 |
WINTER FLOUNDER |
semester |
borrow |
survey |
single |
survey,commerical |
year,semester |
7 |
YELLOWTAIL FLOUNDER |
semester |
borrow |
survey |
semester |
survey,commerical |
year,semester |
8 |
HADDOCK |
semester |
borrow |
survey |
semester |
survey,commerical |
year,semester |
15 |
ATLANTIC HERRING |
quarter |
borrow |
survey |
quarter |
survey,commerical |
year,semester |
8 |
ATLANTIC MACKEREL |
semester |
borrow |
survey |
semester |
survey,commerical |
year,semester |
10 |
SPINY DOGFISH |
year |
borrow |
survey |
single |
|||
WINTER SKATE |
semester |
borrow |
survey |
single |
|||
SILVER HAKE |
semester |
borrow |
survey |
semester |
survey,commerical |
year,semester |
8 |
Length sample availability
The figures shown are generated by the function aggregate_landings from the R package mscatch. Each figure represents the number and frequency of fish lengths taken by port agents in each year, season, market category, gear type (represented as a fleet)