An evaluation of length-frequency sampling procedures and subsequent data analysis for purse seine-caught yellowfin tuna in the Eastern Pacific Ocean

ENGLISH: A two-stage sampling design is used to estimate the variances of the numbers of yellowfin in different age groups caught in the eastern Pacific Ocean. For purse seiners, the primary sampling unit (n) is a brine well containing fish from a month-area stratum; the number of fish lengths (m) measured from each well are the secondary units. The fish cannot be selected at random from the wells because of practical limitations. The effects of different sampling methods and other factors on the reliability and precision of statistics derived from the length-frequency data were therefore examined. Modifications are recommended where necessary. Lengths of fish measured during the unloading of six test wells revealed two forms of inherent size stratification: 1) short-term disruptions of existing pattern of sizes, and 2) transition zones between long-term trends in sizes. To some degree, all wells exhibited cyclic changes in mean size and variance during unloading. In half of the wells, it was observed that size selection by the unloaders induced a change in mean size. As a result of stratification, the sequence of sizes removed from all wells was non-random, regardless of whether a well contained fish from a single set or from more than one set. The number of modal sizes in a well was not related to the number of sets. In an additional well composed of fish from several sets, an experiment on vertical mixing indicated that a representative sample of the contents may be restricted to the bottom half of the well. The contents of the test wells were used to generate 25 simulated wells and to compare the results of three sampling methods applied to them. The methods were: (1) random sampling (also used as a standard), (2) protracted sampling, in which the selection process was extended over a large portion of a well, and (3) measuring fish consecutively during removal from the well. Repeated sampling by each method and different combinations indicated that, because the principal source of size variation occurred among primary units, increasing n was the most effective way to reduce the variance estimates of both the age-group sizes and the total number of fish in the landings. Protracted sampling largely circumvented the effects of size stratification, and its performance was essentially comparable to that of random sampling. Sampling by this method is recommended. Consecutive-fish sampling produced more biased estimates with greater variances. Analysis of the 1988 length-frequency samples indicated that, for age groups that appear most frequently in the catch, a minimum sampling frequency of one primary unit in six for each month-area stratum would reduce the coefficients of variation (CV) of their size estimates to approximately 10 percent or less. Additional stratification of samples by set type, rather than month-area alone, further reduced the CV's of scarce age groups, such as the recruits, and potentially improved their accuracy. The CV's of recruitment estimates for completely-fished cohorts during the 198184 period were in the vicinity of 3 to 8 percent. Recruitment estimates and their variances were also relatively insensitive to changes in the individual quarterly catches and variances, respectively, of which they were composed. SPANISH: Se usa un diseño de muestreo de dos etapas para estimar las varianzas de los números de aletas amari11as en distintos grupos de edad capturados en el Océano Pacifico oriental. Para barcos cerqueros, la unidad primaria de muestreo (n) es una bodega de salmuera que contenía peces de un estrato de mes-área; el numero de ta11as de peces (m) medidas de cada bodega es la unidad secundaria. Limitaciones de carácter practico impiden la selección aleatoria de peces de las bodegas. Por 10 tanto, fueron examinados los efectos de distintos métodos de muestreo y otros factores sobre la confiabilidad y precisión de las estadísticas derivadas de los datos de frecuencia de ta11a. Se recomiendan modificaciones donde sean necesarias. Las ta11as de peces medidas durante la descarga de seis bodegas de prueba revelaron dos formas de estratificación inherente por ta11a: 1) perturbaciones a corto plazo en la pauta de ta11as existente, y 2) zonas de transición entre las tendencias a largo plazo en las ta11as. En cierto grado, todas las bodegas mostraron cambios cíclicos en ta11a media y varianza durante la descarga. En la mitad de las bodegas, se observo que selección por ta11a por los descargadores indujo un cambio en la ta11a media. Como resultado de la estratificación, la secuencia de ta11as sacadas de todas las bodegas no fue aleatoria, sin considerar si una bodega contenía peces de un solo lance 0 de mas de uno. El numero de ta11as modales en una bodega no estaba relacionado al numero de lances. En una bodega adicional compuesta de peces de varios lances, un experimento de mezcla vertical indico que una muestra representativa del contenido podría estar limitada a la mitad inferior de la bodega. Se uso el contenido de las bodegas de prueba para generar 25 bodegas simuladas y comparar los resultados de tres métodos de muestreo aplicados a estas. Los métodos fueron: (1) muestreo aleatorio (usado también como norma), (2) muestreo extendido, en el cual el proceso de selección fue extendido sobre una porción grande de una bodega, y (3) medición consecutiva de peces durante la descarga de la bodega. EI muestreo repetido con cada método y distintas combinaciones de n y m indico que, puesto que la fuente principal de variación de ta11a ocurría entre las unidades primarias, aumentar n fue la manera mas eficaz de reducir las estimaciones de la varianza de las ta11as de los grupos de edad y el numero total de peces en los desembarcos. El muestreo extendido evito mayormente los efectos de la estratificación por ta11a, y su desempeño fue esencialmente comparable a aquel del muestreo aleatorio. Se recomienda muestrear con este método. El muestreo de peces consecutivos produjo estimaciones mas sesgadas con mayores varianzas. Un análisis de las muestras de frecuencia de ta11a de 1988 indico que, para los grupos de edad que aparecen con mayor frecuencia en la captura, una frecuencia de muestreo minima de una unidad primaria de cada seis para cada estrato de mes-área reduciría los coeficientes de variación (CV) de las estimaciones de ta11a correspondientes a aproximadamente 10% 0 menos. Una estratificación adicional de las muestras por tipo de lance, y no solamente mes-área, redujo aun mas los CV de los grupos de edad escasos, tales como los reclutas, y mejoró potencialmente su precisión. Los CV de las estimaciones del reclutamiento para las cohortes completamente pescadas durante 1981-1984 fueron alrededor de 3-8%. Las estimaciones del reclutamiento y sus varianzas fueron también relativamente insensibles a cambios en las capturas de trimestres individuales y las varianzas, respectivamente, de las cuales fueron derivadas. (PDF contains 70 pages)

