Southeastern Bering Sea Research Articles

Modelling the multiple action pathways of projected climate change on the Pacific cod (Gadus macrocephalus) early life stages

Understanding how future ocean conditions will impact early life stages and population recruitment of fishes is critical for adapting fisheries communities to climate change. In this study, we incorporated projected changes in physical and biological ecosystem dynamics from an oceanographic model into a mechanistic individual-based model for larval and juvenile stages of the Pacific cod (Gadus macrocephalus) in the eastern Bering Sea. We particularly investigated the impacts of ocean currents, temperature, prey density, and pCO2 on the hatching success, growth, survival, and spatial distribution of this species during 2021–2100. We evaluated two CO2 emission scenarios: RCP8.5 (high CO2 emissions, low mitigation efforts) and RCP4.5 (medium CO2 emissions and mitigation efforts). We found that the increase in temperature and decrease in prey density were the main drivers of faster growth rates and lower survival through increased starvation by the end of the century. Conversely, pCO2 had negligible impacts, which suggests that this species might be resilient to ocean acidification. The largest effects were observed under the high CO2 emission scenario, while the RCP4.5 projections displayed minimal impacts. We also identified an area with favourable conditions in the southeastern Bering Sea that will likely persist in future decades. This study provides relevant information on the future impacts of climate change on Pacific cod, and our results can be used to implement and inform climate-ready management for this important stock in Alaska.

Modeling the larval growth and survival of Pacific cod (Gadus macrocephalus) in the eastern Bering Sea

The eastern Bering Sea (EBS) is a highly productive ecosystem that supports several important commercial species such as the Pacific cod (Gadus macrocephalus). Climate variability affects the population dynamics of this stock throughout its life stages, especially early life stages, since they are particularly susceptible to environmental changes. In recent decades, warm and cold stanzas (i.e., 3–5 year periods) have been observed in the EBS, and there is evidence that they can modulate the recruitment of this stock, causing important socioeconomic impacts. Using a mechanistic individual-based model, this study investigates the spatial and temporal variability of growth and survival of Pacific cod's early life stages during 2000–2020. We examined changes by year and over space and compared our results with published literature to validate our model. We found that temperature played a key role in modulating the survival of fish larvae, observing an increase in starvation events in warmer years or locations. Periods or areas with low prey density, especially small-bodied copepods, also contributed to increased starvation. The average temperature in the fish habitat was negatively correlated with recruitment estimates from the stock assessment model. Growth was primarily temperature-driven; however, food-limited growth became more frequent when larvae were smaller during cold years. Spatially, we found that the environmental conditions in the southeastern Bering Sea may favor larval survival but reduce growth, and higher mortality may be persistent on the middle and outer shelf. Our model produces results that agree with previous field studies, and it offers a valuable tool to investigate other ecological questions on the impact of the environment on early life stages of fishes.

Biological and environmental covariates of juvenile sockeye salmon distribution and abundance in the southeastern Bering Sea, 2002-2018.

Climate change is altering the distribution and abundance of marine species, especially in Arctic and sub-Arctic regions. In the eastern Bering Sea, home of the world's largest run of sockeye salmon (Oncorhynchus nerka), juvenile sockeye salmon abundance has increased and their migration path shifted north with warming, 2002-2018. The reasons for these changes are poorly understood. For these sockeye salmon, we quantify environmental and biological covariate effects within spatio-temporal species distribution models. Spatio-temporally, with respect to juvenile sockeye salmon densities: (1) sea surface temperature had a nonlinear effect, (2) large copepod, Calanus, a minor prey item, had no effect, (3) age-0 pollock (Gadus chalcogrammus), a major prey item during warm years, had a positive linear effect, and (4) juvenile pink salmon (O. gorbuscha) had a positive linear effect. Temporally, annual biomass of juvenile sockeye salmon was nonlinearly related to sea temperature and positively related to age-0 pollock and juvenile pink salmon abundance. Results indicate that sockeye salmon distributed with and increased in abundance with increases in prey, and reached a threshold for optimal temperatures in the eastern Bering Sea. Changes in population dynamics and distribution of sockeye salmon in response to environmental variability have potential implications for projecting specific future food securities and management of fisheries in Arctic waters.

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.

Seasonal Dynamics of Primary Production in the Southeastern Bering Sea Assessed Using Continuous Temporal and Vertical Dissolved Oxygen and Chlorophyll‐a Measurements

