Distribution pattern of Polar cod (<Emphasis Type="Italic">Boreogadus saida</Emphasis>) larvae and larval fish assemblages in relation to oceanographic parameters in the northern Bering Sea and Chukchi Sea

Polar cod, Boreogadus saida, is a key species in the Arctic Ocean ecosystem. We examined the distribution pattern of B. saida and other fish larvae in relation to oceanographic parameters, including sea surface temperature and salinity (SST and SSS), the mode of temperature and salinity within the water column (F _temp and F _sal), and the temporal duration between the date of sea ice retreat and the date of field surveys (dSRT) in the northern Bering Sea and Chukchi Sea during the summers of 2008 and 2013. Sampling was conducted onboard the T/S Oshoro-Maru using a bongo net for 2 years. At sampling stations, the temperature and salinity were measured using conductivity-temperature-depth profiler casts. We calculated dSRT from satellite derived polar gridded sea ice concentration data. A total of 1186 individuals comprising 7 families and 16 species were collected, with B. saida (35 %) and Ammodytes hexapterus (27 %) dominating the catch in number. Based on the species composition (cluster analysis), the sampling stations were divided into four groups. Pleuronectidae dominated group A, which was characterized by relatively high temperature (SST and F _temp), while B. saida dominated group B, characterized by low temperatures. A. hexapterus dominated group C and D, which had similar temperatures and salinities, but group C had a shorter dSRT than group D. The latter was also differentiated by the presence of Lumpenus sp. B. saida, which were most abundant in regions where temperatures ranged from ?2 to 0.5 °C. These results, the geographical variation of group B sampling stations, suggest there are two separate spawning areas characterized by similar fish communities and related to consistently cold bottom temperatures. Moreover, stations with a higher abundance of smaller sized larval B. saida were characterized by a short dSRT, whereas stations with a lower abundance and a larger size were characterized by a long dSRT.