Detecting and inferring cause of change in an Alaska nearshore marine ecosystem

Vertical Distribution of Rocky Intertidal Organisms Shifts With Sea-Level Variability on the Northeast Pacific Coast

Disentangling the effects of cyclical variability in environmental forcing and long-term climate change on natural communities is a major challenge for ecologists, managers, and policy makers across ecosystems. Here we examined whether the vertical distribution of rocky intertidal taxa has shifted with sea-level variability occurring at multiple temporal scales and/or long-term anthropogenic sea-level rise (SLR). Because of the distinct zonation characteristic of intertidal communities, any shift in tidal dynamics or average sea level is expected to have large impacts on community structure and function. We found that across the Northeast Pacific Coast (NPC), sea level exhibits cyclical seasonal variability, tidal amplitude exhibits ecologically significant variability coherent with the 18.6-year periodicity of lunar declination, and long-term sea-level rise is occurring. Intertidal taxa largely do not exhibit significant vertical distribution shifts coherent with short-term (monthly to annual) sea-level variability but do exhibit taxa-specific vertical distribution shifts coherent with cyclical changes in lunar declination and long-term SLR at decadal timescales. Finally, our results show that responses to cyclical celestial mechanics and SLR vary among taxa, primarily according to their vertical distribution. Long-term SLR is occurring on ecologically relevant scales, but the confounding effects of cyclical celestial mechanics make interpreting shifts in zonation or community structure challenging. Such cyclical dynamics alternatingly amplify and dampen long-term SLR impacts and may modify the impacts of other global change related stressors, such as extreme heat waves and swell events, on intertidal organisms living at the edge of their physiological tolerances. As a result, intertidal communities will likely experience cyclical periods of environmental stress and concomitant nonlinear shifts in structure and function as long-term climate change continues. Our results demonstrate that consistent, large-scale monitoring of marine ecosystems is critical for understanding natural variability in communities and documenting long-term change.

Macronutrient composition of sea otter diet with respect to recolonization, life history, and season in southern Southeast Alaska

The sea otter (Enhydra lutris) population of Southeast Alaska has been growing at a higher rate than other regions along the Pacific coast. While good for the recovery of this endangered species, rapid population growth of this apex predator can create a human-wildlife conflict, negatively impacting commercial and subsistence fishing. Previous foraging studies throughout the sea otter range have shown they will reduce invertebrate prey biomass when recolonizing an area. The goal of this study was to examine and quantify the energy content of sea otter diets through direct foraging observations and prey collection. Our study area, Prince of Wales Island in southern Southeast Alaska, exhibits a gradient of sea otter recolonization, thus providing a natural experiment to test diet change in regions with different recolonization histories. Sea otter prey items were collected in three seasons (spring, summer, and winter) to measure caloric value and lipid and protein content. We observed 3523 sea otter dives during the spring and summer. A majority of the sea otter diet consisted of clams. Sea otters in newly recolonized areas had lower diet diversity, higher energetic intake rates (EIR, kcal/min), and prey had higher energy content (kcal/g). Females with pups had the highest diet diversity and the lowest EIR. Sea otter EIR were higher in the fall and winter vs. spring and summer. Sea cucumber energy and lipid content appeared to correspond with times when sea otters consumed the highest proportion of sea cucumbers. These caloric variations are an important component of understanding ecosystem-level effects sea otters have in the nearshore environment.

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.

Brown bear–sea otter interactions along the Katmai coast: terrestrial and nearshore communities linked by predation

Sea otters were extirpated throughout much of their range by the maritime fur trade in the 18th and 19th centuries, including the coast of Katmai National Park and Preserve in southcentral Alaska. Brown bears are an important component of the Katmai ecosystem where they are the focus of a thriving ecotourism bear-viewing industry as they forage in sedge meadows and dig clams in the extensive tidal flats that exist there. Sea otters began reoccupying Katmai in the 1970s where their use of intertidal clam resources overlapped that of brown bears. By 2008, the Katmai sea otter population had grown to an estimated 7,000 animals and was likely near carrying capacity; however, in 2006–2015, the age-at-death distribution (AADD) of sea otter carcasses collected at Katmai included a higher-than-expected proportion of prime-age animals compared to most other sea otter populations in Alaska. The unusual AADD warranted scientific investigation, particularly because the Katmai population is part of the Threatened southwest sea otter stock. Brown bears in Katmai are known to prey on marine mammals and sea otters, but depredation rates are unknown; thus, we investigated carnivore predation, especially by brown bears, as a potential explanation for abnormally high prime-age otter mortality. We installed camera traps at two island-based marine mammal haulout sites within Katmai to gather direct evidence that brown bears prey on seals and sea otters. Over a period of two summers, we gathered photo evidence of brown bears making 22 attempts to prey on sea otters of which nine (41%) were successful and 12 attempts to prey on harbor seals of which one (8%) was successful. We also developed a population model based on the AADD to determine if the living population is declining, as suggested by the high proportion of prime-age animals in the AADD. We found that the population trend predicted by the modeled AADDs was contradictory to aerial population surveys that indicated the population was not in steep decline but was consistent with otter predation. Future work should focus on the direct and indirect effects these top-level predators have on each other and the coastal community that connects them.

Gene Expression Profiles in Two Razor Clam Populations: Discerning Drivers of Population Status

