Copy this text
Global patterns and inferences of tuna movements and trophodynamics
A global dataset of carbon stable isotope (δ13C) values from yellowfin, bigeye, and albacore tuna muscle tissue (n = 4275) was used to develop a novel tool to infer broad-scale movement and residency patterns of these highly mobile marine predators. This tool was coupled with environmental models and lipid content (C:N ratio) of tuna muscle tissues to examine ocean warming impacts on tuna ecology and bioenergetic condition across Longhurst provinces. Over a 16-year study period (2000–2015), latitudinal gradients in tuna δ13C values were consistent, with values decreasing with increasing latitude. Tuna δ13C values, reflecting modelled global phytoplankton δ13C landscapes (“isoscapes”), were largely related to spatial changes in oxygen concentrations at depth and temporal changes in sea surface temperature. Observed tuna isoscapes (δ13CLScorr), corrected for lipid content and the Suess effect (oceanic changes in CO2 over time), were subtracted from model-predicted baseline isoscapes (Δ13Ctuna-phyto) to infer spatial movement and residency patterns of the different tuna species. Stable isotope niche width was calculated for each Longhurst province using Δ13Ctuna-phyto and baseline-corrected nitrogen isotope (δ15Ntuna-phyto) values to further quantify isotopic variability as evidence of movements across isoscapes. A high degree of movement—defined as the deviation from expected range of Δ13Ctuna-phyto values— was evident in four Longhurst provinces: Guinea current coast, North Atlantic drift, Pacific equatorial divergence, and the North Pacific equatorial counter current. The highest level of population dispersal (variability in Δ13Ctuna-phyto values) was observed in Longhurst provinces within the western and central Pacific Oceans and in the Guinea current coast. While lipid content was low in yellowfin and bigeye, high and variable lipid stores in albacore muscle were consistent with seasonal movements between productive foraging and oligotrophic spawning habitats. Our ability to characterize tuna movement patterns without ambiguity remains challenged by uncertainty in trophic discrimination factors and ecological (e.g. diet variability) processes. However, this study illustrates that model-corrected δ13C values are a valuable, relatively cost-effective tool for identifying potential areas of mixing across management zones, particularly when electronic tagging studies are limited or absent. Stable isotope analyses of tuna tissues can therefore be an additional tool for guiding spatial stock assessments on top predator movement, dispersal patterns, and how they may be altered under a changing climate.
Keyword(s)
Albacore tuna, Bigeye tuna, Yellowfin tuna, Carbon isotope analysis, Movement indicators, Global ocean, Fisheries management
Full Text
File | Pages | Size | Access | |
---|---|---|---|---|
Author's final draft | 48 | 3 Mo | ||
Multimedia component 1. | 1 | 297 Ko | ||
Multimedia component 2 | 2 | 316 Ko | ||
Multimedia component 3 | 2 | 319 Ko | ||
Multimedia component 4 | 2 | 313 Ko | ||
Multimedia component 5 | 2 | 313 Ko | ||
Multimedia component 6 | 1 | 339 Ko | ||
Fig. S1. | - | 332 Ko | ||
Fig. S3. GAMM produced isoscapes of lipid and Suess corrected δ13C values (δ13CLScorr) in (A) yellowfin, (B) bigeye, and (C) albacore tunas sampled between 2000–2007 and 2008–2015. | - | 501 Ko | ||
Fig. S4. Tuna muscle (lipid and Suess-corrected) and modelled phytoplankton (Suess corrected) (A) δ13C values and (B) δ13C separation (Δ13Ctuna-phyto) in relation to distance from the equator. | - | 570 Ko | ||
Publisher's official version | 19 | 3 Mo |