Comparison of two approaches for estimating natural mortality based on longevity

Vetter (1988) noted that her review of the estimation of the instantaneous natural mortality rate (M) was initiated by a discussion among colleagues that identified M as the single most impor ta nt but least well-estimated parameter in fishery models. A lthough much has been accomplished in the inter vening years, M remains one of the most difficult parameters to estimate in fishery stock assessments. A number of novel approaches using tagging and telemetry data provide promise for making reliable direct estimates of M for a given stock (Hearn et al., 1998 ; Frusher and Hoenig, 2001; Hightower et al., 2001; Latour et al., 2003; Pollock et al., 2004). However, such methods are often impracticable and fishery scientists must approximate M by using estimates made for other stocks of the same or similar species or by predicting M from features of the species’ life history (Beverton and Holt, 1959; Beverton, 1963; Alverson and Carney, 1975; Pauly, 1980; Hoenig, 1983; Peterson and Wroblewski, 1984; Roff, 1984; Gunderson and Dygert, 1988; Chen and Watanabe, 1989; Charnov, 1993; Jensen, 1996; Lorenzen, 1996).

Seasonal growth of King George whiting (Sillaginodes punctata) estimated from length-at-age samples of the legal-size harvest

Samples of 11,000 King George whiting (Sillaginodes punctata) from the South Australian commercial and recreational catch, supplemented by research samples, were aged from otoliths. Samples were analyzed from three coastal regions and by sex. Most sampling was undertaken at fish processing plants, from which only fish longer than the legal minimum length were obtained. A left-truncated normal distribution of lengths at monthly age was therefore employed as model likelihood. Mean length-at-monthly-age was described by a generalized von Bertalanffy formula with sinusoidal seasonality. Likelihood standard deviation was modeled to vary allometrically with mean length. A range of related formulas (with 6 to 8 parameters) for seasonal mean length at age were compared. In addition to likelihood ratio tests of relative fit, model selection criteria were a minimum occurrence of high uncertainties (>20% SE), of high correlations (>0.9, >0.95, and >0.99) and of parameter estimates at their biological limits, and we sought a model with a minimum number of parameters. A generalized von Bertalanffy formula with t0 fixed at 0 was chosen. The truncated likelihood alleviated the overestimation bias of mean length at age that would otherwise accrue from catch samples being restricted to legal sizes.

An egg production method for estimating spawning biomass of pelagic fish: Application to the northern anchovy, Engraulis mordax

Fishery scientists engaged in estimating the size of free-swimming populations have never had a technique available to them whereby all the parameters could be estimated from a resource survey and where no parameter values need to be assumed. Recognizing the need for a technique of this kind, the staff of the Coastal Fisheries Resources Division of the Southwest Fisheries Center (SWFC) devised an egg production method for anchovy biomass assessment. Previously, anchovy biomass was estimated by approximate methods derived from a long-time series and anchovy larval abundance, which required about 5 ma of shiptime each year to integrate the area under a seasonal spawning curve. One major assumption used in the larval abundance census method is that there is constant proportionality between larval numbers and spawning biomass. This has now proved to be erroneous. (PDF file contains 105 pages.)

Veränderungen sensorischer und instrumenteller Parameter während der Gefrierlagerung von Ostseedorsch (Gadus morhua) bei unterschiedlichen Temperaturen

Changes in sensory and instrumental quality parameter sand in thawing drip, cooking drip and total drip loss of frozen stored Baltic cod fillets (Gadus morhua) at different storage temperatures were investigated. Cod fillets stored at –20 °C and –30 °C exhibited the lowest drip losses and obtained the highest sensory scores. Drip losses were found to be highest in cod fillets stored at –10°C and in double frozen fillets stored at –20 °C. These two experiments also gave the lowest sensory scores. The texture parameters increased during storage parallel with storage time. The waterbinding capacity was lowest at –10 °C and almost constant at –30 °C. There is a good correlation between the sensory scores for “tough” and the instrumental texture measurement for hardness and chewiness.

Growth of skipjack, Katsuwonus pelamis, and yellowfin, Thunnus albacares, tunas in the Eastern Pacific Ocean, as estimated from tagging data

