Essential omega-3 fatty acids are depleted in sea ice and pelagic algae of the Central Arctic Ocean.

Microalgae are the main source of the omega-3 fatty acids eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), essential for the healthy development of most marine and terrestrial fauna including humans. Inverse correlations of algal EPA and DHA proportions (% of total fatty acids) with temperature have led to suggestions of a warming-induced decline in the global production of these biomolecules and an enhanced importance of high latitude organisms for their provision. The cold Arctic Ocean is a potential hotspot of EPA and DHA production, but consequences of global warming are unknown. Here, we combine a full-seasonal EPA and DHA dataset from the Central Arctic Ocean (CAO), with results from 13 previous field studies and 32 cultured algal strains to examine five potential climate change effects; ice algae loss, community shifts, increase in light, nutrients, and temperature. The algal EPA and DHA proportions were lower in the ice-covered CAO than in warmer peripheral shelf seas, which indicates that the paradigm of an inverse correlation of EPA and DHA proportions with temperature may not hold in the Arctic. We found no systematic differences in the summed EPA and DHA proportions of sea ice versus pelagic algae, and in diatoms versus non-diatoms. Overall, the algal EPA and DHA proportions varied up to four-fold seasonally and 10-fold regionally, pointing to strong light and nutrient limitations in the CAO. Where these limitations ease in a warming Arctic, EPA and DHA proportions are likely to increase alongside increasing primary production, with nutritional benefits for a non-ice-associated food web.

Male reproductive system of the deep-sea acorn worm Quatuoralisia malakhovi (Hemichordata, Enteropneusta, Torquaratoridae) from the Bering Sea.

BACKGROUND: The deep-sea acorn worm Quatuoralisia malakhovi belongs to the phylum Hemichordata, class Enteropneusta, family Torquaratoridae, which was described in 2005. Owing to their epibenthic lifestyle and deep-sea habitat features, torquaratorids differ anatomically from shallow-water acorn worms; however, their morphology and fine structure are poorly studied. We have the opportunity to make three complete detailed series of histological sections of Q. malakhovi and to study the microscopic anatomy, histology and fine structure of the reproductive system of this acorn worm using scanning and transmission electron microscopy. RESULTS: The sexes of Q. malakhovi are separate and indistinguishable externally. The lobed testes occupy the dorsal side of the genital wings and distinctly bulge into the peribranchial cavity by their mature lobes. The central part of the testis is always submerged into the genital wing and opens via a single gonad pore. The monociliary muscle cells stretch along the external wall of the testis and surround the gonad pore, probably taking part in the contraction of the testis lobes for spawning. The germinative epithelium of the testis contains spermatogenic cells at different stages of development and interstitial cells. Yolk cells are not found. Interstitial cells embrace the spermatogonia and spermatogenic columns, providing horizontal compartmentalization of the germinative epithelium, and contain numerous phagosomes with remnants of degenerating spermatogenic cells. The testis wall contains haemal lacunae, which are usually located on the side opposite the gonad pore. We describe the fine structure of spermatogonia, spermatocytes clustered in spermatogenic columns, spermatids, and spermatozoa. Spermatozoa are of the ectaquasperm type and consist of an acorn-shaped head and a flagellum 18-25 µm long. The sperm head includes a beak-shaped acrosomal part, a spherical nucleus and a midpiece containing a ring of 5 or rarely 6 mitochondria. CONCLUSIONS: The male reproductive system and sperm structure of Q. malakhovi, a representative of the family Torquaratoridae, have a number of differences from shallow-water acorn worms; however, the spermatogenesis and sperm structure of Q. malakhovi generally follow the pattern of the other three enteropneust families, and the phylogenetic significance of these deviations should be the subject of further research.

Nemertodermatida-Endosymbionts of Deep-Sea Acorn Worms (Hemichordata, Torquaratoridae).