AbstractThe Bering Sea and other high‐latitude systems are experiencing unprecedented changes related to climate warming and the consequent loss of sea ice. Understanding how such changes influence primary production is a pressing question. Here, we quantify primary production rates in the southeastern Bering Sea over 4 years (2016–2019) at daily to weekly time resolution. A moored high‐resolution Profiling Crawler was used in combination with in situ sampling to collect physical and biological data in the ocean's upper 50 m. We used dissolved oxygen (O2) data to estimate gross and net primary production (NPP), and chlorophyll‐a (Chl‐a), temperature, and irradiance data to model NPP rates. We then assessed seasonal variation in phytoplankton production rates, phytoplankton community growth rates () and explored to what extent summer phytoplankton community growth rates may be influenced by nitrogen limitation. Our analyses revealed that the majority of gross primary production (GPP) and net community production occurs in association with the spring phytoplankton bloom. After the bloom, the water column generally experienced low GPP and net biological carbon consumption. Phytoplankton growth rates were commonly suppressed in late summer due to apparent nitrogen depletion. Using high temporal resolution vertically resolved measurements of O2 and Chl‐a, our analyses provides important insight into how biogeochemical cycles, phytoplankton community growth rates, and carbon available for export vary seasonally in the southeastern Bering Sea.

Meiofauna in the southeastern Bering Sea: community composition and structuring environmental factors

The Bering Sea is the second largest marginal sea in the North Pacific and is one of the areas with highest biological productivity in high-latitude waters. The continental shelf of the Bering Sea hosts large populations of marine mammals and fishery resources. However, the smaller organisms in benthic ecosystems, including meiofauna, have been largely overlooked in this area, despite their potential importance in ecosystem functioning and the resultant biogeochemical cycles. This study analyzed spatial differences in the total abundance and community structure of the metazoan meiofauna at five stations around the Bering Canyon, located at the southeastern margin of the Bering Sea. Their association with environmental factors in sediments was also studied. The results confirmed that the investigated stations had meiofaunal standing stocks that were comparable to those of other Arctic seas. Among the investigated sediment biological and geochemical parameters (total organic carbon, median grain size, prokaryotic cell numbers, etc.), multivariate analyses showed that the C/N of organic matter in sediments was the main factor associated with meiofaunal community structure.

Long-Term Biophysical Observations and Climate Impacts in US Arctic Marine Ecosystems

In 1995, the first of a nearly continuous sequence of biophysical moorings was deployed at a site (M2) on the southeastern Bering Sea shelf. Over the next 15 years, 10 additional mooring sites were initiated. The resultant long-term biophysical mooring array extends over 1,800 km from the southern Bering Sea to the northern Chukchi Sea, covering most of the US Arctic. It provides a full range of oceanographic data for researchers, stakeholders, and managers. In addition, these data sets have been critical for the validation of regional ocean models. The ocean temperature data have quantified regional warming and formed the basis for understanding how warmer temperatures and loss of sea ice are modifying these high-latitude marine ecosystems. Changes observed in the context of observations from the mooring program include delayed spring bloom, low abundances of large crustacean zooplankton and crab species, seabird die-offs, changes in ocean acidification, northward expansion of subarctic fish species, and shifts in the ranges of marine mammal species.

The northern Bering Sea zooplankton community response to variability in sea ice: evidence from a series of warm and cold periods

Recent, unprecedented losses of sea ice have resulted in widespread changes in the northern Bering Sea ecosystem, and this study explores the zooplankton community response. Time-series observations were used to identify zooplankton community changes in the northern (>60°N) Bering Sea (NBS) over a 17 yr period (2002-2018). The overall objective was to determine if the changes in zooplankton populations previously described for the southeastern Bering Sea shelf (<60°N) were also observed in the NBS over alternating warm and cold periods. Particular attention was paid to more recent (2014-2018) years that showed significant losses of sea ice in the NBS (2017/2018) in comparison to a prior warm period (2003-2005) and an intervening cold period (2006-2013). A multivariate framework (redundancy analysis) was used to explore correlations with environmental conditions, and differences in mean abundance across the differing warm and cold periods were tested. The NBS zooplankton community had different responses across each warm and cold period, and the primary driver for the differences in response was sea ice. Redundancy analysis demonstrated that the zooplankton community during the second warm period experienced greater variability compared to the prior warm period. The zooplankton community had higher abundances of small copepods and meroplankton and reduced abundances of Calanus spp. and chaetognaths during the most recent warm period. This suggests that the NBS zooplankton will not be impacted by reduced sea ice when the ice coverage extends south of 60°N, but show community change once a minimum threshold in ice extent and timing of retreat is reached. Shifts in the zooplankton community may have had cascading effects on higher trophic levels that were evident during the latter warm period.

