Variations in eastern North Pacific demersal fish biomass based on the U.S. west coast groundfish bottom trawl survey (2003–2010)

In response to declining biomass of Northeast Pacific groundfish in the late 1990s and to improve the scientific basis for management of the fishery, the Northwest Fisheries Science Center standardized and enhanced their annual bottom trawl survey in 2003. The survey was expanded to include the entire area along the U.S. west coast at depths of 55–1280 m. Coast-wide biomass and species richness significantly decreased during the first eight years (2003–10) of this fishery-independent survey. We observed an overall tendency toward declining biomass for 62 dominant taxa combined (fishery target and nontarget species) and four of seven subgroups (including cartilaginous fish, flatfishes, shelf rockfishes, and other shelf species), despite increasing or variable biomass trends in individual species. These decreases occurred during a period of reduced catch for groundfish along the shelf and upper slope regions relative to historical rates. We used information from multiple stock assessments to aggregate species into three groups: 1) with strong recruitment, 2) without strong recruitment in 1999, and 3) with unknown recruitment level. For each group, we evaluated whether declining biomass was primarily related to depletion (using year as a proxy) or environmental factors (i.e., variation in the Pacific Decadal Oscillation). According to Akaike’s information criterion, changes in aggregate biomass for species with strong recruitment were more closely related to year, whereas those with no strong recruitment were more closely related to climate. The significant decline in biomass for species without strong recruitment confirms that factors other than depletion of the exceptional 1999 year class may be responsible for the observed decrease in biomass along the U.S. west coast.

Length and sex effects on the spatial structure of catches of Pacific halibut (Hippoglossus stenolepis) on longline gear

Field experiments were conducted to test the hypotheses that Pacific halibut (Hippoglossus stenolepis) display small-scale spatial structure within longline catches, relative to other species and empty hooks, or within-species based on sex or length. Sequential hook-by-hook inventories, along with length and sex data, were taken at thirty-one survey stations. Two-dimensional spatial statistics were used to test for 1) aggregation, defined as the clustering of individuals within a given demographic of size or sex over small intervals of distance; and 2) segregation, defined as the sequential occurrence of individuals within a given demographic of size or sex, uninterrupted by other observations, irrespective of the distance between individuals. Statistically significant structure was detected within catches that is more commonly associated with fish length than sex. Significant spatial structuring occurred at 60% of all stations tested. Significant aggregation of halibut of legal length for commercial retention (≥82 cm) was detected at 44% of stations and aggregation of sublegal-size halibut was detected at 11%. Maleand female-based aggregations were observed at 22% and 11% of stations, respectively. Significant segregation of females was observed at 20% of stations, male segregation occurred at 8% of stations, and segregation by size at 16% of stations. Understanding small-scale spatial structure within longline catches may help us interpret changes in survey and commercial catch data. If structure is generated by behavior, then observed size-at-age or relative sex-ratios may be biased relative to underlying distributions. Although physical processes such as gape limitation should remain stable over the time, dynamic processes may be spatially and temporally variabl

Prediction and confirmation of seasonal migration of Pacific sardine (Sardinops sagax) in the California Current Ecosystem

During the last century, the population of Pacific sardine (Sardinops sagax) in the California Current Ecosystem has exhibited large fluctuations in abundance and migration behavior. From approximately 1900 to 1940, the abundance of sardine reached 3.6 million metric tons and the “northern stock” migrated from offshore of California in the spring to the coastal areas near Oregon, Washington, and Vancouver Island in the summer. In the 1940s, the sardine stock collapsed and the few remaining sardine schools concentrated in the coastal region off southern California, year-round, for the next 50 years. The stock gradually recovered in the late 1980s and resumed its seasonal migration between regions off southern California and Canada. Recently, a model was developed which predicts the potential habitat for the northern stock of Pacific sardine and its seasonal dynamics. The habitat predictions were successfully validated using data from sardine surveys using the daily egg production method; scientific trawl surveys off the Columbia River mouth; and commercial sardine landings off Oregon, Washington, and Vancouver Island. Here, the predictions of the potential habitat and seasonal migration of the northern stock of sardine are validated using data from “acoustic–trawl” surveys of the entire west coast of the United States during the spring and summer of 2008. The estimates of sardine biomass and lengths from the two surveys are not significantly different between spring and summer, indicating that they are representative of the entire stock. The results also confirm that the model of potential sardine habitat can be used to optimally apply survey effort and thus minimize random and systematic sampling error in the biomass estimates. Furthermore, the acoustic–trawl survey data are useful to estimate concurrently the distributions and abundances of other pelagic fishes.

Effect of environmental conditions on the distribution of Pacific mackerel (Scomber japonicus) larvae in the California Current system

