Echoes of the 2013–2015 Marine Heat Wave in the Eastern Bering Sea and Consequent Biological Responses

1. Introduction

An unprecedented large-scale warm water anomaly emerged in the NE Pacific in late 2013 (Freeland, 2014; Freeland, 2015; Freeland and Ross, 2019). It was first reported by Freeland (2014) based on (a) Argo data along Line P and (b) sea surface temperature (SST) data, both from January 2014, when SST anomaly in the Gulf of Alaska exceeded 4 standard deviations (Figure 1). As Freeland (2014, p. 58) noted: “Four standard deviations away from the mean state is HUGE. … Again we see huge deviations from normal conditions with temperature in January 2014 4.4 standard deviations above the mean, salinity 3 standard deviations below the mean and as a consequence of both low salinity and high temperature, we see extremely low surface density by 4.4 standard deviations.”

During 2014, this anomaly dubbed “The Blob” (Bond et al., 2015) shifted east and north, and by June 2014 reached the Bering Sea (Figure 2) (Freeland, 2015, Fig. 6-1).

The Blob lasted through 2016 in the upper 0-100-m layer yet lingered through 2017-2018 in the subsurface 100-200-m layer as documented by Freeland and Ross (2019), who also noted a possible return in late 2018 of a warm anomaly (Figure 3), which actually occurred in 2019 as reported by Scannell et al. (2020). During the first years, The Blob’s manifestation in the surface layer diminished, while its manifestation in the subsurface markedly increased as noted by Freeland and Ross (2019), Zhi et al. (2019), and Holser et al. (2022), among others (Figure 3).

The Blob’s impact on the marine ecology of the Northeast Pacific was dramatic in terms of magnitude, duration, and spatial extent. Ecosystem implications of The Blob are well documented by now in the California Current (Cavole et al., 2016) and around the Gulf of Alaska (GOA; Suryan et al., 2021). Less is known about physical and biological manifestations of The Blob in the Bering Sea, particularly in the Eastern Bering Sea (EBS) and Northern Bering Sea (NBS), largely because of the uncertain timeline of The Blob’s propagation from the open NE Pacific to the Bering Sea. The present review is an attempt to address this issue.

Given the regional circulation pattern that connects the GOA with EBS via the Alaskan Stream, Aleutian Current, and Aleutian North Slope Current (Figure 4), one would expect to see surface waters affected by the Blob in the open NE Pacific to be advected by the above currents into the SE Bering Sea and farther downstream, along the EBS Shelf, all the way to the Bering Strait. In the first part of this paper, we review available evidence of The Blob’s propagation to the Bering Sea and eventually to the Bering Strait. We also review publicly available data on SST, sea ice cover, and cold pool on the EBS Shelf, as well as other physical factors that directly or indirectly affect marine biota, particularly fish. In the second part of this paper, we review ecological consequences of The Blob’s propagation to and across the Bering Sea, with a strong focus on the EBS Shelf (EBSS) and NBS Shelf (NBSS).

2. Data and Methods

We reviewed English-language literature on The Blob and its physical and biological ramifications in the Northeast Pacific, from the Gulf of Alaska to the Bering Sea and Bering Strait. While our main focus was on the peer-review papers published in academic journals, we also searched for data and data analyses in the grey literature, particularly in regular ecosystem assessments published by the NOAA. Since this review is squarely focused on The Blob that emerged in 2013, our time scope was limited to the last decade, from 2013 till now. Rare exceptions were made for papers describing similar anomalies in the past; a comparison of such anomalies with The Blob can be quite useful.

Numerous biological consequences of The Blob were reported in the literature, including The Blob’s effects on zooplankton and ecologically significant and commercially important species of fish and invertebrates, their life stages, and species-specific and ontogeny-associated spatial shifts in distribution and abundance (e.g., Mueter et al., 2021; Suryan et al., 2021). Biological events linked (hypothetically or conclusively) to The Blob’s passage include phytoplankton blooms such as coccolithophore blooms (Ladd et al., 2018; Matson et al., 2019), and toxic algal blooms (harmful algal blooms, HABs) (Peterson et al., 2016). In turn, toxic algal blooms cause mass die-offs of sea birds and marine mammals (e.g., whales).

