A-SCALA: an age-structured statistical catch-at-length analysis for assessing tuna stocks in the eastern tropical Pacific Ocean

English: We describe an age-structured statistical catch-at-length analysis (A-SCALA) based on the MULTIFAN-CL model of Fournier et al. (1998). The analysis is applied independently to both the yellowfin and the bigeye tuna populations of the eastern Pacific Ocean (EPO). We model the populations from 1975 to 1999, based on quarterly time steps. Only a single stock for each species is assumed for each analysis, but multiple fisheries that are spatially separate are modeled to allow for spatial differences in catchability and selectivity. The analysis allows for error in the effort-fishing mortality relationship, temporal trends in catchability, temporal variation in recruitment, relationships between the environment and recruitment and between the environment and catchability, and differences in selectivity and catchability among fisheries. The model is fit to total catch data and proportional catch-at-length data conditioned on effort. The A-SCALA method is a statistical approach, and therefore recognizes that the data collected from the fishery do not perfectly represent the population. Also, there is uncertainty in our knowledge about the dynamics of the system and uncertainty about how the observed data relate to the real population. The use of likelihood functions allow us to model the uncertainty in the data collected from the population, and the inclusion of estimable process error allows us to model the uncertainties in the dynamics of the system. The statistical approach allows for the calculation of confidence intervals and the testing of hypotheses. We use a Bayesian version of the maximum likelihood framework that includes distributional constraints on temporal variation in recruitment, the effort-fishing mortality relationship, and catchability. Curvature penalties for selectivity parameters and penalties on extreme fishing mortality rates are also included in the objective function. The mode of the joint posterior distribution is used as an estimate of the model parameters. Confidence intervals are calculated using the normal approximation method. It should be noted that the estimation method includes constraints and priors and therefore the confidence intervals are different from traditionally calculated confidence intervals. Management reference points are calculated, and forward projections are carried out to provide advice for making management decisions for the yellowfin and bigeye populations. Spanish: Describimos un análisis estadístico de captura a talla estructurado por edad, A-SCALA (del inglés age-structured statistical catch-at-length analysis), basado en el modelo MULTIFAN- CL de Fournier et al. (1998). Se aplica el análisis independientemente a las poblaciones de atunes aleta amarilla y patudo del Océano Pacífico oriental (OPO). Modelamos las poblaciones de 1975 a 1999, en pasos trimestrales. Se supone solamente una sola población para cada especie para cada análisis, pero se modelan pesquerías múltiples espacialmente separadas para tomar en cuenta diferencias espaciales en la capturabilidad y selectividad. El análisis toma en cuenta error en la relación esfuerzo-mortalidad por pesca, tendencias temporales en la capturabilidad, variación temporal en el reclutamiento, relaciones entre el medio ambiente y el reclutamiento y entre el medio ambiente y la capturabilidad, y diferencias en selectividad y capturabilidad entre pesquerías. Se ajusta el modelo a datos de captura total y a datos de captura a talla proporcional condicionados sobre esfuerzo. El método A-SCALA es un enfoque estadístico, y reconoce por lo tanto que los datos obtenidos de la pesca no representan la población perfectamente. Además, hay incertidumbre en nuestros conocimientos de la dinámica del sistema e incertidumbre sobre la relación entre los datos observados y la población real. El uso de funciones de verosimilitud nos permite modelar la incertidumbre en los datos obtenidos de la población, y la inclusión de un error de proceso estimable nos permite modelar las incertidumbres en la dinámica del sistema. El enfoque estadístico permite calcular intervalos de confianza y comprobar hipótesis. Usamos una versión bayesiana del marco de verosimilitud máxima que incluye constreñimientos distribucionales sobre la variación temporal en el reclutamiento, la relación esfuerzo-mortalidad por pesca, y la capturabilidad. Se incluyen también en la función objetivo penalidades por curvatura para los parámetros de selectividad y penalidades por tasas extremas de mortalidad por pesca. Se usa la moda de la distribución posterior conjunta como estimación de los parámetros del modelo. Se calculan los intervalos de confianza usando el método de aproximación normal. Cabe destacar que el método de estimación incluye constreñimientos y distribuciones previas y por lo tanto los intervalos de confianza son diferentes de los intervalos de confianza calculados de forma tradicional. Se calculan puntos de referencia para el ordenamiento, y se realizan proyecciones a futuro para asesorar la toma de decisiones para el ordenamiento de las poblaciones de aleta amarilla y patudo.