We modeled the probability of capturing Pacif ic mackerel (Scomber japonicus) larvae as a function of environmental variables for the Southern California Bight (SCB) most years from 1951 through 2008 and Mexican waters offshore of Baja California from 1951 through 1984. The model exhibited acceptable fit, as indicated by the area under a receiver-operating-characteristic curve of 0.80 but was inconsistent with the zero catches that occurred frequently in the 2000s. Two types of spawners overlapped spatially within the survey area: those that exhibited peak spawning during April in the SCB at about 15.5°C and a smaller group that exhibited peak spawning in August near Punta Eugenia, Mexico, at 20°C or greater. The SCB generally had greater zooplankton than Mexican waters but less appropriate (lower) geostrophic f lows. Mexican waters generally exhibited greater predicted habitat quality than the SCB in cold years. Predicted quality of the habitat in the SCB was greater from the 1980s to 2008 than in the earlier years of the survey primarily because temperatures and geostrophic flows were more appropriate for larvae. However, stock size the previous year had a larger effect on predictions than any environmental variable, indicating that larval Pacific mackerel did not fully occupy the suitable habitat during most years.

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

Postrelease survival, vertical and horizontal movements, and thermal habitats of five species of pelagic sharks in the central Pacific Ocean

From 2001 to 2006, 71 pop-up satellite archival tags (PSATs) were deployed on five species of pelagic shark (blue shark [Prionace glauca]; shortfin mako [Isurus oxyrinchus]; silky shark [Carcharhinus falciformis]; oceanic whitetip shark [C. longimanus]; and bigeye thresher [Alopias superciliosus]) in the central Pacific Ocean to determine species-specific movement patterns and survival rates after release from longline fishing gear. Only a single postrelease mortality could be unequivocally documented: a male blue shark which succumbed seven days after release. Meta-analysis of published reports and the current study (n=78 reporting PSATs) indicated that the summary effect of postrelease mortality for blue sharks was 15% (95% CI, 8.5–25.1%) and suggested that catch-and-release in longline fisheries can be a viable management tool to protect parental biomass in shark populations. Pelagic sharks displayed species-specific depth and temperature ranges, although with significant individual temporal and spatial variability in vertical movement patterns, which were also punctuated by stochastic events (e.g., El Niño-Southern Oscillation). Pelagic species can be separated into three broad groups based on daytime temperature preferences by using the unweighted pair-group method with arithmetic averaging clustering on a Kolmogorov-Smirnov Dmax distance matrix: 1) epipelagic species (silky and oceanic whitetip sharks), which spent >95% of their time at temperatures within 2°C of sea surface temperature; 2) mesopelagic-I species (blue sharks and shortfin makos, which spent 95% of their time at temperatures from 9.7° to 26.9°C and from 9.4° to 25.0°C, respectively; and 3) mesopelagic-II species (bigeye threshers), which spent 95% of their time at temperatures from 6.7° to 21.2°C. Distinct thermal niche partitioning based on body size and latitude was also evident within epipelagic species.

Genetic variation between outer-coastal and fjord populations of Pacific herring (Clupea pallasii) in the eastern Gulf of Alaska

Pacific herring (Clupea pallasii) from the Gulf of Alaska were screened for temporal and spatial genetic variation with 15 microsatellite loci. Thirteen collections were examined in this study: 11 from Southeast Alaska and 2 from Prince William Sound, Alaska. Although FST values were low, a neighbor-joining tree based on genetic distance, homogeneity, and FST values revealed that collectively, the Berners Bay and Lynn Canal (interior) collections were genetically distinct from Sitka Sound and Prince of Wales Island (outer-coastal) collections. Temporal genetic variation within regions (among three years of Berners Bay spawners and between the two Sitka Sound spawners) was zero, whereas 0.05% was attributable to genetic variation between Berners Bay and Sitka Sound. This divergence may be attributable to environmental differences between interior archipelago waters and outer-coast habitats, such as differences in temperature and salinity. Early spring collections of nonspawning Lynn Canal herring were nearly genetically identical to collections of spawning herring in Berners Bay two months later—an indication that Berners Bay spawners over-winter in Lynn Canal. Southeast Alaskan herring (collectively) were significantly different from those in Prince William Sound. This study illustrates that adequate sample size is needed to detect variation in pelagic fish species with a large effective population size, and microsatellite markers may be useful in detecting low-level genetic divergence in Pacific herring in the Gulf of Alaska.

Otolith morphometrics and population structure of Pacific sardine (Sardinops sagax) along the west coast of North America