3. Formation, Evolution, and Propagation of The Blob to the Bering Sea

Results of our review are arranged along-stream, presenting evidence of The Blob and documenting a time line of The Blob’s propagation from the open Northeast Pacific into the Bering Sea, then along the EBS Shelf to the Bering Strait. Naturally, the downstream narrative is ordered chronologically.

The Blob’s sudden emergence at Ocean Station Papa (OSP) in late 2013 was analyzed in detail by Freeland and Ross (2019) who noted (ibid., p. 5): “Fig. 5 shows the very sudden onset of The Blob at the end of 2013 and that it was fully present in January 2014. This agrees with the sudden onset seen in the SST maps (Fig. 1).” Freeland and Ross (2019) also noted that “a large salinity anomaly occurs at the same time and this anomaly is not density-compensating, rather there is a large fresh anomaly with salinity departing from the mean by 3.0 standard deviations. The result is that a very large density anomaly occurs with the surface waters becoming less dense by 4.3 standard deviations. … The result was the creation of a substantial barrier to wind-induced mixing and as the mixed-layer deepening season progressed through early 2014, with low storm activity, we saw anomalously shallow mixed layers.” (ibid., p. 6). In absolute terms, the near-surface salinity in January 2014 was 32.45 psu, a drop of 0.2 in a year from the mean surface salinity of 34.65 psu in 2007-2013. The stratification barrier notwithstanding, The Blob gradually extended to the subsurface over 2014-2016 as documented by Freeland and Ross (2019). Temperature and salinity data collected in 2014-2017 by instrumented northern elephant seals in the NE Pacific were analyzed by Holser et al. (2022) to provide evidence of The Blob extending eventually to 1000 m depth. The vertical extent of The Blob is essential with regard to its propagation to the Bering Sea via the Aleutian Straits where the shallow depths of some passes will constrain water transport as discussed below; for more information see the latest geomorphological study of the Aleutian passes by Zimmerman and Prescott (2021).

Zhi et al. (2019) have shown that The Blob’s appearance as a temperature anomaly was preceded and actually preconditioned by a large negative salinity anomaly in the central North Pacific that originated in 2010-2011 in the western Pacific. This salinity anomaly played a key role during 2012-2013 in the observed formation of a shallow mixed layer in the NE Pacific and emergence of a large SST anomaly (The Blob) by trapping heat in the shallow mixed layer. The shallow mixed layer anomaly was more dependent on the local subsurface salinity anomaly in the 100-150 m layer, which was much stronger than the surface salinity anomaly. Zhi et al. (2019) have shown that the salinity anomaly in the 100-150 m layer originated west of the dateline, where it was freshened, subducted, and advected eastward.

The propagation path of The Blob across the North Pacific tentatively traced by Zhi et al. (2019) lies near the long-term mean path of the Polar Front (PF) mapped by Belkin et al. (2002) from >200 hydrographic sections (Figure 5); therefore, The Blob’s propagation to the northeast might have been facilitated by a geostrophic current associated with the PF.

