No abstract

Considerable uncertainty remains over how increasing atmospheric CO2 and anthropogenic climate changes are affecting open‐ocean marine ecosystems from phytoplankton to top predators. Biological time series data are thus urgently needed for the world's oceans. Here, we use the carbon stable isotope composition of tuna to provide a first insight into the existence of global trends in complex ecosystem dynamics and changes in the oceanic carbon cycle. From 2000 to 2015, considerable declines in δ13C values of 0.8‰–2.5‰ were observed across three tuna species sampled globally, with more substantial changes in the Pacific Ocean compared to the Atlantic and Indian Oceans. Tuna recorded not only the Suess effect, that is, fossil fuel‐derived and isotopically light carbon being incorporated into marine ecosystems, but also recorded profound changes at the base of marine food webs. We suggest a global shift in phytoplankton community structure, for example, a reduction in 13C‐rich phytoplankton such as diatoms, and/or a change in phytoplankton physiology during this period, although this does not rule out other concomitant changes at higher levels in the food webs. Our study establishes tuna δ13C values as a candidate essential ocean variable to assess complex ecosystem responses to climate change at regional to global scales and over decadal timescales. Finally, this time series will be invaluable in calibrating and validating global earth system models to project changes in marine biota.

Global anthropogenic mercury (Hg) emissions to the atmosphere since industrialization are widely considered to be responsible for a significant increase in surface ocean Hg concentrations. Still unclear is how those inputs are converted into toxic methylmercury (MeHg) then transferred and biomagnified in oceanic food webs. We used a unique long-term and continuous dataset to explore the temporal Hg trend and variability of three tropical tuna species (yellowfin, bigeye, and skipjack) from the southwestern Pacific Ocean between 2001 and 2018 (n = 590). Temporal trends of muscle nitrogen (δ15N) and carbon (δ13C) stable isotope ratios, amino acid (AA) δ15N values and oceanographic variables were also investigated to examine the potential influence of trophic, biogeochemical and physical processes on the temporal variability of tuna Hg concentrations. For the three species, we detected significant inter-annual variability but no significant long-term trend for Hg concentrations. Inter-annual variability was related to the variability in tuna sampled lengths among years and to tuna muscle δ15N and δ13C values. Complementary AA-and model-estimated phytoplankton δ15N values suggested the influence of baseline processes with enhanced tuna Hg concentrations observed when dinitrogen fixers prevail, possibly fuelling baseline Hg methylation and/or MeHg bioavailability at the base of the food web. Our results show that MeHg trends in top predators do not necessary capture the increasing Hg concentrations in surface waters suspected at the global oceanic scale due to the complex and variable processes governing Hg deposition, methylation, bioavailability and biomagnification. This illustrates the need for long-term standardized monitoring programs of marine biota worldwide.Please note that this is an author-produced PDF of an article accepted for publication following peer review. The definitive publisher-authenticated version is available on the publisher Web site.

No abstract

