Spatial Distribution and Temporal Patterns of Cassin's Auklet Foraging and Their Euphausiid Prey in a Variable Ocean Environment

Temporal changes in the diet composition and trophic level of walleye pollock (Gadus chalcogrammus) inhabiting the middle-eastern coast of Korea

The influences of oceanographic changes on diet composition and trophic level for pollock (Gadus chalcogrammus) inhabiting the East Sea off the Korean coast were examined based on stomach content and stable isotope analyses during 2016 and 2017. The diets of pollock consisted mainly of benthic crustaceans (particularly carid shrimps and euphausiids) and cephalopods, with a predominance of teleosts in the diets of larger individuals in deeper habitats. In 2016, amphipods, carid shrimps and cephalopods featured strongly in pollock diets, and the contribution of amphipods decreased in the diets of larger individuals and deeper depths. In 2017, euphausiids dominated at shallower depths, whereas the contributions of carid shrimps and teleosts increased in deeper habitats. Body-size-related differences in carbon stable isotope (δ13C) values were present in both 2016 and 2017, but size-related differences in nitrogen stable isotope (δ15N) values were only observed in 2017. The increased contribution of euphausiids during 2017 resulted in a distinct decrease in the trophic level of pollock compared to co-occurring higher trophic level predators, which can be linked to changes in habitat water temperature. Combined stomach contents and isotopic analyses provide a more comprehensive understanding of how fish diets and trophic levels fluctuate with changes in the type and abundance of prey resources in response to environmental changes.

Submesoscale coupling of krill and whales revealed by aggregative Lagrangian coherent structures

In the marine environment, dynamic physical processes shape biological productivity and predator-prey interactions across multiple scales. Identifying pathways of physical-biological coupling is fundamental to understand the functioning of marine ecosystems yet it is challenging because the interactions are difficult to measure. We examined submesoscale (less than 100 km) surface current features using remote sensing techniques alongside ship-based surveys of krill and baleen whale distributions in the California Current System. We found that aggregative surface current features, represented by Lagrangian coherent structures (LCS) integrated over temporal scales between 2 and 10 days, were associated with increased (a) krill density (up to 2.6 times more dense), (b) baleen whale presence (up to 8.3 times more likely) and (c) subsurface seawater density (at depths up to 10 m). The link between physical oceanography, krill density and krill-predator distributions suggests that LCS are important features that drive the flux of energy and nutrients across trophic levels. Our results may help inform dynamic management strategies aimed at reducing large whales ship strikes and help assess the potential impacts of environmental change on this critical ecosystem.

Factors affecting the seasonal distribution and biomass of E. pacifica and T. spinifera along the Pacific coast of Canada: A spatiotemporal modelling approach

Euphausiids are a keystone species in coastal food webs due to their high lipid content and seasonally high biomass. Understanding the habitat and environmental drivers that lead to areas of high biomass, or ‘hotspots’, and their seasonal persistence, will support the identification of important foraging regions for mid- and upper- trophic level predators. We quantify the distribution of hotspots of the two dominant species of euphausiid in the north-east Pacific Ocean: Euphausia pacifica and Thysanoessa spinifera, as well as euphausiid larvae (mixed species). The Canadian coast encompasses the northern California Current Ecosystem and the transition zone to the Alaska current, and is a highly productive region for fisheries, marine mammals, and seabirds. We used spatiotemporal modelling to predict the distribution of these three euphausiid groups in relation to geomorphic and environmental variables during the important spring-summer months (April through September) when euphausiid biomass is highest. We quantified the area, intensity, and persistence of biomass hotspots across months according to specific oceanographic ecosections developed for marine spatial planning purposes. Persistent hotspots of both adult species were predicted to occur along the 200 m depth contour of the continental slope; however, differences were predicted on the shallower Dixon shelf, which was a key area for T. spinifera, and within the Juan de Fuca Eddy system where E. pacifica hotspots occurred. The continental slope along the west coast of Vancouver Island was the only persistent hotspot region common between both adult species and euphausiid larvae. Larval distribution was more correlated with T. spinifera than E. pacifica biomass. Hotspots of adults were more persistent across months than hotspots of euphausiid larvae, which were seasonally patchy. The persistence of biomass hotspots of forage species through periods of low overall biomass could maintain trophic connectivity through perturbation events and increase ecosystem resilience to climate change.

Modeling predator and prey hotspots: Management implications of baleen whale co-occurrence with krill in Central California

As global ocean-bound commerce increases, managing human activities has become important in reducing conflict with threatened wildlife. This study investigates environmental factors determining abundance and distribution of blue whales (Balaenoptera musculus), humpback whales (Megaptera novaeangliae) and their prey (Euphausia pacifica and Thysanoessa spinifera) in central California. We provide insights into environmental drivers of the ecology and distribution of these species, model whale distributions and determine coincident hotspots of whales and their prey that will help decrease human threats to whales and protect critical feeding habitat. We developed separate predictive models of whale abundances (using negative binomial regression on count data) and krill abundance (using a two-part hurdlemodel combining logistic and negative binomial regressions) over a 14 year period (2004-2017). Variables included in situ surface and midwater oceanographic measures (temperature, salinity, and fluorescence), basin-scale climate indices, and bathymetric- and distance-related data. Predictions were applied to 1 km2 cells spanning the study area for May, June, July, and September during each of the 14 years of surveys to identify persistent distribution patterns. Both whales and krill were found to consistently use the northeast region of Cordell Bank, the Farallon Escarpment, and the shelf-break waters. The main identified blue whale hotspots were also krill hotspots, while co-occurrence was more limited and varied seasonally for humpback whales and krill. These results are valuable in identifying patterns in important areas of ecological interaction to assist management of whales. Areas north of Cordell Bank are of particular management concern since they overlap with the end of the San Francisco Bay northern shipping lane. Our findings can help decrease threats to whales, particularly in important foraging areas, by supporting implementation of vessel management and informing potential conflicts with other human uses.