Some statistical aspects of the long-term gill net monitoring programme for pike Esox lucius in Windermere (English Lake District)

For more than 55 years, data have been collected on the population of pike Esox lucius in Windermere, first by the Freshwater Biological Association (FBA) and, since 1989, by the Institute of Freshwater Ecology (IFE) of the NERC Centre for Ecology and Hydrology. The aim of this article is to explore some methodological and statistical issues associated with the precision of pike gill net catches and catch-per-unit-effort (CPUE) data, further to those examined by Bagenal (1972) and especially in the light of the current deployment within the Windermere long-term sampling programme. Specifically, consideration is given to the precision of catch estimates from gill netting, including the effects of sampling different locations, the effectiveness of sampling for distinguishing between years, and the effects of changing fishing effort.

Temporal and spatial variations in the quality of water in El Gergal Reservior, Seville, Spain

It is often difficult to define ‘water quality’ with any degree of precision. One approach is that suggested by Battarbee (1997) and is based on the extent to which individual lakes have changed compared with their natural ‘baseline’ status. Defining the base-line status of artificial lakes and reservoirs however, is, very difficult. In ecological terms, the definition of quality must include some consideration of their functional characteristics and the extent to which these characteristics are self-sustaining. The challenge of managing lakes in a sustainable way is particularly acute in semi-arid, Mediterranean countries. Here the quality of the water is strongly influenced by the unpredictability of the rainfall as well as year-to-year variations in the seasonal averages. Wise management requires profound knowledge of how these systems function. Thus a holistic approach must be adopted and the factors influencing the seasonal dynamics of the lakes quantified over a range of spatial and temporal scales. In this article, the authors describe some of the ways in which both long-term and short-term changes in the weather have influenced the seasonal and spatial dynamics of phytoplankton in El Gergal, a water supply reservoir situated in the south of Spain. The quality of the water stored in this reservoir is typically very good but surface blooms of algae commonly appear during warm, calm periods when the water level is low. El Gergal reservoir is managed by the Empresa Municipal de Abastecimiento y Saneamiento (EMASESA) and supplies water for domestic, commercial and industrial use to an area which includes the city of Seville and twelve of its surrounding towns (ca. 1.3 million inhabitants). El Gergal is the last of two reservoirs in a chain of four situated in the Rivera de Huelva basin, a tributary of the Guadalquivir river. It was commissioned by EMASESA in 1979 and since then the company has monitored its main limnological parameters on, at least, a monthly basis and used this information to improve the management of the reservoir. As a consequence of these intensive studies the physical, chemical and biological information acquired during this period makes the El Gergal database one of the most complete in Spain. In this article the authors focus on three ‘weather-related’ effects that have had a significant impact on the composition and distribution of phytoplankton in El Gergal: (i) the changes associated with severe droughts; (ii) the spatial variations produced by short-term changes in the weather; (iii) the impact of water transfers on the seasonal dynamics of the dinoflagellate Ceratium.

