Support for Integrated Ecosystem Assessments of NOAA’s National Estuarine Research Reserves System (NERRS), Volume II: Assessment of Ecological Condition and Stressor Impacts in Subtidal Waters of the North Carolina NERRS

A study was conducted to assess the status of ecological condition and potential human-health risks in subtidal estuarine waters throughout the North Carolina National Estuarine Research Reserve System (NERRS) (Currituck Sound, Rachel Carson, Masonboro Island, and Zeke’s Island). Field work was conducted in September 2006 and incorporated multiple indicators of ecosystem condition including measures of water quality (dissolved oxygen, salinity, temperature, pH, nutrients and chlorophyll, suspended solids), sediment quality (granulometry, organic matter content, chemical contaminant concentrations), biological condition (diversity and abundances of benthic fauna, fish contaminant levels and pathologies), and human dimensions (fish-tissue contaminant levels relative to human-health consumption limits, various aesthetic properties). A probabilistic sampling design permitted statistical estimation of the spatial extent of degraded versus non-degraded condition across these estuaries relative to specified threshold levels of the various indicators (where possible). With some exceptions, the status of these reserves appeared to be in relatively good to fair ecological condition overall, with the majority of the area (about 54%) having various water quality, sediment quality, and biological (benthic) condition indicators rated in the healthy to intermediate range of corresponding guideline thresholds. Only three stations, representing 10.5% of the area, had one or more of these indicators rated as poor/degraded in all three categories. While such a conclusion is encouraging from a coastal management perspective, it should be viewed with some caution. For example, although co-occurrences of adverse biological and abiotic environmental conditions were limited, at least one indicator of ecological condition rated in the poor/degraded range was observed over a broader area (35.5%) represented by 11 of the 30 stations sampled. In addition, the fish-tissue contaminant data were not included in these overall spatial estimates; however, the majority of samples (77% of fish that were analyzed, from 79%, of stations where fish were caught) contained inorganic arsenic above the consumption limits for human cancer risks, though most likely derived from natural sources. Similarly, aesthetic indicators are not reflected in these spatial estimates of ecological condition, though there was evidence of noxious odors in sediments at many of the stations. Such symptoms reflect a growing realization that North Carolina estuaries are under multiple pressures from a variety of natural and human influences. These data also suggest that, while the current status of overall ecological condition appears to be good to fair, long-term monitoring is warranted to track potential changes in the future. This study establishes an important baseline of overall ecological condition within NC NERRS that can be used to evaluate any such future changes and to trigger appropriate management actions in this rapidly evolving coastal environment. (PDF contains 76 pages)

Simulation of Tail Weight Distributions in Biological Year 1986–2006 Landings of Brown Shrimp, Farfantepenaeus aztecus, from the Northern Gulf of Mexico Fishery

Size distribution within re- ported landings is an important aspect of northern Gulf of Mexico penaeid shrimp stock assessments. It reflects shrimp population characteristics such as numerical abundance of various sizes, age structure, and vital rates (e.g. recruitment, growth, and mortality), as well as effects of fishing, fishing power, fishing practices, sampling, size-grading, etc. The usual measure of shrimp size in archived landings data is count (C) the number of shrimp tails (abdomen or edible portion) per pound (0.4536 kg). Shrimp are marketed and landings reported in pounds within tail count categories. Statistically, these count categories are count class intervals or bins with upper and lower limits expressed in C. Count categories vary in width, overlap, and frequency of occurrence within the landings. The upper and lower limits of most count class intervals can be transformed to lower and upper limits (respectively) of class intervals expressed in pounds per shrimp tail, w, the reciprocal of C (i.e. w = 1/C). Age based stock assessments have relied on various algorithms to estimate numbers of shrimp from pounds landed within count categories. These algorithms required un- derlying explicit or implicit assumptions about the distribution of C or w. However, no attempts were made to assess the actual distribution of C or w. Therefore, validity of the algorithms and assumptions could not be determined. When different algorithms were applied to landings within the same size categories, they produced different estimates of numbers of shrimp. This paper demonstrates a method of simulating the distribution of w in reported biological year landings of shrimp. We used, as examples, landings of brown shrimp, Farfantepenaeus aztecus, from the northern Gulf of Mexico fishery in biological years 1986–2006. Brown shrimp biological year, Ti, is defined as beginning on 1 May of the same calendar year as Ti and ending on 30 April of the next calendar year, where subscript i is the place marker for biological year. Biological year landings encompass most if not all of the brown shrimp life cycle and life span. Simulated distributions of w reflect all factors influencing sizes of brown shrimp in the landings within a given biological year. Our method does not require a priori assumptions about the parent distributions of w or C, and it takes into account the variability in width, overlap, and frequency of occurrence of count categories within the landings. Simulated biological year distributions of w can be transformed to equivalent distributions of C. Our method may be useful in future testing of previously applied algorithms and development of new estimators based on statistical estimation theory and the underlying distribution of w or C. We also examine some applications of biological year distributions of w, and additional variables derived from them.

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.

Demersal groundfish densities in trawlable and untrawlable habitats off Washington: implications for the estimation of habitat bias in trawl surveys

Demersal groundfish densities were estimated by conducting a visual strip-transect survey via manned submersible on the continental shelf off Cape Flattery, Washington. The purpose of this study was to evaluate the statistical sampling power of the submersible survey as a tool to discriminate density differences between trawlable and untrawlable habitats. A geophysical map of the study area was prepared with side-scan sonar imagery, multibeam bathymetry data, and known locations of historical NMFS trawl survey events. Submersible transects were completed at randomly selected dive sites located in each habitat type. Significant differences in density between habitats were observed for lingcod (Ophiodon elongatus), yelloweye rockfish (Sebastes ruberrimus), and tiger rockfish (S. nigrocinctus) individually, and for “all rockfish” and “all flatfish” in the aggregate. Flatfish were more than ten times as abundant in the trawlable habitat samples than in the untrawlable samples, whereas rockfish as a group were over three times as abundant in the untrawlable habitat samples. Guidelines for sample sizes and implications for the estimation of the continental shelf trawl-survey habitat-bias are considered. We demonstrate an approach that can be used to establish sample size guidelines for future work by illustrating the interplay between statistical sampling power and 1) habitat specific-density differences, 2) variance of density differences, and 3) the proportion of untrawlable area in a habitat.

Methods of sensory evaluation of quality and application of statistical methods in sensory evaluation problems with special reference to fishery products

The quality of raw and processed fishery products depend on several factors like physiological conditions at the time of capture, morphological differences, rigor mortis, species, rate of icing and subsequent storage conditions. Sensory evaluation is still the most reliable method for evaluation of the freshness of raw processed fishery products. Sophisticated methods like Intelectron fish tester, cell fragility technique and chemical and bacteriological methods like estimation of trimethylamine, hypoxanthine, carbonyl compounds, volatile acid and total bacterial count have no doubt been developed for accessing the spoilage in fish products.