Effects of fishing on growth traits: a simulation analysis

Abstract—Fisheries often target individuals based on size. Size-selective fishing can create selection differentials on life-history traits and, when those traits have a genetic basis, may cause evolution. The evolution of life history traits affects potential yield and sustainability of fishing, and it is therefore an issue for fishery management. Yet fishery managers usually disregard the possibility of evolution, because little guidance is available to predict evolutionary consequences of management strategies. We attempt to provide some generic guidance. We develop an individual-based model of a population with overlapping generations and continuous reproduction. We simulate model populations under size-selective fishing to generate and quantify selection differentials on growth. The analysis comprises a variety of common life-history and fishery characteristics: variability in growth, correlation between von Bertalanffy growth parameters (K and L∞), maturity rate, natural mortality rate (M), M/K ratio, duration of spawning season, fishing mortality rate (F), maximum size limit, slope of selectivity curve, age at 50% selectivity, and duration of fishing season. We found that each characteristic affected the magnitude of selection differentials. The most vulnerable stocks were those with a short spawning or fishing season. Under almost all life-history and fishery characteristics examined, selection differentials created by realistic fishing mortality rates are considerable.

Simulation of Tail Weight Distributions in Biological Year 1986–2006 Landings of Brown Shrimp, Farfantepenaeus aztecus, from the Northern Gulf of Mexico Fishery

Size distribution within re- ported landings is an important aspect of northern Gulf of Mexico penaeid shrimp stock assessments. It reflects shrimp population characteristics such as numerical abundance of various sizes, age structure, and vital rates (e.g. recruitment, growth, and mortality), as well as effects of fishing, fishing power, fishing practices, sampling, size-grading, etc. The usual measure of shrimp size in archived landings data is count (C) the number of shrimp tails (abdomen or edible portion) per pound (0.4536 kg). Shrimp are marketed and landings reported in pounds within tail count categories. Statistically, these count categories are count class intervals or bins with upper and lower limits expressed in C. Count categories vary in width, overlap, and frequency of occurrence within the landings. The upper and lower limits of most count class intervals can be transformed to lower and upper limits (respectively) of class intervals expressed in pounds per shrimp tail, w, the reciprocal of C (i.e. w = 1/C). Age based stock assessments have relied on various algorithms to estimate numbers of shrimp from pounds landed within count categories. These algorithms required un- derlying explicit or implicit assumptions about the distribution of C or w. However, no attempts were made to assess the actual distribution of C or w. Therefore, validity of the algorithms and assumptions could not be determined. When different algorithms were applied to landings within the same size categories, they produced different estimates of numbers of shrimp. This paper demonstrates a method of simulating the distribution of w in reported biological year landings of shrimp. We used, as examples, landings of brown shrimp, Farfantepenaeus aztecus, from the northern Gulf of Mexico fishery in biological years 1986–2006. Brown shrimp biological year, Ti, is defined as beginning on 1 May of the same calendar year as Ti and ending on 30 April of the next calendar year, where subscript i is the place marker for biological year. Biological year landings encompass most if not all of the brown shrimp life cycle and life span. Simulated distributions of w reflect all factors influencing sizes of brown shrimp in the landings within a given biological year. Our method does not require a priori assumptions about the parent distributions of w or C, and it takes into account the variability in width, overlap, and frequency of occurrence of count categories within the landings. Simulated biological year distributions of w can be transformed to equivalent distributions of C. Our method may be useful in future testing of previously applied algorithms and development of new estimators based on statistical estimation theory and the underlying distribution of w or C. We also examine some applications of biological year distributions of w, and additional variables derived from them.