Revealing the extent of sea otter impacts on bivalve prey through multi-trophic monitoring and mechanistic models.

Sea otters are apex predators that can exert considerable influence over the nearshore communities they occupy. Since facing near extinction in the early 1900s, sea otters are making a remarkable recovery in Southeast Alaska, particularly in Glacier Bay, the largest protected tidewater glacier fjord in the world. The expansion of sea otters across Glacier Bay offers both a challenge to monitoring and stewardship and an unprecedented opportunity to study the top-down effect of a novel apex predator across a diverse and productive ecosystem. Our goal was to integrate monitoring data across trophic levels, space, and time to quantify and map the predator-prey interaction between sea otters and butter clams Saxidomus gigantea, one of the dominant large bivalves in Glacier Bay and a favoured prey of sea otters. We developed a spatially-referenced mechanistic differential equation model of butter clam dynamics that combined both environmental drivers of local population growth and estimates of otter abundance from aerial survey data. We embedded this model in a Bayesian statistical framework and fit it to clam survey data from 43 intertidal and subtidal sites across Glacier Bay. Prior to substantial sea otter expansion, we found that butter clam density was structured by an environmental gradient driven by distance from glacier (represented by latitude) and a quadratic effect of current speed. Estimates of sea otter attack rate revealed spatial heterogeneity in sea otter impacts and a negative relationship with local shoreline complexity. Sea otter exploitation of productive butter clam habitat substantially reduced the abundance and altered the distribution of butter clams across Glacier Bay, with potential cascading consequences for nearshore community structure and function. Spatial variation in estimated sea otter predation processes further suggests that community context and local environmental conditions mediate the top-down influence of sea otters on a given prey. Overall, our framework provides high-resolution insights about the interaction among components of this food web and could be applied to a variety of other systems involving invasive species, epidemiology or migration.

Ecosystem response persists after a prolonged marine heatwave.

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.

Keystone predators govern the pathway and pace of climate impacts in a subarctic marine ecosystem.

Predator loss and climate change are hallmarks of the Anthropocene yet their interactive effects are largely unknown. Here, we show that massive calcareous reefs, built slowly by the alga Clathromorphum nereostratum over centuries to millennia, are now declining because of the emerging interplay between these two processes. Such reefs, the structural base of Aleutian kelp forests, are rapidly eroding because of overgrazing by herbivores. Historical reconstructions and experiments reveal that overgrazing was initiated by the loss of sea otters, Enhydra lutris (which gave rise to herbivores capable of causing bioerosion), and then accelerated with ocean warming and acidification (which increased per capita lethal grazing by 34 to 60% compared with preindustrial times). Thus, keystone predators can mediate the ways in which climate effects emerge in nature and the pace with which they alter ecosystems.

Size, growth, and density data for shallow-water sea urchins from Mexico to the Aleutian Islands, Alaska, 1956-2016.

Size, growth, and density have been studied for North American Pacific coast sea urchins Strongylocentrotus purpuratus, S. droebachiensis, S. polyacanthus, Mesocentrotus (Strongylocentrotus) franciscanus, Lytechinus pictus, Centrostephanus coronatus, and Arbacia stellata by various workers at diverse sites and for varying lengths of time from 1956 to present. Numerous peer-reviewed publications have used some of these data but some data have appeared only in graduate theses or the gray literature. There also are data that have never appeared outside original data sheets. Motivation for studies has included fisheries management and environmental monitoring of sewer and power plant outfalls as well as changes associated with disease epidemics. Studies also have focused on kelp restoration, community effects of sea otters, basic sea urchin biology, and monitoring. The data sets presented here are a historical record of size, density, and growth for a common group of marine invertebrates in intertidal and nearshore environments that can be used to test hypotheses concerning future changes associated with fisheries practices, shifts of predator distributions, climate and ecosystem changes, and ocean acidification along the Pacific Coast of North America and islands of the north Pacific. No copyright restrictions apply. Please credit this paper when using the data.