Using the empirical Bayes method to estimate and evaluate bycatch rates of seabirds from individual fishing vessels

Minimizing bycatch of seabirds is a major goal of the U.S. National Marine Fisheries Service. In Alaska waters, the bycatch (i.e., inadvertent catches) of seabirds has been an incidental result of demersal groundfish longline fishery operations. Notably, the endangered short-tailed albatross (Phoebastria albatrus) has been taken in this groundfish fishery. Bycatch rates of seabirds from individual vessels may be of particular interest because vessels with high bycatch rates may not be functioning effectively with seabird avoidance gears, and there may be a need for suggestions on how to use these avoidance gears more effectively. Therefore, bycatch estimates are usually made on an individual vessel basis and then summed to obtain the total estimate for the entire fleet.

A catch-free stock assessment model with application to goliath grouper (Epinephelus itajara) off southern Florida

Many modern stock assessment methods provide the machinery for determining the status of a stock in relation to certain reference points and for estimating how quickly a stock can be rebuilt. However, these methods typically require catch data, which are not always available. We introduce a model-based framework for estimating reference points, stock status, and recovery times in situations where catch data and other measures of absolute abundance are unavailable. The specif ic estimator developed is essentially an age-structured production model recast in terms relative to pre-exploitation levels. A Bayesian estimation scheme is adopted to allow the incorporation of pertinent auxiliary information such as might be obtained from meta-analyses of similar stocks or anecdotal observations. The approach is applied to the population of goliath grouper (Epinephelus itajara) off southern Florida, for which there are three indices of relative abundance but no reliable catch data. The results confirm anecdotal accounts of a marked decline in abundance during the 1980s followed by a substantial increase after the harvest of goliath grouper was banned in 1990. The ban appears to have reduced fishing pressure to between 10% and 50% of the levels observed during the 1980s. Nevertheless, the predicted fishing mortality rate under the ban appears to remain substantial, perhaps owing to illegal harvest and depth-related release mortality. As a result, the base model predicts that there is less than a 40% chance that the spawning biomass will recover to a level that would produce a 50% spawning potential ratio.

The relationship of shell dimensions and shell volume to live weight and soft tissue weight in the mangrove clam, Polymesoda erosa (Solander, 1786) from northern Australia

Shell dimensions (length, height, width) and shell volume were evaluated as estimators of growth for Polymesoda erosa in northern Australia. Each parameter was a good estimator when applied to live weight (r2 values of 76-96 percent), but not to soft tissue weight (wet, dry, or ash-free dry weight) (r2 values of 13-32 percent). The b value for shell volume to weight relationship of clams collected during the dry season (June to October) was signifi cantly different than for those collected in the wet season (February to April).

Using poststratification to improve abundance estimates from multispecies surveys: a study of juvenile flatfishes

Population assessments seldom incorporate habitat information or use previously observed distributions of fish density. Because habitat affects the spatial distribution of fish density and overall abundance, the use of habitat information and previous estimates of fish density can produce more precise and less biased population estimates. In this study, we describe how poststratification can be applied as an unbiased estimator to data sets that were collected under a probability sampling design, typical of many multispecies trawl surveys. With data from a multispecies survey of juvenile flatfish, we show how poststratification can be applied to a data set that was not collected under a probability sampling design, where both the precision and the bias are unknown. For each of four species, three estimates of total abundance were compared: 1) unstratified; 2) poststratified by habitat; and 3) poststratified by habitat and fish density (high fish density and low fish density) in nearby years. Poststratification by habitat gave more precise and (or) less design-biased estimates than an unstratified estimator for all species in all years. Poststratification by habitat and fish density produced the most precise and representative estimates when the sample size in the high fish-density and low fish-density strata were sufficient (in this study, n≥20 in the high fish-density stratum, n≥9 in the low fish-density stratum). Because of the complexities of statistically testing the annual stratified data, we compared three indices of abundance for determining statistically significant changes in annual abundance. Each of the indices closely approximated the annual differences of the poststratified estimates. Selection of the most appropriate index was dependent upon the species’ density distribution within habitat and the sample size in the different habitat areas. The methods used in this study are particularly useful for estimating individual species abundance from multispecies surveys and for retrospective st

Seabed classification for trawlability determined with a multibeam echo sounder on Snakehead Bank in the Gulf of Alaska

Rockfishes (Sebastes spp.) tend to aggregate near rocky, cobble, or generally rugged areas that are difficult to survey with bottom trawls, and evidence indicates that assemblages of rockfish species may differ between areas accessible to trawling and those areas that are not. Consequently, it is important to determine grounds that are trawlable or untrawlable so that the areas where trawl survey results should be applied are accurately identified. To this end, we used multibeam echosounder data to generate metrics that describe the seafloor: backscatter strength at normal and oblique incidence angles, the variation of the angle-dependent backscatter strength within 10° of normal incidence, the scintillation of the acoustic intensity scattered from the seafloor, and the seafloor rugosity. We used these metrics to develop a binary classification scheme to estimate where the seafloor is expected to be trawlable. The multibeam echosounder data were verified through analyses of video and still images collected with a stereo drop camera and a remotely operated vehicle in a study at Snakehead Bank, ~100 km south of Kodiak Island in the Gulf of Alaska. Comparisons of different combinations of metrics derived from the multibeam data indicated that the oblique-incidence backscatter strength was the most accurate estimator of trawlability at Snakehead Bank and that the addition of other metrics provided only marginal improvements. If successful on a wider scale in the Gulf of Alaska, this acoustic remote-sensing technique, or a similar one, could help improve the accuracy of rockfish stock assessments.