Worm-like endosymbionts were found in the hepatic region of the digestive tract of the deep-sea acorn worm Quatuoralisia malakhovi Ezhova et Lukinykh, 2022 (family Torquaratoridae) from the Bering Sea. The symbionts were assigned to the taxon Nemertodermatida on the basis of histological examination. Torquaratoridae are similar in feeding type to holothuroids, which have also been found to have Xenacoelomorpha endosymbionts.

Digestive System Anatomy and Feeding Mechanism of Quatuoralisia malakhovi (Hemichordata, Torquaratoridae).

The digestive system was anatomically studied in the deep-sea enteropneust Quatuoralisia mala-khovi. It was shown that lateral collar lips are twisted in such a way that they form a ciliary groove that leads to an internal channel, through which collected detritus particles are transferred to peripheral pharyngeal channels. The size of the selected particles ranges from 1-6 to 100-200 µm, which corresponds to feeding on the remains of planktonic diatoms. A fecal cord was observed to act as an anchor that holds the heavily watered jelly-like body of Torquaratoridae at the sea floor during feeding.

Shifts in the composition and distribution of Pacific Arctic larval fish assemblages in response to rapid ecosystem change.

The Pacific Arctic marine ecosystem has undergone rapid changes in recent years due to ocean warming, sea ice loss, and increased northward transport of Pacific-origin waters into the Arctic. These climate-mediated changes have been linked to range shifts of juvenile and adult subarctic (boreal) and Arctic fish populations, though it is unclear whether distributional changes are also occurring during the early life stages. We analyzed larval fish abundance and distribution data sampled in late summer from 2010 to 2019 in two interconnected Pacific Arctic ecosystems: the northern Bering Sea and Chukchi Sea, to determine whether recent warming and loss of sea ice has restricted habitat for Arctic species and altered larval fish assemblage composition from Arctic- to boreal-associated taxa. Multivariate analyses revealed the presence of three distinct multi-species assemblages across all years: (1) a boreal assemblage dominated by yellowfin sole (Limanda aspera), capelin (Mallotus catervarius), and walleye pollock (Gadus chalcogrammus); (2) an Arctic assemblage composed of Arctic cod (Boreogadus saida) and other common Arctic species; and (3) a mixed assemblage composed of the dominant species from the other two assemblages. We found that the wind- and current-driven northward advection of warmer, subarctic waters and the unprecedented low-ice conditions observed in the northern Bering and Chukchi seas beginning in 2017 and persisting into 2018 and 2019 have precipitated community-wide shifts, with the boreal larval fish assemblage expanding northward and offshore and the Arctic assemblage retreating poleward. We conclude that Arctic warming is most significantly driving changes in abundance at the leading and trailing edges of the Chukchi Sea larval fish community as boreal species increase in abundance and Arctic species decline. Our analyses document how quickly larval fish assemblages respond to environmental change and reveal that the impacts of Arctic borealization on fish community composition spans multiple life stages over large spatial scales.

Joint species distribution modeling reveals a changing prey landscape for North Pacific right whales on the Bering shelf.

The eastern North Pacific right whale (NPRW) is the most endangered population of whale and has been observed north of its core feeding ground in recent years with low sea ice extent. Sea ice and water temperature are important drivers for zooplankton dynamics within the whale's core feeding ground in the southeastern Bering Sea, seasonally forming stable fronts along the shelf that give rise to distinct zooplankton communities. A northward shift in NPRW distribution driven by changing distribution of prey resources could put this species at increased risk of entanglement and vessel strikes. We modeled the abundance of NPRW prey, Calanus glacialis, Neocalanus, and Thysanoessa species, using a dynamic biophysical food web model of nine zooplankton guilds in the Bering shelf zooplankton community during a period of warming (2006-2016). This model is unique among prior zooplankton studies from the region in that it includes density dependence, thereby allowing us to ask whether species interactions influence zooplankton dynamics. Modeling confirmed the importance of sea ice and ocean temperature to zooplankton dynamics in the region. Density-independent growth drove community dynamics, while dependent factors were comparatively minimal. Overall, Calanus responded to environment terms, with the strength and direction of response driven by copepodite stage. Neocalanus and Thysanoessa responses were weaker, likely due to their primary occurrence on the outer shelf. We also modeled the steady-state (equilibrium) abundance of Calanus in conditions with and without wind gusts to test whether advection of outer shelf species might disrupt the steady-state dynamics of Calanus abundance; the results did not support disruption. Given the annual fall sampling design, we interpret our results as follows: low-ice-extent winters induced stronger spring winds and weakened fronts on the shelf, thereby advecting some outer shelf species into the study region; increased development rates in these warm conditions influenced the proportion of C. glacialis copepodite stages over the season. Residual correlation suggests missing drivers, possibly predators, and phytoplankton bloom composition. Given the continued loss of sea ice in the region and projected continued warming, our findings suggest that C. glacialis will move northward, and thus, whales may move northward to continue targeting them.

