Skip Navigation


ICES Journal of Marine Science: Journal du Conseil Advance Access originally published online on November 6, 2006
ICES Journal of Marine Science: Journal du Conseil 2007 64(1):178-191; doi:10.1093/icesjms/fsl015
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
64/1/178    most recent
fsl015v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Ellis, N.
Right arrow Articles by Wang, Y.-G.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Ellis, N.
Right arrow Articles by Wang, Y.-G.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© 2006 International Council for the Exploration of the Sea. Published by Oxford Journals. All rights reserved. For Permissions, please email: journals.permissions@oxfordjournals.org

Effects of fish density distribution and effort distribution on catchability

Nick Ellis1 and You-Gan Wang2

1 CSIRO Division of Marine and Atomospheric Research, PO Box 120, Cleveland, Queensland 4163, Australia
2 CSIRO Mathematical and Information Sciences, Long Pocket Laboratories, 120 Meiers Road, Indooroopilly, Queensland 4068, Australia

Correspondence to Y-G. Wang: tel: +61 7 3214 2700; fax: +61 7 3214 2855; e-mail: you-gan.wang{at}csiro.au

Ellis, N., and Wang, Y-G. 2007. Effects of fish density distribution and effort distribution on catchability – ICES Journal of Marine Science, 64, 178–191.

The effects of fish density distribution and effort distribution on the overall catchability coefficient are examined. Emphasis is also on how aggregation and effort distribution interact to affect overall catch rate [catch per unit effort (cpue)]. In particular, it is proposed to evaluate three indices, the catchability index, the knowledge parameter, and the aggregation index, to describe the effectiveness of targeting and the effects on overall catchability in the stock area. Analytical expressions are provided so that these indices can easily be calculated. The average of the cpue calculated from small units where fishing is random is a better index for measuring the stock abundance. The overall cpue, the ratio of lumped catch and effort, together with the average cpue, can be used to assess the effectiveness of targeting. The proposed methods are applied to the commercial catch and effort data from the Australian northern prawn fishery. The indices are obtained assuming a power law for the effort distribution as an approximation of targeting during the fishing operation. Targeting increased catchability in some areas by 10%, which may have important implications on management advice.

Keywords: aggregation, catchability, catch and effort data, density distribution, effort distribution, stock assessment

Received 31 October 2005; accepted 25 August 2006; advance access publication 6 November 2006.


    Introduction
 Top
 Introduction
 Relative fish density...
 Effect of aggregation
 Effect of effort distribution
 Analysis of catch and...
 Analysis of commercial catch...
 Discussion
 Appendix
 References
 
Catchability is defined as the proportion of available fish in the population that would be caught by a unit of effort. If the fishing effort is randomly distributed in the sense that it is not related to fish abundance, catchability is also the probability of an individual fish being caught by a unit of effort. However, when the distribution of fishing effort depends on fish abundance (attributable to fish aggregation and fisher targeting), the overall catchability for the stock will be affected, depending on the degree of aggregation and targeting (Francis et al., 2003). One way to explain this is to attribute it to changes in overall catchability for the whole area of the stock. For convenience, we refer to this phenomenon as density-dependent catchability.

The simplest approach to catch and effort data is based on the well-known catch equation of Baranov (1918), which is described by Gulland (1983, p. 105) and Xiao (2005, 2006). The catch equation presumes that fishing effort is randomly distributed over a fishing ground in the sense that it is unrelated to fish abundance. This is clearly not the case for most fisheries because of (i) the aggregation behaviour of fish, and (ii) the targeting behaviour of fishers. In such cases, the traditional catch equation is inappropriate for describing the relationship between catch and effort (Crecco and Overholtz, 1990; Richards and Schnute, 1992; Paloheimo and Chen, 1993), because the overall catchability is affected by the extent of correlation between effort and abundance, resulting in the so-called density-dependent catchability. Paloheimo and Dickie (1964) proposed density-dependent catchabilities for Georges Bank haddock, and Crecco and Overholtz (1990) tested and supported their theory. Shardlow (1993) also found that the catchability increases with abundance in a salmon fishery. Furthermore, Paloheimo and Dickie (1964) argued that density-dependent catchability exists for most demersal and pelagic fisheries. Indeed, Swain and Wade (2003) found that abundance was positively correlated with effort for the snow crab fishery, and Hilborn and Walters (1992, pp. 175–192) discussed the challenges of obtaining abundance indices from spatial catch and effort data in the presence of hyperstability and hyperdepletion. The results from those authors suggest that density-dependent catchability should be incorporated into population models and into the estimation methods, to obtain less biased estimates of effective effort, which are essential for stock assessment models.