An environmentally based growth model that uses finite difference calculus with maximum likelihood method: its application to the brackish water bivalve Corbicula japonica in Lake Abashiri, Japan

We present a growth analysis model that combines large amounts of environmental data with limited amounts of biological data and apply it to Corbicula japonica. The model uses the maximum-likelihood method with the Akaike information criterion, which provides an objective criterion for model selection. An adequate distribution for describing a single cohort is selected from available probability density functions, which are expressed by location and scale parameters. Daily relative increase rates of the location parameter are expressed by a multivariate logistic function with environmental factors for each day and categorical variables indicating animal ages as independent variables. Daily relative increase rates of the scale parameter are expressed by an equation describing the relationship with the daily relative increase rate of the location parameter. Corbicula japonica grows to a modal shell length of 0.7 mm during the first year in Lake Abashiri. Compared with the attain-able maximum size of about 30 mm, the growth of juveniles is extremely slow because their growth is less susceptible to environmental factors until the second winter. The extremely slow growth in Lake Abashiri could be a geographical genetic variation within C. japonica.

Observations on the behaviour of Polyphemus pediculus and functions of its extremities in the process of feeding. [Translation of: Materialy po biologii i gidrologii Volzhskikh vodokhranilishch, pp. 13-16. Moscow, 1965]

The author summarises observations on the behaviour of Polyphemus pediculus and functions of its extremities in the process of feeding. The crustacean Polyphemus pediculus seizes its prey, kills it and pulverises its food with the help of its extremities. Therefore for a study of its feeding method was necessary not only to have been acquainted in detail with the structure of its extremities, but also to have observed their interaction for the accomplishment of the stated functions.

Organization, functions, and achievements of the Inter-American Tropical Tuna Commission

ENGLISH: The Inter-American Tropical Tuna Commission (IATTC) operates under the authority and direction of a Convention originally entered into by the governments of Costa Rica and the United States. The Convention, which came into force in 1950, is open to the adherence by other governments whose nationals participate in the fisheries for tropical tunas in the eastern Pacific Ocean. The member nations of the Commission now are t in addition to Costa Rica and the United States, Canada, France, Japan, Mexico, Nicaragua, and Panama.This report is a description of the organization, functions, and achievements of the Commission. It has been prepared to provide in a convenient format answers to requests for information concerning the Commission. It replaces a similar, earlier report (Carroz, 1965), which is now largely outdated. SPANISH: La Comisión Interamericana del Atún Tropical (CIAT) funciona bajo la autoridad y dirección de un Convenio firmado originalmente por los gobiernos de Costa Rica y los Estados Unidos de America. El Convenio, que entro en vigencia en 1950, se encuentra libre para que otros gobiernos cuyos ciudadanos participen en la pesca de atunes tropicales en el Océano Pacifico oriental se afilien a el. Las naciones miembros de la Comisión, además de Costa Rica y los Estados Unidos, son Cañada, Francia, Japón, México, Nicaragua y Panamá. Este informe es una descripción de la organización, funciones y resultados de la Comisión. Ha sido preparado para suministrar en forma conveniente respuestas a preguntas sobre la Comisión. Reemplaza un informe anterior similar (Carroz 1965), que ya es anticuado en su mayor parte.

Organization, functions, and achievements of the Inter-American Tropical Tuna Commission