Advection of SST anomalies along the PF was studied by Belkin and Shotwell (2012) and Shotwell et al. (2014) from Hadley 1-degree SST climatology. The Polar Front Current associated with the PF acts as a conduit of oceanic anomalies originating in the Northwest Pacific. These anomalies propagate along the PF into the GOA where the PF retroflects and extends westward along the GOA shelf towards the Aleutians and into the Bering Sea. The Polar Front Current carries relatively warm waters into the GOA and to the Aleutians; therefore intensification (relaxation) of the Polar Front Current results in warming (cooling) of the high-latitude Northeast Pacific. Intensity of along-front advection varies on seasonal, interannual, and decadal scales. The unique succession of annual advective events in 1956-1966 was followed by a 20-year relaxation epoch, except for the 1971-1972 cold anomaly that traveled around the GOA and along the Aleutians. In 1986-1987, a major cold anomaly propagated into the GOA. The eastward advection episodes in 1990-1991 (warm and cold) and 1994 (warm) were followed by a major propagation event in 1996- 1997 that culminated in an exceptionally strong warm event in the GOA in 1997, when SST anomaly south of Kodiak exceeded +2.5°C.

The series of selected monthly SSTA maps presented by Freeland (2015) and Freeland and Ross (2019) reveal a northward expansion of The Blob into the GOA and eastward shift of The Blob towards the west coast of North America (Figure 2), where SSTA exceeded +6°C off Southern California (Gentemann et al., 2017). As a result of The Blob’s shift east and north, its eastern and northern peripheries have become subjected to the advection by the poleward eastern limb of the Alaska Gyre and adjoining poleward coastal currents such as the Vancouver Island Coastal Current, the Haida Current, and the Alaska Coastal Current. Therefore, The Blob’s manifestation in coastal SST records can be expected all the way around the GOA, along the Alaska Peninsula and Aleutian Islands. Indeed, the monthly time series of SST from the Prince William Sound shows a robust cold-to-warm-epoch transition in late 2014 (von Biela et al., 2019, Fig. 3) consistent with the occupation of the Sound by The Blob (Figure 6). Based on the same data, one can argue that the cold epoch of 2006-2013 abruptly ended in June 2013 with the arrival of the Blob’s leading edge. Immediately downstream from the Prince William Sound, along the Seward Line, a warm water temperature anomaly of unprecedented magnitude and duration was recorded in 2013-2018, consistent with The Blob’s advection by the Alaska Coastal Current (Figure 7).

Our view of The Blob’s advection by coastal currents around the GOA is shared by Walsh et al. (2018, pp. S41-S42): “The warming was primarily confined to the inner GOA shelf in September 2014, suggesting that heat was advected along-shore within the Alaska Coastal Current. By spring 2015 the shelf was uniformly warm and water remained 1°–2°C warmer than normal through September 2016. … The warmth of the Bering Sea in 2016 was unprecedented in the historical record.”

The monthly SSTA maps in Figure 2 show The Blob appearing in the southern Bering Sea in June 2014 (when SSTA reached 3 standard deviations at 180°W), then moving east (likely advected by the Aleutian North Slope Current) to occupy the Bristol Bay in August 2014.

Based on 20+ years of observations from moorings, satellite-tracked drifters, and hydrographic drifters, Stabeno et al. (2016) concluded that there are two principal pathways for the Pacific waters into the Bering Sea: one (major) via Unimak Pass at 164°W, and another (minor) via Amukta Pass at 172°W, and estimated that “A typical transit type from Unimak Pass to Bering Strait is >13 months, and from Amukta Pass to Bering Strait via the Bering Slope Current is >8 months.” (ibid., p. 13). These estimates are essential to the timeline of The Blob’s propagation along the EBSS as The Blob is carried by ocean currents.

The precise timeline of The Blob’s propagation into the Bering Sea and onto the EBSS is especially important since the cold epoch of 2006-2013 ended almost simultaneously with The Blob’s penetration to the Bering Sea and EBSS. The available time series of physical parameters (SST, sea ice extent, bottom temperature, cold pool extent etc.) can be sometimes interpreted variously. Depending on an interpretation, the cold-to-warm-epoch transition in 2013-2014 was either (a) triggered by The Blob, or (b) initiated independently from The Blob and before The Blob’s arrival to the EBS, and then reinforced by The Blob.

