Assessment of Atlantic Red Drum for 1999: Northern and southern regions

An assessment of the status of the Atlantic stock of red drum is conducted using recreational and commercial data from 1986 through 1998. This assessment updates data and analyses from the 1989, 1991, 1992 and 1995 stock assessments on Atlantic coast red drum (Vaughan and Helser, 1990; Vaughan 1992; 1993; 1996). Since 1981, coastwide recreational catches ranged between 762,300 pounds in 1980 and 2,623,900 pounds in 1984, while commercial landings ranged between 60,900 pounds in 1997 and 422,500 pounds in 1984. In weight of fish caught, Atlantic red drum constitute predominantly a recreational fishery (ranging between 85 and 95% during the 1990s). Commercially, red drum continue to be harvested as part of mixed species fisheries. Using available length-frequency distributions and age-length keys, recreational and commercial catches are converted to catch in numbers at age. Separable and tuned virtual population analyses are conducted on the catch in numbers at age to obtain estimates of fishing mortality rates and population size (including recruitment to age 1). In tum, these estimates of fishing mortality rates combined with estimates of growth (length and weight), sex ratios, sexual maturity and fecundity are used to estimate yield per recruit, escapement to age 4, and static (or equilibrium) spawning potential ratio (static SPR, based on both female biomass and egg production). Three virtual analysis approaches (separable, spreadsheet, and FADAPT) were applied to catch matrices for two time periods (early: 1986-1991, and late: 1992-1998) and two regions (Northern: North Carolina and north, and Southern: South Carolina through east coast of Florida). Additional catch matrices were developed based on different treatments for the catch-and-release recreationally-caught red drum (B2-type). These approaches included assuming 0% mortality (BASEO) versus 10% mortality for B2 fish. For the 10% mortality on B2 fish, sizes were assumed the same as caught fish (BASEl), or positive difference in size distribution between the early period and the later period (DELTA), or intermediate (PROP). Hence, a total of 8 catch matrices were developed (2 regions, and 4 B2 assumptions for 1986-1998) to which the three VPA approaches were applied. The question of when offshore emigration or reduced availability begins (during or after age 3) continues to be a source of bias that tends to result in overestimates of fishing mortality. Additionally, the continued assumption (Vaughan and Helser, 1990; Vaughan 1992; 1993; 1996) of no fishing mortality on adults (ages 6 and older), causes a bias that results in underestimates of fishing mortality for adult ages (0 versus some positive value). Because of emigration and the effect of the slot limit for the later period, a range in relative exploitations of age 3 to age 2 red drum was considered. Tuning indices were developed from the MRFSS, and state indices for use in the spreadsheet and FADAPT VPAs. The SAFMC Red Drum Assessment Group (Appendix A) favored the FADAPT approach with catch matrix based on DELTA and a selectivity for age 3 relative to age 2 of 0.70 for the northern region and 0.87 for the southern region. In the northern region, estimates of static SPR increased from about 1.3% for the period 1987-1991 to approximately 18% (15% and 20%) for the period 1992-1998. For the southern region, estimates of static SPR increased from about 0.5% for the period 1988-1991 to approximately 15% for the period 1992-1998. Population models used in this assessment (specifically yield per recruit and static spawning potential ratio) are based on equilibrium assumptions: because no direct estimates are available as to the current status of the adult stock, model results imply potential longer term, equilibrium effects. Because current status of the adult stock is unknown, a specific rebuilding schedule cannot be determined. However, the duration of a rebuilding schedule should reflect, in part, a measure of the generation time of the fish species under consideration. For a long-lived, but relatively early spawning, species as red drum, mean generation time would be on the order of 15 to 20 years based on age-specific egg production. Maximum age is 50 to 60 years for the northern region, and about 40 years for the southern region. The ASMFC Red Drum Board's first phase recovery goal of increasing %SPR to at least 10% appears to have been met. (PDF contains 79 pages)

Hydroacoustic surveys: a non-destructive approach to monitoring fish distributions at National Marine Sanctuaries

