TUNPOP, a computer simulation model of the yellowfin tuna population and the surface tuna fishery of the Eastern Pacific Ocean

ENGLISH: Mathematical documentation of TUNP0P, an age-structured computer simulation model of the yellowfin tuna population and surface tuna fishery of the eastern Pacific Ocean, is described. Example runs of the model are presented and discussed, and the sensitivity of the model output to changes in various parameters is examined. SPANISH: Se describe la documentación matemática de TUNP0P, un modelo computador de simulación basado en la edad de la población del atún aleta amarilla y de la pesca atunera epipelágíca del Océano Pacífico oriental. Se presentan y se discuten ejemplos de las pasadas del modelo, y se examina la sensibilidad de los resultados de salida con relación a los cambios de varios parámetros. (PDF contains 47 pages.)

TUNPOP: A simulation of the dynamics and structure of the yellowfin tuna stock and surface fishery of the Eastern Pacific Ocean

ENGLISH: Three distinct versions of TUNP0P, an age-structured computer simulation model of the eastern Pacific yellowfin tuna, Thunnus albacores, stock and surface tuna fishery, are used to reveal mechanisms which appear to have a significant effect on the fishery dynamics. Real data on this fishery are used to make deductions on the distribution of the fish and to show how that distribution might influence events in the fishery. The most important result of the paper is that the concept of the eastern Pacific yellowfin tuna stock as a homogeneous unit is inadequate to represent the recent history of the fishery. Inferences are made on the size and distribution of the underlying stock as well as its potential yield to the surface fishery as a result of alterations in the level and distribution of the effort. SPANISH: Se han empleado tres versiones diferentes de TUNP0P, un modelo de simulación de la computadora (basado en la estructura de la edad) de la población y la pesca epipelágica del atún aleta amarilla, Tbunnus albacares, del Pacífico oriental, para revelar los mecanismos que parecen tener un efecto importante en la dinámica pesquera. Se emplean los datos verdaderos de esta pesca para hacer deducciones sobre la distribución de los peces y para mostrar cómo puede influir esta distribución en los eventos de pesca. La conclusión más importante de este estudio es que el concepto de que la población del aleta amarilla del Pacífico oriental es una unidad homogénea, es inadecuado para representar la historia reciente de pesca. Se teoriza sobre la talla y distribución de la población subyacente como también sobre su producción potencial en la pesca epipelágica al cambiar el nivel y distribución del esfuerzo.

Simulation stochastique des structures démographiques. Application aux populations d'albacores, Thunnus albacares

A computer program has been written in order to generate a population of fishes following a Von Bertalanffy growth curve with a random Gaussian variability for birth dates and growth parameters K and L ∞. Standard deviations for these 3 parameters are chosen separately for each run. Fishing and natural mortalities are applied to this population. Using as an input parameters usually taken for yellowfin in the eastern Atlantic, the simulation suggests a standard deviation between 1 and 2 months for the birth dates in this population. It also indicates that increasing levels of fishing mortalities must produce a better agreement between age and length for the larger fish.

Evaluating methods for estimating rare events with zero-heavy data: a simulation model estimating sea turtle bycatch in the pelagic longline fishery

Estimating rare events from zero-heavy data (data with many zero values) is a common challenge in fisheries science and ecology. For example, loggerhead sea turtles (Caretta caretta) and leatherback sea turtles (Dermochelys coriacea) account for less than 1% of total catch in the U.S. Atlantic pelagic longline fishery. Nevertheless, the Southeast Fisheries Science Center (SEFSC) of the National Marine Fisheries Service (NMFS) is charged with assessing the effect of this fishery on these federally protected species. Annual estimates of loggerhead and leatherback bycatch in a fishery can affect fishery management and species conservation decisions. However, current estimates have wide confidence intervals, and their accuracy is unknown. We evaluate 3 estimation methods, each at 2 spatiotemporal scales, in simulations of 5 spatial scenarios representing incidental capture of sea turtles by the U.S. Atlantic pelagic longline fishery. The delta-log normal method of estimating bycatch for calendar quarter and fishing area strata was the least biased estimation method in the spatial scenarios believed to be most realistic. This result supports the current estimation procedure used by the SEFSC.

A low-level geographic information system for coastal zone management with applications to Brunei Darussalam: Part III:: simulation and tracking oil spills

B:RUN is a low-level GIS software designed to help formulate options for the management of the coastal zone of Brunei Darussalam. This contribution presents the oil spill simulation module of B:RUN. This simple module, based largely on wind and sea surface current vector parameters, may be helpful in formulating relevant oil spill contingency plans. It can be easily adapted to other areas, as can the B:RUN software itself.

A simulation model of the dynamics of Peruvian anchoveta (Engraulis ringens)

The article describes the key elements of a model simulating the dynamics of the anchoveta (Engraulis ringens) in the Peruvian upwelling system (4 degrees to 14 degrees South). This model, based on coupled differential equations, is parametrized mainly using empirical data and functional relationships presented in two volumes issued by ICLARM in 1987 and 1989, and may thus be viewed as test of the hypotheses presented therein. Results to date suggest that present knowledge of mechanisms controlling the anchoveta stock is essentially consistent, and sufficient to build a model reflecting essential features of the stock biomass and recruitment dynamics.

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.