Changes in the vertical distribution of age‐0 walleye pollock (Gadus chalcogrammus) during warm and cold years in the southeastern Bering Sea

AbstractIn the last 20 years, the southeastern Bering Sea has shifted its thermal variability to longer‐term (4–6 years) ocean–ecosystem temperature stanzas. Age‐0 walleye pollock (Gadus chalcogrammus) populations respond to thermal changes with horizontal (east–west) shifts in spatial distribution over the continental shelf, though there are limited data on whether thermally mediated vertical shifts in distribution also occur. Vertical shifts may impact predator–prey overlap between age‐0 pollock and their lipid‐rich prey, calanoid copepods and euphausiids, resulting in different feeding conditions that ultimately affect fish body condition prior to winter onset. For this study, we analyzed acoustic backscatter measured during acoustic trawl surveys over the southeastern Bering Sea shelf in cold years (2011, 2012) and warm years (2014 and 2016) to determine the vertical distribution of age‐0 pollock. This study presents evidence that age‐0 pollock changes in vertical distribution were related to changes in ocean annual temperature. Age‐0 pollock went from occurring deeper in the water column during cold periods to being surface‐oriented during warm periods, potentially exacerbating spatial mismatches between pollock and prey. We relate patterns in a vertical position to physical water column properties, feeding, and bioenergetic condition of collected pollock and discuss implications for recruitment success during different thermal oceanographic stanzas.

Delineating yellowfin sole (Limanda aspera) reproduction in the northern Bering Sea provides information across the eastern Bering Sea continental shelf

Yellowfin sole (Limanda aspera) is an abundant, commercially harvested flatfish that ranges across the northern and southeastern Bering Sea continental shelf. In recent years, the summer bottom trawl survey of the southeastern Bering Sea (SEBS) conducted by the National Marine Fisheries Service’s Alaska Fisheries Science Center (AFSC) extended into the northern Bering Sea (NBS) in August. This opportunity allowed for describing yellowfin sole reproductive parameters in the NBS, such as female length and age at maturation, reproductive status, and sex ratio distribution. Estimates of 50% female sexual maturity were 10.11 years (A50; 95% CI: 9.47–10.76 years; n = 209) and 28.47 cm (L50; 95% CI: 27.16–29.68 cm; n = 212). Histology indicated much of the mature population was approaching the end of spawning. Yellowfin sole spawning timing (summer) appears to be synchronous across the NBS and SEBS. There is evidence that part of the yellowfin sole spawning population in the NBS is connected to a recognized spawning migratory group from the SEBS. A higher proportion of females than males was observed, which varied by year and stratum. This was negatively related to warmer bottom temperature and positively related to location (western longitudes and northern latitudes). Yellowfin sole females in the NBS exhibit a similar size and age of maturation, depth of spawning, and sex ratio proportions with those females inhabiting the SEBS. These results provide further information for fisheries managers on the yellowfin sole stock throughout this region. Spatial and temporal aspects of reproduction, however, should be more thoroughly investigated.

Dynamically downscaled projections of ocean acidification for the Bering Sea

A regional ocean biogeochemical model for the Bering Sea is used to dynamically downscale three Earth System Models from the CMIP5 archive under the RCP 8.5 and RCP 4.5 scenarios. These continuous model runs, completed in conjunction with the Alaska Climate Integrated Modeling Project (ACLIM), span the 2006–2100 timeframe and project continued warming, freshening, and ocean acidification (OA) for the Bering Sea shelf region over the 21st Century, with larger magnitude changes in the RCP 8.5 scenario. The downscaled projections suggest that annual average surface seawater aragonite saturation state (Ωarag) for the Bering Sea shelf will decrease by 0.63–0.86 under RCP 8.5 and 0.18–0.43 under RCP 4.5 by 2100. Surface pH values decrease by 0.31–0.35 under RCP 8.5 and 0.07–0.13 under RCP 4.5. Seasonally, Ωarag < 1 conditions start to emerge for ∼2 months per year during winter between 2015 and 2030 under both climate change scenarios. Under RCP 8.5, the duration of these undersaturated conditions grows to ∼5 months per year by 2100, occurring from mid-October through mid-March. Under RCP 4.5, these conditions remain constrained to 2–3 months per year by 2100. In both scenarios, summer months maintain conditions of Ωarag > 1 due to primary productivity, though the maximum in Ωarag is greatly reduced under RCP 8.5. Spatially, the regions of greatest pH and Ωarag decline are the southeastern Bering Sea shelf and the outer shelf domain near the shelf break. Linear trends in carbonate variables between our downscaled simulations and the Earth System Model (ESM) output are comparable and indistinguishable compared to the model spread. However, bottom water trends differ somewhat between the ESM and our downscaled simulations, with the latter more consistently resolving the different shelf domains. The OA information provided by these downscaled simulations can help inform biological sensitivity experiments and longterm strategic planning for marine fisheries management.