The First Discovery of Trematodes (Digenea) in the Deep-Sea Acorn Worms Torquaratoridae (Hemichordata, Enteropneusta).

Trematodes found in the enteropneust hemichordates are described for the first time. Metacercariae have been found in the trunk coelom, in the collar coelom, in the proboscis coelom, and in the glomerulus of the deep-sea torquaratorid Quatuoralisia malakhovi Ezhova et Lukinykh, 2022. This is the first find of parasites in the glomerulus of acorn worms. The taxonomy of the found trematodes is discussed.

Marine biodiversity refugia in a climate-sensitive subarctic shelf.

The subarctic shelf of the Eastern Bering Sea (EBS) is one of the world's most productive marine environments, exposed to drastic climate changes characterized by extreme fluctuations in temperature, sea ice concentration, timing, and duration. These climatic changes elicit profound responses in species distribution, abundance, and community composition. Here, we examined the patterns of alpha and temporal beta diversity of 159 marine taxa (66 vertebrates and 93 invertebrate species) from 29 years (1990-2018) of species observations from the NOAA bottom trawl surveys in the EBS. Based on these data, we identified geographically distinct refugial zones in the northern and southern regions of the middle shelf, defined by high species richness and similarity in community species composition over time. These refugial zones harbor higher frequencies of occurrence for representative taxa relative to the regions outside of refugia. We also explored the primary environmental factors structuring marine biodiversity distributions, which underpinned the importance of the winter sea ice concentration to alpha and temporal beta diversity. The spatial biodiversity distributions between high and low winter sea ice regimes highlighted contrasting signals. In particular, the latter showed elevated species richness compared to the former. Further, the temporal beta diversity between the high and low winter sea ice periods underpinned an overall increase in the compositional similarity of marine communities in the EBS. Despite these spatiotemporal differences in biodiversity distributions, the identified refugia represent safe havens of marine biodiversity in the EBS. Distinguishing these areas can help facilitate conservation and management efforts under accelerated and ongoing climatic changes.

Feeding Characteristics of Deep-Sea Acorn Worm (Hemichordata, Enteropneusta, Torquaratoridae) from the Bering Sea.

Deep-sea hemichordates Torquaratoridae gen. sp. reach high abundance up to 12 spec. m-2 at the depths of 1830-2130 m on the slope of the Volcanologists Massif in the south-western part of the Bering Sea, dominating in the benthic community at these depths. Their abundance exceeds by two orders the values recorded earlier. In order to clarify this phenomenon, we examined the gut contents of Torquaratoridae gen. sp. The detritus particles and frustules of planktonic diatoms Thalassiosira, Coscinodiscus, Actinocyclus, Chaetoceros, Neodenticula, and Grammatophora were the most common in the gut, as well as the remains of skeletons of benthic invertebrates with little admixture of mineral particles. According to obtained data, Torquaratoridae gen. sp. are mobile deposit feeders with high selectivity to fresh phytodetritus, able to compete with holothurians occupying similar trophic niche. Unusually high abundance of acorn worms is apparently related to high organic matter flux to the seafloor as a result of spring phytoplankton bloom in the surface water layer of the Bering Sea.

Tracking long-distance migration of marine fishes using compound-specific stable isotope analysis of amino acids.

