Natal Origin and Trans‐Pacific Migration of Pacific Bluefin Tuna ( Thunnus orientalis ) Inferred From Otolith Oxygen and Carbon Isotope Analysis

Geochemical chronologies were constructed from otoliths of adult Pacific bluefin tuna (PBT) to investigate the timing of age-specific egress of juveniles from coastal nurseries in the East China Sea or Sea of Japan to offshore waters of the Pacific Ocean. Element:Ca chronologies were developed for otolith Li, Mg, Mn, Zn, Sr, and Ba, and our assessment focused on the section of the otolith corresponding to the age-0 to age-1 + interval. Next, we applied a common time-series approach to geochemical profiles to identify divergences presumably linked to inshore-offshore migrations. Conspicuous geochemical shifts were detected during the juvenile interval for Mg:Ca, Mn:Ca, and Sr:Ca that were indicative of coastal-offshore transitions or egress generally occurring for individuals approximately 4–6 mo. old, with later departures (6 mo. or older) linked to overwintering being more limited. Changepoints in otolith Ba:Ca profiles were most common in the early age-1 period (ca. 12–16 mo.) and appear associated with entry into upwelling areas such as the California Current Large Marine Ecosystem following trans-Pacific migrations. Natal origin of PBT was also predicted using the early life portion of geochemical profile in relation to a baseline sample comprised of age-0 PBT from the two primary spawning areas in the East China Sea and Sea of Japan. Mixed-stock analysis indicated that the majority (66%) of adult PBT in our sample originated from the East China Sea, but individuals of Sea of Japan origin were also detected in the Ryukyu Archipelago.

Natal origin of subadult (age-1) Pacific bluefin tuna (PBT, Thunnus orientalis) from the California Current Large Marine Ecosystem (CCLME) was determined using natural tracers in ear stones (otoliths). Age-0 PBT collected from the two known spawning areas in the western Pacific Ocean (East China Sea, Sea of Japan) were used to establish baseline signatures from otolith cores over 4 years (2014–2017) based on a suite of trace elements (Li, Mg, Mn, Sr, Zn and Ba). Distinct chemical signatures existed in the otolith cores of age-0 PBT collected from the two spawning areas, with overall classification accuracy ranging 73–93% by year. Subadult PBT collected in the CCLME over the following 4 years (2015–2018) were then age-class matched to baselines using mixed-stock analysis. Natal origin of trans-Pacific migrants in the CCLME ranged 43–78% from the East China Sea and 22–57% from the Sea of Japan, highlighting the importance of both spawning areas for PBT in the CCLME. This study provides the first estimates on the natal origin of subadult PBT in this ecosystem using otolith chemistry and expands upon the application of these natural tracers for population connectivity studies for this species.

Juvenile Pacific bluefin tuna (PBT, Thunnus orientalis) are known to migrate from western Pacific spawning grounds to their eastern Pacific nursery and feeding grounds in the California Current Large Marine Ecosystem (CCLME), but the timing, durations, and fraction of the population that makes these migrations need to be better understood for improved management. To complement recent work focused on stable isotope and radiotracer approaches (“tracer toolbox”; Madigan et al., 2014) we explored the suitability of combining longitudinal analyses of otolith microstructure and trace elemental composition in age ∼1–2 PBT (n = 24, 66–76 cm curved fork length) for inferring the arrival of individuals in the CCLME. Element:Ca ratios in transverse otolith sections (9–12 rows, triplicate ablations from primordium to edge, ø50 μm) were quantified for eight elements: Li, Mg, Mn, Co, Cu, Zn, Sr, and Ba, which was followed by microstructure analysis to provide age estimates corresponding to each ablation spot. Age estimates from otoliths ranged from 328 to 498 d post-hatch. The combined elemental signatures of four elements (Ba, Mg, Co, Cu) showed a significant increase at the otolith edge in approximately half of the individuals (30–60 d before catch). Given the different oceanographic properties of oligotrophic open Pacific vs. high nutrient, upwelling CCLME waters, this signal is consistent with the entry of the fish into the CCLME, which was estimated to occur primarily in July after a transoceanic migration of ∼1.5–2.0 months. Our approach comprises a useful addition to the available tracer toolbox and can provide additional and complementary understanding of trans-Pacific migration patterns in PBT.