A sudden collapse in distribution of Pacific sardine (Sardinops sagax) off southwestern Australia enables an objective re-assessment of biomass estimates

Fishery-independent estimates of spawning biomass (BSP) of the Pacific sardine (Sardinops sagax) on the south and lower west coasts of Western Australia (WA) were obtained periodically between 1991 and 1999 by using the daily egg production method (DEPM). Ichthyoplankton data collected during these surveys, specifically the presence or absence of S. sagax eggs, were used to investigate trends in the spawning area of S. sagax within each of four regions. The expectation was that trends in BSP and spawning area were positively related. With the DEPM model, estimates of BSP will change proportionally with spawning area if all other variables remain constant. The proportion of positive stations (PPS), i.e., stations with nonzero egg counts — an objective estimator of spawning area — was high for all south coast regions during the early 1990s (a period when the estimated BSP was also high) and then decreased after the mid-1990s. There was a decrease in PPS from the mid-1990s to 1999. The particularly low estimates in 1999 followed a severe epidemic mass mortality of S. sagax throughout their range across southern Australia. Deviations from the expected relationship between BSP and PPS were used to identify uncertainty around estimates of BSP. Because estimation of spawning area is subject to less sampling bias than estimation of BSP, the deviation in the relation between the two provides an objective basis for adjusting some estimates of the latter. Such an approach is particularly useful for fisheries management purposes when sampling problems are suspected to be present. The analysis of PPS undertaken from the same set of samples from which the DEPM estimate is derived will help provide information for stock assessments and for the management of purse-seine fisheries.

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.

Applications in adaptive cluster sampling of Gulf of Alaska rockfish

Adaptive cluster sampling (ACS) has been the subject of many publications about sampling aggregated populations. Choosing the criterion value that invokes ACS remains problematic. We address this problem using data from a June 1999 ACS survey for rockfish, specifically for Pacific ocean perch (Sebastes alutus), and for shortraker (S. borealis) and rougheye (S. aleutianus) rockfish combined. Our hypotheses were that ACS would outperform simple random sampling (SRS) for S. alutus and would be more applicable for S. alutus than for S. borealis and S. aleutianus combined because populations of S. alutus are thought to be more aggregated. Three alternatives for choosing a criterion value were investigated. We chose the strategy that yielded the lowest criterion value and simulated the higher criterion values with the data after the survey. Systematic random sampling was conducted across the whole area to determine the lowest criterion value, and then a new systematic random sample was taken with adaptive sampling around each tow that exceeded the fixed criterion value. ACS yielded gains in precision (SE) over SRS. Bootstrapping showed that the distribution of an ACS estimator is approximately normal, whereas the SRS sampling distribution is skewed and bimodal. Simulation showed that a higher criterion value results in substantially less adaptive sampling with little tradeoff in precision. When time-efficiency was examined, ACS quickly added more samples, but sampling edge units caused this efficiency to be lessened, and the gain in efficiency did not measurably affect our conclusions. ACS for S. alutus should be incorporated with a fixed criterion value equal to the top quartile of previously collected survey data. The second hypothesis was confirmed because ACS did not prove to be more effective for S. borealis-S. aleutianus. Overall, our ACS results were not as optimistic as those previously published in the literature, and indicate the need for further study of this sampling method.

Age, growth, mortality, and distribution of pinfish (Lagodon rhomboides) in Tampa Bay and adjacent Gulf of Mexico waters

Life-history dynamics of pinfish (Lagodon rhomboides) were examined from data derived from random station surveys conducted in Tampa Bay and adjacent Gulf of Mexico waters during 1993–97. In addition, patterns in spatial distribution and abundance in Gulf of Mexico waters were investigated. Ages determined from whole otoliths ranged from 0 to 7 years, and von Bertalanffy growth models for males and females were not significantly different. Von Bertalanffy growth model parameters were L∞=219.9 mm SL, k =0.33/yr, and t0 =–1.10 years for all fish combined. High gonadosomatic indices during October–December indicated that some spawning may occur in Tampa Bay. Estimated lengths at 50% maturity were 132 mm SL for males and 131 mm SL for females. Total instantaneous mortality rates derived from the Chapman-Robson estimator ranged from 0.88 to 1.08/yr, and natural mortality was estimated to be 0.78/yr. In Gulf of Mexico waters, pinfish catch rates declined with increasing depth, and most pinfish were caught in <17 m of water. Length distributions showed that pinfish segregate by size with increasing depth.