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.

A model of the pelagic ecossystem in the eastern tropical Pacific Ocean

English: Recent calls for a more holistic approach to fisheries management have motivated development of trophic mass-balance models of ecosystems that underlie fisheries production. We developed a model hypothesis of the pelagic ecosystem in the eastern tropical Pacific Ocean (ETP) to gain insight into the relationships among the various species in the system and to explore the ecological implications of alternative methods of harvesting tunas. We represented the biomasses of and fluxes between the principal elements in the ecosystem with Ecopath, and examined the ecosystem's dynamic, time-series behavior with Ecosim. We parameterized the model for 38 species or groups of species, and described the sources, justifications, assumptions, and revisions of our estimates of the various parameters, diet relations, fisheries landings, and fisheries discards in the model. We conducted sensitivity analyses with an intermediate version of the model, for both the Ecopath mass-balance and the dynamic trajectories predicted by Ecosim. The analysis showed that changes in the basic parameters for two components at middle trophic levels, Cephalopods and Auxis spp., exert the greatest influence on the system. When the Cephalopod Q/B and Auxis spp. P/B were altered from their initial values and the model was rebalanced, the trends of the biomass trajectories predicted by Ecosim were not sensitive, but the scaling was sensitive for several components. We described the review process the model was subjected to, which included reviews by the IATTC Purse-seine Bycatch Working Group and by a working group supported by the National Center for Ecological Analysis and Synthesis. We fitted the model to historical time series of catches per unit of effort and mortality rates for yellowfin and bigeye tunas in simulations that incorporated historical fishing effort and a climate driver to represent the effect of El Niño-Southern Oscillation-scale variation on the system. The model was designed to evaluate the possible ecological implications of fishing for tunas in various ways. We recognize that a model cannot possibly represent all the complexity of a pelagic ocean ecosystem, but we believe that the ETP model provides insight into the structure and function of the pelagic ETP. Spanish: Llamamientos recientes hacia un enfoque más holístico al ordenamiento de la pesca han motivado el desarrollo de modelos tróficos de balance de masas de los ecosistemas que sostienen la producción pesquera. Desarrollamos una hipótesis modelo del ecosistema pelágico en el Océano Pacífico oriental tropical (POT) con miras a mejorar los conocimientos de las relaciones entre las distintas especies en el sistema y explorar las implicaciones ecológicas de métodos alternativos de capturar atunes. Con Ecopath representamos las biomasas de los elementos principales en el ecosistema, y los flujos entre los mismos, y con Ecosim examinamos el comportamiento dinámico del ecosistema con el tiempo. Parametrizamos el modelo para 38 especies o grupos de especies (denominados “componentes” del modelo), y describimos las fuentes, justificaciones, supuestos, y revisiones de nuestras estimaciones de los distintos parámetros, relaciones basadas en dieta, capturas retenidas de las pesquerías, y descartes de las mismas en el modelo. Realizamos análisis de sensibilidad con una versión intermedia del modelo, para el balance de masas de Ecopath y las trayectorias dinámicas predichas por Ecosim también. El análisis demostró que cambios en los parámetros básicos para dos componentes en niveles tróficos medianos, Cefalópodos y Auxis spp., ejercieron la mayor influencia sobre el sistema. Cuando se alteraron el Q/B de los Cefalópodos y el P/B de los Auxis spp. de sus valores iniciales y se balanceó el modelo de nuevo, las tendencias de las trayectorias de la biomasa predichas por Ecosim no fueron sensibles, pero la escala fue sensible para varios componentes. Describimos el proceso de revisión al que fue sujeto el modelo, inclusive revisiones por el Grupo de Trabajo sobre Captura Incidental de la CIAT y un grupo de trabajo apoyado por el Centro Nacional para Síntesis y Análisis Ecológicos. Ajustamos el modelo a series de tiempo históricas de capturas por unidad de esfuerzo y tasas de mortalidad de atunes aleta amarilla y patudo en simulaciones que incorporaron esfuerzo de pesca histórico e impulsos climáticos para representar el efecto de variaciones a escala de El Niño-Oscilación del Sur sobre el sistema. El modelo fue diseñado para evaluar las posibles implicaciones ecológicas de la pesca atunera de varias formas. Reconocemos la imposibilidad de que el modelo represente toda la complejidad de un ecosistema oceánico pelágico, pero creemos que el modelo del POT mejora los conocimientos de la estructura y función del POT pelágico.