Alaskan Seabird Die-Offs

Prior to 2015, seabird die-offs in Alaskan waters (Northern Gulf of Alaska, eastern Bering Sea, eastern Chukchi Sea) were rare, typically occurred in mid-winter, and were linked to epizootic disease events (Bodenstein et al., 2015) or above-average ocean temperatures associated with strong El Nino-Southern Oscillation events. An exception was late summer of 1997, a year with unusually warm waters in the southeastern Bering Sea, when possibly as many as 10% of the several million short-tailed shearwaters (Ardenna tenuitortris) present in the area died of starvation because their principal prey, euphausiids, were scarce or hard to find in a widespread coccolithophore bloom (Baduini et al., 2001). Since 2015, such mass mortality events have been annual occurrences. In 2015, between 470,000 and 1.03 million common murres (Uria aalge) were estimated to have died in the Gulf of Alaska due to anomalous ocean conditions and the impacts on forage fish associated with the 2014–2016 Pacific marine heatwave (Piatt et al., 2020; Arimitsu et al., 2021). Another large die-off event occurred during late summer of 2019 (Figure 1), when over 10,000 carcasses of short-tailed shearwaters washed onto beaches of Bristol Bay and the Alaska Peninsula in the southeastern Bering Sea. Total mortality during that summer was likely much greater than reported given the region’s expansive coastline and sparse human population.

Decreased lipid storage in juvenile Bering Sea crabs (Chionoecetes spp.) in a warm (2014) compared to a cold (2012) year on the southeastern Bering Sea

Decreased lipid storage in juvenile Bering Sea crabs (Chionoecetes spp.) in a warm (2014) compared to a cold (2012) year on the southeastern Bering Sea

The potential utility of otolith microchemistry as an indicator of nursery origins in Pacific halibut (Hippoglossus stenolepis)

Identifying the nursery origins of individual Pacific halibut (Hippoglossus stenolepis) is critical to a refined understanding of regional productivity, the effects of environmental forcing on population dynamics, and the spatial impacts of mortality and bycatch. Otolith microchemistry may provide a tool for identifying fish origins, but a precursor to using microchemistry-based assignment models is to evaluate the extent to which there exists detectable spatial structure in these natural tags. Here, we examined assignment accuracy as a function of spatial scale using canonical discriminant function analyses (DFA) applied to microchemical data from age-2 Pacific halibut in the western Gulf of Alaska (GOA) and southeastern Bering Sea (SEBS). Element:calcium ratios were assayed for fourteen trace elements using inductively coupled plasma mass spectrometry; δ13C and δ18O were determined via isotope ratio mass spectrometry. Substantial assignment success (∼75–90 %) was observed at spatial scales that are consistent with stock management. Elemental signatures were defined primarily by δ18O, δ13C, and 88Sr:48Ca, with a minor contribution from 55Mn:48Ca. Individuals were most commonly mis-assigned to adjacent locations. However, discontinuities were observed within the western GOA that suggest that elemental signatures do not vary along strictly longitudinal and latitudinal clines. These results highlight the need to exercise caution when attempting to use otolith microchemistry to assign fish to their origins when the baseline elemental data are unable to resolve missing nursery sources from the locations that are included in the discrimination models.

Long-term dynamics for stocks of sablefish Anoplopoma fimbria in the western Bering Sea and prospects for their fishery