Comparison of habitat-based indices of abundance with fishery-independent biomass estimates from bottom trawl surveys

Rockfish species are notoriously difficult to sample with multispecies bottom trawl survey methods. Typically, biomass estimates have high coefficients of variation and can fluctuate outside the bounds of biological reality from year to year. This variation may be due in part to their patchy distribution related to very specific habitat preferences. We successfully modeled the distribution of five commercially important and abundant rockf ish species. A two-stage modeling method (modeling both presence-absence and abundance) and a collection of important habitat variables were used to predict bottom trawl survey catch per unit of effort. The resulting models explained between 22% and 66% of the variation in rockfish distribution. The models were largely driven by depth, local slope, bottom temperature, abundance of coral and sponge, and measures of water column productivity (i.e., phytoplankton and zooplankton). A year-effect in the models was back-transformed and used as an index of the time series of abundance. The abundance index trajectories of three of five species were similar to the existing estimates of their biomass. In the majority of cases the habitat-based indices exhibited less interannual variability and similar precision when compared with stratified survey-based biomass estimates. These indices may provide for stock assessment models a more stable alternative to current biomass estimates produced by the multispecies bottom trawl survey in the Gulf of Alaska.

Age validation, growth, mortality, and demographic modeling of spotted gully shark (Triakis megalopterus) from the southeast coast of South Africa

This study documents validation of vertebral band-pair formation in spotted gully shark (Triakis megalopterus) with the use of fluorochrome injection and tagging of captive and wild sharks over a 21-year period. Growth and mortality rates of T. megalopterus were also estimated and a demographic analysis of the species was conducted. Of the 23 OTC (oxytetracycline) -marked vertebrae examined (12 from captive and 11 from wild sharks), seven vertebrae (three from captive and four from wild sharks) exhibited chelation of the OTC and fluoresced under ultraviolet light. It was concluded that a single opaque and translucent band pair was deposited annually up to at least 25 years of age, the maximum age recorded. Reader precision was assessed by using an index of average percent error calculated at 5%. No significant differences were found between male and female growth patterns (P>0.05), and von Bertalanffy growth model parameters for combined sexes were estimated to be L∞=1711.07 mm TL, k=0.11/yr and t0=–2.43 yr (n=86). Natural mortality was estimated at 0.17/yr. Age at maturity was estimated at 11 years for males and 15 years for females. Results of the demographic analysis showed that the population, in the absence of fishing mortality, was stable and not significantly different from zero and particularly sensitive to overfishing. At the current age at first capture and natural mortality rate, the fishing mortality rate required to result in negative population growth was low at F>0.004/ yr. Elasticity analysis revealed that juvenile survival was the principal factor in explaining variability in population growth rate.

Using experiments and expert judgment to model catchability of Pacific rockfishes in trawl surveys, with application to bocaccio (Sebastes paucispinis) off British Columbia