ENGLISH: The average linear growth rate of skipjack in the eastern Pacific is less than 1 mm per day except for fish 375 to 424 mm in length at release. The growth rate shows a decrease with increasing length and increasing time at liberty. The growth rate of fish in the length range of about 43 to 57 cm is apparently more rapid in the eastern Pacific than in the western Pacific. Dsing data for the northeastern and southeastern Pacific combined, K and ~ were estimated to be 0.658 (on an annual basis) and 885 mm, respectively, by the ungrouped method and 0.829 and 846 mm, respectively, by the grouped method. Sensitivity analyses have shown however, that the estimates of these parameters are poorly determined by the sum of squares method used to derive them. Estimates of K and ~ for the eastern Pacific tend to be lower and higher, respectively, than those for the western Pacific. The average linear growth rate of yellowfin in the eastern Pacific is a little less than 1 mm per day for fish between about 25 and 100 cm in length at release. The growth appears to be most rapid in Area 2 (Revillagigedo Islands) and slowest in Areas 1 (Baja California), 5 (Central America- Colombia), and 6 (Ecuador-Peru). There is considerable variation in the growth rates of individual fish. The growth does not show a decrease with increasing length or increasing time at liberty so realistic estimates of the parameters of the von Bertalanffy or other similar equations cannot be calculated from these data. If realistic estimates of these parameters are to be secured larger fish must be tagged and released or many more long-term returns from fish to about 100 cm in length at release must be obtained. The growth patterns for the eastern Pacific, central Pacific and eastern Atlantic found by most other investigators differ from one another and from those found in the present study. Some of these differences may be real and others may be due to deficiencies in the data or the methods of analysis. Estimates obtained from tagging data are believed to be realistic provided the tags do not inhibit the growth of the fish. It appears that the growth rates of single- and double-tagged fish are the same; this indicates, though not unequivocally, that the tags do not inhibit the growth. SPANISH: La tasa media de crecimiento lineal del barrilete en el Pacífico oriental es inferior a lmm/día, excepto en el caso de peces de entre 375y 424mm de longitud de liberación. La tasa de crecimiento disminuye a medida que aumenta la longitud y el tiempo en libertad. La tasa de crecimiento de peces de entre unos 43 y 57 cm de longitud parece ser mayor en el Pacífico oriental que en el occidental. A partir de datos del Pacífico nororiental y suroriental combinados, se estimaron K y loo en 0.658 (anual) y 885mm, respectivamente, usando el método no agrupado, y 0.829 y 846mm, respectivamente, usando el método agrupado. Sin embargo, los análisis de sensitividad han demostrado que el método de suma de cuadrados utilizado para derivar las estimaciones de estos parámetros las determina con poca precisión. Las estimaciones de K y loo para el Pacífico oriental suelen ser inferiores y superiores, respectivamente, a los del Pacífico occidental. La tasa media de crecimiento lineal del aleta amarilla en el Pacífico oriental es ligeramente inferior a lmm/día para los peces de entre unos 25y 100cmde longitud de liberación. El crecimiento parece ser más rápido en el Area 2(Islas Revillagigedo),y más lento en las Areas 1(Baja California), 5 (Centroamérica-Colombia), y 6 (Ecuador-Perú). Las tasas de crecimiento de peces individuales varían considerablemente. El crecimiento no muestra una disminuciónconun aumento en la longitud o en el tiempo en libertad, y por consecuencia no se se pueden calcular estimaciones realistas de los parámetros de la ecuación de von Bertalanffy u otras ecuaciones similares a partir de estos datos. Para obtener estimaciones realistas de estos parámetros sería necesario marcar peces mayores u obtener muchas más devoluciones a largo plazo de marcas de peces de unos 100cm de longitud de liberación. Los patrones de crecimiento correspondientes al Pacífico oriental, Pacífico central, y Atlántico oriental descubiertos por la mayoría de los investigadores son diferentes entre síy también de los del presente estudio. Es posibleque algunas de estas diferencias sean verdaderas, mientras que otras se deban a faltas en los datos on en los métodos analíticos utilizados. Se considera que las estimaciones obtenidas a partir de los datos de marcado son realistas, suponiendo siempre que las marcas no impidan el crecimiento de los peces. Parece ser que las tasas de crecimiento de peces con una marca y con dos son idénticas, lo cual indica, aunque sin certeza total, que las marcas no ejercen tal efecto. (PDF contains 76 pages.)

Growth and age composition of northern bluefin tuna, Thunnus thynnus, caught in the Eastern Pacific Ocean, as estimated from length-frequency data, with comments on trans-Pacific migrations

