Abstract
Understanding mechanisms behind variability in early life survival of marine fishes through modeling efforts can improve predictive capabilities for recruitment success under changing climate conditions. Walleye pollock (Theragra chalcogramma) support the largest single-species commercial fishery in the United States and represent an ecologically important component of the Bering Sea ecosystem. Variability in walleye pollock growth and survival is structured in part by climate-driven bottom-up control of zooplankton composition. We used two modeling approaches, informed by observations, to understand the roles of prey quality, prey composition, and water temperature on juvenile walleye pollock growth: (1) a bioenergetics model that included local predator and prey energy densities, and (2) an individual-based model that included a mechanistic feeding component dependent on larval development and behavior, local prey densities and size, and physical oceanographic conditions. Prey composition in late-summer shifted from predominantly smaller copepod species in the warmer 2005 season to larger species in the cooler 2010 season, reflecting differences in zooplankton composition between years. In 2010, the main prey of juvenile walleye pollock were more abundant, had greater biomass, and higher mean energy density, resulting in better growth conditions. Moreover, spatial patterns in prey composition and water temperature lead to areas of enhanced growth, or growth 'hot spots', for juvenile walleye pollock and survival may be enhanced when fish overlap with these areas. This study provides evidence that a spatial mismatch between juvenile walleye pollock and growth 'hot spots' in 2005 contributed to poor recruitment while a higher degree of overlap in 2010 resulted in improved recruitment. Our results indicate that climate-driven changes in prey quality and composition can impact growth of juvenile walleye pollock, potentially severely affecting recruitment variability.