Sablefish is an endemic species of the North Pacific. Its range extends from California Peninsula, along the Pacific coast of the US and Canada to Aleutian Islands and further, along the Pacific coast of Kamchatka and the Kuriles to the central part of Honshu Island. They dwell also in the Bering Sea and southeastern Okhotsk Sea. Sablefish are the most abundant in the southeastern Bering Sea and in the Gulf of Alaska, that is conditioned by favorable conditions for their larvae and juveniles. In the Asian part of the range, the environments are generally more severe, and reproduction of sablefish is rather risky. Following to the results of modern genetic studies, the sablefish stocks are distinguished by high genetic homogeneity that suggests a common population with the main spawning grounds in the southeastern Bering Sea, at the Pacific coasts of Aleutian Islands, in the Gulf of Alaska, and at the coasts of British Columbia, Washington, Oregon and California. Dynamics of the sablefish biomass is considered on the data of bottom and midwater trawl surveys conducted by TINRO in 2003–2020, fishery statistics, and accessible data of NOAA (USA). Sharp increasing of the biomass and annual catches is noted both in the eastern and western Bering Sea in the last few years because of appearance of several strong year-classes. Western Bering Sea stock depends on migration of recruits from the common spawning grounds in the southeastern Bering Sea. For the western Bering Sea, two main ways of such migration are possible: i) active migration of juveniles with benthic habitat; and ii) passive transfer of pelagic larvae and early juveniles across the Bering Sea through the system of surface currents. The latter mechanism supports the sablefish recruitment in the bays of the western Bering Sea and, to a lesser extent, at the eastern coast of Kamchatka. Sablefish in the West Bering Sea fishery zone were caught in 2010–2020 mostly as by-catch for trawling and longline fishery (93 %), other 7 % were landed by specialized longline fishery. The basic points for managing the sablefish fishery in the West Bering Sea zone are defined. About 400 t of sablefish is permissible to catch annually in the West Bering Sea fishery zone in conditions of modern high stock of this species. This value includes 100–120 t that will inevitably be caught as by-catch and the rest of 280–300 t is a foreseeable resource for organization of specialized fishery.

The Influence of a Deep‐Water Intrusion on the Distribution of Chrysaora melanaster in the Southeastern Bering Sea

AbstractScyphozoan jellyfish are conspicuous components of marine ecosystems, which during a bloom, can impact food web structure and economically important fisheries. Jellyfish biomass in the southeastern Bering Sea (SEBS) is primarily composed of Chrysaora melanaster and has varied widely over the past four decades, yet the underlying causes of these biomass fluctuations remain unclear. The present study investigated the spatial and seasonal dynamics of C. melanaster along the Alaska Peninsula in the SEBS using an adaptive resolution imaging sonar system and nets in June–July and September 2018. The abundance of C. melanaster was high in coastal waters near the Alaska Peninsula, with peak densities occurring east of Unimak Island in both sampling periods. The current pattern revealed by an acoustic Doppler current profiler showed that cold, deep water from Bering Canyon flowed onto the shelf east of Unimak Pass. Differences in the strength of this cross‐isobath, deep‐water intrusion were observed during the two sampling periods. This flow may have influenced the distribution of C. melanaster by constraining the mixed coastal water near of Unimak Island in June‐July and by transporting medusae westward, into deeper water in September, after the deep‐water intrusion weakened. The present study provides a new mechanistic understanding of how regional‐scale ocean variability in the SEBS affects C. melanaster distribution.

Climate warming and the loss of sea ice: the impact of sea-ice variability on the southeastern Bering Sea pelagic ecosystem

AbstractWe investigated relationships among three metrics of sea-ice cover in eight regions of the eastern Bering Sea and the abundance of Calanus copepods, jellyfish medusae, and year-class strength of walleye pollock (Gadus chalcogrammus). In summer, Calanus spp. were more abundant over the middle shelf when sea ice lingered late into spring, and, to a lesser extent, when February sea-ice cover was heavy. Between 1982 and 1999, there were no significant (p ≤ 0.05) relationships between the amount or timing of sea-ice cover and pollock recruitment. However, between 2000 and 2015, pollock year-class strength was positively correlated with sea ice in the outer and middle shelves, with 17 of 24 regressions significant. Pollock year-class strength was best predicted by days with sea-ice cover after February. Pollock recruitment was positively influenced by copepod numbers, particularly in the middle shelf, with r2 values from 0.36 to 0.47. We hypothesize that the Calanus spp. present in the southeastern Bering Sea are primarily Calanus glacialis that have been advected south in association with sea ice. None of our sea-ice metrics explained the variance in jellyfish biomass. Jellyfish biomass in our study area in the pollock age-0 year was not correlated with pollock recruitment 3 years later.

Distinctive spring phytoplankton bloom in the Bering Strait in 2018: A year of historically minimum sea ice extent