The time series of abundance indices for many groundfish populations, as determined from trawl surveys, are often imprecise and short, causing stock assessment estimates of abundance to be imprecise. To improve precision, prior probability distributions (priors) have been developed for parameters in stock assessment models by using meta-analysis, expert judgment on catchability, and empirically based modeling. This article presents a synthetic approach for formulating priors for rockfish trawl survey catchability (qgross). A multivariate prior for qgross for different surveys is formulated by using 1) a correction factor for bias in estimating fish density between trawlable and untrawlable areas, 2) expert judgment on trawl net catchability, 3) observations from trawl survey experiments, and 4) data on the fraction of population biomass in each of the areas surveyed. The method is illustrated by using bocaccio (Sebastes paucipinis) in British Columbia. Results indicate that expert judgment can be updated markedly by observing the catch-rate ratio from different trawl gears in the same areas. The marginal priors for qgross are consistent with empirical estimates obtained by fitting a stock assessment model to the survey data under a noninformative prior for qgross. Despite high prior uncertainty (prior coefficients of variation ≥0.8) and high prior correlation between qgross, the prior for qgross still enhances the precision of key stock assessment quantities.

Use of stereo camera systems for assessment of rockfish abundance in untrawlable areas and for recording pollock behavior during midwater trawls

We describe the application of two types of stereo camera systems in fisheries research, including the design, calibration, analysis techniques, and precision of the data obtained with these systems. The first is a stereo video system deployed by using a quick-responding winch with a live feed to provide species- and size- composition data adequate to produce acoustically based biomass estimates of rockfish. This system was tested on the eastern Bering Sea slope where rockfish were measured. Rockfish sizes were similar to those sampled with a bottom trawl and the relative error in multiple measurements of the same rockfish in multiple still-frame images was small. Measurement errors of up to 5.5% were found on a calibration target of known size. The second system consisted of a pair of still-image digital cameras mounted inside a midwater trawl. Processing of the stereo images allowed fish length, fish orientation in relation to the camera platform, and relative distance of the fish to the trawl netting to be determined. The video system was useful for surveying fish in Alaska, but it could also be used broadly in other situations where it is difficult to obtain species-composition or size-composition information. Likewise, the still-image system could be used for fisheries research to obtain data on size, position, and orientation of fish.

Advances in methods for determining fecundity: application of the new methods to some marine fishes

Estimation of individual egg production (realized fecundity) is a key step either to understand the stock and recruit relationship or to carry out fisheries-independent assessment of spawning stock biomass using egg production methods. Many fish are highly fecund and their ovaries may weigh over a kilogram; therefore the work time can be consuming and require large quantities of toxic fixative. Recently it has been shown for Atlantic cod (Gadus morhua) that image analysis can automate fecundity determination using a power equation that links follicles per gram ovary to the mean vitellogenic follicular diameter (the autodiametric method). In this article we demonstrate the precision of the autodiametric method applied to a range of species with different spawning strategies during maturation and spawning. A new method using a solid displacement pipette to remove quantitative fecundity samples (25, 50, 100, and 200 milligram [mg]) is evaluated, as are the underlying assumptions to effectively fix and subsample the ovary. Finally, we demonstrate the interpretation of dispersed formaldehyde-fixed ovarian samples (whole mounts) to assess the presence of atretic and postovulatory follicles to replace labor intensive histology. These results can be used to estimate down regulation (production of atretic follicles) of fecundity during maturation.

Incorporating fishery observer data into an integrated catch-at-age and multiyear tagging model for estimating mortality rates and abundance

Tagging experiments are a useful tool in fisheries for estimating mortality rates and abundance of fish. Unfortunately, nonreporting of recovered tags is a common problem in commercial fisheries which, if unaccounted for, can render these estimates meaningless. Observers are often employed to monitor a portion of the catches as a means of estimating reporting rates. In our study, observer data were incorporated into an integrated model for multiyear tagging and catch data to provide joint estimates of mortality rates (natural and f ishing), abundance, and reporting rates. Simulations were used to explore model performance under a range of scenarios (e.g., different parameter values, parameter constraints, and numbers of release and recapture years). Overall, results indicated that all parameters can be estimated with reasonable accuracy, but that fishing mortality, reporting rates, and abundance can be estimated with much higher precision than natural mortality. An example of how the model can be applied to provide guidance on experimental design for a large-scale tagging study is presented. Such guidance can contribute to the successful and cost-effective management of tagging programs for commercial fisheries.

