The vulnerability of sharks, skates, and rays to ocean deoxygenation: Physiological mechanisms, behavioral responses, and ecological impacts.

Levels of dissolved oxygen in open ocean and coastal waters are decreasing (ocean deoxygenation), with poorly understood effects on marine megafauna. All of the more than 1000 species of elasmobranchs (sharks, skates, and rays) are obligate water breathers, with a variety of life-history strategies and oxygen requirements. This review demonstrates that although many elasmobranchs typically avoid hypoxic water, they also appear capable of withstanding mild to moderate hypoxia with changes in activity, ventilatory responses, alterations to circulatory and hematological parameters, and morphological alterations to gill structures. However, such strategies may be insufficient to withstand severe, progressive, or prolonged hypoxia or anoxia where anaerobic metabolic pathways may be used for limited periods. As water temperatures increase with climate warming, ectothermic elasmobranchs will exhibit elevated metabolic rates and are likely to be less able to tolerate the effects of even mild hypoxia associated with deoxygenation. As a result, sustained hypoxic conditions in warmer coastal or surface-pelagic waters are likely to lead to shifts in elasmobranch distributions. Mass mortalities of elasmobranchs linked directly to deoxygenation have only rarely been observed but are likely underreported. One key concern is how reductions in habitat volume as a result of expanding hypoxia resulting from deoxygenation will influence interactions between elasmobranchs and industrial fisheries. Catch per unit of effort of threatened pelagic sharks by longline fisheries, for instance, has been shown to be higher above oxygen minimum zones compared to adjacent, normoxic regions, and attributed to vertical habitat compression of sharks overlapping with increased fishing effort. How a compound stressor such as marine heatwaves alters vulnerability to deoxygenation remains an open question. With over a third of elasmobranch species listed as endangered, a priority for conservation and management now lies in understanding and mitigating ocean deoxygenation effects in addition to population declines already occurring from overfishing.

Shark conservation requires mortality-limiting regulations amid global change.

Despite anti-finning laws aimed at conserving sharks, Worm et al. have revealed that global shark mortality rates have surprisingly risen over the past decade, driven in large part by increased demand for meat. Here, we discuss the importance of this study, underscoring the need for broader regulations addressing overall shark mortality amid threats from global change.

Highly active fish in low oxygen environments: vertical movements and behavioural responses of bigeye and yellowfin tunas to oxygen minimum zones in the eastern Pacific Ocean.

Oxygen minimum zones in the open ocean are predicted to significantly increase in volume over the coming decades as a result of anthropogenic climatic warming. The resulting reduction in dissolved oxygen (DO) in the pelagic realm is likely to have detrimental impacts on water-breathing organisms, particularly those with higher metabolic rates, such as billfish, tunas, and sharks. However, little is known about how free-living fish respond to low DO environments, and therefore, the effect increasing OMZs will have cannot be predicted reliably. Here, we compare the responses of two active predators (bigeye tuna Thunnus obesus and yellowfin tuna Thunnus albacares) to DO at depth throughout the eastern Pacific Ocean. Using time-series data from 267 tagged tunas (59,910 days) and 3D maps of modelled DO, we find that yellowfin tuna respond to low DO at depth by spending more time in shallower, more oxygenated waters. By contrast, bigeye tuna, which forage at deeper depths well below the thermocline, show fewer changes in their use of the water column. However, we find that bigeye tuna increased the frequency of brief upward vertical excursions they performed by four times when DO at depth was lower, but with no concomitant significant difference in temperature, suggesting that this behaviour is driven in part by the need to re-oxygenate following time spent in hypoxic waters. These findings suggest that increasing OMZs will impact the behaviour of these commercially important species, and it is therefore likely that other water-breathing predators with higher metabolic rates will face similar pressures. A more comprehensive understanding of the effect of shoaling OMZs on pelagic fish vertical habitat use, which may increase their vulnerability to surface fisheries, will be important to obtain if these effects are to be mitigated by future management actions. Supplementary Information: The online version contains supplementary material available at 10.1007/s00227-023-04366-2.

Linking vertical movements of large pelagic predators with distribution patterns of biomass in the open ocean.

