Statistical inference about the relative efficiency of a new survey protocol, based on paired-tow survey calibration data

Paired-tow calibration studies provide information on changes in survey catchability that may occur because of some necessary change in protocols (e.g., change in vessel or vessel gear) in a fish stock survey. This information is important to ensure the continuity of annual time-series of survey indices of stock size that provide the basis for fish stock assessments. There are several statistical models used to analyze the paired-catch data from calibration studies. Our main contributions are results from simulation experiments designed to measure the accuracy of statistical inferences derived from some of these models. Our results show that a model commonly used to analyze calibration data can provide unreliable statistical results when there is between-tow spatial variation in the stock densities at each paired-tow site. However, a generalized linear mixed-effects model gave very reliable results over a wide range of spatial variations in densities and we recommend it for the analysis of paired-tow survey calibration data. This conclusion also applies if there is between-tow variation in catchability.

Perspectives from agriculture: advances in livestock breeding - implications for aquaculture genetics

In this paper we present livestock breeding developments that could be taken into consideration in the genetic improvement of farmed aquaculture species, especially in freshwater fish. Firstly, the current breeding objective in aquatic species has focused almost exclusively on the improvement of body weight at harvest or on growth related traits. This is unlikely to be sufficient to meet the future needs of the aquaculture industry. To meet future demands breeding programs will most likely have to include additional traits, such as fitness related ones (survival, disease resistance), feed efficiency, or flesh quality, rather than only growth performance. In order to select for a multi-trait breeding objective, genetic variation in traits of interest and the genetic relationships among them need to be estimated. In addition, economic values for these traits will be required. Generally, there is a paucity of data on variable and fixed production costs in aquaculture, and this could be a major constraint in the further expansion of the breeding objectives. Secondly, genetic evaluation systems using the restricted maximum likelihood method (REML) and best linear unbiased prediction (BLUP) in a framework of mixed model methodology could be widely adopted to replace the more commonly used method of mass selection based on phenotypic performance. The BLUP method increases the accuracy of selection and also allows the management of inbreeding and estimation of genetic trends. BLUP is an improvement over the classic selection index approach, which was used in the success story of the genetically improved farmed tilapia (GIFT) in the Philippines, with genetic gains from 10 to 20 per cent per generation of selection. In parallel with BLUP, optimal genetic contribution theory can be applied to maximize genetic gain while constraining inbreeding in the long run in selection programs. Thirdly, by using advanced statistical methods, genetic selection can be carried out not only at the nucleus level but also in lower tiers of the pyramid breeding structure. Large scale across population genetic evaluation through genetic connectedness using cryopreserved sperm enables the comparison and ranking of genetic merit of all animals across populations, countries or years, and thus the genetically superior brood stock can be identified and widely used and exchanged to increase the rate of genetic progress in the population as a whole. It is concluded that sound genetic programs need to be established for aquaculture species. In addition to being very effective, fully pedigreed breeding programs would also enable the exploration of possibilities of integrating molecular markers (e.g., genetic tagging using DNA fingerprinting, marker (gene) assisted selection) and reproductive technologies such as in-vitro fertilization using cryopreserved spermatozoa.

Identification and estimation of size from the beaks of 18 species of cephalopods from the Pacific Ocean

A method of identifying the beaks and estimating body weight and mantle length of 18 species of cephalopods from the Pacific Ocean is presented. Twenty specimens were selected from each of the following cephalopod species: Symplectoteuthis oualaniensis, Dosidicus gigas, Ommastrephes bartramii, S. luminosa, Todarodes pacificus, Nototodarus hawaiiensis, Ornithoteuthis volalilis, Hyaloteuthis pelagica, Onychoteuthis banksii, Pterygioteuthis giardi, Abraliopsis affinis, A. felis, Liocranchia reinhardti, Leachia danae, Histioteuthis heteropsis, H. dofleini, Gonalus onyx, and Loligo opalescens. Dimensions measured on the upper and lower beak are converted to ratios and compared individually among the species using an analysis of variance procedure with Tukey's omega and Duncan's multiple range tests. Significant differences (P =0.05) observed among the species' beak ratio means and structural characteristics are used to construct artificial keys for the upper and lower beaks of the 18 species. Upper and lower beak dimensions are used as independent variables in a linear regression model with mantle length and body weight (log transformed). (PDF file contains 56 pages.)

Population dynamics of the skipjack tuna (Katsuwonus pelamis) of the Eastern Pacific Ocean