The detection of Fukushima-derived radionuclides in Pacific bluefin tuna (PBFT) that crossed the Pacific Ocean to the California Current Large Marine Ecosystem (CCLME) in 2011 presented the potential to use radiocesium as a tracer in highly migratory species. This tracer requires that all western Pacific Ocean emigrants acquire the (134)Cs signal, a radioisotope undetectable in Pacific biota prior to the Fukushima accident in 2011. We tested the efficacy of the radiocesium tracer by measuring (134)Cs and (137)Cs in PBFT (n = 50) caught in the CCLME in 2012, more than a year after the Fukushima accident. All small PBFT (n = 28; recent migrants from Japan) had (134)Cs (0.7 ± 0.2 Bq kg(-1)) and elevated (137)Cs (2.0 ± 0.5 Bq kg(-1)) in their white muscle tissue. Most larger, older fish (n = 22) had no (134)Cs and only background levels of (137)Cs, showing that one year in the CCLME is sufficient for (134)Cs and (137)Cs values in PBFT to reach pre-Fukushima levels. Radiocesium concentrations in 2012 PBFT were less than half those from 2011 and well below safety guidelines for public health. Detection of (134)Cs in all recent migrant PBFT supports the use of radiocesium as a tracer in migratory animals in 2012.

Contribution rates of different spawning and feeding grounds to adult Pacific bluefin tuna (Thunnus orientalis) in the northwestern Pacific Ocean

Plastics in the Pacific: Assessing risk from ocean debris for marine birds in the California Current Large Marine Ecosystem

: While specifically focusing on trophic interactions affecting habitat utilization and foraging pattern of California sea lions (CSL) in the California Current Large Marine Ecosystem (CCLME), the long-term goal of our modeling approach is to better understand and characterize biological hotspots (i.e., the aggregation of multiple marine organisms over multiple trophic levels) off the U.S. west coast and in other regions where similar fully-coupled ecosystem models may be implemented. As such, our research represents a major step towards a predictive model that can provide fundamental knowledge about: (1) the spatial and temporal distribution of key marine organisms over multiple trophic levels, and (2) natural and anthropogenic variability in ecosystem structure and trophic interactions. The main research objective is to quantify habitat utilization and trophic interactions in the CCLME by considering patterns of covariability between environmental variables (e.g., temperature, primary production) and foraging patterns and success of middle (forage fish) and higher (sea lions) trophic level organisms. Since our numerical experiments are designed to isolate patterns of variability on seasonal to interannual timescales during normal and extreme years, we focus our analysis on identifying shifts in habitat utilization (e.g., shelf vs. offshore foraging) in the CCLME. We will also explore which features and environmental properties control foraging success in different sub-regions of the CCLME (e.g., onset and duration of upwelling season on the shelf, eddy variability).

DNA metabarcoding is an important tool for molecular ecology. However, its effectiveness hinges on the quality of reference sequence databases and classification parameters employed. Here we evaluate the performance of MiFish 12S taxonomic assignments using a case study of California Current Large Marine Ecosystem fishes to determine best practices for metabarcoding. Specifically, we use a taxonomy cross-validation by identity framework to compare classification performance between a global database comprised of all available sequences and a curated database that only includes sequences of fishes from the California Current Large Marine Ecosystem. We demonstrate that the regional database provides higher assignment accuracy than the comprehensive global database. We also document a tradeoff between accuracy and misclassification across a range of taxonomic cutoff scores, highlighting the importance of parameter selection for taxonomic classification. Furthermore, we compared assignment accuracy with and without the inclusion of additionally generated reference sequences. To this end, we sequenced tissue from 597species using the MiFish 12S primers, adding 252species to GenBank's existing 550 California Current Large Marine Ecosystem fish sequences. We then compared species and reads identified from seawater environmental DNA samples using global databases with and without our generated references, and the regional database. The addition of new references allowed for the identification of 16 additional native taxa representing 17.0% of total reads from eDNA samples, including species with vast ecological and economic value. Together these results demonstrate the importance of comprehensive and curated reference databases for effective metabarcoding and the need for locus-specific validation efforts.