Walters (2003) has warned that "simple, nonspatial ratio (cpue) estimates should not be used". An alternative approach is to apply generalized linear models (GLMs) to take account of the many factors that affect catchability (Bishop et al., 2000; Campbell, 2004). The objective is to obtain a closer relationship between average catch per unit effort (cpue) that reflects the stock status for stock management. It is therefore essential to define an average cpue which, by taking account of how the allocation of effort responds to aggregation, is proportional to stock abundance. The constant of proportionality should then be independent of time and effort.

In Australian prawn fisheries, fishers use their knowledge of aggregation behaviour, together with their experience and modern searching technology, to enhance catch rates (). The extent of fisher targeting in general depends on the aggregation level and the cost of finding high densities of fish. Therefore, the overall catchability will depend on the distribution of fishing effort in relation to the distribution of the target species. In fact, Salthaug and Aanes (2003) showed that catchability of two migratory stocks was strongly related to the concentration of the fleet. Here we examine the effects of both fish distribution (aggregation) and effort distribution on catchability in a way similar to that of Swain and Sinclair (1994) and Swain and Morin (1996) and describe effort distribution using a parameter that measures the knowledge of fish density distribution. A change in catchability is a joint effect of the aggregation of fish and the knowledge of fishers. We also demonstrate our concepts by analysing the commercial catch and effort data for Australia's northern prawn fishery.


    Relative fish density distribution
 Top
 Introduction
 Relative fish density...
 Effect of aggregation
 Effect of effort distribution
 Analysis of catch and...
 Analysis of commercial catch...
 Discussion
 Appendix
 References
 
Following Swain and Sinclair (1994), assume that fish density at geographic position x is given by d(x)=Nf(x), where N is the abundance and f(x)>0 is given by a probability density function. We call f(x) the relative fish density distribution. (For a table of all symbols used herein, the reader is referred to Table 1.)


View this table:
[in this window]
[in a new window]

 
Table 1. Symbols used throughout the paper.

 
Let e(x) be the density of effort at position x, such that the total effort over area A is E={int}Ae(x)dx. Then the catch density c(x) is given by


Formula 015M1

(1)
Quantity p is a constant independent of x called the instantaneous catchability. It is the proportion of fish at x that would be caught if subjected to a unit of effort density. Its units depend on the units of effort; for trawl effort, where the units are swept area (and effort density is dimensionless), p is the probability that a fish subjected to a single tow will be caught by the trawl gear. In contrast, for a hook fishery, where the units of effort are, for example, hook-minutes, p is the proportion of fish over a local area that would be caught by one hook-minute per square metre.

In general, let a be a unit of effort. The (spatial) catchability q(x) at position x is defined as


Formula 015M2

(2)
It is the expected proportion of the population caught by a unit of effort at x. It depends on position through the density of the fish.

The average catchability q is the average of q(x) weighted by effort:


Formula 015M3

(3)
It is the proportion of the population caught per unit effort. The average catchability combines the distributions of both fish and effort. If fishing is assumed to be uniform over the entire area of fish stock A, we find that q=pa/A, which differs from the instantaneous catchability by a constant factor. Quantity q is the analogue of p at the scale of the fishery. For trawl effort, if we define the unit of effort as the singly swept-out area of the entire fishery, so that a=A, then, numerically, q=p. For targeted trawling, one would expect q/p>1.

In general q/p, which we will refer to as the catchability index {kappa}, may be used to quantify the effect of fish aggregation and fisher knowledge (or targeting) on the stock. It is usually much easier to estimate {kappa} than to estimate either q or p separately from catch and effort data.


    Effect of aggregation
 Top
 Introduction
 Relative fish density...
 Effect of aggregation
 Effect of effort distribution
 Analysis of catch and...
 Analysis of commercial catch...
 Discussion
 Appendix
 References
 
In the appendix, following Swain and Sinclair (1994), we consider one-dimensional examples of the function f(x) that include a spread parameter {sigma} to quantify the degree of fish aggregation. We quantify the degree of aggregation in a way that is not dependent on such a form, so define a function G(r), the proportion of fish distributed in areas with fish density <r. This is given by


Formula 015M4

(4)

We can interpret G(r) as the cumulative density function of a random variable R. The actual fish densities in the sea can be thought of as samples from the probability density of R. Equivalently, the relative fish densities, f(x), can be thought of as realizations of R/N. We show that the variance of R/N provides a measure of the degree of fish aggregation. After some algebra, one can show that the nth moment of R/N is given by


Formula 015M5

(5)
When f(x) is the standard normal density in one dimension, we find, using the identity Var[X]{equiv}E[X2]–E[X]2, that