Many predator species make regular excursions from near-surface waters to the twilight (200 to 1,000 m) and midnight (1,000 to 3,000 m) zones of the deep pelagic ocean. While the occurrence of significant vertical movements into the deep ocean has evolved independently across taxonomic groups, the functional role(s) and ecological significance of these movements remain poorly understood. Here, we integrate results from satellite tagging efforts with model predictions of deep prey layers in the North Atlantic Ocean to determine whether prey distributions are correlated with vertical habitat use across 12 species of predators. Using 3D movement data for 344 individuals who traversed nearly 1.5 million km of pelagic ocean in [Formula: see text]42,000 d, we found that nearly every tagged predator frequented the twilight zone and many made regular trips to the midnight zone. Using a predictive model, we found clear alignment of predator depth use with the expected location of deep pelagic prey for at least half of the predator species. We compared high-resolution predator data with shipboard acoustics and selected representative matches that highlight the opportunities and challenges in the analysis and synthesis of these data. While not all observed behavior was consistent with estimated prey availability at depth, our results suggest that deep pelagic biomass likely has high ecological value for a suite of commercially important predators in the open ocean. Careful consideration of the disruption to ecosystem services provided by pelagic food webs is needed before the potential costs and benefits of proceeding with extractive activities in the deep ocean can be evaluated.

A vision for incorporating human mobility in the study of human-wildlife interactions.

As human activities increasingly shape land- and seascapes, understanding human-wildlife interactions is imperative for preserving biodiversity. Habitats are impacted not only by static modifications, such as roads, buildings and other infrastructure, but also by the dynamic movement of people and their vehicles occurring over shorter time scales. Although there is increasing realization that both components of human activity substantially affect wildlife, capturing more dynamic processes in ecological studies has proved challenging. Here we propose a conceptual framework for developing a 'dynamic human footprint' that explicitly incorporates human mobility, providing a key link between anthropogenic stressors and ecological impacts across spatiotemporal scales. Specifically, the dynamic human footprint integrates a range of metrics to fully acknowledge the time-varying nature of human activities and to enable scale-appropriate assessments of their impacts on wildlife behaviour, demography and distributions. We review existing terrestrial and marine human-mobility data products and provide a roadmap for how these could be integrated and extended to enable more comprehensive analyses of human impacts on biodiversity in the Anthropocene.

Direct measurement of cruising and burst swimming speeds of the shortfin mako shark (Isurus oxyrinchus) with estimates of field metabolic rate.

The shortfin mako shark is a large-bodied pursuit predator thought to be capable of the highest swimming speeds of any elasmobranch and potentially one of the highest energetic demands of any marine fish. Nonetheless, few direct speed measurements have been reported for this species. Here, animal-borne bio-loggers attached to two mako sharks were used to provide direct measurements of swimming speeds, kinematics and thermal physiology. Mean sustained (cruising) speed was 0.90 m s-1 (±0.07 s.d.) with a mean tail-beat frequency (TBF) of 0.51 Hz (±0.16 s.d.). The maximum burst speed recorded was 5.02 m s-1 (TBFmax = 3.65 Hz) from a 2 m long female. Burst swimming was sustained for 14 s (mean speed = 2.38 m s-1 ), leading to a 0.24°C increase in white muscle temperature in the 12.5 min after the burst. Routine field metabolic rate was estimated at 185.2 mg O2 kg-1 h-1 (at 18°C ambient temperature). Gliding behaviour (zero TBF) was more frequently observed after periods of high activity, especially after capture when internal (white muscle) temperature approached 21°C (ambient temperature: 18.3°C), indicating gliding probably functions as an energy recovery mechanism and limits further metabolic heat production. The results show shortfin mako sharks generally cruise at speeds similar to other endothermic fish - but faster than ectothermic sharks - with the maximum recorded burst speed being among the highest so far directly measured among sharks, tunas and billfishes. This newly recorded high-oxygen-demand performance of mako sharks suggests it may be particularly vulnerable to habitat loss due to climate-driven ocean deoxygenation.

Building use-inspired species distribution models: Using multiple data types to examine and improve model performance.

Species distribution models (SDMs) are becoming an important tool for marine conservation and management. Yet while there is an increasing diversity and volume of marine biodiversity data for training SDMs, little practical guidance is available on how to leverage distinct data types to build robust models. We explored the effect of different data types on the fit, performance and predictive ability of SDMs by comparing models trained with four data types for a heavily exploited pelagic fish, the blue shark (Prionace glauca), in the Northwest Atlantic: two fishery dependent (conventional mark-recapture tags, fisheries observer records) and two fishery independent (satellite-linked electronic tags, pop-up archival tags). We found that all four data types can result in robust models, but differences among spatial predictions highlighted the need to consider ecological realism in model selection and interpretation regardless of data type. Differences among models were primarily attributed to biases in how each data type, and the associated representation of absences, sampled the environment and summarized the resulting species distributions. Outputs from model ensembles and a model trained on all pooled data both proved effective for combining inferences across data types and provided more ecologically realistic predictions than individual models. Our results provide valuable guidance for practitioners developing SDMs. With increasing access to diverse data sources, future work should further develop truly integrative modeling approaches that can explicitly leverage the strengths of individual data types while statistically accounting for limitations, such as sampling biases.