Acoustic data are invaluable information sources for characterizing the distribution and abundance of midtrophic-level organisms (micronekton). These organisms play a pivotal role in the ecosystem as prey of top predators and as predators of low-trophic-level organisms. Although shipboard-ADCP (acoustic Doppler current profiler) acoustic backscatter signal intensity cannot provide an absolute biomass estimate, it may be a useful proxy to investigate variability in the distribution and relative density of micronekton. This study used acoustic recordings data spread across 19 years (1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017) from 54 ADCP cruises in New Caledonia's subtropical EEZ (exclusive economic zone) to assess seasonal and interannual variabilities and spatial distribution of micronekton. The dataset was composed of two different ADCPs: 150 kHz for the first period, followed by 75 kHz for more recent years. We examined the 20-120 m averaged scattering layer. Using the few cruises with concurrent EK60 measurements, we proposed that the backscatter from the ADCPs and 70 kHz EK60 were sufficiently closely linked to allow the use of the backscatter signal from the ADCPs in a combined dataset over the full time series. We then designed a GAMM (generalized additive mixed model) model that takes into account the two ADCP devices as well as temporal variability. After accounting for the effect of the devices, we showed that the acoustic signal was mainly driven by diel vertical migration, season, year, and ENSO (El Niño-Southern Oscillation). In a second step, a consensus model between two statistical approaches (GAMM and SVM) (support vector machine) was constructed, linking the nighttime 20-120 m backscatter to the oceanographic and geographic environment. This model showed that sea surface temperature was the main factor driving backscatter variability in the EEZ, with intensified backscatter during the austral summer (December to May) in the northern part of the EEZ. We showed that acoustic density differed significantly, spatially and temporally from micronekton biomass predicted for the same period by the SEAPODYM-MTL (mid-trophic level) ecosystem model. The seasonal cycle given by ADCP data lagged behind the SEAPODYM-MTL Please note that this is an author-produced PDF of an article accepted for publication following peer review. The definitive publisher-authenticated version is available on the publisher Web site.seasonal cycle by around three months. Reasons to explain these differences and further needs in observation and modeling were explored in the discussion. In addition to providing new insights for micronekton dynamics in this EEZ (i.e., the science needed for ecosystem-based fisheries management), the data should help improve our ability to model this key trophic component. Highlights► 20 years of ADCP acoustic data give insights on spatiotemporal micronekton dynamics. ► Sea-surface temperature is the main factor driving the backsca...

The western Baltic Sea cod (WBC) stock is at historically low levels, mainly attributed to high fishing pressure and low recruitment. Stable stock assessment metrics suggested recovery potential, given appropriate fisheries management measures. However, changing environmental conditions violate stability assumptions, may negatively affect WBC, and challenge the resource management. The present study explored 42 years of changes in WBC biological parameters. WBC body condition gradually decreased over the last decades for juveniles and adults, with a rapid decrease in recent years when a single cohort dominated the overfished stock. The hepato-somatic index and the muscle weight decreased by 50% and 10% in the last 10 years, respectively, suggesting severely decreasing energy reserves and productivity. The changes in energy reserves were associated with changes in environmental conditions (increase in bottom water temperature, expansion of hypoxic areas during late summer/autumn), and changes in diet composition (less herring). A key bottleneck is the warming and longer-lasting summer period when WBC, trapped between warmed shallow waters and hypoxic deeper waters, have to mobilize energy reserves to account for reduced feeding opportunities and thermal stress. Our results suggest that stock recovery is unlikely to happen by fisheries management alone if environmental trajectories remain unchanged.

No abstract

Understanding the ecological mechanisms underpinning distribution patterns is vital in managing populations of mobile marine species. This study is a first step towards an integrated description of the habitats and spatial distributions of marine predators in the Natural Park of the Coral Sea, one of the world’s largest marine-protected areas at about 1.3 million km2, covering the entirety of New Caledonia’s pelagic waters. The study aims to quantify the benefit of including a proxy for prey abundance in predator niche modelling, relative to other marine physical variables. Spatial distributions and relationships with environmental data were analysed using catch per unit of effort data for three fish species (albacore tuna, yellowfin tuna and dolphinfish), sightings collected from aerial surveys for three cetacean guilds (Delphininae, Globicephalinae and Ziphiidae) and foraging locations identified from bio-tracking for three seabird species (wedge-tailed shearwater, Tahiti petrel and red-footed booby). Predator distributions were modelled as a function of a static covariate (bathymetry), oceanographic covariates (sea surface temperature, chlorophyll-a concentration and 20 °C-isotherm depth) and an acoustically derived micronekton preyscape covariate. While distributions were mostly linked to bathymetry for seabirds, and chlorophyll and temperature for fish and cetaceans, acoustically derived prey abundance proxies slightly improved distribution models for all fishes and seabirds except the Tahiti petrel, but not for the cetaceans. Predicted spatial distributions showed that pelagic habitats occupied by predator fishes did not spatially overlap. Finally, predicted habitats and the use of the preyscapes in predator habitat modelling were discussed.