ENGLISH: The spawning of Pacific northern bluefin tuna, Thunnus thynnus, takes place only in the western Pacific Ocean (WPO), but substantial numbers of the juveniles migrate to the eastern Pacific Ocean (EPO), where they remain for several months, or longer, and the.n return to the WPO. Lengthfrequency and tagging data show that many bluefin arrive in the EPO as 1-and 2-year olds, and remain there for one or two fishing seasons before returning to the WPO. The proportion of the fish which make the west-to-east migration varies among years. The numbers of 1-, 2-, 3-, 4, and >4 –year olds in the catches of the EPO are estimated for most years of the 1952-1991 period. SPANISH: EI desove del atun aleta azul del norte del Pacifico, Thunnus thynnus, ocurre solamente en el Océano Pacifico occidental (WPO), pero números substanciales de los juveniles migran al Océano Pacifico oriental (OPO), donde permanecen unos meses, 0 mas, antes de regresar al WPO. Datos de marcado y frecuencia de talla indican que muchos aletas azules llegan al OPO a 1 o 2 anos de edad, y permanecen alIi una 0 dos temporadas de pesca antes de regresar al WPO. La proporcion de los peces que migra del oeste al este varia entre anos. Se estima el numero de peces de 1, 2, 3, 4, Y>4 anos de edad en las capturas del OPO para la mayoria de los anos del periodo de 1952-1991. (PDF contains 40 pages.)

Determination and modelling of water quality parameter in Abuja metropolis

This research work involves the determination and modelling of water parameter such as pH, temperature, turbidity, chloride, hardness. The result of the analysis was used as important operating variables to generate a model equation of pH, hardness, temperature, turbidity and chloride. The values obtained from the model equation were compared with those from experiment. On an average bases the values were close. These parameters can be used to monitor the extent of pollution of pond water and to monitor stress and diseases of fish. The experimental data of pH was in the range of 6.7 to 6.9 while the modelled result was also between 6.7 to 7.0. The turbidity experimental value was close to the modelled value also. The chloride value for the experimental data was in the range of 25.32 to 35.0. The total hardness value ranges between 4.5 to 65.1 mg/l while the modelled result ranges between 11.025 to 68.402 mg/l. The result was within the acceptable limit of world health organization standard on water quality parameter.

Reference-growth rate – a simple and handy parameter summarizing the influence of environmental conditions

Abstract Environmental changes may have an impact on life conditions of the fish, e.g. food supply for the fish. The prevailing environmental conditions apply evenly to all age groups of one stock. Small fish have high growth rates, whereas large fish grow with low rates. But, it can be shown on the basis of the von Bertalanffy-growth model that it is sufficient to know only the growth rate of one single age group to compute the growth rates of all other age groups. The growth rate of a reference fish GRF (e.g. a fish with a body mass of 1 kg) was introduced as a reference growth describing the current food condition of all age groups of the stock. As an example a time series of the reference-growth rate of the northern cod stock (NAFO, 3K) was computed for the time span 1979 to 1999. For the northern cod stock it can be observed that environmental conditions caused growth rates below the long-term mean for seven years in a row. After a prolonged hunger period the fish stock collapsed in 1992 also by the impact of fisheries - and this was probably not a coincidence. Now, with the reference-growth rate GRF a simple and handy parameter was found to summarize the influence of the environmental conditions on growth and other derived models and therefore makes it easier to compute the influence of environmental changes within stock assessment. Zusammenfassung Veränderungen der Umwelt können Auswirkungen auf die Lebensbedingungen der Fische haben, z. B. auf das Nahrungsangebot der Fische. Die vorherrschenden Umgebungsbedingungen wirken gleichmäßig auf alle Altersgruppen eines Bestandes, wobei typischer Weise kleineFische hohe Wachstumsraten haben, während die großen Fische mit niedrigen Raten wachsen. Auf der Grundlage des von Bertalanffy-Wachstumsmodells kann gezeigt werden, dass es ausreicht, nur die Wachstumsrate von einer einzigen Altersgruppe zu kennen, um die Wachstumsraten von allen anderen Altersgruppen berechnen zu können. Die Wachstumsrate eines Referenz-Fisches (z.B. eines Fisches mit einer Körpermasse von 1 kg) wurde als Referenz-Wachstum GRF eingeführt, die den aktuellen Zustand des Nahrungsangebots füralle Altersgruppen des Bestandes beschreibt. Als Beispiel wurde einer Zeitreihe der Referenz-Wachstumsraten des nördlichen Kabeljaubestandes (NAFO, 3K) für die Zeitsraum 1979 bis 1999 berechnet. Für diesen Kabeljaubestand war zu beobachten, dass Umgebungsbedingungen für sieben Jahre in Folge Wachstumsraten unter dem langjährigen Mittelwert verursachten. Nach einer längeren Hungerperiode kollabierte dieser Fischbestand im Jahr 1992 auch durch den Einfluß der Fischerei - und dies war sicher kein Zufall. Jetzt, mit der Referenz-Wachstumsrate GRF, ist ein einfacher und handlicher Parameter gefunden, der es gestattet den Einfluss der Umweltbedingungen auf die Wachstumsbedingungen und andere davon abgeleitete Modelle zusammenzufassen. Dies macht es einfach, den Einfluss von Umweltveränderungen innerhalb der Bestandsabschätzungen zu berechnen.