Juvenile survival and movements of two threatened oceanic sharks in the North Atlantic Ocean inferred from tag-recovery data.

Understanding population dynamics, movements, and fishing mortality is critical to establish effective shark conservation measures across international boundaries in the ocean. There are few survival and dispersal estimates of juveniles of oceanic shark species in the North Atlantic despite it being one of the most fished regions in the world. Here we provide estimates of dispersal, survival, and proportion of fishing mortality in the North Atlantic for two threatened oceanic sharks: the blue shark (Prionace glauca) and the shortfin mako shark (Isurus oxyrinchus). Our results are based on multi-event models applied to tag-recovery data of 700 blue sharks and 132 shortfin makos tagged over a decade. A total of 60 blue sharks (8.57% of tagged) and 30 makos (22.73%) were recovered by the longline fishery between 2009 and 2017. Tag-reporting rate (percentage of returned information when a tagged shark was caught) was estimated to be high (0.794 ± 0.232 SE). Mean annual survival, as predicted from the models, was higher for blue shark (0.835 ± 0.040 SE) than for shortfin mako (0.618 ± 0.189 SE). Models predicted that fishing caused more than a half of total mortality in the study area for both species (0.576 ± 0.209), and more than a third of tagged individuals dispersed from the study area permanently (0.359 ± 0.073). Our findings, focused mainly on juveniles from oceanic areas, contribute to a better understanding of shark population dynamics in the North Atlantic and highlight the need for further conservation measures for both blue shark and shortfin mako, such as implementing efficient bycatch mitigation measures and static/dynamic time-area closures in the open ocean.

Circles in the sea: annual courtship "torus" behaviour of basking sharks Cetorhinus maximus identified in the eastern North Atlantic Ocean.

Groups of basking sharks engaged in circling behaviour are rarely observed, and their function remains enigmatic in the absence of detailed observations. Here, underwater and aerial video recordings of multiple circling groups of basking sharks during late summer (August and September 2016-2021) in the eastern North Atlantic Ocean showed groups numbering between 6 and 23 non-feeding individuals of both sexes. Sharks swam slowly in a rotating "torus" (diameter range: 17-39 m), with individuals layered vertically from the surface to a maximum depth of 16 m. Within a torus, sharks engaged in close-following, echelon, close-flank approach or parallel-swimming behaviours. Measured shark total body lengths were 5.4-9.5 m (mean LT : 7.3 m ± 0.9 s.d.; median: 7.2 m, n = 27), overlapping known lengths of sexually mature males and females. Males possessed large claspers with abrasions that were also observed on female pectoral fins. Female body colouration was paler than that of males, similar to colour changes observed during courtship and mating in other shark species. Individuals associated with most other members rapidly (within minutes), indicating toroidal behaviours facilitate multiple interactions. Sharks interacted through fin-fin and fin-body contacts, rolling to expose the ventral surfaces to following sharks, and breaching behaviour. Toruses formed in late summer when feeding aggregations in zooplankton-rich thermal fronts switched to non-feeding following and circling behaviours. Collectively, the observations explain a courtship function for toruses. This study highlights northeast Atlantic coastal waters as a critical habitat supporting courtship reproductive behaviour of endangered basking sharks, the first such habitat identified for this species globally.

Meta-Fish-Lib: A generalised, dynamic DNA reference library pipeline for metabarcoding of fishes.

The accuracy and reliability of DNA metabarcoding analyses depend on the breadth and quality of the reference libraries that underpin them. However, there are limited options available to obtain and curate the huge volumes of sequence data that are available on public repositories such as NCBI and BOLD. Here, we provide a pipeline to download, clean and annotate mitochondrial DNA sequence data for a given list of fish species. Features of this pipeline include (a) support for multiple metabarcode markers; (b) searches on species synonyms and taxonomic name validation; (c) phylogeny assisted quality control for identification and removal of misannotated sequences; (d) automatically generated coverage reports for each new GenBank release update; and (e) citable, versioned DOIs. As an example we provide a ready-to-use curated reference library for the marine and freshwater fishes of the U.K. To augment this reference library for environmental DNA metabarcoding specifically, we generated 241 new MiFish-12S sequences for 88 U.K. marine species, and make available new primer sets useful for sequencing these. This brings the coverage of common U.K. species for the MiFish-12S fragment to 93%, opening new avenues for scaling up fish metabarcoding across wide spatial gradients. The Meta-Fish-Lib reference library and pipeline is hosted at https://github.com/genner-lab/meta-fish-lib.

Climate-driven deoxygenation elevates fishing vulnerability for the ocean's widest ranging shark.

