Using Predicted Patterns of 3D Prey Distribution to Map King Penguin Foraging Habitat

King penguins (Aptenodytes patagonicus) are an iconic Southern Ocean species, but the prey distributions that underpin their at-sea foraging tracks and diving behaviour remain unclear. We conducted simultaneous acoustic surveys off South Georgia and tracking of king penguins breeding ashore there in Austral summer 2017 to gain insight into habitat use and foraging behaviour. Acoustic surveys revealed ubiquitous deep scattering layers (DSLs; acoustically detected layers of fish and other micronekton that inhabit the mesopelagic zone) at c. 500 m and shallower ephemeral fish schools. Based on DNA extracted from penguin faecal samples, these schools were likely comprised of lanternfish (an important component of king penguin diets), icefish (Channichthyidae spp.) and painted noties (Lepidonotothen larseni). Penguins did not dive as deep as DSLs, but their prey-encounter depth-distributions, as revealed by biologging, overlapped at fine scale (10s of m) with depths of acoustically detected fish schools. We used neural networks to predict local scale (10 km) fish echo intensity and depth distribution at penguin dive locations based on environmental correlates, and developed models of habitat use. Habitat modelling revealed that king penguins preferentially foraged at locations predicted to have shallow and dense (high acoustic energy) fish schools associated with shallow and dense DSLs. These associations could be used to predict the distribution of king penguins from other colonies at South Georgia for which no tracking data are available, and to identify areas of potential ecological significance within the South Georgia and the South Sandwich Islands marine protected area.

First tracking of the oceanic spawning migrations of Australasian short-finned eels (Anguilla australis)

Anguillid eel populations have declined dramatically over the last 50 years in many regions of the world, and numerous species are now under threat. A critical life-history phase is migration from freshwater to distant oceans, culminating in a single life-time spawning event. For many anguillids, especially those in the southern hemisphere, mystery still shrouds their oceanic spawning migrations. We investigated the oceanic spawning migrations of the Australasian short-finned eel (Anguilla australis) using pop-up satellite archival tags. Eels were collected from river estuaries (38° S, 142° E) in south-eastern temperate Australia. In 2019, 16 eels were tracked for up to about 5 months, ~ 2620 km from release, and as far north as the tropical Coral Sea (22° S, 155° E) off the north-east coast of Australia. Eels from southern Australia appeared to access deep water off the Australian coast via two main routes: (i) directly east via Bass Strait, or (ii) south-east around Tasmania, which is the shortest route to deep water. Tagged eels exhibited strong diel vertical migrations, alternating between the warm euphotic zone (~ 100–300 m, 15–20 °C) at night and the mesopelagic zone (~ 700–900 m, 6–8 °C) during the day. Marine predators, probably lamnid sharks, tuna, or marine mammals, ended many eel migrations (at least ~ 30%), largely before the eels had left the Australian continental shelf. The long and risky marine migrations of Australasian eels highlight the need for better information on the processes contributing to eel mortality throughout the life cycle, including the impacts of future changes to oceanic currents, predator abundance and direct anthropogenic disturbances.

Episodic subduction patches in the western North Pacific identified from BGC-Argo float data

Subduction associated with mesoscale eddies is an important but difficult-to-observe process that can efficiently export carbon and oxygen to the mesopelagic zone (100–1000 dbar). Using a novel BGC-Argo dataset covering the western North Pacific (20–50∘ N, 120–180∘ E), we identified imprints of episodic subduction using anomalies in dissolved oxygen and spicity, a water mass marker. These subduction patches were present in 4.0 % (288) of the total profiles (7120) between 2008 and 2019, situated mainly in the Kuroshio Extension region between March and August (70.6 %). Roughly 31 % and 42 % of the subduction patches were identified below the annual permanent pycnocline depth (300 m vs. 450 m) in the subpolar and subtropical regions, respectively. Around half (52 %) of these episodic events injected oxygen-enriched waters below the maximum annual permanent thermocline depth (450 dbar), with >20 % occurring deeper than 600 dbar. Subduction patches were detected during winter and spring when mixed layers are deep. The oxygen inventory within these subductions is estimated to be on the order of 64 to 152 g O2/m2. These mesoscale events would markedly increase oxygen ventilation as well as carbon removal in the region, both processes helping to support the nutritional and metabolic demands of mesopelagic organisms. Climate-driven patterns of increasing eddy kinetic energies in this region imply that the magnitude of these processes will grow in the future, meaning that these unexpectedly effective small-scale subduction processes need to be better constrained in global climate and biogeochemical models.