Trends in nitrate-nitrogen, nitrite-nitrogen and phosphorus concentration in Ebonyi river, Nigeria

A baseline study of the limnological integrity of Ebonyi River, a tropical lotic system in south-eastern Nigeria was conducted between September 2006 and February 2008 to assess its potential in enhancing fisheries production for the benefit of the rural poor, who depend on the resources of the river for survival. The parameters measured were nitrate-nitrogen, nitrite-nitrogen and phosphate-phosphorus. Results show that nitrate varied between 40.43mg/L in September 2006 and 1.73mg/L in December 2007, Showing significant difference (P<0.01) among months. The values recorded for nitrites varied between 0.2mg/L in September 2006 and 0.4mg/L in February 2008, showing significant (P<0.01) variation among months. Values recorded for phosphorus was highest (0.05mg/L) in the month of October 2006 while the least mean value (0.32mg/L) was recorded in the month of May 2007 and showed significant (p<0.01) variation in monthly means. It was concluded that the values of the measured parameters falls within tolerable range for enhanced fisheries development in the area. KEYWORDS: Limnology, Tropical, River, Monthly, Mean, Variation

An assessment of mussel mortality caused by a drop in the water level of Lake Kariba

The level of Lake Kariba steadily fell during the period 1 June 1979 to 2 Feb 1980, except for a 2-wk period during Dec when it was allowed to rise slightly. Following this the level was again drawn down in anticipation of the Upper Zambezi flood water reaching the lake. At its highest level in June 1979 the lake was 487.42 m above sea level but by Feb 1980 it had dropped to 484.53 m, a total drop of 2.89 m. This left a considerable area of exposed shoreline and a large number of stranded mussels. This report presents the results of an attempt to estimate the mussel mortality, carried out from 28 Jan to 1 Feb 1980. The study area extended from the Charara river mouth to Andora harbour with a total of 24 stations.

Biological report - fish mortality, River Conder, 21.7.79

Following a large scale fish mortality caused by a considerable spillage of a pentachlbrophenol-based disinfectant from a mushroom farm, a biological survey of the river Conder was undertaken on 23.7.79. A spillage of this same make of disinfectant from this mushroom farm in February 1976 caused a similar large scale fish mortality. On that occasion, the pollution appeared to have little effect on the benthic invertebrates of the river Conder, with no dead organisms found, although the small beck which received the spillage was quite badly affected. This recent spillage of disinfectant involved a considerably greater volume (approximately 900 gallons) than the 1976 discharge and so this report looks at whether invertebrate fauna of the main river were affected on this occasion. Kick samples were taken at a number of sites in the river Conder from the confluence with the polluted beck to Conder Green. A control sample was also taken upstream of the confluence with the polluted beck. All samples were examined on site and live and dead invertebrates noted and identified as accurately as possible.

Estimation of discard mortality of sablefish (Anoplopoma fimbria) in Alaska longline fisheries

Sablefish (Anoplopoma fimbria) are often caught incidentally in longline fisheries and discarded, but the extent of mortality after release is unknown, which creates uncertainty for estimates of total mortality. We analyzed data from 10,427 fish that were tagged in research surveys and recovered in surveys and commercial fisheries up to 19 years later and found a decrease in recapture rates for fish originally captured at shallower depths (210–319 m) during the study, sustaining severe hooking injuries, and sustaining amphipod predation injuries. The overall estimated discard mortality rate was 11.71%. This estimate is based on an assumed survival rate of 96.5% for fish with minor hooking injuries and the observed recapture rates for sablefish at each level of severity of hook injury. This estimate may be lower than what actually occurs in commercial fisheries because fish are likely not handled as carefully as those in our study. Comparing our results with data on the relative occurrence of the severity of hooking injuries in longline fisheries may lead to more accurate accounting of total mortality attributable to fishing and to improved management of this species.