ENGLISH: Catches of skipjack tuna supporting major fisheries in parts of the western, central and eastern Pacific Ocean have increased in recent years; thus, it is important to examine the dynamics of the fishery to determine man's effect on the abundance of the stocks. A general linear hypothesis model was developed to standardize fishing effort to a single vessel size and gear type. Standardized effort was then used to compute an index of abundance which accounts for seasonal variability in the fishing area. The indices of abundance were highly variable from year to year in both the northern and southern areas of the fishery but indicated a generally higher abundance in the south. Data from 438 fish tagged and recovered in the eastern Pacific Ocean were used to compute growth curves. A least-squares technique was used to estimate the parameters of the von Bertalanffy growth function. Two estimates of the parameters were made by analyzing the same data in different ways. For the first set of estimates, K= 0.819 on an annual instantaneous basis and L= 729 mm; for the second, K = 0.431 and L=881. These compared well with estimates derived using the Chapman-Richards growth function, which includes the von Bertalanffy function as a special case. It was concluded that the latter function provided an adequate empirical fit to the skipjack data since the more complicated function did not significantly improve the fit. Tagging data from three cruises involving 8852 releases and 1777 returns were used to compute mortality rates during the time the fish were in the fishery. Two models were used in the analyses. The best estimates of the catchability coefficient (q) in the north and south were 8.4 X 10- 4 and 5.0 X 10- 5 respectively. The other loss rate (X), which included losses due to emigration, natural mortality and mortality due to carrying a tag, was 0.14 on an annual instantaneous basis for both areas. To detect the possible effect of fishing on abundance and total yield, the relation between abundance and effort and between total catch and effort was examined. It was found that at levels of intensity observed in the fishery, fishing does not appear to have had any measurable effect on the stocks. It was concluded therefore that the total catch could probably be increased by substantially increasing total effort beyond the present level, and that the fluctuations in abundance are fishery-independent. The estimates of growth, mortality and fishing effort were used to compute yield-per-recruitment isopleths for skipjack in both the northern and southern areas. For a size at first entry of about 425 mm, the yield per recruitment was calculated at 3 pounds in the north and 1.5 pounds in the south. In both areas it would be possible to increase the yield per recruitment by increasing fishing effort. It was not possible to assess potential production of the skipjack stocks fished in the eastern Pacific, except to note that the fishery had not affected their abundance and that they were certainly under-exploited. It was concluded that the northern and southern stocks could support increased harvests, especially the latter. SPANISH: Las capturas de atún barrilete que sostienen las pesquerías principales de la parte occidental, central y oriental del Océano Pacífico han aumentado en los últimos años; así que es importante examinar la dinámica de la pesquería para determinar el efecto que pueda tener sobre la abundancia de los stocks. Se desarrolló un modelo hipotético, lineal para standardizar el esfuerzo de pesca a un solo tamaño de barco y tipo de arte. Luego se usó el esfuerzo standardizado para computar un índice de la abundancia que pueda dar razón de la variabilidad estacional en el área de pesca. Los índices de la abundancia variaron mucho de un año a otro tanto en el área septentrional como en el área meridional de la pesquería, pero indicaron una abundancia generalmente superior en el sur. Se emplearon los datos de 438 peces marcados y recuperados en el Océano Pacífico oriental para computar las curvas de crecimiento. Una técnica de mínimos cuadrados fue usada para estimar los parámetros de la función de crecimiento de van Bertalanffy. Se hicieron dos estimativos de los parámetros mediante el análisis de los mismos datos, de diferente manera. Para el primer juego de estimativos, K=0.819 sobre una base anual instantánea y L∞=729 mm; para el segundo, K=0.431 y L∞=881. Estos se correlacionaron bien con los estimativos obtenidos usando la función de crecimiento de Chapman-Richards, que incluye la de von Bertalanffy como un caso especial. Se decidió que la última función proveía un ajuste empírico, adecuado a los datos del barrilete, ya que la función más complicada no mejoró significativamente el ajuste. Los datos de marcación de tres cruceros incluyendo 8852 liberaciones y 1777 retornos, fueron usados para computar las tasas de mortalidad durante el tiempo en que los peces estuvieron en la pesquería. Se usaron dos modelos en los análisis. Los mejores estimativos del coeficiente de capturabilidad (q) en el norte y en el sur fueron 8.4 X 10-4 y 5.0 X 10-5 , respectivamente. La otra tasa de pérdida (X), la cual incluyó pérdidas debidas a la emigración, mortalidad natural y mortalidad debida a llevar una marca, fue 0.14 sobre una base anual instantánea para las dos áreas. Con el fin de descubrir el efecto que posiblemente pueda tener la pesca sobre la abundancia y el rendimiento total, se examinó la relación entre la abundancia y el esfuerzo y entre la captura total y el esfuerzo. Se encontró que a los niveles de la intensidad observada en la pesquería, la pesca no parece haber tenido ningún efecto perceptible en los stocks. Por lo tanto se decidió que mediante un aumento substancial del esfuerzo total, más allá del nivel actual, la captura total probablemente podría aumentarse, y que las fluctuaciones de la abundancia son independientes de la pesquería. Los estimativos del crecimiento, mortalidad y esfuerzo de pesca fueron usados para computar las isopletas del rendimiento por recluta del barrilete, tanto en las áreas del norte como del sur. Para una talla de primera entrada de unos 425 mm, el rendimiento por recluta fue calculado en 3 libras en el norte y 1.5 libras en el sur. En ambas áreas sería posible aumentar el rendimiento por recluta mediante un aumento del esfuerzo de pesca. No fue posible determinar la producción potencial de los stocks del barrilete pescado en el Pacífico oriental, excepto para observar que la pesquería no ha afectado su abundancia y que ciertamente se encuentran subexplotados. Se concluyó que los stocks norte y sur pueden soportar un aumento en el rendimiento, especialmente este último. (PDF contains 274 pages.)

Semiempirical down-scaling of GCM output to the local scale for temperature, precipitation, and runoff

EXTRACT (SEE PDF FOR FULL ABSTRACT): An empirically derived multiple linear regression model is used to relate a local-scale dependent variable (either temperature, precipitation, or surface runoff) measured at individual gauging stations to six large-scale independent variables (temperature, precipitation, surface runoff, height to the 500-mbar pressure surface, and the zonal and meridional gradient across this surface). ...The area investigated is the western United States. ... The calibration data set is from 1948 through 1988 and includes data from 268 joint temperature and precipitation stations, 152 streamflow stations (which are converted to runoff data), and 24 gridded 500-mbar pressure height nodes.