The long-distance migrations by marine fishes are difficult to track by field observation. Here, we propose a new method to track such migrations using stable nitrogen isotopic composition at the base of the food web (δ15 NBase ), which can be estimated by using compound-specific isotope analysis. δ15 NBase exclusively reflects the δ15 N of nitrate in the ocean at a regional scale and is not affected by the trophic position of sampled organisms. In other words, δ15 NBase allows for direct comparison of isotope ratios between proxy organisms of the isoscape and the target migratory animal. We initially constructed a δ15 NBase isoscape in the northern North Pacific by bulk and compound-specific isotope analyses of copepods (n = 360 and 24, respectively), and then we determined retrospective δ15 NBase values of spawning chum salmon (Oncorhynchus keta) from their vertebral centra (10 sections from each of two salmon). We then estimated the migration routes of chum salmon during their skeletal growth by using a state-space model. Our isotope tracking method successfully reproduced a known chum salmon migration route between the Okhotsk and Bering seas, and our findings suggest the presence of a new migration route to the Bering Sea Shelf during a later growth stage.

Trait-based climate vulnerability assessments in data-rich systems: An application to eastern Bering Sea fish and invertebrate stocks.

Trait-based climate vulnerability assessments based on expert evaluation have emerged as a rapid tool to assess biological vulnerability when detailed correlative or mechanistic studies are not feasible. Trait-based assessments typically view vulnerability as a combination of sensitivity and exposure to climate change. However, in some locations, a substantial amount of information may exist on system productivity and environmental conditions (both current and projected), with potential disparities in the information available for data-rich and data-poor stocks. Incorporating this level of detailed information poses challenges when conducting, and communicating uncertainty from, rapid vulnerability assessments. We applied a trait-based vulnerability assessment to 36 fish and invertebrate stocks in the eastern Bering Sea (EBS), a data-rich ecosystem. In recent years, the living marine resources of the EBS and Aleutian Islands have supported fisheries worth more than US $1 billion of annual ex-vessel value. Our vulnerability assessment uses projections (to 2039) from three downscaled climate models, and graphically characterizes the variation in climate projections between climate models and between seasons. Bootstrapping was used to characterize uncertainty in specific biological traits and environmental variables, and in the scores for sensitivity, exposure, and vulnerability. The sensitivity of EBS stocks to climate change ranged from "low" to "high," but vulnerability ranged between "low" and "moderate" due to limited exposure to climate change. Comparison with more detailed studies reveals that water temperature is an important variable for projecting climate impacts on stocks such as walleye pollock (Gadus chalcogrammus), and sensitivity analyses revealed that modifying the rule for determining vulnerability increased the vulnerability scores. This study demonstrates the importance of considering several uncertainties (e.g., climate projections, biological, and model structure) when conducting climate vulnerability assessments, and can be extended in future research to consider the vulnerability of user groups dependent on these stocks.

Arctic source for elevated atmospheric mercury (Hg0) in the western Bering Sea in the summer of 2013.

Measurements of gaseous elemental mercury (Hg0) in the marine boundary layer of the western Bering Sea were performed using an automatic mercury analyzer RA 915+ (Ltd. "Lumex", St. Petersburg, Russia) aboard the Russian research vessel Academician M.A. Lavrentev from 3 to 20 August 2013. Hg0 concentrations varied from 0.3 to 2.1ng/m3 (n=4783); the average value (1.1±0.3ng/m3) was lower than both the background range of the Northern Hemisphere (1.5-1.7ng/m3) and average values previously observed in the Bering Sea, and corresponded to the background concentrations of the Southern Hemisphere (1.1-1.3ng/m3). Maximum Hg0 concentrations were observed within air masses that came from the lower troposphere of the central Arctic. Under these conditions, Hg0 ranged between 1.1 and 2.1ng/m3 with an average of 1.5±0.2ng/m3 (n=1183). Except for these periods, Hg0 concentrations during the rest of the study varied from 0.3 to 1.8ng/m3, with an average value of 1.0±0.2ng/m3 (n=3600). Our results support the hypothesis that, in the summer, air masses from the central Arctic Ocean can be an exporter of mercury to lower latitudes. Perhaps the atmospheric transport of elevated concentrations of Hg0 into lower latitudes may have implications for the biologic and economic health of important fisheries, such as the Bering Sea.

