Large copepods as leading indicators of walleye pollock recruitment in the southeastern Bering Sea: Sample-Based and spatio-temporal model (VAST) results

Abstract In the southeastern Bering Sea, years with above-average sea ice extent have been associated with elevated abundances of large lipid-rich copepods in both net tows and age-0 walleye pollock (Gaduschalcogrammus) diets, followed by high overwinter pollock survival and stronger recruitment into the fishery three years later. In this study, we directly compare large copepod abundances of taxa important in age-0 pollock diets (Calanus glacialis/marshallae, Metridia pacifica and Neocalanus spp.) to age-3 pollock abundance from stock assessments for year classes 2002–2015 (excluding 2013). Copepod samples from bongo nets were collected on fishery oceanography surveys in late summer during warm (low ice) and cold (high ice) climate stanzas. In addition to mean estimates of large copepods among stations (sample-based), we implemented a spatial delta-generalized linear mixed model using the Vector Autoregressive Spatio-Temporal (VAST) package to account for variations in survey coverage among years. Positive significant linear relationships were found for both of these estimates of large copepods during the age-0 year and the abundance of these pollock at age-3. Correlations were stronger using VAST model-based indices compared with sample-based indices (R2 = 0.74 versus 0.43, respectively, all copepod taxa combined). A significant relationship was also observed between the bottom cold pool (