Climate-driven expansions of ocean hypoxic zones are predicted to concentrate pelagic fish in oxygenated surface layers, but how expanding hypoxia and fisheries will interact to affect threatened pelagic sharks remains unknown. Here, analysis of satellite-tracked blue sharks and environmental modelling in the eastern tropical Atlantic oxygen minimum zone (OMZ) shows shark maximum dive depths decreased due to combined effects of decreasing dissolved oxygen (DO) at depth, high sea surface temperatures, and increased surface-layer net primary production. Multiple factors associated with climate-driven deoxygenation contributed to blue shark vertical habitat compression, potentially increasing their vulnerability to surface fisheries. Greater intensity of longline fishing effort occurred above the OMZ compared to adjacent waters. Higher shark catches were associated with strong DO gradients, suggesting potential aggregation along suitable DO gradients contributed to habitat compression and higher fishing-induced mortality. Fisheries controls to counteract deoxygenation effects on shark catches will be needed as oceans continue warming.

Spatio-temporal genetic tagging of a cosmopolitan planktivorous shark provides insight to gene flow, temporal variation and site-specific re-encounters.

Migratory movements in response to seasonal resources often influence population structure and dynamics. Yet in mobile marine predators, population genetic consequences of such repetitious behaviour remain inaccessible without comprehensive sampling strategies. Temporal genetic sampling of seasonally recurring aggregations of planktivorous basking sharks, Cetorhinus maximus, in the Northeast Atlantic (NEA) affords an opportunity to resolve individual re-encounters at key sites with population connectivity and patterns of relatedness. Genetic tagging (19 microsatellites) revealed 18% of re-sampled individuals in the NEA demonstrated inter/multi-annual site-specific re-encounters. High genetic connectivity and migration between aggregation sites indicate the Irish Sea as an important movement corridor, with a contemporary effective population estimate (Ne) of 382 (CI = 241-830). We contrast the prevailing view of high gene flow across oceanic regions with evidence of population structure within the NEA, with early-season sharks off southwest Ireland possibly representing genetically distinct migrants. Finally, we found basking sharks surfacing together in the NEA are on average more related than expected by chance, suggesting a genetic consequence of, or a potential mechanism maintaining, site-specific re-encounters. Long-term temporal genetic monitoring is paramount in determining future viability of cosmopolitan marine species, identifying genetic units for conservation management, and for understanding aggregation structure and dynamics.

Optimal searching behaviour generated intrinsically by the central pattern generator for locomotion.

Efficient searching for resources such as food by animals is key to their survival. It has been proposed that diverse animals from insects to sharks and humans adopt searching patterns that resemble a simple Lévy random walk, which is theoretically optimal for 'blind foragers' to locate sparse, patchy resources. To test if such patterns are generated intrinsically, or arise via environmental interactions, we tracked free-moving Drosophila larvae with (and without) blocked synaptic activity in the brain, suboesophageal ganglion (SOG) and sensory neurons. In brain-blocked larvae, we found that extended substrate exploration emerges as multi-scale movement paths similar to truncated Lévy walks. Strikingly, power-law exponents of brain/SOG/sensory-blocked larvae averaged 1.96, close to a theoretical optimum (µ ≅ 2.0) for locating sparse resources. Thus, efficient spatial exploration can emerge from autonomous patterns in neural activity. Our results provide the strongest evidence so far for the intrinsic generation of Lévy-like movement patterns.

Translating Marine Animal Tracking Data into Conservation Policy and Management.

There have been efforts around the globe to track individuals of many marine species and assess their movements and distribution, with the putative goal of supporting their conservation and management. Determining whether, and how, tracking data have been successfully applied to address real-world conservation issues is, however, difficult. Here, we compile a broad range of case studies from diverse marine taxa to show how tracking data have helped inform conservation policy and management, including reductions in fisheries bycatch and vessel strikes, and the design and administration of marine protected areas and important habitats. Using these examples, we highlight pathways through which the past and future investment in collecting animal tracking data might be better used to achieve tangible conservation benefits.

Persistence of environmental DNA in marine systems.

As environmental DNA (eDNA) becomes an increasingly valuable resource for marine ecosystem monitoring, understanding variation in its persistence across contrasting environments is critical. Here, we quantify the breakdown of macrobial eDNA over a spatio-temporal axis of locally extreme conditions, varying from ocean-influenced offshore to urban-inshore, and between winter and summer. We report that eDNA degrades 1.6 times faster in the inshore environment than the offshore environment, but contrary to expectation we find no difference over season. Analysis of environmental covariables show a spatial gradient of salinity and a temporal gradient of pH, with salinity-or the biotic correlates thereof-most important. Based on our estimated inshore eDNA half-life and naturally occurring eDNA concentrations, we estimate that eDNA may be detected for around 48 h, offering potential to collect ecological community data of high local fidelity. We conclude by placing these results in the context of previously published eDNA decay rates.