Simulation of physically mediated variability in prey resources of a larval fish: a three-dimensional NPZ model

A three-dimensional biophysical nutrient–phytoplankton–zooplankton model was used to investigate the spatial and temporal dynamics of food resources for young walleye pollock in the western Gulf of Alaska, to further understanding of recruitment processes for pollock. We modeled nitrogen, phytoplankton, a large herbivorous grazer parameterized as Neocalanus spp. (the biomass dominant copepod in the Gulf), and the 13 stages (egg, naupliar and copepodite) of Pseudocalanus spp. (a major constituent of the diet of pollock) so that the appropriate size class of food for each size of larval pollock was represented. Model results identified an area between the Semidi and Shumagin Islands that may not be suitable as a nursery area early in the year due to low prey abundance. Modeled mesoscale eddies, previously hypothesized to be important for larval pollock retention in Shelikof Strait, contained higher prey concentrations than the surrounding waters when they were cyclonic. This work also help to understand the consistency of pollock spawning in time and space in Shelikof Strait, by examining the timing and location of prey availability which, along with transport, narrows the window for optimal spawning.     

Anomalous transport of walleye pollock larvae linked to ocean and atmospheric patterns in May 1996

Larvae from a large aggregation of walleye pollock spawning in early spring in Shelikof Strait, Gulf of Alaska, are normally transported to the south-west in the vigorous Alaska Coastal Current. In the spring of 1996, anomalous winds resulted in unusually weak transport in the Shelikof Strait sea valley. The main aggregation of larval pollock in the Shelikof region was surveyed four times in 1996 over a period of about 40 days, including finer-scale sampling of the leading south-western edge of the larval distribution. The south-western edge of the larval distribution showed weak transport up the sea valley for a period of about 10 days, corresponding to the observations of currents, after which many larvae were transported over the shelf region to the west. These observations are unique in over 15 years of monitoring larval transport patterns and demonstrate how anomalous weather, and hence current patterns, influence variability in larval transport.     

Observed patches of walleye pollock eggs and larvae in Shelikof Strait, Alaska: their characteristics, formation and persistence

Using observations from 38 ichthyoplankton surveys conducted near Shelikof Strait, Alaska between 1979 and 1992, we characterized the horizontal distribution and spatial patchiness of the early life stages of walleye pollock ( Theragra chalcogtamma ). Lloyd's index of patchiness ranged from 3.9-6.1 for eggs and 3.9–16.2 for larvae. This index was size (age) dependent: low for eggs, high for newly hatched larvae, then decreasing through late larval stage. By the early juvenile stage, patchiness increased as pollock began to school. The percentage of larvae in a patch (defined as the percentage of larvae present at stations where larval counts exceeded the mean by one standard deviation during the given survey) varied greatly (26–92%). Larval distributions were used to deduce physical mechanisms responsible for patches. Three categories of patches were identified: those created by interaction of larvae with time-dependent currents, those in the vicinity of Sutwik Island, and those associated with eddies. Simulation experiments were utilized to examine processes influencing patch formation and the role of larval swimming. Between 5 and 6 weeks after hatching, larvae have swimming abilities that enable them to maintain a patch already created by physical mechanisms.     

Contributions of FOCI research to forecasts of year-class strength of walleye pollock in Shelikof Strait, Alaska

NOAA's Fisheries Oceanography Coordinated Investigations (FOCI) contributes information to help forecast year-class strength of walleye pollock (Theragra chalcogramma ) in the Gulf of Alaska. Quantitative estimates of recruitment are obtained from models of stock assessment and stock projection employing information supplied by FOCI. To generate its information, FOCI convenes specialists in marine biology, physical and fisheries oceanography, meteorology, and statistics to assemble and analyse relevant biological and physical time series with respect to recruitment and processes hypothesized to influence fish survival. Statistical methods encompass linear and nonlinear regression, stochastic simulation modelling, transfer function time series modelling, and tree-modelling regression. The current database consists of 31 years of data, and analyses have identified factors that affect ocean stratification and circulation during spring and summer of the fish's birth year as being important to recruitment. A conceptual model of the recruitment process serves as the framework for a recruitment forecast scheme. A stochastic mathematical simulation model of the conceptual model produces similarities between simulated and observed recruitment time series. FOCI has successfully forecast recruitment observed over the past several years.     

Physical transport of young pollock larvae (Theragra chalcogramma) near Shelikof Strait as inferred from a hydrodynamic model