Ocean acidification is intensifying and hypoxia is projected to expand in the California Current large marine ecosystem as a result of processes associated with the global emission of CO2. Observed changes in the California Current outpace those in many other areas of the ocean, underscoring the pressing need to adopt management approaches that can accommodate uncertainty and the complicated dynamics forced by accelerating change. We argue that changes occurring in the California Current large marine ecosystem provide opportunities and incentives to adopt an integrated, systems-level approach to resource management to preserve existing ecosystem services and forestall abrupt change. Practical options already exist to maximize the benefits of management actions and ameliorate impending change in the California Current, for instance, adding ocean acidification and hypoxia to design criteria for marine protected areas, including consideration of ocean acidification and hypoxia in fisheries management decisions, and fully enforcing existing laws and regulations that govern water quality and land use and development.

Forecasts of sea surface temperature anomalies (SSTAs) provide essential information to stakeholders of marine resources in coastal ecosystems, such as the California Current Large Marine Ecosystem (CCLME), at management-relevant monthly-to-annual time scales. Diagnosing dynamical sources of predictability and the mechanisms differentiating skill among forecasts is required for verification and improvement in operational forecasting systems. Using retrospective forecasts (1982–2020) from a four-member subset of the North American Multi-Model Ensemble (NMME), we evaluate the conditional skill of SSTA forecasts in the CCLME at monthly resolution for lead times up to 10.5 months. Forecasts from ensemble members with relatively small SSTA errors at shorter lead times retain higher skill at longer lead times, with the most substantial and long-lasting increases for forecasts initialized in the fall and early spring. The “best” low-error SSTA forecasts are characterized by increased skill in the prediction of North Pacific atmospheric circulation [sea level pressure (SLP) and 200-hPa geopotential height] the month prior to the evaluation of SSTA errors in the CCLME and exhibit more realistic progressions of anomalous SLP. The Pacific meridional mode (PMM) emerges as a diagnostic of skillful North Pacific atmosphere–ocean coupling, as forecasts that correctly simulate the PMM and its associated SLP variability increase the SSTA prediction skill in the CCLME in the fall through spring. Predictable coupled ocean–atmosphere modes provide a target for enhancing predictability with early detection of the onset of a deterministic progression emerging from stochastic atmospheric variability. Significance Statement Global forecast systems provide near-term climate predictions that inform the management of marine resources, such as those of the California Current Large Marine Ecosystem. In this study, we probe the processes which lead forecasts to succeed or fail at predicting sea surface temperatures in the California Current at seasonal time scales among retrospective forecasts from the North American Multimodel Ensemble. We demonstrate that forecasts which best simulate sea surface temperatures at the earliest lead times sustain advantages in forecast skill and find that correctly simulating extratropical atmospheric circulation increases the predictive skill of sea surface temperatures in the northeast Pacific in the following lead times. Our results offer North Pacific atmospheric circulation as a target for forecast model improvement that would additionally enhance ocean forecasts.

We analyze and compare changes in ocean acidification metrics caused by anthropogenic carbon (Canth) accumulation in the North Pacific Ocean and California Current Large Marine Ecosystem (CCLME). The greatest declines in pH and carbonate mineral saturation state occur near the surface, coincident with the highest Canth concentrations. However, maximal increases in the partial pressure of carbon dioxide (pCO2) and hydrogen ion concentration occur subsurface where Canth values are lower. We attribute dissimilar sensitivities of these metrics to background ocean chemistry, which has naturally high pCO2 and low buffering capacity in subsurface waters due to accumulated byproducts of organic matter respiration, which interacts with Canth. In the CCLME, rising subsurface pCO2 has increased the frequency, duration, and intensity of hypercapnia (pCO2 ≥ 1,000 μatm) on the continental shelf. Our findings suggest that hypercapnia induced by Canth accumulation can co‐occur with hypoxia in the CCLME and is an additional modern stressor for marine organisms.