The broad distribution of Pacific sardine (Sardinops sagax) along the Pacif ic coast of North America makes it difficult for fisheries managers to identify regional stocks of this dominant small pelagic species. An investigation of morphometric characteristics of otoliths of Pacific sardine across most of their range revealed regional differences in populations. In a survey of over 2000 otoliths, all ages (with an emphasis on age-1 recruits) were compared. Principal components analysis, multivariate analysis of variance, and a novel method derived from regression and residuals calculations, termed perimeter-weight profiles (PWPs), revealed otolith similarities and differences. The results of the different approaches to statistical comparisons did not always agree. Sardine otoliths from Mexican waters were generally lighter and more lobate than those from U.S. and Canadian populations. Age-1 otoliths from northern California in 2006–07 tended to be heavier and smoother than those from other areas, including year-class cohorts from southern California. Comparisons of age-groups and year-classes of northern California otoliths with the use of the PWP models indicated signif icant trends in year-to-year patterns. In conjunction with other established indices of population structure, otolith PWPs are a useful tool for identifying local and regional stocks of Pacific sardine and may help distinguish populations of other fish species as well.

Ocean distribution of the American shad (Alosa sapidissima) along the Pacific coast of North America

We examined the incidental catches of American shad (Alosa sapidissima) taken during research cruises and in commercial and recreational landings along the Pacific coast of North America during over 30 years of sampling. Shad, an introduced species, was mainly found over the shallow continental shelf, and largest catches and highest frequency of occurrences were found north of central Oregon, along the coasts of Washington and Vancouver Island, and in California around San Francisco Bay. Migrations to the north off Washington and Vancouver were seen during spring to fall, but we found no evidence for large-scale seasonal migrations to the south during the fall or winter. The average weight of shad increased in deeper water. Sizes were also larger in early years of the study. Most were caught over a wide range of sea surface temperatures (11–17°C) and bottom temperatures (6.4–8.0°C). Abundance of shad on the continental shelf north of 44°N was highly correlated with counts of shad at Bonneville Dam on the Columbia River in the same year. Counts were negatively related to average weights and also negatively correlated with the survival of hatchery coho salmon (Oncorhynchus kisutch), indicating that survival of shad is favored by warm ocean conditions. Examining the catch during research cruises and commercial and recreational landings, we concluded that American shad along the Pacific coast have adapted to the prevailing environmental conditions and undertake only moderate seasonal migrations compared with the long seasonal migrations of shad along the Atlantic coast of North America. We suggest that the large spawning populations in the Columbia River and San Francisco Bay areas explain most of the distributional features along the Pacific coast.

Ontogenetic patterns and temperature-dependent growth rates in early life stages of Pacific cod (Gadus macrocephalus)

Pacific cod (Gadus macrocephalus) is an important component of fisheries and food webs in the North Pacific Ocean and Bering Sea. However, vital rates of early life stages of this species have yet to be described in detail. We determined the thermal sensitivity of growth rates of embryos, preflexion and postflexion larvae, and postsettlement juveniles. Growth rates (length and mass) at each ontogenetic stage were measured in three replicate tanks at four to five temperatures. Nonlinear regression was used to obtain parameters for independent stage-specific growth functions and a unified size- and temperature-dependent growth function. Specific growth rates increased with temperature at all stages and generally decreased with increases in body size. However, these analyses revealed a departure from a strict size-based allometry in growth patterns, as reduced growth rates were observed among preflexion larvae: the reduction in specific growth rate between embryos and free-swimming larvae was greater than expected based on body size differences. Growth reductions in the preflexion larvae appear to be associated with increased metabolic rates and the transition from endogenous to exogenous feeding. In future studies, experiments should be integrated across life transitions to more clearly define intrinsic ontogenetic and size-dependent growth patterns because these are critical for evaluations of spatial and temporal variation in habitat quality.

Feeding ecology of juvenile Pacific salmon (Oncorhynchus spp.) in a northeast Pacific fjord: diet, availability of zooplankton, selectivity for prey, and potential competition for prey resources

We investigated the feeding ecology of juvenile salmon during the critical early life-history stage of transition from shallow to deep marine waters by sampling two stations (190 m and 60 m deep) in a northeast Pacific fjord (Dabob Bay, WA) between May 1985 and October 1987. Four species of Pacific salmon—Oncorhynchus keta (chum) , O. tshawytscha (Chinook), O. gorbuscha (pink), and O. kisutch (coho)—were examined for stomach contents. Diets of these fishes varied temporally, spatially, and between species, but were dominated by insects, euphausiids, and decapod larvae. Zooplankton assemblages and dry weights differed between stations, and less so between years. Salmon often demonstrated strongly positive or negative selection for specific prey types: copepods were far more abundant in the zooplankton than in the diet, whereas Insecta, Araneae, Cephalapoda, Teleostei, and Ctenophora were more abundant in the diet than in the plankton. Overall diet overlap was highest for Chinook and coho salmon (mean=77.9%)—species that seldom were found together. Chum and Chinook salmon were found together the most frequently, but diet overlap was lower (38.8%) and zooplankton biomass was not correlated with their gut fullness (%body weight). Thus, despite occasional occurrences of significant diet overlap between salmon species, our results indicate that interspecific competition among juvenile salmon does not occur in Dabob Bay.