An advective model was used to simulate the drift of larval walleye pollock ( Theragra chalcogramma ) over a 40-day period (late April through early June) near Shelikof Strait, Alaska. This model was used: (i) to assess how much of the observed change in larval positions during that period can be explained by transport at fixed depth; (ii) to demonstrate that observed change can be related to mean large-scale meteorological forcing; and (iii) to investigate accumulation of larvae in specific areas near the coast. Based on availability of larval and circulation data, three years were studied: 1988, 1989 and 1991. Velocity fields generated from a hydrodynamic model driven by winds and runoff were used to advect particles seeded in accordance with observed larval distributions in late April of each year. The modelled larvae were tracked at 40 m depth, corresponding to the mean depth of sampled larvae and the depth of neutrally buoyant drifters employed in field studies. Specific features observed in late May larval surveys were reproduced by the model, such as the accumulation of larvae in a shoal area downstream of the strait. Differences among the modelled years include extensive flushing of larvae to the south-west in 1988 and 1991, vs. limited flushing in 1989. These differences appear related to the mean large-scale atmospheric pressure patterns for April-May of those years.     

A synthesis of biological and physical processes affecting the feeding environment of larval walleye pollock (Theragra chalcogramma) in the eastern Bering Sea

Biological and physical phenomena that affect conditions for larval survival and eventual recruitment differ in the oceanic and shelf regions. In the oceanic region, eddies are a common feature. While their genesis is not well known, eddies have unique biophysical characteristics and occur with such regularity that they likely affect larval survival. High concentrations of larval pollock often are associated with eddies. Some eddies are transported onto the shelf, thereby providing larvae to the Outer Shelf Domain. Advection, rather than local production, dominated the observed springtime increase in chlorophyll (often a correlate of larval food) in the oceanic region. Over two-thirds of the south-eastern shelf, eddies are absent and other phenomena are important. Sea ice is a feature of the shelf region: its interannual variability (time of arrival, persistence, and areal extent) affects developmental rate of larvae, timing of the phytoplankton bloom (and potentially the match/mismatch of larvae and prey), and abundance and distribution of juvenile pollock. In the oceanic region, interannual variation in food for first-feeding pollock larvae is determined by advection; in the shelf region, it is the coupled dynamics of the atmosphere–ice–ocean system.     

Interannual variation in the coastal distribution of a juvenile gadid in the northeast Pacific Ocean: The relevance of wind and effect on recruitment

Drift of propagules occurs within many populations inhabiting flow fields. This affects the number of propagules that rejoin their source population (recruitment) and plays a role in adaptive spatial redistribution. We focus on the cause and consequence of interannual variation in geographic distribution of population density among five cohorts of young‐of‐the‐year (age‐0) juvenile walleye pollock Gadus chalcogrammus in the western Gulf of Alaska (GOA). The coastal GOA is a wind‐driven advective system. Walleye pollock spawn during spring and their eggs and larvae drift southwestward; by late summer, age‐0 juveniles are variously distributed over the shelf. We found that high population densities of age‐0 juveniles (ca. 6 months old) near the southwestward exit of the Alaska Coastal Current from the GOA corresponded with high abundance of larvae from the major spawning area upstream, but did not translate into high abundance at older ages. Further, offshore and upwelling‐favorable winds were associated with the high downstream abundance and presumed export. In contrast, downwelling‐favorable (northeasterly) wind during and shortly after spawning (April–May) was associated with high recruitment at age 1. Finally, we found that recruitment also increased with apparent retention of age‐0 juveniles in favorable habitat upstream near the main spawning area. We hypothesize that wind‐related retention in superior upstream habitat favors recruitment. Our results argue for including wind‐driven transport in future walleye pollock recruitment models. We encourage more work on the juvenile stage of marine fishes aimed at understanding how transport and species‐specific habitat suitability interact to affect population response to large‐scale forcing.     

Comparison of factors affecting recruitment variability of walleye pollock Theragra chalcogramma in the Pacific Ocean and the Sea of Japan off northern Japan