Growth of Skipjack tuna (Katsuwonus pelamis) in the eastern Pacific Ocean as estimated from tagging data

English: Data obtained from tagging experiments initiated during 1953-1958 and 1969-1981 for skipjack tuna from the coastal eastern Pacific Ocean (EPO) are reanalyzed, using the Schnute generalized growth model. The objective is to provide information that can be used to generate a growth transition matrix for use in a length-structured population dynamics model. The analysis includes statistical approaches to include individual variability in growth as a function of length at release and time at liberty, measurement error, and transcription error. The tagging data are divided into northern and southern regions, and the results suggest that growth rates differ between the two regions. The Schnute model provides a significantly better fit to the data than the von Bertalanffy model, a sub-model of the Schnute model, for the northern region, but not for the southern region. Individual variation in growth is best described as a function of time at liberty and as a function of growth increment for the northern and southern regions, respectively. Measurement error is a significant part of the total variation, but the results suggest that there is no bias caused by the measurement error. Additional information, particularly for small and large fish, is needed to produce an adequate growth transition matrix that can be used in a length-structured population dynamics model for skipjack tuna in the EPO. Spanish: Los datos obtenidos de los experimentos de marcado iniciados durante los períodos de 1953- 1958 y de 1969-1981 para el atún barrilete en las costas del Océano Pacífico Oriental (OPO) fueron analizados nuevamente, utilizando el modelo de crecimiento generalizado de Schnute. El objetivo es brindar información que sea útil para producir una matriz sobre la tran-sición de crecimiento que pueda utilizarse en un modelo de dinámica poblacional estructurado por talla. El análisis usa enfoques estadísticos para poder incluir la variabilidad individual del crecimiento como función de la talla de liberación y tiempo en libertad, el error de medición, y el error de transcripción. Los datos de marcado son divididos en regiones norte y sur, y los resultados sugieren que las tasas de crecimiento en las dos regiones son diferentes. En la región norte, pero no en la región sur, el modelo de Schnute se ajusta significativamente mejor a los datos que el modelo von Bertalanffy, un sub-modelo del modelo de Schnute. La mejor descripción de la variación individual en el crecimiento es como una función del tiempo en libertad y como una función del incremento de crecimiento para las regiones norte y sur, respectivamente. El error de medición es una parte significativa de la variación total, pero los resultados sugieren que no existe un sesgo causado por el error de medición. Se necesita información adicional, particularmente para peces pequeños y grandes, para poder producir una matriz de transición de crecimiento adecuada que pueda utilizarse en el modelo de dinámica poblacional estructurado por tallas para el atún barrilete en el OPO.

Distributions and abundances of Pacific sardine (Sardinops sagax) and other pelagic fishes in the California Current Ecosystem during spring 2006, 2008, and 2010, estimated from acoustic–trawl surveys

The abundances and distributions of coastal pelagic fish species in the California Current Ecosystem from San Diego to southern Vancouver Island, were estimated from combined acoustic and trawl surveys conducted in the spring of 2006, 2008, and 2010. Pacific sardine (Sardinops sagax), jack mackerel (Trachurus symmetricus), and Pacific mackerel (Scomber japonicus) were the dominant coastal pelagic fish species, in that order. Northern anchovy (Engraulis mordax) and Pacific herring (Clupea pallasii) were sampled only sporadically and therefore estimates for these species were unreliable. The estimates of sardine biomass compared well with those of the annual assessments and confirmed a declining trajectory of the “northern stock” since 2006. During the sampling period, the biomass of jack mackerel was stable or increasing, and that of Pacific mackerel was low and variable. The uncertainties in these estimates are mostly the result of spatial patchiness which increased from sardine to mackerels to anchovy and herring. Future surveys of coastal pelagic fish species in the California Current Ecosystem should benefit from adaptive sampling based on modeled habitat; increased echosounder and trawl sampling, particularly for the most patchy and nearshore species; and directed-trawl sampling for improved species identification and estimations of their acoustic target stren