Decadal shifts in autumn migration timing by Pacific Arctic beluga whales are related to delayed annual sea ice formation.

Migrations are often influenced by seasonal environmental gradients that are increasingly being altered by climate change. The consequences of rapid changes in Arctic sea ice have the potential to affect migrations of a number of marine species whose timing is temporally matched to seasonal sea ice cover. This topic has not been investigated for Pacific Arctic beluga whales (Delphinapterus leucas) that follow matrilineally maintained autumn migrations in the waters around Alaska and Russia. For the sympatric Eastern Chukchi Sea ('Chukchi') and Eastern Beaufort Sea ('Beaufort') beluga populations, we examined changes in autumn migration timing as related to delayed regional sea ice freeze-up since the 1990s, using two independent data sources (satellite telemetry data and passive acoustics) for both populations. We compared dates of migration between 'early' (1993-2002) and 'late' (2004-2012) tagging periods. During the late tagging period, Chukchi belugas had significantly delayed migrations (by 2 to >4 weeks, depending on location) from the Beaufort and Chukchi seas. Spatial analyses also revealed that departure from Beaufort Sea foraging regions by Chukchi whales was postponed in the late period. Chukchi beluga autumn migration timing occurred significantly later as regional sea ice freeze-up timing became later in the Beaufort, Chukchi, and Bering seas. In contrast, Beaufort belugas did not shift migration timing between periods, nor was migration timing related to freeze-up timing, other than for southward migration at the Bering Strait. Passive acoustic data from 2008 to 2014 provided independent and supplementary support for delayed migration from the Beaufort Sea (4 day yr-1 ) by Chukchi belugas. Here, we report the first phenological study examining beluga whale migrations within the context of their rapidly transforming Pacific Arctic ecosystem, suggesting flexible responses that may enable their persistence yet also complicate predictions of how belugas may fare in the future.

Sources and distribution of sedimentary organic matter along the northern Bering and Chukchi Seas.

In this study, lignin-derived phenols were used to determine the sources and distribution of sedimentary organic matter along the northern Bering Sea and Chukchi Sea of the Arctic Ocean. The lignin parameter syringyl/vanillyl (S/V) and cinnamyl/vanillyl (C/V) ratios are used to indicate vegetation sources; and the ratios of vanillic acid/vanillin, (Ad/Al)v and syringic acid/syringaldehyde, (Ad/Al)s are used as indicators of lignin diagenesis. Results showed the predominance of woody gymnosperm signal at the easternmost location in the northern Bering Sea, a mixture of refractory non-woody angiosperm and fresher gymnosperm tissues in the Chukchi Sea, and signal of fresher woody gymnosperm tissues in the northernmost locations in the Chukchi Sea. The lignin materials showed gradual increase in decomposition stage during transport along the northern Bering Sea. Hydrodynamic sorting process, which is the retention of coarser materials nearshore and transportation of finer particles farther offshore, most probably occurred along the east coast of the northern Bering Sea. In Chukchi Sea, the non-woody angiosperm tissues could have originated from the Canadian Arctic and gymnosperm tissues could be from the Russian Arctic side. The fresher materials in the northernmost Chukchi Sea could have been transported here via the ice-rafting process. Detection of fresh lignin materials and the occurrence of lignin decomposition mean that this region could be sensitive to the impact of climate change.

Distribution of Dall's Porpoise, Phocoenoides dalli, in the North Pacific and Bering Sea, Based on T/S Oshoro Maru 2012 Summer Cruise Data.