Increasing Hydrostatic Pressure Impacts the Prokaryotic Diversity during Emiliania huxleyi Aggregates Degradation

In the dark ocean, the balance between the heterotrophic carbon demand and the supply of sinking carbon through the biological carbon pump remains poorly constrained. In situ tracking of the dynamics of microbial degradation processes occurring on the gravitational sinking particles is still challenging. Our particle sinking simulator system (PASS) intends to mimic as closely as possible the in situ variations in pressure and temperature experienced by gravitational sinking particles. Here, we used the PASS to simultaneously track geochemical and microbial changes that occurred during the sinking through the mesopelagic zone of laboratory-grown Emiliania huxleyi aggregates amended by a natural microbial community sampled at 105 m depth in the North Atlantic Ocean. The impact of pressure on the prokaryotic degradation of POC and dissolution of E. huxleyi-derived calcite was not marked compared to atmospheric pressure. In contrast, using global O2 consumption monitored in real-time inside the high-pressure bottles using planar optodes via a sapphire window, a reduction of respiration rate was recorded in surface-originated community assemblages under increasing pressure conditions. Moreover, using a 16S rRNA metabarcoding survey, we demonstrated a drastic difference in transcriptionally active prokaryotes associated with particles, incubated either at atmospheric pressure or under linearly increasing hydrostatic pressure conditions. The increase in hydrostatic pressure reduced both the phylogenetic diversity and the species richness. The incubation at atmospheric pressure, however, promoted an opportunistic community of “fast” degraders from the surface (Saccharospirillaceae, Hyphomonadaceae, and Pseudoalteromonadaceae), known to be associated with surface phytoplankton blooms. In contrast, the incubation under increasing pressure condition incubations revealed an increase in the particle colonizer families Flavobacteriaceae and Rhodobacteraceae, and also Colwelliaceae, which are known to be adapted to high hydrostatic pressure. Altogether, our results underline the need to perform biodegradation experiments of particles in conditions that mimic pressure and temperature encountered during their sinking along the water column to be ecologically relevant.

Uncertainties of particulate organic carbon concentrations in the mesopelagic zone of the Atlantic ocean

Measurements of particulate organic carbon (POC) in the open ocean provide grounds for estimating oceanic carbon budgets and for modelling carbon cycling. The majority of the published POC measurements have been collected at the sea surface. Thus, POC stocks in the upper layer of the water column are relatively well constrained. However, our understanding of the POC distribution and its dynamics in deeper areas is modest due to insufficient in POC measurements. Moreover, the accuracy of published POC estimates is not always quantified, and neither is it fully understood. In this study, we determined the POC concentrations of samples collected in the upper 500 m during an Atlantic Meridional Transect and described a method for quantifying its experimental uncertainties using duplicate measurements. The analysis revealed that the medians of the total experimental uncertainties associated with our POC concentrations in the productive and mesopelagic zones were 2.5(±1.2) mg/m3 and 2.6(±0.6) mg/m3, respectively. In relative terms, these uncertainties corresponded to ~14% and ~ 35% of POC concentrations, respectively. However, despite our best efforts, we could explain only ~ 21% of the total experimental POC uncertainty. The potential sources of this unexplained portion of uncertainty are discussed.

First Tracking of the Oceanic Spawning Migrations of Australasian Short- Finned Eels (Anguilla Australis)

Across their near-global range, anguillid eel populations have declined dramatically over the last fifty years and many species are now threatened. A critical life history phase is migration from freshwater to distant oceans, culminating in their single life-time spawning event. For many anguillids, especially those in the southern hemisphere, mystery still shrouds their oceanic spawning migrations. We investigated the oceanic spawning migrations of the Australasian short-finned eel (Anguilla australis) using pop-up satellite archival tags. Eels were collected from river estuaries (38°S, 142°E) in south-eastern temperate Australia. In 2019, sixteen eels were tracked for up to about 5 months, ~2620 km from release, and as far north as the tropical Coral Sea (22°S, 155°E) off the north-east coast of Australia. Eels from southern Australia appeared to access deep water off the Australian coast via two main routes: (i) directly east via Bass Strait, or (ii) south-east around Tasmania, which is the shortest route to deep water. Tagged eels exhibited strong diel vertical migrations, alternating between the warm euphotic zone (~100 to 300 m, 15 to 20°C) at night, and the mesopelagic zone (~700 m to 900 m, 6 to 8°C) during the day. Marine predators, probably lamnid sharks, tuna, or marine mammals, ended many eel migrations (at least ~30%), with many eels lost before leaving the Australian continental shelf. The long and risky marine migrations of Australasian eels highlight the need for better information on the processes contributing to eel mortality across the life cycle, including the impacts of future changes to oceanic currents, predator abundance and direct anthropogenic disturbances.

