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.

Prediction of mesophotic coral distributions in the Au‘au Channel, Hawaii

The primary objective of this study was to predict the distribution of mesophotic hard corals in the Au‘au Channel in the Main Hawaiian Islands (MHI). Mesophotic hard corals are light-dependent corals adapted to the low light conditions at approximately 30 to 150 m in depth. Several physical factors potentially influence their spatial distribution, including aragonite saturation, alkalinity, pH, currents, water temperature, hard substrate availability and the availability of light at depth. Mesophotic corals and mesophotic coral ecosystems (MCEs) have increasingly been the subject of scientific study because they are being threatened by a growing number of anthropogenic stressors. They are the focus of this spatial modeling effort because the Hawaiian Islands Humpback Whale National Marine Sanctuary (HIHWNMS) is exploring the expansion of its scope—beyond the protection of the North Pacific Humpback Whale (Megaptera novaeangliae)—to include the conservation and management of these ecosystem components. The present study helps to address this need by examining the distribution of mesophotic corals in the Au‘au Channel region. This area is located between the islands of Maui, Lanai, Molokai and Kahoolawe, and includes parts of the Kealaikahiki, Alalākeiki and Kalohi Channels. It is unique, not only in terms of its geology, but also in terms of its physical oceanography and local weather patterns. Several physical conditions make it an ideal place for mesophotic hard corals, including consistently good water quality and clarity because it is flushed by tidal currents semi-diurnally; it has low amounts of rainfall and sediment run-off from the nearby land; and it is largely protected from seasonally strong wind and wave energy. Combined, these oceanographic and weather conditions create patches of comparatively warm, calm, clear waters that remain relatively stable through time. Freely available Maximum Entropy modeling software (MaxEnt 3.3.3e) was used to create four separate maps of predicted habitat suitability for: (1) all mesophotic hard corals combined, (2) Leptoseris, (3) Montipora and (4) Porites genera. MaxEnt works by analyzing the distribution of environmental variables where species are present, so it can find other areas that meet all of the same environmental constraints. Several steps (Figure 0.1) were required to produce and validate four ensemble predictive models (i.e., models with 10 replicates each). Approximately 2,000 georeferenced records containing information about mesophotic coral occurrence and 34 environmental predictors describing the seafloor’s depth, vertical structure, available light, surface temperature, currents and distance from shoreline at three spatial scales were used to train MaxEnt. Fifty percent of the 1,989 records were randomly chosen and set aside to assess each model replicate’s performance using Receiver Operating Characteristic (ROC), Area Under the Curve (AUC) values. An additional 1,646 records were also randomly chosen and set aside to independently assess the predictive accuracy of the four ensemble models. Suitability thresholds for these models (denoting where corals were predicted to be present/absent) were chosen by finding where the maximum number of correctly predicted presence and absence records intersected on each ROC curve. Permutation importance and jackknife analysis were used to quantify the contribution of each environmental variable to the four ensemble models.

Growth, mortality, and hatchdate distributions of larval and juvenile spotted seatrout (Cynoscion nebulosus) in Florida Bay, Everglades National Park

Life history aspects of larval and, mainly, juvenile spotted seatrout (Cynoscion nebulosus) were studied in Florida Bay, Everglades National Park, Florida. Collections were made in 1994−97, although the majority of juveniles were collected in 1995. The main objective was to obtain life history data to eventually develop a spatially explicit model and provide baseline data to understand how Everglades restoration plans (i.e. increased freshwater flows) could influence spotted seatrout vital rates. Growth of larvae and juveniles (<80 mm SL) was best described by the equation loge standard length = –1.31 + 1.2162 (loge age). Growth in length of juveniles (12–80 mm SL) was best described by the equation standard length = –7.50 + 0.8417 (age). Growth in wet weight of juveniles (15–69 mm SL) was best described by the equation loge wet-weight = –4.44 + 0.0748 (age). There were no significant differences in juvenile growth in length of spotted seatrout in 1995 between three geographical subdivisions of Florida Bay: central, western, and waters adjacent to the Gulf of Mexico. We found a significant difference in wet-weight for one of six cohorts categorized by month of hatchdate in 1995, and a significant difference in length for another cohort. Juveniles (i.e. survivors) used to calculate weekly hatchdate distributions during 1995 had estimated spawning times that were cyclical and protracted, and there was no correlation between spawning and moon phase. Temperature influenced otolith increment widths during certain growth periods in 1995. There was no evidence of a relationship between otolith growth rate and temperature for the first 21 increments. For increments 22–60, otolith growth rates decreased with increasing age and the extent of the decrease depended strongly in a quadratic fashion on the temperature to which the fish was exposed. For temperatures at the lower and higher range, increment growth rates were highest. We suggest that this quadratic relationship might be influenced by an environmental factor other than temperature. There was insufficient information to obtain reliable inferences on the relationship of increment growth rate to salinity.

Assessing the precision of frequency distributions estimated from trawl-survey samples

In trawl surveys a cluster of fish are caught at each station, and fish caught together tend to have more similar characteristics, such as length, age, stomach contents etc., than those in the entire population. When this is the case, the effective sample size for estimates of the frequency distribution of a population characteristic can, therefore, be much smaller than the number of fish sampled during a survey. As examples, it is shown that the effective sample size for estimates of length-frequency distributions generated by trawl surveys conducted in the Barents Sea, off Namibia, and off South Africa is on average approximately one fish per tow. Thus many more fish than necessary are measured at each station (location). One way to increase the effective sample size for these surveys and, hence, increase the precision of the length-frequency estimates, is to reduce tow duration and use the time saved to collect samples at more stations.