Dall's porpoise (Phocoenoides dalli) is a small toothed cetacean, widely inhabiting the North Pacific Ocean and adjacent seas, between about 30 and 62°N; however, only limited studies of its ecology have been made in nearshore areas. A cetacean sighting survey lasting 60 days was conducted during the 2012 summer cruise of the T/S Oshoro Maru (Hokkaido University, Japan) in the North Pacific Ocean and Bering Sea. Based on this data, the distribution of Dall's porpoises and the factors controlling it in the pelagic habitat were investigated. A total of 808 individual Dall's porpoises in 166 groups were sighted during a total of 469.6 hr and 4946.6 nm observations. The cruise consisted of three legs and the average porpoise group size was significantly larger in Leg 1. The sightings were concentrated at water depths of less than 1000 m and near eastern Aleutian passes. Sighting clusters were found on the 200 m isobath of the southeastern Bering Sea continental slope. There was a peak in sightings where the sea surface temperature (SST) was relatively cold, between 5 and 7°C. Although similar track routes were taken in Leg 1 and Leg 3, the number of sightings per unit effort was larger in Leg 1. This difference may have arisen from the significant rise in SST as the season progressed. Relatively large group size found in this study might relate with prey abundance along the Aleutian Islands.

Subtle genetic population structure in Pacific halibut Hippoglossus stenolepis.

Pacific halibut Hippoglossus stenolepis from 10 sampling locations throughout their range were investigated for signs of population structure. Two genetic data sets were created: (1) all individuals (n = 828) at few anonymous microsatellite markers (number of loci = 16); (2) fewer individuals (n = 435) genotyped at anonymous as well as expressed sequence-tag linked microsatellites (number of loci = 61). A combination of multidimensional scaling plots, discriminant analysis of principal components and pairwise differentiation estimates suggested that samples from the Aleutian Islands, particularly the western Aleutian Islands, were genetically distinct from samples collected in other regions. In addition, outlier analyses found that two markers linked to expressed sequence tags may be under directional selection and could explain the differentiation among samples. These results confirm findings from previous research and suggest that population structure may exist within a current management unit (i.e. International Pacific Halibut Commission Regulatory Area 4B).

Body size affects individual winter foraging strategies of thick-billed murres in the Bering Sea.

Foraging and migration often require different energetic and movement strategies. Though not readily apparent, constraints during one phase might influence the foraging strategies observed in another. For marine birds that fly and dive, body size constraints likely present a trade-off between foraging ability and migration as smaller bodies reduce flight costs, whereas larger bodies are advantageous for diving deeper. This study examines individual wintering strategies of deep diving thick-billed murres (Uria lomvia) breeding at three colonies in the south-eastern Bering Sea: St Paul, St George and Bogoslof. These colonies, arranged north to south, are located such that breeding birds forage in a gradient from shelf to deep-water habitats. We used geolocation time-depth recorders and stable isotopes from feathers to determine differences in foraging behaviour and diet of murres during three non-breeding periods, 2008-2011. Body size was quantified by a principal component analysis (wing, culmen, head+bill and tarsus length). A hierarchical cluster analysis identified winter foraging strategies based on individual movement, diving behaviour and diet (inferred from stable isotopes). Structural body size differed by breeding island. Larger birds from St Paul had higher wing loading than smaller birds from St George. Larger birds, mainly from St Paul, dove to deeper depths, spent more time in the Bering Sea, and likely consumed higher trophic-level prey in late winter. Three winter foraging strategies were identified. The main strategy, employed by small birds from all three breeding colonies in the first 2 years, was characterized by high residency areas in the North Pacific south of the Aleutians and nocturnal diving. In contrast, 31% of birds from St Paul remained in the Bering Sea and foraged mainly during the day, apparently feeding on higher trophic-level prey. Throat feather stable isotopes indicated that individuals exhibited flexibility in the use of this colony-specific foraging strategy. The third strategy only occurred in 2010/2011, when birds dove more and deeper, suggesting limited prey resources. Foraging strategies partitioned with respect to annual differences, presumably in response to shifts in distribution of prey, and were linked to body size. The presence of a colony-specific wintering strategy suggests the potential for overwinter survival differences between these populations.