The Japanese Pacific stock (JPS) and the northern Japan Sea stock (JSS) of walleye pollock Theragra chalcogramma are mainly distributed in the Pacific Ocean and the Sea of Japan off northern Japan, respectively. This paper summarizes and compares the factors affecting the recruitment variability of these two stocks. Spawning season is from December to March for both stocks. JPS recruitment has a positive relationship with the water temperature in January and February, whereas that of JSS has a negative relationship with the water temperature in January, February, and April. One possible reason for this is that pollock larvae have an optimum growth temperature of approximately 5 °C in the field. Drift of early life stages also appears to be an important influence on the recruitment of both stocks. Because the current generated by the northwest wind carries eggs of JPS into the main larval nursery ground, JPS recruitment is enhanced in years when the northwest wind is predominant in February. On the other hand, early life stages of JSS are transported into the nursery ground by the Tsushima Warm Current. However, this current also carries early life stages into the Sea of Okhotsk and offshore, resulting in poor JSS recruitment in years when this current is strong in March. In contrast to JPS, the recruitment of which is significantly impacted by cannibalism, young pollock have not been found in the stomachs of adult JSS. Warm temperatures in the Sea of Japan seem to induce the separation of young and adult pollock, and the shape of the stock–recruitment relationship also suggests that cannibalism is not important for JSS. Based on this knowledge, and on the hatch date distributions of larvae and juveniles, we propose mechanisms that can explain the recruitment fluctuations for JPS and JSS pollock.     

Use of juvenile salmon growth and temperature change indices to predict groundfish post age‐0 yr class strengths in the Gulf of Alaska and eastern Bering Sea

Juvenile marine growth (SW1) of salmon and a new temperature change (TC) index were evaluated as ecosystem indicators and predictors for the post age‐0 year class strength (YCS) of groundfish in the Gulf of Alaska (GOA) and eastern Bering Sea (EBS). Our hypothesis was that SW1, as measured on the scales of adult Pacific salmon (Oncorhynchus spp.), is a proxy for ocean productivity on the continental shelf, a rearing area for young salmon and groundfish. Less negative TC index values are the result of a cool late summer followed by a warm spring, conditions favorable for groundfish YCS. In the GOA, SW1 was a positive predictor of age‐1 pollock (Theragra chalcogramma), but not age‐2 sablefish (Anoplopoma fimbria) YCS, indicating that the growth of the Karluk River sockeye salmon that enter Shelikof Strait is a proxy for ocean conditions experienced by age‐0 pollock. Contrary to our hypotheses, the TC index was a negative predictor of GOA pollock YCS; and the SW1 a negative predictor of EBS pollock and cod YCS since the 1980s. Recent fisheries oceanography survey results provide insight into possible mechanisms to support the inverse SW1 and YCS relationship. For the EBS, the TC index was a significant positive predictor for pollock and cod YCS, supporting the hypothesis that a cool late summer followed by a warm spring maximizes the over‐wintering survival of pollock and cod (Gadus macrocephalus), especially since the 1980s. The TC and SW1 index showed value for the assessment of pollock and cod, but not sablefish.     

A review of the distribution, ecology and population dynamics of age-0 walleye pollock in the Gulf of Alaska

This review paper synthesizes published research and unpublished data on the abundance and distribution patterns, ecology and population dynamics of walleye pollock ( Theragra chalcogramma ) during their first year of life (age-0) in the Gulf of Alaska. Distribution patterns have been described using mainly trawl catches, but recently, acoustic methodology has been employed, especially in examining vertical distributions. Although age-0 pollock are found throughout the Gulf, the highest catches occurred west of Kodiak Island. Pollock are pelagic for at least their first 6 months of life but show an ontogenetic increase in depth distribution superimposed on a pronounced diel vertical migration at a larger size. Daily growth rates are variable depending on year, season and area, and growth generally ceases during the winter. The diet of age-0 pollock shifts from mainly copepods in early juveniles to euphausiids by fall, with epibenthic organisms becoming important during the winter months. Feeding occurs mainly at night in surface waters. Age-0 pollock are most frequently associated with gelatinous zooplankton (medusae) and older pollock. Many predators on age-0 pollock have been identified; the most important are arrowtooth flounder ( Atheresthes stomias ), adult pollock, puffins ( Fratercula spp.), murres ( Una spp.), harbor seals ( Phoca vitulina richardsi ) and Steller sea lions ( Eumetopiasjubatus ). Modelling provides some insight into the population dynamics of these juveniles and environmental conditions which interannually affect their survival. These results are discussed relative to the importance of age-0 pollock in the recruitment of this species and to their role in the pelagic ecosystem.     

Observation and analysis of fishery processes: larval pollock at Shelikof Strait, Alaska

