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.

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.

Age validation of Dover sole (Microstomus pacificus) by means of bomb radiocarbon

We used bomb radiocarbon (14C) in this age validation study of Dover sole (Microstomus pacificus). The otoliths of Dover sole, a commercially important fish in the North Pacific, are difficult to age and ages derived from the current break-andburn method were not previously validated. The otoliths used in this study were chosen on the basis of estimated birth year and for the ease of interpreting growth zone patterns. Otolith cores, material representing years 0 through 3, were isolated and analyzed for 14C. Additionally, a small number of otoliths with difficult-to-interpret growth patterns were analyzed for 14C to help determine age interpretation. The measured Dover sole 14C values in easier-to-interpret otoliths were compared with a 14C reference chronology for Pacific halibut (Hippoglossus stenolepis) in the North Pacific. We used an objective statistical analysis where sums of squared residuals between otolith 14C values of Dover sole and the reference chronology were examined. Our statistical analysis also included a procedure where the Dover sole 14C values were standardized to the reference chronology. These procedures allowed an evaluation of aging error. The 14C results indicated that the Dover sole age estimates from the easier-to-interpret otoliths with the break-and-burn method are accurate. This study validated Dover sole ages from 8 to 47 years.

Age, growth, and reproduction of dolphinfish (Coryphaena hippurus) caught off the coast of North Carolina

Age, growth, and reproductive data were obtained from dolphinfish (Coryphaena hippurus, size range: 89 to 1451 mm fork length [FL]) collected between May 2002 and May 2004 off North Carolina. Annual increments from scales (n=541) and daily increments from sagittal otoliths (n=107) were examined; estimated von Bertalanffy parameters were L∞ (asymptotic length)=1299 mm FL and k (growth coefficient)=1.08/yr. Daily growth increments reduced much of the residual error in length-at-age estimates for age-0 dolphinfish; the estimated average growth rate was 3.78 mm/day during the first six months. Size at 50% maturity was slightly smaller for female (460 mm FL) than male (475 mm FL) dolphinfish. Based on monthly length-adjusted gonad weights, peak spawning occurs from April through July off North Carolina; back-calculated hatching dates from age-0 dolphinfish and prior reproductive studies on the east coast of Florida indicate that dolphinfish spawning occurs year round off the U.S. east coast and highest levels range from January through June. No major changes in length-at-age or size-at-maturity have occurred since the early 1960s, even after substantial increases in fishery landings.

Estimates of growth and comparisons of growth rates determined from length- and age-based models for populations of purple wrasse (Notolabrus fucicola)

Growth of a temperate reefa-ssociated fish, the purple wrasse (Notolabrus fucicola), was examined from two sites on the east coast of Tasmania by using age- and length-based models. Models based on the von Bertalanffy growth function, in the standard and a reparameterized form, were constructed by using otolith-derived age estimates. Growth trajectories from tag-recaptures were used to construct length-based growth models derived from the GROTAG model, in turn a reparameterization of the Fabens model. Likelihood ratio tests (LRTs) determined the optimal parameterization of the GROTAG model, including estimators of individual growth variability, seasonal growth, measurement error, and outliers for each data set. Growth models and parameter estimates were compared by bootstrap confidence intervals, LRTs, and randomization tests and plots of bootstrap parameter estimates. The relative merit of these methods for comparing models and parameters was evaluated; LRTs combined with bootstrapping and randomization tests provided the most insight into the relationships between parameter estimates. Significant differences in growth of purple wrasse were found between sites in both length- and age-based models. A significant difference in the peak growth season was found between sites, and a large difference in growth rate between sexes was found at one site with the use of length-based models.

Age determination and growth of the night shark (Carcharhinus signatus) off the northeastern Brazilian coast

Age and growth of the night shark (Carcharhinus signatus) from areas off northeastern Brazil were determined from 317 unstained vertebral sections of 182 males (113–215 cm total length [TL]), 132 females (111.5–234.9 cm) and three individuals of unknown sex (169–242 cm). Although marginal increment (MI) analysis suggests that band formation occurs in the third and fourth trimesters in juveniles, it was inconclusive for adults. Thus, it was assumed that one band is formed annually. Births that occur over a protracted period may be the most important source of bias in MI analysis. An estimated average percent error of 2.4% was found in readings for individuals between two and seventeen years. The von Bertalanffy growth function (VBGF) showed no significant differences between sexes, and the model derived from back-calculated mean length at age best represented growth for the species (L∞=270 cm, K=0.11/yr, t0=–2.71 yr) when compared to the observed mean lengths at age and the Fabens’ method. Length-frequency analysis on 1055 specimens (93–260 cm) was used to verify age determination. Back-calculated size at birth was 66.8 cm and maturity was reached at 180–190 cm (age 8) for males and 200–205 cm (age ten) for females. Age composition, estimated from an age-length key, indicated that juveniles predominate in commercial catches, representing 74.3% of the catch. A growth rate of 25.4 cm/yr was estimated from birth to the first band (i.e. juveniles grow 38% of their birth length during the first year), and a growth rate of 8.55 cm/yr was estimated for eight- to ten-year-old adults.

Age and growth estimates of the winter skate (Leucoraja ocellata) in the western Gulf of Maine

Age and growth estimates for the winter skate (Leucoraja ocellata) were estimated from vertebral band counts on 209 fish ranging in size from 145 to 940 mm total length (TL). An index of average percent error (IAPE) of 5.8% suggests that our aging method represents a precise approach to the age assessment of L. ocellata. Marginal increments were significantly different between months (Kruskal-Wallis P<0.001) and a distinct trend of increasing monthly increment growth began in July. Estimates of von Bertalanffy growth parameters suggest that females attain a slightly larger asymptotic TL (L∞=1374 mm) than males (L∞=1218 mm) and grow more slowly (k=0.059 and 0.074, respectively). The oldest ages obtained for the winter skate were 19 years for males and 18 years for females, which corresponded to total lengths of 932 mm and 940 mm, respectively. The results indicate that the winter skate exhibits the characteristics that have made other elasmobranch populations highly susceptible to exploitation by commercial fisheries.