Under climate change, the Bering–Chukchi continental shelves have experienced drastic sea ice loss, breaking the minimum sea ice extent record in 2018. In this paper, we report a distinctive spring phytoplankton bloom in the Bering Strait based on year-round mooring observations during 2016–2018. A near-bottom fluorescence sensor mounted on a mooring and satellite observations indicate that, despite the early sea ice retreat in spring of 2018 (mid-April), the spring phytoplankton bloom was not observed for more than 20 days after sea ice retreat. In contrast, in 2017 when sea ice retreat did not occur as early as it did in 2018 (early May), the bloom and sea ice retreat occurred in tighter synchrony. Mooring (2016–2018) and satellite observations (2003–2011, 2013–2018) suggest that an open-water bloom rather than an ice-edge bloom occurred in the Bering Strait in the spring of 2018. Early sea ice retreat apparently provided unfavorable conditions for an ice-edge bloom due to weak sunlight and winter-like weather. Mid-winter sea ice loss was promoted by unusual southerly winds, bringing in warm water from the south. For the Bering Strait, the end of April is a reference date for sea ice retreat that distinguishes between an open-water bloom and an ice-edge bloom, which is approximately one month later than that in the southeastern Bering Sea.

Foraminiferal Ecology and Role in Nitrogen Benthic Cycle in the Hypoxic Southeastern Bering Sea

Southeastern Bering Sea is one of the highest surface productivity area in the open ocean due to strong upwelling along the Bering canyon. However, the benthic geochemistry and organisms living in the area have been largely overlooked. In August 2017, surface sediment was sampled from four stations along a transect at depths between 1536 and 103 meters in the Bering canyon with JAMSTEC R/V Mirai. Bottom-water hypoxia was recorded in the two deepest stations (1536 and 536 m). At these stations, the oxygen penetrated down to 5 mm in the sediment due to siltier and much organic-rich sediments in the deeper stations while oxygen penetration was about 20 mm at stations 103 and 197 meters deep with coarse-grained sediment stations. Foraminiferal number of species and abundances were higher in the Unimak pass depression station E2 (197 m). Abundance did not change significantly between stations, suggesting that foraminiferal densities are not affected by the hypoxic conditions but are rather controlled by organic matter and nutrients availability. At the upper bathyal and middle bathyal stations, living foraminiferal communities were in general dominated by Uvigerina peregrina, Nonionella pulchella, Elphidium batialis, Globobulimina pacifica, Reophax spp., and Bolivina spathulata while the shallower stations exhibited large densities of Uvigerina peregrina, Cibicidoides wuellerstorfi, Recurvoidella bradyi, Globocassidulina subglobosa, and Portatrochammina pacifica. More than 50% of the individuals have a potential to accumulate nitrate in their cell (from 3 to 648 mmol/L; which is from 100 to 4000 times larger than the highest concentration measured in pore water). Onboard denitrification measurements confirmed that B. spathulata, N. pulchella and G. pacifica could reduce nitrate through denitrification and foraminiferal denitrification could contribute over 6% to benthic nitrate reduction at the southeast Bering Sea. Although the foraminiferal contributions were smaller than those measured at other hypoxic areas, our study quantitatively revealed the significance of eukaryotic microbes on benthic nitrogen cycles at this area.

Die–offs, reproductive failure, and changing at–sea abundance of murres in the Bering and Chukchi Seas in 2018

Common and thick–billed murres are among the most numerous and widespread seabirds in the northern hemisphere though they appear to be especially susceptible to mass die–off events. During the spring and summer of 2018, the Bering Sea experienced warmer than average sea temperatures following a winter of unprecedented, near complete lack of sea ice. To determine if breeding murres were negatively affected by these warm sea temperatures, surveys at most of the major murre breeding colonies in the eastern Bering and eastern Chukchi seas were conducted during the 2018 breeding season. Nearly all colonies surveyed experienced near-complete reproductive failure. Concurrently, above average levels of murre mortality were observed, primarily at St. Lawrence Island and in the Bering Strait region. The timing of this mortality is somewhat unusual and concerning as it occurred during the breeding season. Based on surveys at sea, overall murre abundance offshore was generally lower in 2018 compared to recent years (2012–2017), particularly in the northern Bering Sea, yet mean density for thick–billed murres increased in the Chukchi Sea. Reproductive failure is rare at monitored murre colonies within the study area and the widespread reproductive failures observed in 2018 were unprecedented over four decades of monitoring, but consistent with the reproductive failure of murres documented in the Gulf of Alaska and southeastern Bering Sea following a large murre die–off in 2015–2016. The combination of large–scale reproductive failure, evidence of elevated levels of mortality, and low abundance of birds offshore suggest that murres experienced severe distress in the Bering and Chukchi sea region during the summer of 2018.