Formula 015M6

(6)
As this is inversely proportional to {sigma}2, and is therefore a measure of fish aggregation, we propose using the variance of R/N for such an aggregation measure in general. In fact, we define a dimensionless aggregation index, thus:


Formula 015M7

(7)
which is the variance of the fish density relative to the mean fish density N/A.


    Effect of effort distribution
 Top
 Introduction
 Relative fish density...
 Effect of aggregation
 Effect of effort distribution
 Analysis of catch and...
 Analysis of commercial catch...
 Discussion
 Appendix
 References
 
We now consider the effect on catchability when the effort distribution varies for a given fish distribution. We assume that the effort distribution is a power function of the fish density, i.e. e(x){propto}(f(x)){gamma}. The constant of proportionality is E/{int}(f(x)){gamma}dx, to obtain a total effort of E. Without loss of generality, we will simply assume that E=1. We will refer to {gamma} as the knowledge parameter. Higher values of {gamma} provide more effort in higher density areas and result in greater catches. The three models discussed by Swain and Sinclair (1994) (appendix) correspond to the cases when {gamma}=0, 1, and {infty}, respectively. The general expression for the catchability index is, from Equations (2) and (3),


Formula 015M8

(8)

For random fishing, {gamma}=0, and the catchability index becomes


Formula 015M9

(9)
where NA is the population within A. If the fish density is uniform within A, then, irrespective of {gamma}, {kappa}{gamma}={kappa}0. This implies that knowledge is ineffective at increasing catchability if there is no aggregation.

The increase in catchability relative to random fishing is given by the ratio


Formula 015M10

(10)
An alternative expression for l({gamma}) can be found in terms of the fish density distribution, using the moment Equation (5); thus,


Formula 015M11

(11)
This form will prove useful in estimating l({gamma}) assuming known {gamma}.

In the case when f(x) is the normal density function, and A={x:|x|≤ z{alpha}{sigma}},


Formula 015M12

(12)
in which z{alpha} is the 100(1–{alpha}/2) percentile value for a standard normal distribution. For example, z0.05=1.96, and A becomes the highest density region containing 95% of the population.

Figure 1 shows how the catchability index changes with the knowledge parameter {gamma}. We plot the catchability index {kappa}{gamma} for two different aggregation models (spread {sigma}=1 and 2) and four different effort models (X=1, 10, {infty}, and z0.05{sigma}). The last case means that fishing is restricted to the most abundant region containing 95% of the population. Catchability always increases with knowledge up to an asymptote that depends on spread, but not on fishing area. This asymptotic value corresponds to fishing at the point of maximum fish density. Catchability decreases as the spatial limits widen. The case {gamma}=0 corresponds to uniform fishing within(–X, X). Note that for an infinite area {kappa}0=0, one must have a finite area to obtain a non-zero catchability for uniform fishing.