Oceanographic drivers of winter habitat use in Cassin's Auklets

Reduced prey abundance and severe weather can lead to a greater risk of mortality for seabirds during the non-breeding winter months. Resource patterns in some regions are shifting and becoming more variable in relation to past conditions, potentially further impacting survival and carryover to the breeding season. As animal tracking technologies and methods to analyze movement data have advanced, it has become increasingly feasible to draw fine-scale inference about how environmental variation affects foraging behavior and habitat use of seabirds during this critical period. Here, we used archival light-sensing tags to evaluate how interannual variation in oceanography affected the winter distribution of Cassin's Auklets from Southeast Farallon Island, California. Thirty-five out of 93 geolocators deployed from 2015 to 2017 were recovered and successfully recorded light-level data, from which geographic positions were estimated. Step-selection functions were applied to identify environmental covariates that best explained winter movement decisions and habitat use, revealing Cassin's Auklets dispersed farther from the colony during a winter with warm SST anomalies, but remained more centralized near the breeding colony during two average winters. Movement patterns were driven by avoidance of areas with higher sea surface temperatures and possible limits of dispersal from the breeding colony, and selection for areas with well-defined mesoscale fronts and cooler surface waters. Through multiple years of tagging and the application of step-selection functions, a robust and widely applied approach for analyzing animal movement in terrestrial species, we show how interannual differences in the movement patterns of a small seabird are driven by oceanographic variability across years. Understanding the winter habitat use of seabirds can help inform changes in population structure and measures of reproductive success, aiding managers in determining potential causes of breeding failures.

Modeling Nonresident Seabird Foraging Distributions to Inform Ocean Zoning in Central California

Seabird aggregations at sea have been shown to be associated with concentrations of prey. Previous research identified Central California as a highly used foraging area for seabirds, with locally breeding seabirds foraging close to their colonies on Southeast Farallon Island. Herein, we focus on nonresident (i.e. non-locally breeding) seabird species off of Central California. We hypothesized that high-use foraging areas for nonresident seabirds would be influenced by oceanographic and bathymetric factors and that spatial and temporal distributions would be similar within planktivorous and generalist foraging guilds but would differ between them. With data collected by the Applied California Current Ecosystem Studies (ACCESS) partnership during cruises between April and October from 2004–2013, we developed generalized linear models to identify high-use foraging areas for each of six nonresident seabird species. The four generalist species are Phoebastria nigripes (black-footed albatross), Ardenna griseus (sooty shearwater), Ardenna creatopus (pink-footed shearwater), and Fulmarus glacialis (northern fulmar). The two planktivorous species are Phalaropus lobatus (red-necked phalarope) and Phalaropus fulicarius (red phalarope). Sea surface temperature was significant for generalist species and sea surface salinity was important for planktivorous species. The distance to the 200-m isobath was significant in five of six models, Pacific Decadal Oscillation with a 3-month lag in four models, and sea surface fluorescence, the distance to Cordell Bank, and depth in three models. We did not find statistically significant differences between distributions of individual seabird species within a foraging guild or between guilds, with the exception of the sooty shearwater. Model results for a multi-use seabird foraging area highlighted the continental shelf break, particularly within the vicinity of Cordell Bank, as the highest use areas as did Marxan prioritization. Our research methods can be implemented elsewhere to identify critical habitat that needs protection as human development pressures continue to expand to the ocean.

Placing Local Aggregations in a Larger-Scale Context: Hierarchical Modeling of Black-Footed Albatross Dispersion

At-sea surveys facilitate the study of the distribution and abundance of marine birds along standardized transects, in relation to changes in the local environmental conditions and large-scale oceanographic forcing. We analyzed the form and the intensity of black-footed albatross (Phoebastria nigripes: BFAL) spatial dispersion off central California, using five years (2004-2008) of vessel-based surveys of seven replicated survey lines. We related BFAL patchiness to local, regional and basin-wide oceanographic variability using two complementary approaches: a hypothesis-based model and an exploratory analysis. The former tested the strength and sign of hypothesized BFAL responses to environmental variability, within a hierarchical atmosphere-ocean context. The latter explored BFAL cross-correlations with atmospheric / oceanographic variables. While albatross dispersion was not significantly explained by the hierarchical model, the exploratory analysis revealed that aggregations were influenced by static (latitude, depth) and dynamic (wind speed, upwelling) environmental variables. Moreover, the largest BFAL patches occurred along the survey lines with the highest densities, and in association with shallow banks. In turn, the highest BFAL densities occurred during periods of negative Pacific Decadal Oscillation index values and low atmospheric pressure. The exploratory analyses suggest that BFAL dispersion is influenced by basin-wide, regional-scale and local environmental variability. Furthermore, the hypothesis-based model highlights that BFAL do not respond to oceanographic variability in a hierarchical fashion. Instead, their distributions shift more strongly in response to large-scale ocean-atmosphere forcing. Thus, interpreting local changes in BFAL abundance and dispersion requires considering diverse environmental forcing operating at multiple scales.