Age and growth of spiny dogfish (Squalus acanthias) in the Gulf of Alaska: analysis of alternative growth models

Ten growth models were fitted to age and growth data for spiny dogfish (Squalus acanthias) in the Gulf of Alaska. Previous studies of spiny dogfish growth have all fitted the t0 formulation of the von Bertalanffy model without examination of alternative models. Among the alternatives, we present a new two-phase von Bertalanffy growth model formulation with a logistically scaled k parameter and which estimates L0. A total of 1602 dogfish were aged from opportunistic collections with longline, rod and reel, set net, and trawling gear in the eastern and central Gulf of Alaska between 2004 and 2007. Ages were estimated from the median band count of three independent readings of the second dorsal spine plus the estimated number of worn bands for worn spines. Owing to a lack of small dogfish in the samples, lengths at age of small individuals were back-calculated from a subsample of 153 dogfish with unworn spines. The von Bertalanffy, two-parameter von Bertalanffy, two-phase von Bertalanffy, Gompertz, two-parameter Gompertz, and logistic models were fitted to length-at-age data for each sex separately, both with and without back-calculated lengths at age. The two-phase von Bertalanffy growth model produced the statistically best fit for both sexes of Gulf of Alaska spiny dogfish, resulting in L∞ = 87.2 and 102.5 cm and k= 0.106 and 0.058 for males and females, respectively.

Defining plausible migration rates based on historical tagging data: a Bayesian mark-recapture model applied to English sole (Parophrys vetulus)

A generalized Bayesian population dynamics model was developed for analysis of historical mark-recapture studies. The Bayesian approach builds upon existing maximum likelihood methods and is useful when substantial uncertainties exist in the data or little information is available about auxiliary parameters such as tag loss and reporting rates. Movement rates are obtained through Markov-chain Monte-Carlo (MCMC) simulation, which are suitable for use as input in subsequent stock assessment analysis. The mark-recapture model was applied to English sole (Parophrys vetulus) off the west coast of the United States and Canada and migration rates were estimated to be 2% per month to the north and 4% per month to the south. These posterior parameter distributions and the Bayesian framework for comparing hypotheses can guide fishery scientists in structuring the spatial and temporal complexity of future analyses of this kind. This approach could be easily generalized for application to other species and more data-rich fishery analyses.

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.

Interdecadal change in growth of sablefish (Anoplopoma fimbria) in the northeast Pacific Ocean

Errors in growth estimates can affect drastically the spawner-perrecruit threshold used to recommend quotas for commercial fish catches. Growth parameters for sablefish (Anoplopoma fimbria) in Alaska have not been updated for stock assessment purposes for more than 20 years, although aging of sablefish has continued. In this study, length-stratified data (1981–93 data from the annual longline survey conducted cooperatively by the Fisheries Agency of Japan and the Alaska Fisheries Science Center of the National Marine Fisheries Service) were updated and corrected for discovered sampling bias. In addition, more recent, randomly collected samples (1996–2004 data from the annual longline survey conducted by the Alaska Fisheries Science Center) were analyzed and new length-at-age and weight-at-age parameters were estimated. Results were similar between this analysis with length-at-age data from 1981 to 2004 and analysis with updated longline survey data through 2010; therefore, we used our initial results from analysis done with data through 2004. We found that, because of a stratified sampling scheme, growth estimates of sablefish were overestimated with the older data (1981–93), and growth parameters used in the Alaskan sablefish assessment model were, thus, too large. In addition, a comparison of the bias-corrected 1981–93 data and the 1996–2004 data showed that, in more recent years, sablefish grew larger and growth differed among regions. The updated growth information improves the fit of the data to the sablefish stock assessment model with biologically reasonable results. These findings indicate that when the updated growth data (1996–2004) are used in the existing sablefish assessment model, estimates of fishing mortality increase slightly and estimates of female spawning biomass decrease slightly. This study provides evidence of the importance of periodically revisiting biological parameter estimates, especially as data accumulate, because the addition of more recent data often will be more biologically realistic. In addition, it exemplifies the importance of correcting biases from sampling that may contribute to erroneous parameter estimates.