Measurements of total scattering spectra from bocaccio (Sebastes paucispinis)

Marine sportfishing in southern California is a huge industry with annual revenues totaling many billions of dollars. However, the stocks of lingcod and six rockfish species have been declared overfished by the Pacific Fisheries Management Council. As part of a multifaceted fisheries management plan, marine conservation areas, covering many million square nautical miles, have been mandated. To monitor the recovery of the rockfish stocks in these areas, scientists are faced with the following challenges: 1) multiple species of rockfish exist in these areas; 2) the species reside near or on the bottom at depths of 80 to 300 m; and 3) they are low in numerical density. To meet these challenges, multifrequency echosounders, multibeam sonar, and cameras mounted on remotely operated vehicles are frequently used (Reynolds et al., 2001). The accuracy and precision of these echosounder results are largely dependent upon the accuracy of the species classification and target strength estimation (MacLennan and Simmonds, 1992).

Behavioral ontogeny in larvae and early juveniles of the giant trevally (Caranx ignobilis) (Pisces: Carangidae)

Behavior of young (8−18 mm SL) giant trevally (Caranx ignobilis), a large coral-reef−associated predator, was observed in the laboratory and the ocean. Size was a better predictor of swimming speed and endurance than was age. Critical speed increased with size from 12 to 40 cm/s at 2.7 cm/s for each mm increase in size. Mean scaled critical speed was 19 body lengths/s and was not size related. Swimming speed in the ocean was 4 to 20 cm/s (about half of critical speed) and varied among areas, but within each area, it increased at 2 cm/s for each mm increase in size. Swimming endurance in the laboratory increased from 5 to 40 km at 5 km for each mm increase in size. Vertical distribution changed ontogenetically: larvae swam shallower, but more variably, and then deeper with growth. Two-thirds of individuals swam directionally with no ontogenetic increase in orientation precision. Larvae swam offshore off open coasts, but not in a bay. In situ observations of C. ignobilis feeding, interacting with pelagic animals, and reacting to reefs are reported. Manusc

The effect of autotrawl systems on the performance of a survey trawl

Three aspects of a survey bottom trawl performance—1) trawl geometry (i.e., net spread, door spread, and headrope height); 2) footrope distance off-bottom; and 3) bridle distance off-bottom—were compared among hauls by using either of two autotrawl systems (equal tension and net symmetry) and hauls conducted with towing cables of equal length and locked winches. The effects of environmental conditions, vessel heave, crabbing (i.e., the difference between vessel heading and actual vessel course over ground), and bottom current on trawl performance with three trawling modes were investigated. Means and standard deviations of trawl geometry measures were not significantly different between autotrawl and locked-winch systems. Bottom trawls performed better with either autotrawl system as compared to trawling with locked winches by reducing the variance and increasing the symmetry of the footrope contact with the bottom. The equal tension autotrawl system was most effective in counteracting effects of environmental conditions on footrope bottom contact. Footrope bottom contact was most inf luenced by environmental conditions during tows with locked winches. Both of the autotrawl systems also reduced the variance and increased the symmetry of bridle bottom contact. Autotrawl systems proved to be effective in decreasing the effects of environmental factors on some aspects of trawl performance and, as a result, have the potential to reduce among-haul variance in catchability of survey trawls. Therefore, by incorporating an autotrawl system into standard survey procedures, precision of survey estimates of relative abundance

Using poststratification to improve abundance estimates from multispecies surveys: a study of juvenile flatfishes