Conditions affecting distributions of larval walleye pollock ( Theragra chalcogramma ) were examined at Shelikof Strait, Alaska, during springtime, 1986 and 1987. Abundance and distribution of larval pollock southwest of the Strait's southern entrance was determined with oblique plankton tows taken each year in May. Infrared images of sea surface temperature patterns were derived from AVHRR scenes obtained by NOAA satellites during each April and May. Pattern displacements between 24-hour-interval images were used to estimate surface motion. Each spring, measurements were taken by remote weather stations and ships, and a nearsurface current meter record was obtained during 1987. Treated as quasi-synoptic, spatial relations between sets of surface temperature, surface flow, and larval pollock distributions show coincidences between submesoscale physical and biological features. The highest larval abundances occurred as patches within a cold plume (1986) and an eddy (1987). These confirm that physical features can retain larval pollock on the continental shelf. Observations are examined for evidence of physical and biological events that jointly can cause such coincidences and foster alternatives for survival during transport to nursery grounds. Explanations for presence of cohorts observed within the 1987 eddy are given in terms of spatial and temporal relationships evident between spawning and hatching areas, hatch date distributions, meanders, eddy generation and movement, background flow, and advection times. The observations, analyses, and results are consistent with the concept of a coupled, fluctuating biophysical process that can emulate variations in larval abundance and provide a multiplicity of system pathways for early-life stages representations.     

Eastern Bering Sea pollock recruitment,abundance, distribution and approach to fishery management

Eastern Bering Sea pollock have two distinctly different stable spawning grounds—along the shelf and in the eastern and central Aleutian Islands between 400 and 500 m water columns. Pollock spawning behavior supports the hypothesis that the shelf and deepwater “basin” spawning pollock are completely independent reproductive stocks. Deepwater pollock inhabit the shelf and, once mature at age 5–6 years, migrate from the shelf onto the continental slope into the Zhemchug, Pribilof, and Bering canyons by the end of winter. Bering Sea pollock recruitment and year class abundance have high annual variability, but there are no clear relationships between pollock year class strength and water temperature, ice distribution or survival on early ontogenesis stages (eggs and larvae). Young-of-the-year fish survival varies dramatically during winter supporting the hypothesis that the Bering Sea pollock recruitment and strength of year class have high annual variability depending on young-of-the-year fish survival during winter. The annual change of physical oceanography condition, productivity and species composition of zooplankton community are associated with great differences in pollock seasonal migrations and distribution, reproduction, survival of recruits at early stages of development and finally with abundance of year classes and total biomass. Implementation of ecosystem-based fishery management most important for application of pollock research both of Russian national program and on base of International Agreements.     

Geographical variations in infection by larval Anisakis simplex and Contracaecum osculatum (Nematoda, Anisakidae) in walleye pollock Theragra chalcogramma stocks off Hokkaido, Japan

ABSTRACT: Stocks of walleye pollock Theragra chalcogramma collected from: (i) the Sea of Japan (off Rebun Island and Kumaishi); (ii) the Pacific coast (off Shikabe and eastern Hokkaido); and (iii) Nemuro Strait off Hokkaido, northern Japan, were examined for anisakid nematodes during December 1999 to February 2000, and the prevalence and abundance of Anisakis simplex and Contracaecum osculatum larvae were compared among the various sampling sites for fish of the same size and age. Anisakis simplex was generally more abundant than C. osculatum . Infection by A. simplex varied between the aforementioned stocks of walleye pollock as well as within stocks, whereby fish from off Rebun Island and Nemuro Strait were infected the most, followed by those from off the Pacific coast and Kumaishi. Infection by C. osculatum differed between the host stocks, and C. osculatum was the most abundant among the fish from Nemuro Strait. The infection variations seemed to be due to differences in host growth rate, host feeding habit, and the distribution of marine mammal final hosts. The results indicate that these two larval nematodes are useful biological indicators for the population study of walleye pollock in Japanese waters.     

The importance of episodic weather events to the ecosystem of the Bering Sea shelf