The Bering Sea had a long memory of The Blob’s passage. The Blob’s residence time in the EBS exceeded at least twice the advection time over which the Bering Slope Current and EBS shelf currents would carry The Blob to the Bering Strait. Indeed, as estimated by Stabeno et al. (2016), a transit time from Unimak Pass and Amukta Pass to the Bering Strait varies between >13 and >8 months, respectively; for a crude estimate, roughly a year. Yet The Blob persisted in the EBS for at least two and a half years (from mid-2014 through 2016) as evidenced by Figure 8.

The enhanced ocean-atmosphere heat flux and atmospheric circulation over the Bering Sea during these years resulted in an unprecedented loss of seasonal sea ice (Figure 9 and Figure 10).

The time series of annual mean anomalies of SST (SSTA) and bottom temperature (BTA) in the EBSS presented by Baker (2021) reveal a rapid cold-to-warm-epoch transition in 2014-2016, with both SSTA and BTA peaking in 2016 at about 3°C and 2°C, respectively (Figure 11).

Baker (2021) documented a stark contrast between surveys in 2010 and 2017-2018: Not only water temperature had risen sharply but the cold pool which occupied a large area in 2010 had all but disappeared in 2018 (Figure 12).

Thorson (2019) constructed a series of annual maps (1982-2017) of bottom temperature in 200 nodes of a quasi-regular grid covering the EBS (Figure 13). These maps show large cold pools during the cold epoch of 2006-2013, then abrupt shrinking of the cold pool in 2014-2015, and a nearly complete disappearance of the cold pool in 2016. Duffy-Anderson et al. (2019) have published yet another map of the cold pool in 2018 (Figure 14) that shows its complete disappearance on the EBSS. The maps of bottom temperature on the EBSS published by Duffy-Anderson et al. (2019), Thorson (2019), and Baker (2021) make evident the disappearance of the cold pool in 2016, its partial recovery in 2017, and disappearance again in 2018.

Johnson et al. (2022) conducted the first long-term monitoring study using acoustics and ancillary oceanographic data (bottom temperature and regional seasonal sea ice (SSI) data) to study the relationship between biological scattering communities and variations of the EBS cold pool across two climate regimes, cold (2006-2013) and warm (2014-2018). They concluded (ibid., p. 210): “In the fall of 2013, the EBS encountered a regime shift that was characterized by less SSI, early receding SSI, increased bottom temperatures, and a zooplankton community structure shift. … The most recent warm regime years have resulted in little to no cold pool formation along the EBS shelf, and while changes have occurred in zooplankton communities between regimes, cold pool dynamics during warm regimes are yet to be fully understood, and perhaps their impact on secondary producer communities may not be felt immediately and could have longer term effects.”

Changes of the entire Bering Shelf over seasonal to century-long time scales have been recently assessed by Danielson et al. (2020), who assembled and analyzed a new temperature-salinity climatology of the Bering and Chukchi Sea shelves from a half-century (1966-2018) hydrographic data base. The Bering Shelf is dominated by decadal-scale variability that precludes detection of a water column trend. However, SST shows a warming trend of 0.22 plus/minus 0.10°C/decade. Heat fluxes in 1979-2018 exhibit no record-length trend over the Bering Sea shelf.

The Bering Sea experienced record-high SST in 2014 (Stabeno et al., 2017), which persisted into 2018 (Thoman et al., 2020) and 2019 (Stabeno and Bell, 2019). The upper ocean (0-300 m) integrated heat content anomalies in the EBS were correlated with SST anomalies (Walsh et al., 2018).

Physical fronts have long been recognized as critically important structural and biophysical component of the Bering Sea ecosystem (Belkin and Cornillon, 2005; Hunt et al., 2008; Belkin, 2016), especially over the broad EBSS, where a few fronts of different physical nature co-exist, waxing and waning as the season progresses (Figure 15). These fronts, in particular their characteristic temperature ranges, must have been affected by a major SST anomaly associated with The Blob’s passage. Also, the recent dramatic reduction of the seasonal ice cover must have affected those fronts whose dynamics is related to the seasonal ice edge (marginal ice zone).