The science of fisheries acoustics and its applicability to resource management have evolved over the past several decades. This document provides a basic description of fisheries acoustics and recommendations on using this technology for research and monitoring of fish distributions and habitats within sanctuaries. It also describes recent efforts aimed at applying fisheries acoustics to Gray’s Reef National Marine Sanctuary (GRNMS) (Figure 1). Historically, methods to assess the underwater environment have included net trawls, diver censuses, hook and line, video, sonar and other techniques deployed in a variety of ways. Fisheries acoustics, using active sonar, relies on the physics of sound traveling through water to quantify the distribution of biota in the water column. By sending a signal of a given frequency through the water column and recording the time of travel and the strength of the reflected signal, it is possible to determine the size and location of fish and estimate biomass from the acoustic backscatter. As a fisheries assessment tool, active hydroacoustics technology is an efficient, non-intrusive method of mapping the water column at a very fine spatial and temporal resolution. It provides a practical alternative to bottom and mid-water trawls, which are not allowed at GRNMS. Passive acoustics, which uses underwater hydrophones to record man-made and natural sounds such as fish spawning calls and sounds produced by marine mammals for communication and echolocation, can provide a useful, complementary survey tool. This report primarily deals with active acoustics, although the integration of active and passive acoustics is addressed as well. (PDF contains 32 pages)

Interpreting spotted dolphin age distributions

Previous work has determined the age distribution from a sample of spotted dolphins (Stenella attenuata) killed in the eastern Pacific tuna purse-seine fishery. In this paper we examine the usefulness of this age distribution for estimating natural mortality rates. The observed age distribution has a deficiency of individuals from 5-15 years and cannot represent a stable age distribution. Sampling bias and errors in age interpretation are examined as possible causes of the "dip" in the observed age structure. Natural mortality rates are estimated for the 15+ age classes based on the assumption that these are sampled representatively. The resulting annual survival rate

Some aspects of the temperature, oxygen and nutrient distributions in Monterey Bay, California. Annual Report, Part 1, 1973

Seasonal variations in temperature, dissolved oxygen, and nutrients in the nearshore areas and in the canyon area of Monterey Bay, California during 1971-1972 were similar~ During upwelling periods, however, water in the nearshore areas was higher in temperature and oxygen and lower in nutrients than water in the canyon area~ This was caused by upwelled water moving north and south of the canyon into counterclockwise and clockwise flow in the northern and southern ends of the bay respectively. The water was heated by insolation and depleted of its nutrients by photosynthesis during this movement. The residence time of water in the nearshore northern and southern bay during upwelling is estimated to be 3 to 8 days, and this fits well into the above circulation pattern and average measured current velocities of 10 to 15 cm/sec~ There is sorne evidence that this circulation pattern and the estimated residence time may be also valid for on-upwelling periods. Upwelling apparently occurred in Monterey Submarine Canyon at rates of 0.4 to 2.9 m/day and was stronger in 1971 than 1972. (PDF contains 107 pages)

Nutrient distributions for Georges Bank and adjacent waters in 1979

In this report we describe the temporal and spatial distributions of inorganic nutrients over Georges Bank and in adjacent waters and discuss major features with respect to tbe nutrient environments of pbytoplankton. Nitrate and orthophosphorus were rapidly depleted from the surface layer of much of the study area in spring, but major differences were found between the shallow areas on Georges Bank and the surrounding stratified waters. In the "well-mixed" area of Georges Bank, the depletion encompassed the entire water column and ammonium became the dominant form of inorganic nitrogen throughout. Dissolved silicon was depleted slowly over central Georges Bank, reaching a minimum concentration in September while orthophosphorus gradually increased during the summer. The nutrient environment of phytoplankton over central Georges Bank may be described as vertically uniform but temporally changing in the relative availability of the various nutrients. In areas that undergo stratification (e.g., the central Gulf of Maine), a quasi-steady state was established as the surface water layer formed, consisting of declining nutrient gradients from below the euphotic layer to the top of the water column. These intergrading nutrient environments are relatively stable through time. Destratification reintroduced nutrients to depleted areas beginning in October; however, dissolved silicon was again depleted over shallow Georges Bank in late autumn though nitrate remained abundant. Slope water has been found to enter the bottom layer of the Gulf of Maine via the Northeast Channel. High nutrient concentrations observed in the bottom water of the Northeast Channel are consistent with this mechanism being the nutrient source for the Gulf of Maine. (PDF file contains 40 pages.)

Geographical distributions of effort and catches of tunas by purse-seine vessels in the Eastern Pacific Ocean during 1965-1998