Climate variability on decadal time scales is generally recognized to influence high‐latitude marine populations. Our recent work in studying air–sea interactions in the Bering Sea suggests that interannual to decadal climate variability is important through its modulation of the frequencies and magnitudes of weather events on intraseasonal time scales. We hypothesize that it is these weather events that directly impact the marine ecosystem of the Bering Sea shelf. The linkages between the event‐scale weather and the ecosystem are illustrated with three examples: walleye pollock (Theragra chalcogramma), Tanner crabs (Chionoecetes bairdi), and coccolithophorid phytoplankton (Emiliania huxleyi). We hypothesize that the strong recruitment of walleye pollock that occurred in 1978, 1982, and 1996 can be attributed in part due to the seasonably strong storms that occurred in the early summer of those years. These storms caused greater than normal mixing of nutrients into the euphotic zone which presumably led to sustained primary productivity after the spring bloom and, possibly, enhanced prey concentrations for pollock larvae and their competitors. Recruitment of Tanner crab was particularly strong for the 1981 and 1984 year‐classes. These years had periods of prominent east wind anomalies along the Alaska Peninsula during the previous winter. Such winds promote flow through Unimak Pass, and hence an enhanced flux of nutrient‐rich water onto the shelf. This mechanism may have ultimately resulted in favorable feeding conditions for Tanner crab larvae. Finally, an unprecedented coccolithophorid bloom occurred over the Bering Sea shelf in the summer of 1997. This summer featured lighter winds and greater insolation than usual after a spring that included a very strong May storm. This combination brought about a warm, nutrient‐poor upper mixed layer by mid‐summer. This provided a competitive advantage for coccolithophorid phytoplankton in 1997 and to a lesser extent in 1998. Unusually high concentrations of coccolithophores persisted for the following two years although physical environmental conditions did not remain favorable. While slow variations in the overall aspects of the physical environment may be important for setting the stage, we propose that the significant multi‐year adjustments in the marine ecosystem of the Bering Sea shelf are more directly caused by major air–sea interaction events on intraseasonal time scales.     

Interannual variability of the early life history of walleye pollock near Shelikof Strait as inferred from a spatially explicit, individual-based model

A coupled biophysical model is used to hindcast the early life history of a population of walleye pollock ( Theragra chalcogramma ), to assess possible physical causes of interannual variability in recruitment. The modelling approach combines a primitive equation, rigid'lid hydrodynamic model with a probabilistic, individual-based biological model of growth, development, and mortality. Individuals are tracked through space using daily velocity fields generated from the hydrodynamic model, along with self-directed vertical migrations appropriate to each life stage in the biological model. The hydrodynamic model is driven with wind and runoff time series appropriate to each year. Biological model output compares favourably with observed spatial distributions for specific years. Lloyd's index of patchiness, calculated from model output, was similar to values calculated from field data. Five noncontiguous years were chosen for hindcasts to span a wide range of meteorological conditions (winds, runoff) and recruitment success. Interannual comparisons suggest that two years of above average recruitment (1978 and 1988), and one year of below average recruitment (1991), experienced flow fields which carried many individuals into the Alaskan Stream. At the same time, the vigorous flow fields generated in each of these years carried some individuals onto the shelf area to the south-west of the spawning site. A year with low runoff and weak winds (1989) exhibited weak circulation, with extended retention of larvae near the spawning site. A year with high runoff (1987) was notable for the strength and frequency of mesoscale eddy activity. Eddies appear capable of both enhancing patchiness of early larvae (through retention) and dissipating patchiness of juveniles (through mesoscale mixing). Larvae retained in an eddy feature exhibit a narrower range of sizes than the population outside that feature.     

Cohort survival patterns of walleye pollock, Theragra chalcogramma, in Shelikof Strait, Alaska: a critical factor analysis

A series of age-specific life tables for walleye pollock ( Theragra chalcogramma ) in the western Gulf of Alaska was compiled for the 1980-91 year classes. The life tables were utilized to perform an exploratory key factor analysis to examine the timing of critical periods in the recruitment process, evidence of density-dependence at different stages and trends in mortality rates. Early larval mortality was significantly correlated with generational mortality (In recruits/spawning bio-mass), but patterns in juvenile mortality also were similar to generational mortality and in some years were clearly dominant in determining the fate of a cohort. Density-dependent mortality, based on the correlation between mortality and initial abundance, was indicated only for the late larval to early juvenile stage. Time trends were marginally significant for juvenile mortality. It is speculated that the observed increase in juvenile mortality is associated with increasing abundance of arrowtooth flounder. Weaknesses in the data base are discussed; these along with the short time series involved make our conclusions tentative and subject to further study. We hypothesize that pollock recruitment levels can be established at any life stage depending on sufficient supply from prior stages, a type of dynamics which can be termed supply dependent multiple life stage control.     