ENGLISH: The Inter-American Tropical Tuna Commission (IATTC) operates under the authority and direction of a Convention originally entered into by the governments of Costa Rica and the United states. The Convention, which came into force in 1950, is open to the adherence by other governments whose nationals participate in the fisheries for tropical tunas in the eastern Pacific Ocean. The member nations of the Commission now are France. Japan, Nicaragua. Panama, and the United States. This report is a description of the organization, functions. and achievements of the Commission. It has been prepared to provide in a convenient format answers to requests for information concerning the Commission. It replaces similar, earlier reports (Carroz, 1965; Spec. Rep., 1), which are now largely outdated. SPANISH: La Comisión Interamericana del Atún Tropical (CIAT) funciona bajo la autoridad y dirección de un Convenio firmado originalmente por los gobiernos de Costa Rica y los Estados Unidos de America. El Convenio, Que entre en vigencia en 1950, se encuentra libre para Que otros gobiernos cuyos ciudadanos participan en la pesca de los atunes tropicales en el Océano Pacifico oriental se afilien a 61. Las naciones actuales que son miembros de la Comisión son: Francia, Japón, Nicaragua, Panamá y los Estados Unidos. Este informe es una descripci6n de la organización, funciones y resultados de la Comisión. Se ha preparado para suministrar en forma conveniente, informaci6n sobre la Comisión. Este informe renueva otros anteriores similares (Carroz, 1965; CIAT, Inf. Esp., 1) que en su mayor parte determinados.

Organization, functions, and achievements of the Inter-American Tropical Tuna Commission

This report is a description of the organization, functions, and achievements of the IATTC. It has been prepared to provide, in a convenient format, answers to requests for information concerning the IATTC. It replaces similar, earlier reports (Carroz, 1965; IATTC Spec. Rep., 1 and 5), which are now largely outdated. In order to make each section of the report independent of the others, some aspects of the IATTC are described in more than one section. For example, work on the early life history of tunas financed by the Overseas Fishery Cooperation Foundation of Japan is mentioned in the subsection entitled Finance, the subsection entitled Biology of tunas and billfishes, and the section entitled RELATIONS WITH OTHER ORGANIZATIONS. Due to space constraints, however, it is not possible to describe the IATTC's activities in detail in this report. Additional information is available in publications of the IATTC, listed in Appendix 6, and in its web site, www.iattc.org. Many abbreviations are used in this report. The names of the organizations or the terms are written out the first time they are used, and, for convenience, they are also listed in the Glossary.

Interactions of age-dependent mortality and selectivity functions in age-based stock assessment models

The natural mortality rate (M) of fish varies with size and age, although it is often assumed to be constant in stock assessments. Misspecification of M may bias important assessment quantities. We simulated fishery data, using an age-based population model, and then conducted stock assessments on the simulated data. Results were compared to known values. Misspecification of M had a negligible effect on the estimation of relative stock depletion; however, misspecification of M had a large effect on the estimation of parameters describing the stock recruitment relationship, age-specific selectivity, and catchability. If high M occurs in juvenile and old fish, but is misspecified in the assessment model, virgin biomass and catchability are often poorly estimated. In addition, stock recruitment relationships are often very difficult to estimate, and steepness values are commonly estimated at the upper bound (1.0) and overfishing limits tend to be biased low. Natural mortality can be estimated in assessment models if M is constant across ages or if selectivity is asymptotic. However if M is higher in old fish and selectivity is dome-shaped, M and the selectivity cannot both be adequately estimated because of strong interactions between M and selectivity.

Integrating methods for determining length-at-age to improve growth estimates for two large scombrids

Fish growth is commonly estimated from length-at-age data obtained from otoliths. There are several techniques for estimating length-at-age from otoliths including 1) direct observed counts of annual increments; 2) age adjustment based on a categorization of otolith margins; 3) age adjustment based on known periods of spawning and annuli formation; 4) back-calculation to all annuli, and 5) back-calculation to the last annulus only. In this study we compared growth estimates (von Bertalanffy growth functions) obtained from the above five methods for estimating length-at-age from otoliths for two large scombrids: narrow-barred Spanish mackerel (Scomberomorus commerson) and broad-barred king mackerel (Scomberomorus semifasciatus). Likelihood ratio tests revealed that the largest differences in growth occurred between the back-calculation methods and the observed and adjusted methods for both species of mackerel. The pattern, however, was more pronounced for S. commerson than for S. semifasciatus, because of the pronounced effect of gear selectivity demonstrated for S. commerson. We propose a method of substituting length-at-age data from observed or adjusted methods with back-calculated length-at-age data to provide more appropriate estimates of population growth than those obtained with the individual methods alone, particularly when faster growing young fish are disproportionately selected for. Substitution of observed or adjusted length-at-age data with back-calculated length-at-age data provided more realistic estimates of length for younger ages than observed or adjusted methods as well as more realistic estimates of mean maximum length than those derived from backcalculation methods alone.