ENGLISH: This report presents fine-scale spatial summaries of annual catch and effort information compiled by the IATTC staff. These data summaries are presented in a graphical format, and display only information collected from purse-seine vessels fishing in the eastern Pacific Ocean (EPO) during 1965-1998. Data collected from baitboat and longline vessels fishing in the EPO are not considered in this report. Equivalent data from Japanese longline vessels fishing in the EPO are presented by Uosaki and Bayliff(1999) for 1988-1992 and by' references cited in that publication for 1956-1987. SPANISH: Este informe presenta resúmenes espaciales a escala fina de información sobre captura y esfuerzo anuales compilada por el personal de la CIAT. Se presentan los resúmenes en formato gráfico, e incluyen solamente información tomada de barcos cerqueros pescando en el Océano Pacifico oriental (OPO) durante 1965-1998. No se consideran en el informe datos provenientes de barcos de camada y palangreros que pescan en el OPO. Se presentan datos equivalentes de barcos palangreros japoneses pescando en el OPO en Uosaki y Bayliff (1999) para 1988-1992 yen las referencias citadas en dicha publicación para 1956-1987. (PDF contains 104 pages.)

Contribution of haematological factors to the distributions and estimation of Eustrongyliides africanus larvae densities in Clarias gariepinus and Clarias anguillaris from Bida floodplain of Nigeria

The contributions of hematological factors to the distribution and estimations of Eustrongylides africanus larvae densities in Clarias gariepinus and C. anguillaris of Bida floodplain of Nigeria were documented for the first time. The hematological factors making the most important contributions to the distributions of E. africanus larvae infections in clarias species are mean corpuscular haemoglobin concentration (MCHC), mean corpuscular haemoglobin (MCH), mean corpuscular volume (MCV) and neutrophils count, in descending order of magnitude; having the manifestations for the months of January, March, September, and December of the year being closely related. Five haematological factors (neutrophils, lymphocytes and eosinophils counts; MCH and MCV) having positive or negative correlation coefficient (r) between 0.50 and 0.85 contributed to the estimated of E.africanus larvae densities in the wild population of Clarias species

Mathematical models for the prediction of temperature distributions resulting from the discharge of heated water into large bodies of water

Mathematical models for heated water outfalls were developed for three flow regions. Near the source, the subsurface discharge into a stratified ambient water issuing from a row of buoyant jets was solved with the jet interference included in the analysis. The analysis of the flow zone close to and at intermediate distances from a surface buoyant jet was developed for the two-dimensional and axisymmetric cases. Far away from the source, a passive dispersion model was solved for a two dimensional situation taking into consideration the effects of shear current and vertical changes in diffusivity. A significant result from the surface buoyant jet analysis is the ability to predict the onset and location of an internal hydraulic jump. Prediction can be made simply from the knowledge of the source Froude number and a dimensionless surface exchange coefficient. Parametric computer programs of the above models are also developed as a part of this study. This report was submitted in fulfillment of Contract No. 14-12-570 under the sponsorship of the Federal Water Quality Administration.

Distributions and abundances of Pacific sardine (Sardinops sagax) and other pelagic fishes in the California Current Ecosystem during spring 2006, 2008, and 2010, estimated from acoustic–trawl surveys

The abundances and distributions of coastal pelagic fish species in the California Current Ecosystem from San Diego to southern Vancouver Island, were estimated from combined acoustic and trawl surveys conducted in the spring of 2006, 2008, and 2010. Pacific sardine (Sardinops sagax), jack mackerel (Trachurus symmetricus), and Pacific mackerel (Scomber japonicus) were the dominant coastal pelagic fish species, in that order. Northern anchovy (Engraulis mordax) and Pacific herring (Clupea pallasii) were sampled only sporadically and therefore estimates for these species were unreliable. The estimates of sardine biomass compared well with those of the annual assessments and confirmed a declining trajectory of the “northern stock” since 2006. During the sampling period, the biomass of jack mackerel was stable or increasing, and that of Pacific mackerel was low and variable. The uncertainties in these estimates are mostly the result of spatial patchiness which increased from sardine to mackerels to anchovy and herring. Future surveys of coastal pelagic fish species in the California Current Ecosystem should benefit from adaptive sampling based on modeled habitat; increased echosounder and trawl sampling, particularly for the most patchy and nearshore species; and directed-trawl sampling for improved species identification and estimations of their acoustic target stren