Population assessments seldom incorporate habitat information or use previously observed distributions of fish density. Because habitat affects the spatial distribution of fish density and overall abundance, the use of habitat information and previous estimates of fish density can produce more precise and less biased population estimates. In this study, we describe how poststratification can be applied as an unbiased estimator to data sets that were collected under a probability sampling design, typical of many multispecies trawl surveys. With data from a multispecies survey of juvenile flatfish, we show how poststratification can be applied to a data set that was not collected under a probability sampling design, where both the precision and the bias are unknown. For each of four species, three estimates of total abundance were compared: 1) unstratified; 2) poststratified by habitat; and 3) poststratified by habitat and fish density (high fish density and low fish density) in nearby years. Poststratification by habitat gave more precise and (or) less design-biased estimates than an unstratified estimator for all species in all years. Poststratification by habitat and fish density produced the most precise and representative estimates when the sample size in the high fish-density and low fish-density strata were sufficient (in this study, n≥20 in the high fish-density stratum, n≥9 in the low fish-density stratum). Because of the complexities of statistically testing the annual stratified data, we compared three indices of abundance for determining statistically significant changes in annual abundance. Each of the indices closely approximated the annual differences of the poststratified estimates. Selection of the most appropriate index was dependent upon the species’ density distribution within habitat and the sample size in the different habitat areas. The methods used in this study are particularly useful for estimating individual species abundance from multispecies surveys and for retrospective st

Effect of type of otolith and preparation technique on age estimation of larval and juvenile spot (Leiostomus xanthurus)

Otoliths of larval and juvenile fish provide a record of age, size, growth, and development (Campana and Neilson, 1985; Thorrold and Hare, 2002). However, determining the time of first increment formation in otoliths (Campana, 2001) and assessing the accuracy (deviation from real age) and precision (repeatability of increment counts from the same otolith) of increment counts are prerequisites for using otoliths to study the life history of fish (Campana and Moksness, 1991). For most fish species, first increment deposition occurs either at hatching, a day after hatching, or after first feeding and yolksac absorption (Jones, 1986; Thorrold and Hare, 2002). Increment deposition before hatching also occurs (Barkmann and Beck, 1976; Radtke and Dean, 1982). If first increment deposition does not occur at hatching, the standard procedure is to add a predetermined number to increment counts to estimate fish age (Campana and Neilson, 1985).

Relative pleopod length as an indicator of size at sexual maturity in slipper (Scyllarides squammosus) and spiny Hawaiian (Panulirus marginatus) lobsters

Body size at gonadal maturity is described for females of the slipper lobster (Scyllarides squammosus) (Scyllaridae) and the endemic Hawaiian spiny lobster (Panulirus marginatus) (Palinuridae) based on microscopic examination of histological preparations of ovaries. These data are used to validate several morphological metrics (relative exopodite length, ovigerous condition) of functional sexual maturity. Relative exopodite length (“pleopod length”) produced consistent estimates of size at maturity when evaluated with a newly derived statistical application for estimating size at the morphometric maturation point (MMP) for the population, identified as the midpoint of a sigmoid function spanning the estimated boundaries of overlap between the largest immature and smallest adult animals. Estimates of the MMP were related to matched (same-year) characterizations of sexual maturity based on ovigerous condition — a more conventional measure of functional maturity previously used to characterize maturity for the two lobster species. Both measures of functional maturity were similar for the respective species and were within 5% and 2% of one another for slipper and spiny lobster, respectively. The precision observed for two shipboard collection series of pleopod-length data indicated that the method is reliable and not dependent on specialized expertise. Precision of maturity estimates for S. squammosus with the pleopod-length metric was similar to that for P. marginatus with any of the other measures (including conventional evidence of ovigerous condition) and greatly exceeded the precision of estimates for S. squammosus based on ovigerous condition alone. The two measures of functional maturity averaged within 8% of the estimated size at gonadal maturity for the respective species. Appendage-to-body size proportions, such as the pleopod length metric, hold great promise, particularly for species of slipper lobsters like S. squammosus for which there exist no other reliable conventional morphological measures of sexual maturity. Morphometric proportions also should be included among the factors evaluated when assessing size at sexual maturity in spiny lobster stocks; previously, these proportions have been obtained routinely only for brachyuran crabs within the Crustacea.