Figure 0151
View larger version (10K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 1. Catchability index {kappa}{gamma} vs. {gamma} for two different aggregation models (spread {sigma}= 1 and 2) and four different effort models (X=1, 10, {infty}, and z0.05{sigma}). The units of {kappa}{gamma} are such that a=1. The fish distribution f(x) follows a normal distribution with mean zero and standard deviation {sigma}, and the effort distribution is proportional to [f(x)]{gamma} but truncated to be within (–X, X).

 

    Analysis of catch and effort data
 Top
 Introduction
 Relative fish density...
 Effect of aggregation
 Effect of effort distribution
 Analysis of catch and...
 Analysis of commercial catch...
 Discussion
 Appendix
 References
 
In the analysis of catch and effort data, we restrict attention to the part of the stock lying within the fishery area A, and treat it as a single stock confined to that area, so precluding migration into or out of the area. This is a nominal fishing area, because parts of it may be unfished or only rarely fished for some times. Abundance N is the same as NA, the abundance within A, and f(x) is the density restricted to A.

Let (ci,t, ei,t) be the catch and effort data collected at time t and in grid (subarea) i of area {Delta}A. If the abundance in grid i is ni,t, we have, from (3),


Formula 015M13

(13)
We can estimate this by plugging in estimates for C and N. We can estimate ni,t in fished grids, because from (1), E[ci,t]=pni,t ei,t/{Delta}A, so that ni,t=ui,t{Delta}A/p, where ui,t is the cpue in grid i at time t. In unfished grids, we assume that ni,t equals the average over fished cells. We then have N=utA/p, where ut is the mean of ui,t. For C, the expected total catch, the estimator is simply the observed total catch. Therefore, a plug-in estimator for {kappa} is


Formula 015M14

(14)
where ût is the overall cpue. Note that ût can be thought of as the estimated cpue under a sampling scheme weighted by effort (because û={sum}ue/{sum}e), whereas ut is the estimated cpue under a uniform sampling scheme. The abundance in unfished grids is actually likely to be less than in fished grids, so N is over-biased and Formula underbiased. Later, in Discussion, we suggest how this limitation might be overcome. As, from the previous paragraph, the fish density satisfies ni,t/{Delta}A=E[ui,t]/p, the empirical density distribution of R at time t takes the set of discrete values of ui,t/{sum}k uk,t at ui,t/p. Therefore, E[R{gamma}] can be estimated by


Formula 015M15

(15)
which, using Equation (11), leads to an estimator of l({gamma}) as


Formula 015M16

(16)
The instantaneous catchability p, which links cpue to abundance, does not depend on effort distribution. We can therefore estimate it under the assumption of uniform fishing using the equation l(0)=Aut/Np=1; then, p is simply the ratio of mean cpue to mean fish density N/A. This agrees with the relationship we obtained previously from Equation (3) for uniform fishing: pa=qA.

The quantity l({gamma}) measures the effect of aggregation on catchability when the knowledge parameter about the stock is fixed as {gamma}. Substituting for p in Equation (16) gives


Formula 015M17

(17)
We can use the first and second moment estimates of R from Equation (5) to estimate the aggregation index {lambda}, thus,


Formula 015M18

(18)
because Var[R]=E[R2]–E[R]2.

The knowledge parameter is estimated by regressing log effort against log cpue and taking the coefficient of log cpue. Actually we use a GLM in which effort eijk follows a Poisson distribution with extra variation, and cpue uijk is used as a covariate, so is assumed known. (It might seem more logical to treat effort as the covariate, and cpue as the response, but then the coefficient to be estimated would be 1/{gamma}, which is unbounded around {gamma}= 0, the expected range of {gamma}.)

The details of the model are


Formula 015M19

(19)
where eijk is the effort in grid i for year j and month k, uijk the cpue, {phi} the dispersion parameter, and {xi} an exponent to be estimated by inspection of the residuals. The dispersion parameter is estimated from the Pearson chi-squared statistic.

Estimates Formula of the knowledge parameter can be plugged into expression (17) to give a second estimate, l(Formula), for the catchability index {kappa}{gamma}.


    Analysis of commercial catch and effort data from Australia's northern prawn fishery
 Top
 Introduction
 Relative fish density...
 Effect of aggregation
 Effect of effort distribution
 Analysis of catch and...
 Analysis of commercial catch...
 Discussion
 Appendix
 References
 
We now consider the commercial catch and effort data from the tiger prawn fisheries in the Gulf of Carpentaria (1970–1996), and we apply the theory of fish aggregation and effort allocation to see how the average cpue can be influenced. The Gulf is divided into large-scale management areas, shown in Figure 2. We focus on the four most important tiger prawn fishery areas: Groote, Vanderlins, Mornington, and Weipa. The finest practical scale on which to estimate fish density is the grid level (6 nautical mile squares), which are the small squares on Figure 2. For a temporal scale, we use month (i.e. we sum catch and effort within each month). This is because prawn distribution (and thus effort) can change markedly from one month to the next.


Figure 0152
View larger version (69K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 2. The Gulf of Carpentaria, Australia, showing large-scale management areas and small-scale grids of 6 nautical miles.

 
We also need to decide on the scale of the stock. Although the large-scale areas are in fact subdivided into smaller management areas, such areas are probably too small to contain a stock, and they have too few grids for reliable estimation. We therefore work on the large-scale areas (or amalgamations thereof): we assume that these areas completely contain the prawn stock whose catchability we are estimating.

The pattern of effort in Mornington shows a marked decline in September and October, suggesting that the fleet moves to neighbouring Vanderlins. We therefore carry out catchability analysis on amalgamated regions as well as on the basic areas. Such regions are SW Gulf (Mornington and Vanderlins), W Gulf (Groote, Mornington, and Vanderlins), and all four areas combined.

In 1980/1981, there was a major expansion of the fishery, so many new 6 nautical mile grids were visited. We restrict attention to the most important months for the tiger prawn fishery, namely August, September, and October. Before the seasonal closures were introduced in 1988, there was year-round effort, but this effort was still concentrated in the months August–November when the prawns are at their largest.

Catchability index {kappa}
Because the northern prawn fishery is a trawl fishery, with effort measured in swept area, we define for convenience our units of effort such that a/A=1. We have computed the catchability index, Formula=ût/ut, for all three months within each year for each region. Grids with only a single day's effort have been excluded. Table 2 shows how many grids are included for each region, month, and year. The results are shown in Figure 3. The dotted horizontal line at {kappa}=1 indicates no density-dependent effect.


Figure 0153
View larger version (38K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 3. Catchability index estimates, Formula, for each month, year, and region in units where a=A. Also shown (dotted) is Formula for the whole tiger prawn season.

 


View this table:
[in this window]
[in a new window]

 
Table 2. Number of grids in a region for which total effort is at least 2 d.

 
Some years are missing because of the lack of grids (Table 2). The dotted line is the index for all months combined. Groote, Vanderlins, and Mornington appear to have {kappa}>1 (i.e. a density-dependent effect) since the period 1980–1985. This is quite clear when the regions are combined (W Gulf). For Mornington, there is a reduction in effort in October, leading to a more variable {kappa} across months. This becomes less variable across months and years when we combine Mornington with Vanderlins (SW Gulf); we consistently get {kappa}{approx}1.1 since 1987. Groote has roughly the same mean, but greater variation across years. The estimate from the three months combined tends to be more consistent. Weipa has no consistent density-dependent effect.

Aggregation index {lambda}
Estimates of {lambda} from Equation (18) are shown in Figure 4 for each region and year. The estimates are on a logarithmic scale. Three curves are given per panel, one for each month. With the exception of Weipa, all regions show a steady decline in aggregation since 1980. As for Mornington, the estimates for October are highly variable.


Figure 0154
View larger version (32K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 4. Aggregation index estimates, Formula, for each month, year, and region. Also shown (dotted) is Formula for the whole tiger prawn season.

 
To explain how these estimates arise, Figure 5 shows details of the abundance distribution in two contrasting years for the W Gulf region. In the left-hand plots we have the relative fish density d(x)/(N/A) [or Af(x)] plotted against relative area x. The range under each curve equals the total area fished. In 1995, this area was about half that in 1985. The departure of the fish density from the average value, as indicated by the dashed line, is generally larger in 1985, which suggests that that year had greater aggregation. This is confirmed in the right-hand plots, which show the empirical density of R, as introduced before Equation (15), after gaussian smoothing. The variance is higher for 1985, so the distribution is indeed more aggregated. There is also slightly greater aggregation in August of both years.


Figure 0155
View larger version (16K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 5. (Left panel) Relative fish density d(x)/(N/A) plotted against fished area x, sorted by decreasing density, for each month in 1985 and 1995. (Right panel) Probability density of relative fish density for the same years. The area under the curve within a fish-density interval denotes the proportion of fish in areas with densities inside that interval.

 
However, this apparent decrease in aggregation with year may be largely an artefact of the shrinkage of the fishery, rather than an actual change in fish distribution. This is because the aggregation parameter, which is derived from the empirical distribution of R, is obtained on fished grids only. The fish density in unfished grids is likely to be less than in fished grids, so the fish distribution is more aggregated than it appears from the fished grids. This downward bias will be greater for the later years with fewer fished grids (Table 2). We revisit the issue of bias later, in Discussion.

Knowledge parameter {gamma}
Using a range of values for {xi} in the GLM (19), {xi}= 2 produces the flattest variance–mean relationship in the Pearson residuals. Grids with a single day's effort can have huge influence on the estimates (Pollock et al., 1997), so we restricted the analysis to grids with ≥2 d effort. The analysis has been carried out using the SAS procedure GENMOD.

Figure 6 shows point estimates of {gamma}jk{equiv}{alpha}j1k1+{delta}jk1. The dotted line at {gamma}=0 corresponds to no targeting. There appears to have been definite targeting since the mid-1980s in Groote, Mornington, and Vanderlins, but not in Weipa. The pattern tends to be fairly consistent across months for Groote and Vanderlins, and the consistency is improved for Vanderlins and Mornington when they are combined (SW Gulf).


Figure 0156
View larger version (28K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 6. Knowledge parameter estimates, Formula, for each month, year, and region.

 
To assess the accuracy of the estimates, we show standard errors for one of the months (October in Weipa, August elsewhere) in Figure 7. The choice of October for Weipa is based on there being more available grids. It is clear that, before the mid-1980s, there was little evidence of targeting anywhere. Although the statistical power is low, this may actually be a real effect owing to the absence of GPS technology. On the other hand, the positive {gamma} for areas in the W Gulf appears a significant effect. The dispersion parameter estimates range from 0.7 to 1.4, the larger values corresponding to "All areas" and Weipa. This further confirms that the relationship between cpue and effort differs between Weipa and the W Gulf.


Figure 0157
View larger version (23K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 7. Estimates of {gamma} in August for each region and each year, with standard errors shown by bars. For Weipa, the {gamma} estimate is for October.

 
The standard errors are governed by the number of grids fished in the region (Table 2). There was a lull in effort in Mornington and Vanderlins around 1975, so fewer grids were visited. Also the level of effort in Weipa was lower than in the other regions for most years. The fewer grids in earlier years may also be due to lack of reporting, because it was not compulsory to record grid-level position in the logbooks.

Using the estimated knowledge parameters for each month, year, and region, we evaluated the alternative estimates of the catchability indices, l(Formula), using Equation (17). Figure 8 shows the scatterplots of l(Formula) against Formula for each month, year, and region. They appear to agree, indicating that the power law for the effort distribution is a reasonable approximation of targeting in fishing.


Figure 0158
View larger version (15K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 8. Scatterplots of l(Formula) against Formula for each month, year, and region; the units of Formula are such that a=A. Different months are denoted by a different symbol.

 

    Discussion
 Top
 Introduction
 Relative fish density...
 Effect of aggregation
 Effect of effort distribution
 Analysis of catch and...
 Analysis of commercial catch...
 Discussion
 Appendix
 References
 
We have demonstrated how a combination of stock clumping and fisher knowledge influences the relationship between catchability and stock size. The overall cpue based on lumped catch and effort over a stock's area is a good index of overall catchability, but the average of the cpues calculated from small units where the fishing is random is a better index for measuring stock abundance. In highly complex fisheries where there are numerous innovations in technology, it can be quite easy (it is almost inevitable) to overlook confounding effects even with a GLM approach with many parameters. Significant coefficients can still be misleading because of confounding (Bishop, 2006).

We further investigated the effect of aggregation and targeting on catchability, following Swain and Sinclair (1994), and suggested three indices to measure the extent of aggregation, targeting, and changes in catchability. Simple estimation methods are also proposed. As catch and effort data are commonly collected in fisheries, the proposed methods can be applied easily to those data to quantify possible effects on catchability.

The various indices presented here are based only on grids fished, because cpue information is not available for unfished grids. However, the population in reality extends over both fished and unfished grids, so the indices may be biased. To base indices on all grids requires assumptions about the cpue in unfished grids. The estimate of {kappa} is unbiased if mean cpue for unfished grids equals ut, and the estimate of {lambda} is unbiased if the fish density distribution is the same in both fished and unfished areas. Because fish density is likely to be less in unfished areas, both {kappa} and {lambda} are probably underestimated. The same estimates of {gamma} would be obtained if the unfished grids were included with cpue set to zero. As fish density is likely to be somewhat greater than zero in those grids, {gamma} is also probably underestimated; the same applies to l({gamma}). The effect of ignoring the unfished grids is consequently to underestimate all four indices. It is therefore of interest to further investigate the effect of these "missing" grids. Other factors, such as the interference competition effect reported by Swain and Wade (2003), may distort the validity of cpue as an abundance index. However, for the northern prawn fishery there is no evidence to suggest that interference competition is important.

This work has relevance to the estimation of fishing power. In studies by Robins et al. (1998), using regression models, and Bishop et al. (2000), using generalized estimating equations, the effect on fishing power of vessel characteristics such as gear size and the presence of GPS plotters during the transition period around 1990 was estimated. Those authors took account of density-dependent catchability by allowing catch to be proportional to a constant power of effort; both works reported the estimate as 1.07, implying that cpue is higher in regions of greater effort. We have confirmed this overall result, but we have also shown that the density-dependence varies both temporally and spatially. Our method filters out the effects attributable to intrinsic vessel characteristics, manifest through the instantaneous catchability p, and captures the effects attributable to improved search capability, manifest through {kappa}. This allows us to build more informed models of fishing power, producing more reliable estimates and sounder advice to fisheries management ().


    Appendix
 Top
 Introduction
 Relative fish density...
 Effect of aggregation
 Effect of effort distribution
 Analysis of catch and...
 Analysis of commercial catch...
 Discussion
 Appendix
 References
 
Comparison with Swain and Sinclair's (1994) findings
Swain and Sinclair (1994) considered three models for fishing effort.

  1. Uniform fishing, or equivalently fishing is random over the area (–xT, xT) where the fish density exceeds some threshold.
  2. Fishing effort at position x is proportional to the density at x for xisin(–{infty},{infty}).
  3. Fishing is only at the location where the fish density is maximum, i.e. at x=0.
In model (i), it is assumed that fishers have no knowledge about the distribution of the fish stock (except that it exceeds some threshold fish density within a finite area). The corresponding catchability can therefore be regarded as a benchmark. In model (ii), the fishers have some knowledge about the distribution of the stock. Model (iii) is the extreme case in which the fishers know exactly where the hotspot is. We would expect {kappa} to be smallest for model (i) and largest for model (iii). Model (ii) is probably the most realistic case.

Assuming f(x)={sigma}–1{phi}(x/{sigma}), the normal density function with a variance {sigma}2, Swain and Sinclair (1994) obtained the corresponding catchability indices for the three effort models as


Formula 015MA1

(A1)
Note that in (i) the fishing area (–xT, xT) depends on both relative fish density distribution and total abundance N. This is why Swain and Sinclair (1994) concluded that spatially uniform changes in abundance also influence catchability to some extent. In fact, such an effect is due to changes in distribution of fishing effort.

We now consider how the catchability changes when the fish density distribution changes for a given fishing pattern e(x). The e(x) is assumed to be

  1. a constant over a fixed area (–X, X), where X> 0,
  2. e(x)={phi}(x),
  3. e(x)={phi}(x)/(2{Phi}(x)–1) for xisin(–X, X).
Unlike case (i) considered by Swain and Sinclair (1994), the fishing area here is fixed in (a). Case (b) is slightly different from case (ii) of Swain and Sinclair (1994), in which the effort pattern changes as fish density changes. In the cases here, fixing the fishing effort function and allowing the fish density distribution to change, q shows the effect of aggregation on catchability.

If we assume that f(x)={sigma}–1{phi}(x/{sigma}), we have the corresponding catchability indices for the above three fishing effort functions:


Formula 015MA2

(A2)
Here, {Phi} (·) is the standard normal cumulative density function. Similar algebra leads to


Formula 015MA3

(A3)

We demonstrate how catchability varies with aggregation in Figure A1. Figure A1(a) shows the fish density distribution for a range of values of spread {sigma}. More highly aggregated distributions have smaller values of {sigma}. The effort pattern in (b) is identical to the distribution labelled 1.0. In (a) and (c), we let X=1, as shown by the dashed line. Figure A1(b) shows the catchability index plotted against spread for the three effort patterns. The catchability always decreases with spread. Case (a) is not very sensitive to {sigma} for small spreads until {sigma}{approx}X, because almost the entire population lies inside X. On the other hand, targeted patterns are more sensitive to {sigma}. The catchability is higher for (c) than for (a) because, although both patterns cover the same area, the effort pattern (c) exploits the aggregation more. Similarly, {kappa}c is higher than {kappa}b because the effort used in (b) outside X is used more efficiently in (c) inside X. The catchability of (a) relative to (b) depends on X; for sufficiently large X, {kappa}a is less than {kappa}b, and for sufficiently small X, {kappa}a is greater than {kappa}b.


Figure 0151
View larger version (7K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure A1. (a) Fish density distribution for a range of values of spread {sigma} over the half-range x>0. Overlaid (dashed line) are the effort patterns (a)–(c). (b) catchability index {kappa} vs. spread for each effort pattern. The units of {kappa} are such that a=1.

 
The use of a one-dimensional function f(x) may appear overly restrictive. However, note that the results for f(x) would hold for any related distribution in which the x-axis was segmented and arbitrarily rearranged. One such distribution is that where the x values are sorted on density. The function f(x) is therefore a prototype for many actual spatial arrangements, so the particular form is not as restrictive and artificial as it may at first appear. Moreover, the x variable of the sorted form of f(x) can also be interpreted as cumulative area, so f(x) can represent two-dimensional distributions too.


    Acknowledgements
 
This research project was partly supported by the Fisheries Research and Development Corporation of Australia. We thank Janet Bishop and two referees for their helpful suggestions and comments on an earlier draft.


    References
 Top
 Introduction
 Relative fish density...
 Effect of aggregation
 Effect of effort distribution
 Analysis of catch and...
 Analysis of commercial catch...
 Discussion
 Appendix
 References
 

    Baranov T. I. (1918) On the question of the biological basis of fisheries. Nauchnyi Issledovatelskii Ikhtiologicheskii Institut Isvestia 1:81–128.

    Bishop J. (2006) Standardizing fishery-dependent catch and effort data in complex fisheries with technology change. Reviews in Fish Biology and Fisheries 16:21–38.[CrossRef]

    Bishop J., Die D. J., Wang Y-G. (2000) A generalized estimating equations approach for analysis of the impact of new technology on a trawl fishery. Australian and New Zealand Journal of Statistics 42:159–177.

    Campbell R. A. (2004) CPUE standardisation and the construction of indices of stock abundance in a spatially varying fishery using general linear models. Fisheries Research 70:209–227.[CrossRef]

    Crecco V. and Overholtz W. J. (1990) Causes of density-dependent catchability for Georges Bank haddock Melanogrammus aeglefinus. Canadian Journal of Fisheries and Aquatic Sciences 47:385–394.

    Francis R. I. C. C., Hurst R. J., Renwick J. A. (2003) Quantifying annual variation in catchability for commercial and research fishing. Fishery Bulletin US 101:293–304.

    Gulland J. A. (1983) Fish Stock Assessment(John Wiley and Sons, Chichester).

    Hilborn R. and Walters C. J. (1992) Quantitative Fisheries Stock Assessment: Choice, Dynamics and Uncertainty(Chapman and Hall, New York).

    Paloheimo L. E. and Chen Y. (1993) Estimation of effective effort from catch-at-age data. Canadian Journal of Fisheries and Aquatic Sciences 50:2421–2428.

    Paloheimo J. E. and Dickie L. E. (1964) Abundance and fishing success. Rapports et Procés-Verbaux des Réunions du Conseil International pour l'Exploration de la Mer 155:152–163.

    Pollock K. H., Hoenig J. M., Jones C. M., Robson D. S., Greene C. J. (1997) Catch rate estimation for roving and access point surveys. North American Journal of Fisheries Management 17:11–19.

    Richards L. J. and Schnute J. T. (1992) Statistical models for estimating CPUE from catch and effort data. Canadian Journal of Fisheries and Aquatic Sciences 49:1315–1327.

    Robins C. M., Wang Y-G., Die D. J. (1998) The impact of global positioning systems and plotters on fishing power in the northern prawn fishery, Australia. Canadian Journal of Fisheries and Aquatic Sciences 55:1645–1651.[CrossRef]

    Salthaug A. and Aanes S. (2003) Catchability and the spatial distribution of fishing vessels. Canadian Journal of Fisheries and Aquatic Sciences 60:259–268.[CrossRef]

    Shardlow T. F. (1993) Components analysis of a density-dependent catchability coefficient in a salmon hook and line fishery. Canadian Journal of Fisheries and Aquatic Sciences 50:513–520.

    Somers I. F. and Wang Y-G. (1997) A simulation model for evaluating seasonal closures in Australia's multispecies northern prawn fishery. North American Journal of Fisheries Management 17:114–130.[CrossRef]

    Swain D. P. and Morin R. (1996) Relationships between geographic distribution and abundance of American plaice (Hippoglossoides platessoides) in the southern Gulf of St Lawrence. Canadian Journal of Fisheries and Aquatic Sciences 53:106–119.[CrossRef]

    Swain D. P. and Sinclair A. F. (1994) Fish distribution and catchability: what is the appropriate measure of distribution? Canadian Journal of Fisheries and Aquatic Sciences 51:1046–1054.

    Swain D. P. and Wade E. J. (2003) Spatial distribution of catch and effort in a fishery for snow crab (Chionoecetes opilio): tests of predictions of the ideal free distribution. Canadian Journal of Fisheries and Aquatic Sciences 60:897–909.[CrossRef]

    Walters C. (2003) Folly and fantasy in the analysis of spatial catch rate data. Canadian Journal of Fisheries and Aquatic Sciences 60:1433–1436.[CrossRef]

    Wang Y-G. and Die D. (1996) Stock-recruitment relationships of the tiger prawns P. esculentus and P. semisulcatus in the Australian northern prawn fishery. Marine and Freshwater Research 47:87–95.[CrossRef][Web of Science]

    Xiao Y. (2005) Catch equations: restoring the missing terms in the nominally generalized Baranov catch equation. Ecological Modelling 181:535–556.[CrossRef]

    Xiao Y. (2006) Catch equations: calculating the instantaneous rate of fishing mortality from catch and back. Ecological Modelling 193:225–252.[CrossRef]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
ICES J. Mar. Sci.Home page
J. J. Poos, F. J. Quirijns, and A. D. Rijnsdorp
Spatial segregation among fishing vessels in a multispecies fishery
ICES J. Mar. Sci., August 12, 2009; (2009) fsp219v1.
[Abstract] [Full Text] [PDF]


Home page
ICES J. Mar. Sci.Home page
S. Mayfield, R. McGarvey, I. J. Carlson, and C. Dixon
Integrating commercial and research surveys to estimate the harvestable biomass, and establish a quota, for an "unexploited" abalone population
ICES J. Mar. Sci., October 1, 2008; 65(7): 1122 - 1130.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
64/1/178    most recent
fsl015v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Ellis, N.
Right arrow Articles by Wang, Y.-G.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Ellis, N.
Right arrow Articles by Wang, Y.-G.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?