Population structure, long-term connectivity, and effective size of mutton snapper (Lutjanus analis) in the Caribbean Sea and Florida Keys

Genetic structure and average long-term connectivity and effective size of mutton snapper (Lutjanus analis) sampled from offshore localities in the U.S. Caribbean and the Florida Keys were assessed by using nuclear-encoded microsatellites and a fragment of mitochondrial DNA. No significant differences in allele, genotype (microsatellites), or haplotype (mtDNA) distributions were detected; tests of selective neutrality (mtDNA) were nonsignificant after Bonferroni correction. Heuristic estimates of average long-term rate of migration (proportion of migrant individuals/generation) between geographically adjacent localities varied from 0.0033 to 0.0054, indicating that local subpopulations could respond independently of environmental perturbations. Estimates of average longterm effective population sizes varied from 341 to 1066 and differed significantly among several of the localities. These results indicate that over time larval drift and interregional adult movement may not be sufficient to maintain population sustainability across the region and that there may be different demographic stocks at some of the localities studied. The estimate of long-term effective population size at the locality offshore of St. Croix was below the minimum threshold size considered necessary to maintain the equilibrium between the loss of adaptive genetic variance from genetic drift and its replacement by mutation. Genetic variability in mutton snapper likely is maintained at the intraregional level by aggregate spawning and random mating of local populations. This feature is perhaps ironic in that aggregate spawning also renders mutton snapper especially vulnerable to overexploitation.

Parameterizing probabilities for estimating age-composition distributions for mixture models

When estimating parameters that constitute a discrete probability distribution {pj}, it is difficult to determine how constraints should be made to guarantee that the estimated parameters { pˆj} constitute a probability distribution (i.e., pˆj>0, Σ pˆj =1). For age distributions estimated from mixtures of length-at-age distributions, the EM (expectationmaximization) algorithm (Hasselblad, 1966; Hoenig and Heisey, 1987; Kimura and Chikuni, 1987), restricted least squares (Clark, 1981), and weak quasisolutions (Troynikov, 2004) have all been used. Each of these methods appears to guarantee that the estimated distribution will be a true probability distribution with all categories greater than or equal to zero and with individual probabilities that sum to one. In addition, all these methods appear to provide a theoretical basis for solutions that will be either maximum-likelihood estimates or at least convergent to a probability distribut

Comparison of average larval fish vertical distributions among species exhibiting different transport pathways on the southeast United States continental shelf

Water currents are vertically structured in many marine systems and as a result, vertical movements by fish larvae and zooplankton affect horizontal transport (Power, 1984). In estuaries, the vertical movements of larvae with tidal periods can result in their retention or ingress (Fortier and Leggett, 1983; Rijnsdorp et al., 1985; Cronin and Forward, 1986; Forward et al., 1999). On the continental shelf, the vertical movements of organisms interact daily and ontogenetically with depth-varying currents to affect horizontal transport (Pillar et al., 1989; Barange and Pillar, 1992; Cowen et al., 1993, 2000; Batchelder et al., 2002).

Year-class formation in Pacific herring (Clupea pallasi) estimated from spawning-date distributions of juveniles in San Francisco Bay, California

Inter and intra-annual variation in year-class strength was analyzed for San Francisco Bay Pacific herring (Clupea pallasi) by using otoliths of juveniles. Juvenile herring were collected from March through June in 1999 and 2000 and otoliths from subsamples of these collections were aged by daily otolith increment analysis. The composition of the year classes in 1999 and 2000 were determined by back-calculating the birth date distribution for surviving juvenile herring. In 2000, 729% more juveniles were captured than in 1999, even though an estimated 12% fewer eggs were spawned in 2000. Spawning-date distributions show that survival for the 2000 year class was exceptionally good for a short (approximately 1 month) period of spawning, resulting in a large abundance of juvenile recruits. Analysis of age at size shows that growth rate increased significantly as the spawning season progressed both in 1999 and 2000. However, only in 2000 were the bulk of surviving juveniles a product of the fast growth period. In the two years examined, year-class strength was not predicted by the estimated number of eggs spawned, but rather appeared to depend on survival of eggs or larvae (or both) through the juvenile stage. Fast growth through the larval stage may have little effect on year-class strength if mortality during the egg stage is high and few larvae are available.

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.