Analysis of permanent magnets as elasmobranch bycatch reduction devices in hook-and-line and longline trials

Previous studies indicate that elasmobranch fishes (sharks, skates and rays) detect the Earth’s geomagnetic field by indirect magnetoreception through electromagnetic induction, using their ampullae of Lorenzini. Applying this concept, we evaluated the capture of elasmobranchs in the presence of permanent magnets in hook-and-line and inshore longline fishing experiments. Hooks with neodymium-iron-boron magnets significantly reduced the capture of elasmobranchs overall in comparison with control and procedural control hooks in the hook-and-line experiment. Catches of Atlantic sharpnose shark (Rhizoprionodon terraenovae) and smooth dogfish (Mustelus canis) were signif icantly reduced with magnetic hook-and-line treatments, whereas catches of spiny dogfish (Squalus acanthias) and clearnose skate (Raja eglanteria) were not. Longline hooks with barium-ferrite magnets significantly reduced total elasmobranch capture when compared with control hooks. In the longline study, capture of blacktip sharks (Carcharhinus limbatus) and southern stingrays (Dasyatis americana) was reduced on magnetic hooks, whereas capture of sandbar shark (Carcharhinus plumbeus) was not affected. Teleosts, such as red drum (Sciaenops ocellatus), Atlantic croaker (Micropogonias undulatus), oyster toadfish (Opsanus tau), black sea bass (Centropristis striata), and the bluefish (Pomatomas saltatrix), showed no hook preference in either hook-and-line or longline studies. These results indicate that permanent magnets, although eliciting species-specific capture trends, warrant further investigation in commercial longline and recreational fisheries, where bycatch mortality is a leading contributor to declines in elasmobranch populations.

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.

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.

Distribution, growth, and mortality of sailfish (Istiophorus platypterus) larvae in the northern Gulf of Mexico

Ichthyoplankton surveys were conducted in shelf and slope waters of the northern Gulf of Mexico during the months of May–September in 2005 and 2006 to investigate the potential role of this region as spawning and nursery habitat of sailfish (Istiophorus platypterus). During the two-year study, 2426 sailfish larvae were collected, ranging in size from 2.0 to 24.3 mm standard length. Mean density for all neuston net collections (n=288) combined was 1.5 sailfish per 1000 m2, and maximum density was observed within frontal features created by hydrodynamic convergence (2.3 sailfish per 1000 m2). Sagittal otoliths were extracted from 1330 larvae, and otolith microstructure analysis indicated that the sailfish ranged in age from 4 to 24 days after hatching (mean=10.5 d, standard deviation [SD]=3.2 d). Instantaneous growth coefficients (g) among survey periods (n=5) ranged from 0.113 to 0.127, and growth peaked during July 2005 collections when density within frontal features was highest. Daily instantaneous mortality rates (Z) ranged from 0.228 to 0.381, and Z was indexed to instantaneous weight-specific growth (G) to assess stage-specific production potential of larval cohorts. Ratios of G to Z were greater than 1.0 for all but one cohort examined, indicating that cohorts were gaining biomass during the majority of months investigated. Stage-specific production potential, in combination with catch rates and densities of larvae, indicates that the Gulf of Mexico likely represents important spawning and nursery habitat for sailfish.

Environmental characterization of seasonal trends and foraging habitat of bottlenose dolphins (Tursiops truncatus) in northern Gulf of Mexico bays

A description of the foraging habitat of a cetacean species is critical for conservation and effective management. We used a fine-scale microhabitat approach to examine patterns in bottlenose dolphin (Tursiops truncatus) foraging distribution in relation to dissolved oxygen, turbidity, salinity, water depth, water temperature, and distance from shore measurements in a highly turbid estuary on the northern Gulf of Mexico. In general, environmental variation in the Barataria Basin marine environment comprises three primary axes of variability (i.e., factors: temperature and dissolved oxygen, salinity and turbidity, and distance and depth) that represent seasonal, spatial-seasonal, and spatial scales, respectively. Foraging sites were differentiated from nonforaging sites by significant differences among group size, temperature, turbidity, and season. Habitat selection analysis on individual variables indicated that foraging was more frequently observed in waters 4–6 m deep, 200–500 m from shore, and at salinity values of around 20 psu. This fine-scale and multivariate approach represents a useful method of exploring the complexity, gradation, and detail of the relationships between environmental variables and the foraging distribution patterns of bottlenose dolphin.