With rapidly changing marine ecosystems, shifts in abundance and distribution are being documented for a variety of intertidal species. We examined two adjacent populations of Pacific razor clams (Siliqua patula) in lower Cook Inlet, Alaska. One population (east) supported a sport and personal use fishery, but this has been closed since 2015 due to declines in abundance, and the second population (west) continues to support commercial and sport fisheries. We used gene expression to investigate potential causes of the east side decline, comparing razor clam physiological responses between east and west Cook Inlet. The target gene profile used was developed for razor clam populations in Alaska based on physiological responses to environmental stressors. In this study, we identified no differences of gene expression between east and west populations, leading to two potential conclusions: (1) differences in factors capable of influencing physiology exist between the east and west and are sufficient to influence razor clam populations but are not detected by the genes in our panel, or (2) physiological processes do not account for the differences in abundance, and other factors such as predation or changes in habitat may be impacting the east Cook Inlet population.

Recursive Bayesian computation facilitates adaptive optimal design in ecological studies

Optimal design procedures provide a framework to leverage the learning generated by ecological models to flexibly and efficiently deploy future monitoring efforts. At the same time, Bayesian hierarchical models have become widespread in ecology and offer a rich set of tools for ecological learning and inference. However, coupling these methods with an optimal design framework can become computationally intractable. Recursive Bayesian computation offers a way to substantially reduce this computational burden, making optimal design accessible for modern Bayesian ecological models. We demonstrate the application of so-called prior-proposal recursive Bayes to optimal design using a simulated data binary regression and the real-world example of monitoring and modeling sea otters in Glacier Bay, Alaska. These examples highlight the computational gains offered by recursive Bayesian methods and the tighter fusion of monitoring and science that those computational gains enable.

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.

Spatio-temporal patterns of northern gannet abundance in a migratory and wintering area

In marine ecosystems, seabird populations are well monitored, thus allowing their use as indicators of system fluctuations at multiple spatio-temporal scales. Population abundance estimates are essential features of any conservation and management measures and initiatives. Population abundances can be used to delineate the distribution range and foraging grounds of species during both breeding and non-breeding periods, with multi-annual monitoring allowing for the inspection of the temporal variability within key marine areas. Taking advantage of long-term monitoring schemes, we examined the annual abundance patterns of the northern gannet Morus bassanus in its southern European migratory flyway. Here, the presence of a topographical feature (i.e. a submarine canyon system) could influence physical processes (e.g. upwelling, alongshore currents, and riverine inputs), thus oceanographically differentiating canyon and shelf ecosystems within a spatially restricted marine area. We assessed seasonal and long-term trend fluctuations of monthly northern gannet abundance using Generalized Additive Mixed Models, yielding only a strong seasonal effect. Moreover, we jointly tested the effect of the phenology and the spatial domain (canyon versus shelf) on northern gannet abundance using Generalised Linear Mixed Models, accounting for the excess of zeros. Northern gannet abundance was higher during the pre-winter and post-winter migration, corresponding to its southward and northward movements in the NE Atlantic, respectively. The effect of the canyon was only detected during the post-winter migration. By age, adult northern gannets were more abundant in some periods of the year (from October to March, including the pre-winter migration, the wintering and the post-winter migration), whereas juveniles and immatures were especially abundant during the pre-winter period. Our results could provide quantitative baseline information of northern gannet spatio-temporal patterns in a migratory and wintering area important for both research and conservation efforts, given the high conservation value of the area within the Natura 2000 network.

Physiological and gene transcription assays to assess responses of mussels to environmental changes

Coastal regions worldwide face increasing management concerns due to natural and anthropogenic forces that have the potential to significantly degrade nearshore marine resources. The goal of our study was to develop and test a monitoring strategy for nearshore marine ecosystems in remote areas that are not readily accessible for sampling. Mussel species have been used extensively to assess ecosystem vulnerability to multiple, interacting stressors. We sampled bay mussels (Mytilus trossulus) in 2015 and 2016 from six intertidal sites in Lake Clark and Katmai National Parks and Preserves, in south-central Alaska. Reference ranges for physiological assays and gene transcription were determined for use in future assessment efforts. Both techniques identified differences among sites, suggesting influences of both large-scale and local environmental factors and underscoring the value of this combined approach to ecosystem health monitoring.

Lactation and resource limitation affect stress responses, thyroid hormones, immune function, and antioxidant capacity of sea otters (Enhydra lutris)

Lactation is the most energetically demanding stage of reproduction in female mammals. Increased energetic allocation toward current reproduction may result in fitness costs, although the mechanisms underlying these trade-offs are not well understood. Trade-offs during lactation may include reduced energetic allocation to cellular maintenance, immune response, and survival and may be influenced by resource limitation. As the smallest marine mammal, sea otters (Enhydra lutris) have the highest mass-specific metabolic rate necessitating substantial energetic requirements for survival. To provide the increased energy needed for lactation, female sea otters significantly increase foraging effort, especially during late-lactation. Caloric insufficiency during lactation is reflected in the high numbers of maternal deaths due to End-Lactation Syndrome in the California subpopulation. We investigated the effects of lactation and resource limitation on maternal stress responses, metabolic regulation, immune function, and antioxidant capacity in two subspecies of wild sea otters (northern: E. l. nereis and southern: E. l. kenyoni) within the California, Washington, and Alaska subpopulations. Lactation and resource limitation were associated with reduced glucocorticoid responses to acute capture stress. Corticosterone release was lower in lactating otters. Cortisol release was lower under resource limitation and suppression during lactation was only evident under resource limitation. Lactation and resource limitation were associated with alterations in thyroid hormones. Immune responses and total antioxidant capacity were not reduced by lactation or resource limitation. Southern sea otters exhibited higher concentrations of antioxidants, immunoglobulins, and thyroid hormones than northern sea otters. These data provide evidence for allocation trade-offs during reproduction and in response to nutrient limitation but suggest self-maintenance of immune function and antioxidant defenses despite energetic constraints. Income-breeding strategists may be especially vulnerable to the consequences of stress and modulation of thyroid function when food resources are insufficient to support successful reproduction and may come at a cost to survival, and thereby influence population trends.