Beverton and Holt equations: spreadsheet functions and uncertainty

This contribution illustrates how modern spreadsheets aid the calculation and visualization of yield models and how the effects of uncertainties may be incorporated using Monte Carlo simulation. It is argued that analogous approaches can be implemented for other assessment models of simple to medium complexity justifying wider use of spreadsheets in fisheries analysis and training.

Maximum likelihood estimation of mortality and growth with individual variability from multiple length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock and derive the underlying length distribution of the population and the catch when there is individual variability in the von Bertalanffy growth parameter L∞. The model is flexible enough to accommodate 1) any recruitment pattern as a function of both time and length, 2) length-specific selectivity, and 3) varying fishing effort over time. The maximum likelihood method gives consistent estimates, provided the underlying distribution for individual variation in growth is correctly specified. Simulation results indicate that our method is reasonably robust to violations in the assumptions. The method is applied to tiger prawn data (Penaeus semisulcatus) to obtain estimates of natural and fishing mortality.

Problems in estimating the response functions of natural climatic recording systems [abstract]

EXTRACT (SEE PDF FOR FULL ABSTRACT): High-resolution proxy records of climate, such as varves, ice cores, and tree-rings, provide the opportunity for reconstructing climate on a year-by-year basis. In order to do so it is necessary to approximate the complex nonlinear response function of the natural recording system using linear statistical models. Three problems with this approach were discussed, and possible solutions were suggested. Examples were given from a reconstruction of Santa Barbara precipitation based on tree-ring records from Santa Barbara County.

Fisheries databases for assessing habitat functions of Cape Romain NWR