Episodic subduction patches in the western North Pacific identified from BGC-Argo float Data

Subduction associated with mesoscale eddies is an important but difficult to observe process that can efficiently export carbon and oxygen to the mesopelagic zone (100–1000 db). Using a novel BGC-Argo dataset covering the western North Pacific (20–50° N, 120–180° E), we identified imprints of episodic subduction using anomalies in dissolved oxygen and spicity, a water mass marker. These subduction patches were present in 4.0 % (288) of the total profiles (7,120) between 2008 and 2019, situated mainly in the Kuroshio Extension region between March and August (70.6 %). Unlike eddy subduction processes observed at higher latitudes, roughly half (52 %) of these episodic events injected carbon- and oxygen-enriched waters below the annual permanent thermocline depth (450 db), with > 20 % occurring deeper than 600 db. Export rates within these subductions are estimated to be on the order of 85–159 mg C m−2 day−1 and 175 to 417 mg O2 m−2 day−1. These mesoscale events would markedly increase carbon removal above that due to biological gravitational settling as well as oxygen ventilation in the region, both helping to support the nutritional and metabolic demands of mesopelagic organisms. Climate-driven patterns of increasing eddy kinetic energies in this region imply that the magnitude of these processes will grow in the future, meaning that these unexpectedly effective small-scale subduction processes need to be better constrained in global climate and biogeochemical models.

Drivers of Spatial Distributions of Basking Shark (Cetorhinus maximus) in the Southwest Pacific

Basking sharks (Cetorhinus maximus) were widely reported throughout New Zealand waters. Once commonly observed, and sometimes in large numbers, basking sharks are now infrequently reported. Basking shark observations are known to be highly variable across years, and their distribution and occurrence have been shown to be influenced by environmental predictors such as thermal fronts, chl-a concentration, and the abundance of prey (zooplankton). Little is known of basking sharks in the South Pacific and more information on distribution, habitat use, and migratory patterns is required to better understand the species’ regional ecology. Here, we used bootstrapped Habitat Suitability Models [HSM, ensembled from Boosted Regression Tree (BRT) and Random Forest (RF) models] to determine the drivers of basking shark distribution, predict habitat suitability and estimated uncertainty in the South Pacific for the first time. High−resolution environmental (1 km2 grid resolution) and biotic data, including inferred prey species, and all available basking shark records across New Zealand’s Exclusive Economic Zone (EEZ) were included in the ensemble HSMs. The most influential driver of modeled basking shark distribution was vertical flux of particulate organic matter at the seabed, which may indicate higher levels of primary production in the surface ocean and higher prey density in the mesopelagic zone and at the seafloor. The BRT and RF models had good predictive power (AUC and TSS > 0.7) and both models performed similarly with low variability in the model fit metrics. Areas of high basking shark habitat suitability included the east and west coasts of the South Island, Puysegur Ridge, and Auckland Island slope. The outputs produced here could be incorporated into future management framework for assessing threat and conservation needs (e.g., spatially explicit risk assessment) for this regionally protected species, as well as providing guidance for future research efforts (e.g., areas of interest for sampling).

Uncertainties of particulate organic carbon concentrations in the mesopelagic zone of the Atlantic ocean

Measurements of particulate organic carbon (POC) in the open ocean provide grounds for estimating oceanic carbon budgets and for modelling carbon cycling. The majority of the published POC measurements have been collected at the sea surface. Thus, POC stocks in the upper layer of the water column are relatively well constrained. However, our understanding of the POC distribution and its dynamics in deeper areas is modest due to insufficient in POC measurements. Moreover, the accuracy of published POC estimates is not always quantified, and neither is it fully understood. In this study, we determined the POC concentrations of samples collected in the upper 500 m during an Atlantic Meridional Transect and described a method for quantifying its experimental uncertainties using duplicate measurements. The analysis revealed that the medians of the total experimental uncertainties associated with our POC concentrations in the productive and mesopelagic zones were 2.5(±1.2) mg/m3 and 2.6(±0.6) mg/m3, respectively. In relative terms, these uncertainties corresponded to ~14% and ~ 35% of POC concentrations, respectively. However, despite our best efforts, we could explain only ~ 21% of the total experimental POC uncertainty. The potential sources of this unexplained portion of uncertainty are discussed.