4. Ecosystem Responses to The Blob in the Eastern Bering Sea

Detection of population-level responses perturbed by The Blob following its entry into the Bering Sea in July 2014 largely depends on the availability of decadal or longer time series, on the life stage most vulnerable to perturbation, and on the duration of any lag between the perturbed and monitored life stages. Time-series records that pre-date the onset of The Blob in the Bering Sea by at least a few years provide context for evaluating any ensuing changes. Detection of biotic responses may be evident within the same season, or delayed if monitored life stages are months or years older than the age of an organism's initial response. Knowledge of such lags informs expectations of when detectable responses should be evident had a perturbation occurred, and the time scale for when they should subside.

In this section we evaluate published records of population-related time series for monitored biota in the eastern Bering Sea that pre-date 2014 by several years, and for which expected lags between the initial effects of The Blob and its expected subsequent appearance in the monitoring records are well established. Most of these records are summarized in the 2022 Ecosystem Status Report, Eastern Bering Sea, edited by Siddon (2022). We begin with consideration of primary productivity, and progress through successively higher trophic levels.

Primary Productivity

Trends in primary productivity from 2003-2022 have been evaluated from ocean color data provided by the satellite-based Moderate Resolution Imaging Spectroradiometer (MODIS) and averaged over spring (April – June) for each of 8 subregions in the Bering Sea (Figure 16; Nielsen et al. 2022). Spring chlorophyll-a concentrations were consistently below the long-term median after 2014 in the off-shelf (>200 m depth) subregion (Figure 16), which includes most of the highly productive Bering Sea Green Belt (Springer et al. 1996). Similar trends are not evident in any of the other subregions.

Copepods and Euphausiids

The abundance of small (<2 mm) and large (>2 mm) copepods, and euphausiids smaller than 15 mm has been estimated from 1998-2022 along the 70 m isobath in May within the middle domain (50-100 m depth) of the SE Bering Sea, and from 1996-2022 within the middle domain in August- September (Figure 17). Abundances of large copepods sampled in May were low from 2002-2006, a period of relatively warm temperatures and low sea ice cover (Figure 9 and Figure 11), then increased by factors of up to 100 during the subsequent colder years of 2007-2014 (Figure 17A). Abundances then declined after 2014 to numbers typical of the 2002-2006 period, with a rebound in 2017. By August-September, large copepods were generally abundant from 2000-2014 except for 2006, 2007 and 2013, with a generally declining trend after 2014 (Figure 17B). In contrast, a generally declining trend of small copepod abundances in the middle domain from 2004-2013 reversed after 2013 and remained elevated through at least 2020 (Figure 17). Also, euphausiid abundances were generally stable during these periods, apart from isolated peaks in May 2017 and August-September 2003.

Jellyfish

Late-summer biomass of jellyfish was exceptionally high in 2014 and 2016 in the upper 25 m layer in the southeastern Bering Sea, returning to pre-2014 levels by 2018 (Figure 18). These increases were not evident during the previous warm period of 2002-2006. In contrast, jellyfish biomass in the northern Bering Sea appeared unrelated to warm- or cold-water periods from 2004-2022. Jellyfish biomass was not estimated in the southeastern Bering Sea in 2020 or odd-numbered years between 2013-2021 inclusive, or in the northeastern Bering Sea in 2008 and 2020.

Alaska Pollock and Pacific Cod