Using experiments and expert judgment to model catchability of Pacific rockfishes in trawl surveys, with application to bocaccio (Sebastes paucispinis) off British Columbia

The time series of abundance indices for many groundfish populations, as determined from trawl surveys, are often imprecise and short, causing stock assessment estimates of abundance to be imprecise. To improve precision, prior probability distributions (priors) have been developed for parameters in stock assessment models by using meta-analysis, expert judgment on catchability, and empirically based modeling. This article presents a synthetic approach for formulating priors for rockfish trawl survey catchability (qgross). A multivariate prior for qgross for different surveys is formulated by using 1) a correction factor for bias in estimating fish density between trawlable and untrawlable areas, 2) expert judgment on trawl net catchability, 3) observations from trawl survey experiments, and 4) data on the fraction of population biomass in each of the areas surveyed. The method is illustrated by using bocaccio (Sebastes paucipinis) in British Columbia. Results indicate that expert judgment can be updated markedly by observing the catch-rate ratio from different trawl gears in the same areas. The marginal priors for qgross are consistent with empirical estimates obtained by fitting a stock assessment model to the survey data under a noninformative prior for qgross. Despite high prior uncertainty (prior coefficients of variation ≥0.8) and high prior correlation between qgross, the prior for qgross still enhances the precision of key stock assessment quantities.

Biomass and reproduction of Pacific sardine (Sardinops sagax) off the Pacific northwestern United States, 2003–2005

The Pacific sardine (Sardinops sagax) is distributed along the west coast of North America from Baja California to British Columbia. This article presents estimates of biomass, spawning biomass, and related biological parameters based on four trawl-ichthyoplankton surveys conducted during July 2003 –March 2005 off Oregon and Washington. The trawl-based biomass estimates, serving as relative abundance, were 198,600 t (coefficient of variation [CV] = 0.51) in July 2003, 20,100 t (0.8) in March 2004, 77,900 t (0.34) in July 2004, and 30,100 t (0.72) in March 2005 over an area close to 200,000 km2. The biomass estimates, high in July and low in March, are a strong indication of migration in and out of this area. Sardine spawn in July off the Pacific Northwest (PNW) coast and none of the sampled fish had spawned in March. The estimated spawning biomass for July 2003 and July 2004 was 39,184 t (0.57) and 84,120 t (0.93), respectively. The average active female sardine in the PNW spawned every 20–40 days compared to every 6–8 days off California. The spawning habitat was located in the southeastern area off the PNW coast, a shift from the northwest area off the PNW coast in the 1990s. Egg production in off the PNW for 2003–04 was lower than that off California and that in the 1990s. Because the biomass of Pacific sardine off the PNW appears to be supported heavily by migratory fish from California, the sustainability of the local PNW population relies on the stability of the population off California, and on local oceanographic conditions for local residence.

Accuracy of sex determination for northeastern Pacific Ocean thornyheads (Sebastolobus altivelis and S. alascanus)

Determining the sex of thornyheads (Sebastolobus alascanus and S. altivelis) can be difficult under field conditions. We assessed our ability to correctly assign sex in the field by comparing results from field observations to results obtained in the laboratory through both macroscopic and microscopic examination of gonads. Sex of longspine thornyheads was more difficult to determine than that of shortspine thornyheads and correct determination of sex was signif icantly related to size. By restricting the minimum size of thornyheads to 18 cm for macroscopic determination of sex we reduced the number of fish with misidentified sex by approximately 65%.

A modeling approach to identify optimal habitat and suitable fishing grounds for neon flying squid (Ommastrephes bartramii) in the Northwest Pacific Ocean

We developed a habitat suitability index (HSI) model to understand and identify the optimal habitat and potential fishing grounds for neon f lying squid (Ommastrephes bartramii) in the Northwest Pacific Ocean. Remote sensing data, including sea surface temperature, sea surface salinity, sea surface height, and chlorophyll-a concentrations, as well as fishery data from Chinese mainland squid f leets in the main fishing ground (150–165°E longitude) from August to October, from 1999 to 2004, were used. The HSI model was validated by using fishery data from 2005. The arithmetic mean modeling with three of the environmental variables—sea surface temperature, sea surface height anomaly, and chlorophyll- a concentrations—was defined as the most parsimonious HSI model. In 2005, monthly HSI values >0.6 coincided with productive fishing grounds and high fishing effort from August to October. This result implies that the model can reliably predict potential f ishing grounds for O. bartramii. Because spatially explicit fisheries and environmental data are becoming readily available, it is feasible to develop a dynamic, near real-time habitat model for improving the process of identifying potential fishing areas for and optimal habitats of neon flying squid.