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.

NMFS Cooperative SharkTagging Program, 1962-93: An Atlas of SharkTag and Recapture Data

The National Marine Fisheries Service (NMFS) Cooperative Shark Tagging Program (CSTP) is part of continuing research directed to the study of the biology of large Atlantic sharks. The CSTP was initiated in 1962 at the Sandy Hook Laboratory in New Jersey under the Department of Interior's U.S. Fish and Wildlife Service (USFWS). During the late 1950's and early 1960's, sharks were considered a liability to the economy of resort communities, of little or no commercial value, and a detriment to fishermen in areas where sharks might damage expensive fishing gear or reduce catches of more commercially valuable species.

Evaluating the quality of bycatch data and bycatch estimates among disparate fisheries

In 2006, the National Marine Fisheries Service, NOAA, initiated development of a national bycatch report that would provide bycatch estimates for U.S. commercial fisheries at the fishery and species levels for fishes, marine mammals, sea turtles, and seabirds. As part of this project, the need to quantify the relative quality of available bycatch data and estimation methods was identified. Working collaboratively with fisheries managers and scientists across the nation, a system of evaluation was developed. Herein we describe the development of this system (the “tier system”), its components, and its application. We also discuss the value of the tier system in allowing fisheries managers to identify research needs and efficiently allocate limited resources toward those areas that will result in the greatest improvement to bycatch data and estimation quality.

Characterization of the U.S. Gulf of Mexico and South Atlantic penaeid and rock shrimp fisheries based on observer data

In July 2007, a mandatory Federal observer program was implemented to characterize the U.S. Gulf of Mexico penaeid shrimp (Farfantepenaeus aztecus, F. duorarum, and Litopenaeus setiferus) fishery. In June 2008, the program expanded to include the South Atlantic penaeid and rock shrimp, Sicyonia spp., fisheries. Data collected from 10,206 tows during 5,197 sea days of observations were analyzed by geographical area and target species. The majority of tows (~70%) sampled were off the coasts of Texas and Louisiana. Based on total hours towed, the highest concentrated effort occurred off South Texas and southwestern Florida. Gear information, such as net characteristics, bycatch reduction devices, and turtle excluder devices were fairly consistent among areas and target species. By species categories, finfish comprised the majority (≥57%) of the catch composition in the Gulf of Mexico and South Atlantic penaeid shrimp fisheries, while in the South Atlantic rock shrimp fishery the largest component (41%) was rock shrimp. Bycatch to shrimp ratios were lower than reported in previous studies for the Gulf of Mexico penaeid shrimp fishery. These decreased ratios may be attributed to several factors, notably decreased shrimp effort and higher shrimp catch per unit of effort (CPUE) in recent years. CPUE density surface plots for several species of interest illustrated spatial differences in distribution. Hot Spot Analyses for shrimp (penaeid and rock) and bycatch species identified areas with significant clustering of high or low CPUE values. Spatial and temporal distribution of protected species interactions were documented.

Estimation of bycatch in shrimp trawl fisheries: a comparison of estimation methods using field data and simulated data

Bycatch, or the incidental catch of nontarget organisms during fi shing operations, is a major issue in U.S. shrimp trawl fisheries. Because bycatch is typically discarded at sea, total bycatch is usually estimated by extrapolating from an observed bycatch sample to the entire fleet with either mean-per-unit or ratio estimators. Using both field observations of commercial shrimp trawlers and computer simulations, I compared five methods for generating bycatch estimates that were used in past studies, a mean-per-unit estimator and four forms of the ratio estimator, respectively: 1) the mean fish catch per unit of effort, where unit effort was a proxy for sample size, 2) the mean of the individual fish to shrimp ratios, 3) the ratio of mean fish catch to mean shrimp catch, 4) the mean of the ratios of fish catch per time fished (a variable measure of effort), and 5) the ratio of mean fish catch per mean time fished. For field data, different methods used to estimate bycatch of Atlantic croaker, spot, and weakfish yielded extremely different results, with no discernible pattern in the estimates by method, geographic region, or species. Simulated fishing fleets were used to compare bycatch estimated by the fi ve methods with “actual” (simulated) bycatch. Simulations were conducted by using both normal and delta lognormal distributions of fish and shrimp and employed a range of values for several parameters, including mean catches of fish and shrimp, variability in the catches of fish and shrimp, variability in fishing effort, number of observations, and correlations between fish and shrimp catches. Results indicated that only the mean per unit estimators provided statistically unbiased estimates, while all other methods overestimated bycatch. The mean of the individual fish to shrimp ratios, the method used in the South Atlantic Bight before the 1990s, gave the most biased estimates. Because of the statistically significant two- and 3-way interactions among parameters, it is unlikely that estimates generated by one method can be converted or corrected to estimates made by another method: therefore bycatch estimates obtained with different methods should not be compared directly.

Estimated ages of red porgy (Pagrus pagrus) from fishery-dependent and fishery-independent data and a comparison of growth parameters

The red porgy, Pagrus pagrus, is an important reef fish in several offshore fisheries along the southeastern United States. We examined samples from North Carolina through southeast Florida from recreational (headboat) and commercial (hook and line) fisheries, as well as samples from a fishery-independent source. Red porgy attain a maximum age of at least 18 years and 733 mm total length. The weight-length relationship is represented by the ln-ln transformed equation: W = 8.85 × 10–6(L)3.06, where W = whole weight in grams, and L = total length in mm. The von Bertalanffy growth equation fitted to the most recent, back-calculated lengths from all the samples is Lt = 644(1 – e –0.15(t + 0.76)). Our study revealed a difference in mean length at age of red porgy from the three sources. Red porgy in fishery-independent collections were smaller at age than specimens examined from fishery-dependent sources. The difference in length-at-age may be related to gear selectivity and have important consequences in the assessment of fish stocks.

Integration of submersible transect data and high-resolution multibeam sonar imagery for a habitat-based groundfish assessment of Heceta Bank, Oregon

In the face of dramatic declines in groundfish populations and a lack of sufficient stock assessment information, a need has arisen for new methods of assessing groundfish populations. We describe the integration of seafloor transect data gathered by a manned submersible with high-resolution sonar imagery to produce a habitat-based stock assessment system for groundfish. The data sets used in this study were collected from Heceta Bank, Oregon, and were derived from 42 submersible dives (1988–90) and a multibeam sonar survey (1998). The submersible habitat survey investigated seafloor topography and groundfish abundance along 30-minute transects over six predetermined stations and found a statistical relationship between habitat variability and groundfish distribution and abundance. These transects were analyzed in a geographic information system (GIS) by using dynamic segmentation to display changes in habitat along the transects. We used the submersible data to extrapolate fish abundance within uniform habitat patches over broad areas of the bank by means of a habitat classification based on the sonar imagery. After applying a navigation correction to the submersible-based habitat segments, a good correlation with major boundaries on the backscatter and topographic boundaries on the imagery were apparent. Extrapolation of the extent of uniform habitats was made in the vicinity of the dive stations and a preliminary stock assessment of several species of demersal fish was calculated. Such a habitat-based approach will allow researchers to characterize marine communities over large areas of the seafloor.