As the largest fishery by volume in the United States, Alaska pollock (Gadus chalcogrammus) in the eastern Bering Sea have been intensively monitored over the last 50 years. Pollock abundance is evaluated by one of the most advanced stock assessments available, which provides estimates of abundance-at-age dating from the mid-1960s (Ianelli et al. 2021). These results show that pollock recruitment is highly variable, with episodic high recruitment events every 4-5 years sustaining the population over time. The two lowest periods of recruitment to age 1 juveniles since 1990, leading by one year, are associated with first the 2002-2006 period of warm SST and low sea ice cover, and second with entry of The Blob into the eastern Bering Sea (Figure 19A). This association suggests that unusually warm surface waters may have contributed to the low recruitments during these two periods. Recruitment to age 1 of Pacific cod (Gadus morhua) responded similarly to the 2002-2006 period of warm SST and low sea ice cover and to The Blob (Figure 19B). Pacific cod also support an important fishery in the eastern Bering Sea.

Pacific Halibut

Pacific halibut (Hippoglossus stenolepis) biomass in the eastern Bering Sea underwent two declines since 1994, both of which were associated with warmer water periods. The first of these declines persisted from 2000-2005, approximately contemporaneous with the 2002-2006 warm water period, and the second began after 2015, following entry of The Blob into the Bering Sea (Figure 20). These biomass estimates are based on weight per unit effort of fish greater than 32 inches (81.6 cm) in fork length caught in fishery-independent setline surveys. The annual mean value of 7.91 lbs/skate in 2022 implies a decline by a factor of nearly 3 from the peak value of 21.4 lbs/skate in 2000.

Red, Snow and Tanner Crabs

Red, Snow and Tanner crabs (Paralithodes camtschaticus, Chionoecetes opilio and C. bairdi, respectively) have all supported major fisheries in the eastern Bering Sea and are all now in a state of decline (Figure 21). Declines of red king crabs began after 2014, becoming severe after 2017. Declines of snow crab males also began after 2014, but females temporarily recovered after 2016 from a decline that began in 2011, reaching their most recent peak in 2018, followed by another decline through 2022. Tanner crabs declined from recent peaks in 2013 (females) or 2014 (males), with females recovering to the time-series mean by 2021.

Salmon

Responses of Pacific salmon to The Blob varied by species. Abundance of juvenile (age 2) Chinook salmon (Oncorhynchus tshawytscha) declined almost monotonically after 2015, but juvenile (ages 2-4) sockeye salmon (O. nerka) and juvenile pink salmon (O. gorbuscha) were relatively abundant during the 2002-2006 water period and after the arrival of The Blob in 2014 in the southeastern Bering Sea (Figure 22). In contrast, trends in the abundances of juvenile chum (O. keta) did not appear related to periods of relatively warm or cool waters in the eastern Bering Sea (Fig. 68 in Murphy et al. 2022).

Seabirds

Reproductive success of five species of piscivorous seabirds on St. George and St. Paul Islands in the southeastern Bering Sea declined nearly in concert following The Blob’s arrival. These five species include common and thick-billed murres (Uria aalge and U. lomvia), black- and red-legged kittiwakes (Rissa tridactyla and R. brevirostris), and red-faced cormorants (Urile urile). Reproductive success of common and thick-billed murres on St. George Island, along with black- and red-legged kittiwakes, and red-faced cormorants on both St. George and St. Paul Islands, declined sharply in 2015, while that of common and thick-billed murres on St. Paul Island declined sharply in 2016 (Figure 23). Reproductive success of all these species recovered by 2019 on both islands.

Range Extensions

Warming bottom waters have led to northwestward expansions or shifts of the distribution of benthic and pelagic fish and shellfish, including commercially important species such as Alaska pollock, Pacific cod, king crab and tanner crab (Landeira et al. 2018, Thorson et al. 2019, Baker 2021, Levine et al. 2023). The cold pool near-disappearance has a major implication for walleye pollock and Pacific cod because these species avoid the cool pool; therefore, they are now expected to colonize the newly available areas of the EBSS that were previously occupied by the cold pool. Baker (2021) documented dramatic spatial shifts of two ecologically important Arctic gadids (polar cod and saffron cod) and two commercially important sub-Arctic gadids (walleye pollock and Pacific cod): While Arctic gadids responded negatively to the warm epoch in the EBSS, the sub-Arctic gadids responded positively to the same changes; moreover, the sub-Arctic gadids now moved massively to the NBSS.

