Some of the longest and most comprehensive marine ecosystem monitoring programs were established in the Gulf of Alaska following the environmental disaster of the Exxon Valdez oil spill over 30 years ago. These monitoring programs have been successful in assessing recovery from oil spill impacts, and their continuation decades later has now provided an unparalleled assessment of ecosystem responses to another newly emerging global threat, marine heatwaves. The 2014–2016 northeast Pacific marine heatwave (PMH) in the Gulf of Alaska was the longest lasting heatwave globally over the past decade, with some cooling, but also continued warm conditions through 2019. Our analysis of 187 time series from primary production to commercial fisheries and nearshore intertidal to offshore oceanic domains demonstrate abrupt changes across trophic levels, with many responses persisting up to at least 5 years after the onset of the heatwave. Furthermore, our suite of metrics showed novel community-level groupings relative to at least a decade prior to the heatwave. Given anticipated increases in marine heatwaves under current climate projections, it remains uncertain when or if the Gulf of Alaska ecosystem will return to a pre-PMH state.

Although climate‐induced shifts in fish distribution have been widely reported at the population level, studies that account for ontogenetic shifts and subregional differences when assessing responses are rare.In this study, groundfish distributional changes in depth, latitude, and longitude were assessed at different size classes by species within nine subregions. We examined large, quality‐controlled datasets of depth‐stratified‐random bottom trawl surveys conducted during summer in three large regions—the Gulf of Alaska and the west coasts of Canada and the United States—over the period 1996–2015, a time period punctuated by a marine “heat wave.” Temporal biases in bottom temperature were minimized by subdividing each region into three subregions, each with short‐duration surveys. Near‐bottom temperatures, weighted by stratum area, were unsynchronized across subregions and exhibited varying subregional interannual variability. The weighted mean bottom depths in the subregions also vary largely among subregions. The centroids (centers of gravity) of groundfish distribution were weighted with catch per unit effort and stratum area for 10 commercially important groundfish species by size class and subregion. Our multivariate analyses showed that there were significant differences in aggregate fish movement responses to warm temperatures across subregions but not among species or sizes. Groundfish demonstrated poleward responses to warming temperatures only in a few subregions and moved shallower or deeper to seek colder waters. The temperature responses of groundfish depended on where they were. Under global warming, groundfish may form geographically distinct thermal ecoregions along the northeast Pacific shelf. Shallow‐depth species exhibited greatly different distributional responses to temperature changes across subregions while deep‐depth species of different subregions tend to have relatively similar temperature responses. Future climate studies would benefit by considering fish distributions on small subregional scales.