Life-history traits of walleye pollock, Theragra chalcogramma, in the northeastern Japan Sea during early to mid 1990s

I examined the age, growth, maturity, mortality, and body condition of walleye pollock, Theragra chalcogramma, in the northeastern Japan Sea (northern Japan Sea population) and evaluated their resilience to exploitation. Walleye pollock were collected in pre-spawning (October 1991-1995) and post-spawning (April 1990-1996) seasons. Estimated ages ranged from 3 to 18 years for both sexes. A von Bertalanffy growth model showed that females had longer asymptotic fork length (460 mm) than males (425 mm). Fifty percent of females and males were mature at 348 mm (4.6 years) and 322 mm (3.9 years), respectively. The instantaneous natural mortality rate was estimated to be 0.22. These life-history traits in the northern Japan Sea population were compared to those in the Bering Sea, the Gulf of Alaska, and the Japan Pacific populations. As a result, female walleye pollock in this population matured at small body sizes, grew rapidly toward small maximum sizes, and had short reproductive lifespans with low size-specific fecundity and poor body condition. Low prey availability and habitat temperatures are considered as a possible mechanism for the small maximum sizes in this population. The potential rate of population increase of both the northern Japan Sea population and other pollock populations tended to be lower than other exploited populations of non-viviparous marine fishes, suggesting potentially lower resilience to exploitation in this population and walleye pollock populations in general.     

Spatial and temporal patterns of walleye pollock (Theragra chalcogramma) spawning in the eastern Bering Sea inferred from egg and larval distributions

Walleye pollock Theragra chalcogramma (pollock hereafter) is a key ecological and economic species in the eastern Bering Sea, yet detailed synthesis of the spatial and temporal patterns of pollock ichthyoplankton in this important region is lacking. This knowledge gap is particularly severe considering that egg and larval distribution are essential to reconstructing spawning locations and early life stages drift pathways. We used 19 yr of ichthyoplankton collections to determine the spatial and temporal patterns of egg and larval distribution. Generalized additive models (GAMs) identified two primary temporal pulses of pollock eggs, the first occurring from 20 February to 31 March and the second from 20 April to 20 May; larvae showed similar, but slightly lagged, pulses. Based on generalized cross‐validation and information theory, a GAM model that allowed for different seasonal patterns in egg density within three unique areas outperformed a GAM that assumed a single fixed seasonal pattern across the entire eastern Bering Sea. This ‘area‐dependent’ GAM predicted the highest densities of eggs (i.e., potential spawning locations) in three major areas of the eastern Bering Sea: near Bogoslof Island (February–April), north of Unimak Island and the Alaska Peninsula (March–April), and around the Pribilof Islands (April–August). Unique temporal patterns of egg density were observed for each area, suggesting that pollock spawning may be more spatially and temporally complex than previously assumed. Moreover, this work provides a valuable baseline of pollock spawning to which future changes, such as those resulting from climate variability, may be compared.     

Linking Northeast Pacific recruitment synchrony to environmental variability

We investigated the hypothesis that synchronous recruitment is due to a shared susceptibility to environmental processes using stock–recruitment residuals for 52 marine fish stocks within three Northeast Pacific large marine ecosystems: the Eastern Bering Sea and Aleutian Islands, Gulf of Alaska, and California Current. There was moderate coherence in exceptionally strong and weak year‐classes and correlations across stocks. Based on evidence of synchrony from these analyses, we used Bayesian hierarchical models to relate recruitment to environmental covariates for groups of stocks that may be similarly influenced by environmental processes based on their life histories. There were consistent relationships among stocks to the covariates, especially within the Gulf of Alaska and California Current. The best Gulf of Alaska model included Northeast Pacific sea surface height as a predictor of recruitment, and was particularly strong for stocks dependent on cross‐shelf transport during the larval phase for recruitment. In the California Current the best‐fit model included San Francisco coastal sea level height as a predictor, with higher recruitment for many stocks corresponding to anomalously high sea level the year before spawning and low sea level the year of spawning. The best Eastern Bering Sea and Aleutian Islands model included several environmental variables as covariates and there was some consistent response across stocks to these variables. Future research may be able to utilize these across‐stock environmental influences, in conjunction with an understanding of ecological processes important across early life history stages, to improve identification of environmental drivers of recruitment.