5. Discussion

The multiple anomalous responses among the taxa we surveyed above in our results section suggest that The Blob caused a major disruption of the marine food web in the eastern Bering Sea. These anomalous responses spanned at least four and probably more trophic levels. While the taxa we surveyed do not represent the entire marine food web, they do represent major components of the trophic levels they occupy. Considered individually, there are many factors that, alone or in combination, might account for the population-level responses of the taxa we surveyed. However, taken together, these responses appear to be broadly consistent with a decrease in primary production that consequently reduced the efficiency of energy transfer across trophic levels (i.e., trophic transfer efficiency, or TTF), resulting in declines in abundance, biomass or reproductive success in higher trophic levels.

Evidence for reduction of TTF following entry of The Blob into the Bering Sea begins with the response of primary productivity. Satellite measurements of chlorophyll-a in spring reflect phytoplankton standing stocks which, based on the reasonable assumption that zooplankton grazing on phytoplankton is negligible through the initial stages of the spring phytoplankton bloom, implies a corresponding intensity of phytoplankton production. Off-shelf chlorophyll-a measurements during spring were above average through 2014, with the 2014 measurements completed before The Blob entered the Bering Sea in July of that year. Subsequent off-shelf measurements have confirmed a general decline (Figure 16). These declines imply reduced primary productivity in the Bering Sea Green Belt, the most productive region of the eastern Bering Sea (Springer et al. 1996), which contributes to primary production on the outer continental shelf of the eastern Bering Sea through advection onto the shelf (Hunt et al. 2008). Considered in isolation, this is tenuous support for concluding that primary productivity in the eastern Bering Sea declined substantially following entry of The Blob, but the simultaneous response of zooplankton lends additional support.

The concurrent declines of large copepods and increases of small copepods beginning in 2014 suggest lower primary productivity as production switched from an ice-algal dominated production scheme to one that is dominated by open-water production (Mueter et al. 2021 and references therein). In any case, this transition lengthens the food chain by increasing the mean number of trophic transfers to higher trophic level consumers such as fish and seabirds. Even slight shifts in TTE can lead to large reductions at higher trophic levels (e.g., Eddy et al. 2021). Such reductions are exacerbated by increased biomass of jellyfish (Figure 18), which diverts additional energy from most higher-order consumers.

Reduced energy flow to higher-order consumers is reflected in the reduced recruitment to age 1 of Alaska pollock and Pacific cod (Figure 19), declining biomass of Pacific halibut (Figure 20), declining abundance of red, snow, and tanner crabs and of juvenile Chinook salmon (Figure 21 and Figure 22), and reduced reproductive success of common and thick-billed murres on St. George Island, and of black- and red-legged kittiwakes, and red-faced cormorants (all piscivorous seabirds) on St. George and St, Paul Islands (Figure 23), following entry of The Blob into the eastern Bering Sea. The delayed decline by one year of murres on St. Paul Island may reflect reduced effects of The Blob near St. Paul Island in 2014 compared with those at St. George Island, the latter island being closer to the continental shelf break. Many but not all of these taxa had similar responses to the prior warm-water period of 2002-2006 in the eastern Bering Sea.

Not all of the higher-order consumer responses we surveyed declined. Juvenile sockeye salmon abundances increased considerably during warm years, especially after entry of The Blob into the southeastern Bering Sea, when subsequent sockeye fishery yields nearly doubled compared with the long-term mean (see Fig. 70 in Cunningham et al. 2022). This likely resulted from concurrent survival increases of Alaska pollock in their first year of life. Juvenile sockeye salmon consumed much higher proportions of age-0 Alaska pollock during the 2002 – 2006 warm period and after entry of The Blob into the Bering Sea than they did during the intervening cold period (Yasumiishi et al. 2019, 2020). High consumption rates of age-0 Alaska pollock by pink and chum salmon, and by other forage fish during these warm periods (Yasumiishi et al. 2020) may account at least in part for the relatively high abundance of juvenile pink salmon (Figure 22), the absence of trends in juvenile chum salmon abundance, and the poor recruitment of Alaska pollock to age 2 (Figure 19).