The US Fish and Wildlife Service Cape Romain National Wildlife Refuge (CRNWR) and the Center for Coastal Environmental Health and Biomolecular Research (CCEHBR) at Charleston are interested in assessing the status of our coastal resources in light of increased coastal development and recreational use. Through an Interagency Agreement (FWS #1448-40181-00-H-001), an ecological characterization was undertaken to describe the status of and potential impacts to resources at CRNWR. This report describes historic fisheries-independent or non-commercial data relevant to CRNWR that can be used to evaluate the role of the Refuge as habitat for nearshore and offshore fish species. The purpose of this document is two-fold, first to give resource managers an understanding of fisheries data that have been collected over the years and, second, to illustrate how these data can be applied to address specific management issues. This report provides an overview of historic fisheries data collected along the southeast coast, as well as basic summaries of that data relevant to CRNWR, indicating how these data can be used to address specific questions of interest to Refuge managers and biologists.

National Centers for Coastal Ocean Science (NCCOS) research highlights in the Chesapeake Bay

The Chesapeake Bay is the largest estuary in the United States. It is a unique and valuable national treasure because of its ecological, recreational, economic and cultural benefits. The problems facing the Bay are well known and extensively documented, and are largely related to human uses of the watershed and resources within the Bay. Over the past several decades as the origins of the Chesapeake’s problems became clear, citizens groups and Federal, State, and local governments have entered into agreements and worked together to restore the Bay’s productivity and ecological health. In May 2010, President Barack Obama signed Executive Order number 13508 that tasked a team of Federal agencies to develop a way forward in the protection and restoration of the Chesapeake watershed. Success of both State and Federal efforts will depend on having relevant, sound information regarding the ecology and function of the system as the basis of management and decision making. In response to the executive order, the National Oceanic and Atmospheric Administration’s National Centers for Coastal Ocean Science (NCCOS) has compiled an overview of its research in Chesapeake Bay watershed. NCCOS has a long history of Chesapeake Bay research, investigating the causes and consequences of changes throughout the watershed’s ecosystems. This document presents a cross section of research results that have advanced the understanding of the structure and function of the Chesapeake and enabled the accurate and timely prediction of events with the potential to impact both human communities and ecosystems. There are three main focus areas: changes in land use patterns in the watershed and the related impacts on contaminant and pathogen distribution and concentrations; nutrient inputs and algal bloom events; and habitat use and life history patterns of species in the watershed. Land use changes in the Chesapeake Bay watershed have dramatically changed how the system functions. A comparison of several subsystems within the Bay drainages has shown that water quality is directly related to land use and how the land use affects ecosystem health of the rivers and streams that enter the Chesapeake Bay. Across the Chesapeake as a whole, the rivers that drain developed areas, such as the Potomac and James rivers, tend to have much more highly contaminated sediments than does the mainstem of the Bay itself. In addition to what might be considered traditional contaminants, such as hydrocarbons, new contaminants are appearing in measurable amounts. At fourteen sites studied in the Bay, thirteen different pharmaceuticals were detected. The impact of pharmaceuticals on organisms and the people who eat them is still unknown. The effects of water borne infections on people and marine life are known, however, and the exposure to certain bacteria is a significant health risk. A model is now available that predicts the likelihood of occurrence of a strain of bacteria known as Vibrio vulnificus throughout Bay waters.

Modeling red sea urchin (Strongylocentrotus franciscanus) growth using six growth functions

The growth of red sea urchins (Strongylocentrotus franciscanus) was modeled by using tag-recapture data from northern California. Red sea urchins (n=211) ranging in test diameter from 7 to 131 mm were examined for changes in size over one year. We used the function Jt+1 = Jt + f(Jt) to model growth, in which Jt is the jaw size (mm) at tagging, and Jt+1 is the jaw size one year later. The function f(Jt), represents one of six deterministic models: logistic dose response, Gaussian, Tanaka, Ricker, Richards, and von Bertalanffy with 3, 3, 3, 2, 3, and 2 minimization parameters, respectively. We found that three measures of goodness of fi t ranked the models similarly, in the order given. The results from these six models indicate that red sea urchins are slow growing animals (mean of 7.2 ±1.3 years to enter the fishery). We show that poor model selection or data from a limited range of urchin sizes (or both) produces erroneous growth parameter estimates and years-to-fishery estimates. Individual variation in growth dominated spatial variation at shallow and deep sites (F=0.246, n=199, P=0.62). We summarize the six models using a composite growth curve of jaw size, J, as a function of time, t: J = A(B – e–Ct) + Dt, in which each model is distinguished by the constants A, B, C, and D. We suggest that this composite model has the flexibility of the other six models and could be broadly applied. Given the robustness of our results regarding the number of years to enter the fishery, this information could be incorporated into future fishery management plans for red sea urchins in northern California.

A model for assessing the likelihood of self-sustaining populations resulting from commercial production of triploid Suminoe oysters (Crassostrea ariakensis) in Chesapeake Bay

Culture of a non-native species, such as the Suminoe oyster (Crassostrea ariakensis), could offset the harvest of the declining native eastern oyster (Crassostrea virginica) fishery in Chesapeake Bay. Because of possible ecological impacts from introducing a fertile non-native species, introduction of sterile triploid oysters has been proposed. However, recent data show that a small percentage of triploid individuals progressively revert toward diploidy, introducing the possibility that Suminoe oysters might establish self-sustaining populations. To assess the risk of Suminoe oyster populations becoming established in Chesapeake Bay, a demographic population model was developed. Parameters modeled were salinity, stocking density, reversion rate, reproductive potential, natural and harvest-induced mortality, growth rates, and effects of various management strategies, including harvest strategies. The probability of a Suminoe oyster population becoming self-sustaining decreased in the model when oysters are grown at low salinity sites, certainty of harvest is high, mini-mum shell length-at-harvest is small, and stocking density is low. From the results of the model, we suggest adopting the proposed management strategies shown by the model to decrease the probability of a Suminoe oyster population becoming self-sustaining. Policy makers and fishery managers can use the model to predict potential outcomes of policy decisions, supporting the ability to make science-based policy decisions about the proposed introduction of triploid Suminoe oysters into the Chesapeake Bay.