Juvenile North Pacific Albacore tuna ( Thunnus alalunga ) support commercial and recreational fisheries in the California Current Large Marine Ecosystem (CCLME), where they forage during summer and fall. The distributions of the commercial and recreational fisheries and estimates of forage availability have varied substantially over the past century. Time‐series quantifying Albacore diet can help link forage composition to variability in Albacore abundance and distribution and, consequently, their availability to fishers. Previous diet studies in the CCLME are of relatively short duration, and long‐term variability in Albacore diet remains poorly understood. We describe the diets of juvenile Albacore from three regions in the CCLME from 2007 to 2019 and use classification and regression tree analysis to explore environmental drivers of variability. Important prey include Northern Anchovy ( Engraulis mordax ), rockfishes ( Sebastes spp.), Boreal Clubhook Squid ( Onychoteuthis borealijaponica ), euphausiids (Order: Euphausiidae), and amphipods (Order: Amphipoda), each contributing >5% mean proportional abundance. Most prey items were short lived species or young‐of‐the‐year smaller than 10 cm. Diet variability was related to environmental conditions over the first 6 months of the year (PDO, sea surface temperature, and NPGO) and conditions concurrent with Albacore capture (region and surface nitrate flux). We describe foraging flexibility over regional and annual scales associated with these environmental influences. Continuous, long‐term studies offer the opportunity to identify flexibility in Albacore foraging behavior and begin to make a predictive link between environmental conditions early in the year and Albacore foraging during summer and fall.

Abstract. The California Current Large Marine Ecosystem (CCLME), a temperate marine region dominated by episodic upwelling, is predicted to experience rapid environmental change in the future due to ocean acidification. The aragonite saturation state within the California Current System is predicted to decrease in the future with near-permanent undersaturation conditions expected by the year 2050. Thus, the CCLME is a critical region to study due to the rapid rate of environmental change that resident organisms will experience and because of the economic and societal value of this coastal region. Recent efforts by a research consortium – the Ocean Margin Ecosystems Group for Acidification Studies (OMEGAS) – has begun to characterize a portion of the CCLME; both describing the spatial mosaic of pH in coastal waters and examining the responses of key calcification-dependent benthic marine organisms to natural variation in pH and to changes in carbonate chemistry that are expected in the coming decades. In this review, we present the OMEGAS strategy of co-locating sensors and oceanographic observations with biological studies on benthic marine invertebrates, specifically measurements of functional traits such as calcification-related processes and genetic variation in populations that are locally adapted to conditions in a particular region of the coast. Highlighted in this contribution are (1) the OMEGAS sensor network that spans the west coast of the US from central Oregon to southern California, (2) initial findings of the carbonate chemistry amongst the OMEGAS study sites, and (3) an overview of the biological data that describes the acclimatization and the adaptation capacity of key benthic marine invertebrates within the CCLME.

Indicators of pelagic forage community shifts in the California Current Large Marine Ecosystem, 1998–2016

Pelagic marine predators face unprecedented challenges and uncertain futures. Overexploitation and climate variability impact the abundance and distribution of top predators in ocean ecosystems. Improved understanding of ecological patterns, evolutionary constraints and ecosystem function is critical for preventing extinctions, loss of biodiversity and disruption of ecosystem services. Recent advances in electronic tagging techniques have provided the capacity to observe the movements and long-distance migrations of animals in relation to ocean processes across a range of ecological scales. Tagging of Pacific Predators, a field programme of the Census of Marine Life, deployed 4,306 tags on 23 species in the North Pacific Ocean, resulting in a tracking data set of unprecedented scale and species diversity that covers 265,386 tracking days from 2000 to 2009. Here we report migration pathways, link ocean features to multispecies hotspots and illustrate niche partitioning within and among congener guilds. Our results indicate that the California Current large marine ecosystem and the North Pacific transition zone attract and retain a diverse assemblage of marine vertebrates. Within the California Current large marine ecosystem, several predator guilds seasonally undertake north-south migrations that may be driven by oceanic processes, species-specific thermal tolerances and shifts in prey distributions. We identify critical habitats across multinational boundaries and show that top predators exploit their environment in predictable ways, providing the foundation for spatial management of large marine ecosystems.