EconMult

EconMult is a general fleet model to be used in fisheries modelling. EconMult has been developed since 1991 as a part of the Multispecies management programme by the Norwegian Research Council at the Norwegian College of Fishery Science (University of Tromsø, Norway).

Model resolution and key variables

EconMult is a discrete time simulation model where the fleet activity is controlled by two variables: Number of vessels (v) (within each fleet segment) and Number of fishing days d) (within each time interval). The biomass units (x) are exogenous variables. The model resolution is determined by four structural variables: Number of fleet groups j), Number of targeted species n), Number of biomass units i) (which may be more than one within each targeted species) and Period length (time step in the simulation). The number of vessels and fishing days therefore are presented in fleet (columns)-targeted species (rows) matrices, while the biomass units is presented in a column vector (X):

V={\begin{pmatrix}v_{1,1}&\cdots &v_{1,n}\\\vdots &\ddots &\vdots \\v_{j,1}&\cdots &v_{j,n}\end{pmatrix}}

$D={\begin{pmatrix}d_{1,1}&\cdots &d_{1,n}\\\vdots &\ddots &\vdots \\d_{j,1}&\cdots &d_{j,n}\end{pmatrix}}$

$X={\begin{pmatrix}x_{1}\\\vdots \\x_{i}\end{pmatrix}}$

$i\geq n>0$
$j>0\,\!$

Catch production

A fishery is in EconMult defined as a unique fleet/targeted species combination. The total catch within each fishery may include all biomass units defined in the model. Each biomass unit vessel catch is computed by Cobb–Douglas production function, applying two input variables: number of fishing days d) and biomass unit x). In the matrix below each column represents a fleet group and each row a targeted species so each element in the matrix is a fishery and gives the vessel catch of the biomass units represented. The biomass units represent all the targeted species. Each catch is represented a Cobb–Douglas production equation as shown in the vessel yield matrix (Y):

$Y_{V}={\begin{pmatrix}{\begin{pmatrix}q_{1,1,1}d_{1,1}v_{1,1}^{\alpha _{1,1,1}-1}x_{1}^{\beta _{1,1,1}}\\\vdots \\q_{1,1,i}d_{1,1}v_{1,1}^{\alpha _{1,1,i}-1}x_{i}^{\beta _{1,1,i}}\end{pmatrix}}&\cdots &{\begin{pmatrix}q_{1,n,1}d_{1,n}v_{1,n}^{\alpha _{1,n,1}-1}x_{1}^{\beta _{1,n,1}}\\\vdots \\q_{1,n,i}d_{1,n}v_{1,n}^{\alpha _{1,n,i}-1}x_{i}^{\beta _{1,n,i}}\end{pmatrix}}\\\vdots &\ddots &\vdots \\{\begin{pmatrix}q_{j,1,1}d_{j,1}v_{j,1}^{\alpha _{j,1,1}-1}x_{1}^{\beta _{j,1,1}}\\\vdots \\q_{j,1,i}d_{j,1}v_{j,1}^{\alpha _{j,1,i}-1}x_{i}^{\beta _{j,1,i}}\end{pmatrix}}&\cdots &{\begin{pmatrix}q_{j,n,1}d_{j,n}v_{j,n}^{\alpha _{j,n,1}-1}x_{1}^{\beta _{j,n,1}}\\\vdots \\q_{j,n,i}d_{j,n}v_{j,n}^{\alpha _{j,n,i}-1}x_{i}^{\beta _{j,n,i}}\end{pmatrix}}\end{pmatrix}}$

The corresponding fleet catches are

$Y={\begin{pmatrix}{\begin{pmatrix}q_{1,1,1}d_{1,1}v_{1,1}^{\alpha _{1,1,1}}x_{1}^{\beta _{1,1,1}}\\\vdots \\q_{1,1,i}d_{1,1}v_{1,1}^{\alpha _{1,1,i}}x_{i}^{\beta _{1,1,i}}\end{pmatrix}}&\cdots &{\begin{pmatrix}q_{1,n,1}d_{1,n}v_{1,n}^{\alpha _{1,n,1}}x_{1}^{\beta _{1,n,1}}\\\vdots \\q_{1,n,i}d_{1,n}v_{1,n}^{\alpha _{1,n,i}}x_{i}^{\beta _{1,n,i}}\end{pmatrix}}\\\vdots &\ddots &\vdots \\{\begin{pmatrix}q_{j,1,1}d_{j,1}v_{j,1}^{\alpha _{j,1,1}}x_{1}^{\beta _{j,1,1}}\\\vdots \\q_{j,1,i}d_{j,1}v_{j,1}^{\alpha _{j,1,i}}x_{i}^{\beta _{j,1,i}}\end{pmatrix}}&\cdots &{\begin{pmatrix}q_{j,n,1}d_{j,n}v_{j,n}^{\alpha _{j,n,1}}x_{1}^{\beta _{j,n,1}}\\\vdots \\q_{j,n,i}d_{j,n}v_{j,n}^{\alpha _{j,n,i}}x_{i}^{\beta _{j,n,i}}\end{pmatrix}}\end{pmatrix}}$

α β and q are parameters, the first two known as output elasticities of effort and biomass respectively, q is often referred to as the catchability coefficient. All the three parameters have the same dimension as the matrix above, e.g. the catchability coefficient:

$Q={\begin{pmatrix}{\begin{pmatrix}q_{1,1,1}\\\vdots \\q_{1,1,i}\end{pmatrix}}&\cdots &{\begin{pmatrix}q_{1,n,1}\\\vdots \\q_{1,n,i}\end{pmatrix}}\\\vdots &\ddots &\vdots \\{\begin{pmatrix}q_{j,1,1}\\\vdots \\q_{j,1,i}\end{pmatrix}}&\cdots &{\begin{pmatrix}q_{j,n,1}\\\vdots \\q_{j,n,i}\end{pmatrix}}\end{pmatrix}}$

gollark: I was lazy, so it reads (can read) slightly outside of the LTS buffer.

gollark: I think `lts_length` is technically capable of causing segfaults in some situations too.

gollark: ```c#define REP(x) x x x x#define T1 "a"#define T2 REP(T1)#define T3 REP(T2)#define T4 REP(T3)#define T5 REP(T4)#define T6 REP(T5)#define T7 REP(T6)```This was part of my failed plan to implement hyperoperations in the preprocessor.

gollark: I said length TERMINATED.

gollark: ```cchar * lts_to_nullterminated(const char * apioform) { uint32_t len = lts_length(apioform); char * bee = malloc(len + 1); memcpy(bee, apioform, len); bee[len] = 0; // for purposes only return bee;}```Gaze upon my wondrous code.

Downloads

Mathematica packages related to EconMult:

Download EconMult.m and PopulationGrowth.m

External links

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

EconMult

Model resolution and key variables

Catch production

See also

Downloads

External links