Entry of The Blob into the Bering Sea also appears to have affected spatial distributions of multiple taxa. Warming bottom waters have led to northwestward expansions or shifts of the distribution of benthic and pelagic fish and shellfish, including commercially important species such as Alaska pollock, Pacific cod, king crab and tanner crab (Baker 2021, Landeira et al. 2018, Thorson et al. 2019, Levine et al. 2023). The Blob’s entry into the Bering Sea also seems to contribute to the proliferation of toxic or otherwise noxious microalgae. The harmful algal blooms (HABs, such as those causing paralytic shellfish poisoning) appear more frequently and spread more widely throughout the eastern Bering Sea as SST rises (Natsuike et al. 2017).

The population- and community-level responses to The Blob we survey here for the eastern Bering Sea differ substantially from comparable responses in the Gulf of Alaska. The effects of The Blob in the Gulf of Alaska were more intense and persistent, and the population- and community-level responses differed markedly (Suryan et al. 2021). In the Gulf of Alaska, microzooplankton biomass decreased in spring but increased in fall, and large copepods increased in response to lower chlorophyll during The Blob years. Also during The Blob years, capelin and herring (both forage fishes) declined in the Gulf of Alaska. Seabirds also declined in the Gulf of Alaska during The Blob years, as they did in the eastern Bering Sea. Differences in the responses of these taxa between the Gulf of Alaska and the eastern Bering Sea likely result at least in part from the differences in the ecological structure of these two marine ecosystems. Biological productivity in the Bering Sea is strongly affected by sea ice extent and phenology, whereas the Gulf of Alaska is ice-free. While fishery yields increased with warming in subpolar waters owing to extension of the annual production cycle (Sherman et al. 2011, 2013), it is not clear whether this effect will be more important than the loss of production associated with loss of sea ice as is occurring in the Bering Sea.

6. Summary and Conclusions

We reviewed various physical drivers and biological responses to the 2013-2015 marine heat wave in the NE Pacific. The anomaly dubbed “The Blob” emerged in 2013-2014 and persisted through 2015-2016, with some signs of its lingering through 2017-2018 and a possible return in 2019. A large salinity anomaly in the surface and subsurface layers preceded The Blob and contributed to its formation by enhancing stratification and trapping heat in the shallow mixed layer. The Blob’s successive appearances in the NE Pacific are suggestive of its advection by coastal currents around the Gulf of Alaska, along the Aleutians, and into the Bering Sea, where it was transported by the Bering Slope Current and shelf currents to the Bering Strait. During the Blob’s development in 2013-2014, advection along the Polar Front likely contributed to the Blob’s propagation to the Gulf of Alaska, along the Aleutians, and into the EBS. Owing to its extreme persistence and exceptional magnitude, The Blob and attendant physical changes impacted, sometimes dramatically, various aspects of the EBS ecosystem. Comparison of time-series of population responses across trophic levels suggests that The Blob lowered primary production during spring, increasing small copepod production and jellyfish, and reducing energy transfer efficiently to higher trophic levels. The multi-year duration of The Blob likely preconditioned the EBS for the record low seasonal sea ice extent and disappearance of the cold pool in 2016 and 2018 that profoundly affected zooplankton, invertebrates, fishes, sea birds, and marine mammals. While the Bering Sea’s water temperature, seasonal sea ice, and cold pool seem to have returned to the long-term mean state in 2022, it remains to be seen whether the Bering Sea ecosystem will return to its previous state in the short term. Likely alternatives include either irreversible changes or hysteresis recovery.