Assessing trends in the density of Atlantic croaker (Micropogonias undulatus): a comparison of passive acoustic and trawl methods

Using data collected simultaneously from a trawl and a hydrophone, we found that temporal and spatial trends in densities of juvenile Atlantic croaker (Micropogonias undulatus) in the Neuse River estuary in North Carolina can be identified by monitoring their sound production. Multivariate analysis of covariance (MA NCOVA) revealed that catch per unit of effort (CPUE) of Atlantic croaker had a significant relationship with the dependent variables of sound level and peak frequency of Atlantic croaker calls. Tests of between-subject correspondence failed to detect relationships between CPUE and either of the call parameters, but statistical power was low. Williamson’s index of spatial overlap indicated that call detection rate (expressed by a 0–3 calling index) was correlated in time and space with Atlantic croaker CPUE. The correspondence between acoustic parameters and trawl catch rates varied by month and by habitat. In general, the calling index had a higher degree of overlap with this species’ density than did the received sound level of their calls. Classification and regression tree analysis identified calling index as the strongest correlate of CPUE. Passive acoustics has the potential to be an inexpensive means of identifying spatial and temporal trends in abundance for soniferous fish species.

Prediction of discard mortality for Alaskan crabs after exposure to freezing temperatures, based on a reflex impairment index

Millions of crabs are sorted and discarded in freezing conditions each year in Alaskan fisheries for Tanner crab (Chionoecetes bairdi) and snow crab (C. opilio). However, cold exposures vary widely over the fishing season and among different vessels, and mortalities are difficult to estimate. A shipboard experiment was conducted to determine whether simple behavioral observations can be used to evaluate crab condition after low-temperature exposures. Crabs were systematically subjected to cold in seven different exposure treatments. They were then tested for righting behavior and six different ref lex actions and held to monitor mortality. Crabs lost limbs, showed ref lex impairment, and died in direct proportion to increases in cold exposure. Righting behavior was a poor predictor of mortality, whereas reflex impairment (scored as the sum of reflex actions that were lost) was an excellent predictor. This composite index could be measured quickly and easily in hand, and logistic regression revealed that the relationship between reflex impairment and mortality correctly predicted 80.0% of the mortality and survival for C. bairdi, and 79.4% for C. opilio. These relationships provide substantial improvements over earlier approaches to mortality estimation and were independent of crab size and exposure temperature.

Modified sorting technique to mitigate the collateral mortality of trawled school prawns (Metapenaeus macleayi)

The potential for changes to onboard handling practices in order to improve the fate of juvenile school prawns (Metapenaeus macleayi) discarded during trawling were investigated in two Australian rivers (Clarence and Hunter) by comparing a purpose-built, water-filled sorting tray against a conventional dry tray across various conditions, including the range of typical delays before the start of sorting the catch (2 min vs. 15 min). Juvenile school prawns (n= 5760), caught during 32 and 16 deployments in each river, were caged and sacrificed at four times: immediately (T0), and at 24 (T24), 72 (T72), and 120 (T12 0) hours after having been discarded. In both rivers, most mortalities occurred between T0 and T24 and, after adjusting for control deaths (<12%), were greatest for the 15-min conventional treatment (up to 41% at T120). Mixed-effects logistic models revealed that in addition to the sampling time, method of sorting, and delay in sorting, the weight of the catch, salinity, and percentage cloud cover were significant predictors of mortality. Although trawling caused some mortalities and comparable stress (measured as L -lactate) in all school prawns, use of the water tray lessened the negative impacts of some of the above factors across both the 2-min and 15-min delays in sorting so that the overall discard mortality was reduced by more than a third. When used in conjunction with selective trawls, widespread application of the water tray should help to improve the sustainability of trawling for school prawns.