Ichthyophonus-infected walleye pollock Theragra chalcogramma (Pallas) in the eastern Bering Sea: a potential reservoir of infections in the North Pacific.

In 2003, the Alaska walleye pollock industry reported product quality issues attributed to an unspecified parasite in fish muscle. Using molecular and histological methods, we identified the parasite in Bering Sea pollock as Ichthyophonus. Infected pollock were identified throughout the study area, and prevalence was greater in adults than in juveniles. This study not only provides the first documented report of Ichthyophonus in any fish species captured in the Bering Sea, but also reveals that the parasite has been present in this region for nearly 20 years and is not a recent introduction. Sequence analysis of 18S rDNA from Ichthyophonus in pollock revealed that consensus sequences were identical to published parasite sequences from Pacific herring and Yukon River Chinook salmon. Results from this study suggest potential for Ichthyophonus exposures from infected pollock via two trophic pathways; feeding on whole fish as prey and scavenging on industry-discharged offal. Considering the notable Ichthyophonus levels in pollock, the low host specificity of the parasite and the role of this host as a central prey item in the Bering Sea, pollock likely serve as a key Ichthyophonus reservoir for other susceptible hosts in the North Pacific.

Plankton community changes during the last 124 000 years in the subarctic Bering Sea derived from sedimentary ancient DNA.

Current global warming results in rising sea-water temperatures, and the loss of sea ice in Arctic and subarctic oceans impacts the community composition of primary producers with cascading effects on the food web and potentially on carbon export rates. This study analyzes metagenomic shotgun and diatom rbcL amplicon sequencing data from sedimentary ancient DNA of the subarctic western Bering Sea that records phyto- and zooplankton community changes over the last glacial-interglacial cycles, including the last interglacial period (Eemian). Our data show that interglacial and glacial plankton communities differ, with distinct Eemian and Holocene plankton communities. The generally warm Holocene period is dominated by picosized cyanobacteria and bacteria-feeding heterotrophic protists, while the Eemian period is dominated by eukaryotic picosized chlorophytes and Triparmaceae. By contrast, the glacial period is characterized by microsized phototrophic protists, including sea ice-associated diatoms in the family Bacillariaceae and co-occurring diatom-feeding crustaceous zooplankton. Our deep-time record of plankton community changes reveals a long-term decrease in phytoplankton cell size coeval with increasing temperatures, resembling community changes in the currently warming Bering Sea. The phytoplankton community in the warmer-than-present Eemian period is distinct from modern communities and limits the use of the Eemian as an analog for future climate scenarios. However, under enhanced future warming, the expected shift toward the dominance of small-sized phytoplankton and heterotrophic protists might result in an increased productivity, whereas the community's potential of carbon export will be decreased, thereby weakening the subarctic Bering Sea's function as an effective carbon sink.

Hydrological characteristics of the Bering Sea in the summer of 2019.

Based on the data of three CTD sections in the northern, northeastern and western Bering Sea of 2019 voyage of Chinese National Arctic Research Expedition (CHINARE), this paper analyzes and studies the hydrological characteristics of the water mass distribution, layered structure, and cline characteristics in different sea areas of the Bering Sea. The results indicate that the hydrological characteristics of the Bering Sea in the summer of 2019 are different from those in the past and that the water mass is warming in many locations. The maximum water temperature reaches 11.13 °C, and the maximum thickness of the warm water is about 32 m. The water mass composition and characteristics of the north-northeast-west sections are significantly different: the BL section has the highest salinity, while the BS section has the lowest salinity, and both the lowest temperature and the largest temperature variation appear in the BL section. The stratification characteristics in all sea areas are noteworthy. In the deep-water seas, there are three types of water masses: upper water (BSW), middle water (BIW) and deep water (BDW) from top to bottom, while two main water masses appear in the shelf waters with the Alaska Coastal Water (ACW) overlies the Bering Sea Shelf Water (BSW). Along the Bering Sea Slope Current (BSC), the water mass is essentially steady. Statically unstable hydrological inversion structure appears near the bottom of the three stations at the northern end of the BL section.