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.

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.

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.

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.

Global collision-risk hotspots of marine traffic and the world's largest fish, the whale shark.

Marine traffic is increasing globally yet collisions with endangered megafauna such as whales, sea turtles, and planktivorous sharks go largely undetected or unreported. Collisions leading to mortality can have population-level consequences for endangered species. Hence, identifying simultaneous space use of megafauna and shipping throughout ranges may reveal as-yet-unknown spatial targets requiring conservation. However, global studies tracking megafauna and shipping occurrences are lacking. Here we combine satellite-tracked movements of the whale shark, Rhincodon typus, and vessel activity to show that 92% of sharks' horizontal space use and nearly 50% of vertical space use overlap with persistent large vessel (>300 gross tons) traffic. Collision-risk estimates correlated with reported whale shark mortality from ship strikes, indicating higher mortality in areas with greatest overlap. Hotspots of potential collision risk were evident in all major oceans, predominantly from overlap with cargo and tanker vessels, and were concentrated in gulf regions, where dense traffic co-occurred with seasonal shark movements. Nearly a third of whale shark hotspots overlapped with the highest collision-risk areas, with the last known locations of tracked sharks coinciding with busier shipping routes more often than expected. Depth-recording tags provided evidence for sinking, likely dead, whale sharks, suggesting substantial "cryptic" lethal ship strikes are possible, which could explain why whale shark population declines continue despite international protection and low fishing-induced mortality. Mitigation measures to reduce ship-strike risk should be considered to conserve this species and other ocean giants that are likely experiencing similar impacts from growing global vessel traffic.

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.

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.

Global spatial risk assessment of sharks under the footprint of fisheries.

Effective ocean management and the conservation of highly migratory species depend on resolving the overlap between animal movements and distributions, and fishing effort. However, this information is lacking at a global scale. Here we show, using a big-data approach that combines satellite-tracked movements of pelagic sharks and global fishing fleets, that 24% of the mean monthly space used by sharks falls under the footprint of pelagic longline fisheries. Space-use hotspots of commercially valuable sharks and of internationally protected species had the highest overlap with longlines (up to 76% and 64%, respectively), and were also associated with significant increases in fishing effort. We conclude that pelagic sharks have limited spatial refuge from current levels of fishing effort in marine areas beyond national jurisdictions (the high seas). Our results demonstrate an urgent need for conservation and management measures at high-seas hotspots of shark space use, and highlight the potential of simultaneous satellite surveillance of megafauna and fishers as a tool for near-real-time, dynamic management.

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.

Oceanic adults, coastal juveniles: tracking the habitat use of whale sharks off the Pacific coast of Mexico.

Eight whale sharks tagged with pop-up satellite archival tags off the Gulf of California, Mexico, were tracked for periods of 14-134 days. Five of these sharks were adults, with four females visually assessed to be pregnant. At least for the periods they were tracked, juveniles remained in the Gulf of California while adults moved offshore into the eastern Pacific Ocean. We propose that parturition occurs in these offshore waters. Excluding two juveniles that remained in the shallow tagging area for the duration of tracking, all sharks spent 65 ± 20.7% (SD) of their time near the surface, even over deep water, often in association with frontal zones characterized by cool-water upwelling. While these six sharks all made dives into the meso- or bathypelagic zones, with two sharks reaching the maximum depth recordable by the tags (1285.8 m), time spent at these depths represented a small proportion of the overall tracks. Most deep dives (72.7%) took place during the day, particularly during the early morning and late afternoon. Pronounced habitat differences by ontogenetic stage suggest that adult whale sharks are less likely to frequent coastal waters after the onset of maturity.

Acoustic telemetry and network analysis reveal the space use of multiple reef predators and enhance marine protected area design.

Marine protected areas (MPAs) are commonly employed to protect ecosystems from threats like overfishing. Ideally, MPA design should incorporate movement data from multiple target species to ensure sufficient habitat is protected. We used long-term acoustic telemetry and network analysis to determine the fine-scale space use of five shark and one turtle species at a remote atoll in the Seychelles, Indian Ocean, and evaluate the efficacy of a proposed MPA. Results revealed strong, species-specific habitat use in both sharks and turtles, with corresponding variation in MPA use. Defining the MPA's boundary from the edge of the reef flat at low tide instead of the beach at high tide (the current best in Seychelles) significantly increased the MPA's coverage of predator movements by an average of 34%. Informed by these results, the larger MPA was adopted by the Seychelles government, demonstrating how telemetry data can improve shark spatial conservation by affecting policy directly.

DNA barcoding identifies a cosmopolitan diet in the ocean sunfish.

The ocean sunfish (Mola mola) is the world's heaviest bony fish reaching a body mass of up to 2.3 tonnes. However, the prey M. mola consumes to fuel this prodigious growth remains poorly known. Sunfish were thought to be obligate gelatinous plankton feeders, but recent studies suggest a more generalist diet. In this study, through molecular barcoding and for the first time, the diet of sunfish in the north-east Atlantic Ocean was characterised. Overall, DNA from the diet content of 57 individuals was successfully amplified, identifying 41 different prey items. Sunfish fed mainly on crustaceans and teleosts, with cnidarians comprising only 16% of the consumed prey. Although no adult fishes were sampled, we found evidence for an ontogenetic shift in the diet, with smaller individuals feeding mainly on small crustaceans and teleost fish, whereas the diet of larger fish included more cnidarian species. Our results confirm that smaller sunfish feed predominantly on benthic and on coastal pelagic species, whereas larger fish depend on pelagic prey. Therefore, sunfish is a generalist predator with a greater diversity of links in coastal food webs than previously realised. Its removal as fisheries' bycatch may have wider reaching ecological consequences, potentially disrupting coastal trophic interactions.

Ocean-wide tracking of pelagic sharks reveals extent of overlap with longline fishing hotspots.

Overfishing is arguably the greatest ecological threat facing the oceans, yet catches of many highly migratory fishes including oceanic sharks remain largely unregulated with poor monitoring and data reporting. Oceanic shark conservation is hampered by basic knowledge gaps about where sharks aggregate across population ranges and precisely where they overlap with fishers. Using satellite tracking data from six shark species across the North Atlantic, we show that pelagic sharks occupy predictable habitat hotspots of high space use. Movement modeling showed sharks preferred habitats characterized by strong sea surface-temperature gradients (fronts) over other available habitats. However, simultaneous Global Positioning System (GPS) tracking of the entire Spanish and Portuguese longline-vessel fishing fleets show an 80% overlap of fished areas with hotspots, potentially increasing shark susceptibility to fishing exploitation. Regions of high overlap between oceanic tagged sharks and longliners included the North Atlantic Current/Labrador Current convergence zone and the Mid-Atlantic Ridge southwest of the Azores. In these main regions, and subareas within them, shark/vessel co-occurrence was spatially and temporally persistent between years, highlighting how broadly the fishing exploitation efficiently "tracks" oceanic sharks within their space-use hotspots year-round. Given this intense focus of longliners on shark hotspots, our study argues the need for international catch limits for pelagic sharks and identifies a future role of combining fine-scale fish and vessel telemetry to inform the ocean-scale management of fisheries.

Repeated, long-distance migrations by a philopatric predator targeting highly contrasting ecosystems.

Long-distance movements of animals are an important driver of population spatial dynamics and determine the extent of overlap with area-focused human activities, such as fishing. Despite global concerns of declining shark populations, a major limitation in assessments of population trends or spatial management options is the lack of information on their long-term migratory behaviour. For a large marine predator, the tiger shark Galeocerdo cuvier, we show from individuals satellite-tracked for multiple years (up to 1101 days) that adult males undertake annually repeated, round-trip migrations of over 7,500 km in the northwest Atlantic. Notably, these migrations occurred between the highly disparate ecosystems of Caribbean coral reef regions in winter and high latitude oceanic areas in summer, with strong, repeated philopatry to specific overwintering insular habitat. Partial migration also occurred, with smaller, immature individuals displaying reduced migration propensity. Foraging may be a putative motivation for these oceanic migrations, with summer behaviour showing higher path tortuosity at the oceanic range extremes. The predictable migratory patterns and use of highly divergent ecosystems shown by male tiger sharks appear broadly similar to migrations seen in birds, reptiles and mammals, and highlight opportunities for dynamic spatial management and conservation measures of highly mobile sharks.

Hierarchical random walks in trace fossils and the origin of optimal search behavior.

Efficient searching is crucial for timely location of food and other resources. Recent studies show that diverse living animals use a theoretically optimal scale-free random search for sparse resources known as a Lévy walk, but little is known of the origins and evolution of foraging behavior and the search strategies of extinct organisms. Here, using simulations of self-avoiding trace fossil trails, we show that randomly introduced strophotaxis (U-turns)--initiated by obstructions such as self-trail avoidance or innate cueing--leads to random looping patterns with clustering across increasing scales that is consistent with the presence of Lévy walks. This predicts that optimal Lévy searches may emerge from simple behaviors observed in fossil trails. We then analyzed fossilized trails of benthic marine organisms by using a novel path analysis technique and find the first evidence, to our knowledge, of Lévy-like search strategies in extinct animals. Our results show that simple search behaviors of extinct animals in heterogeneous environments give rise to hierarchically nested Brownian walk clusters that converge to optimal Lévy patterns. Primary productivity collapse and large-scale food scarcity characterizing mass extinctions evident in the fossil record may have triggered adaptation of optimal Lévy-like searches. The findings suggest that Lévy-like behavior has been used by foragers since at least the Eocene but may have a more ancient origin, which might explain recent widespread observations of such patterns among modern taxa.

Optimal foraging strategies: Lévy walks balance searching and patch exploitation under a very broad range of conditions.

While evidence for optimal random search patterns, known as Lévy walks, in empirical movement data is mounting for a growing list of taxa spanning motile cells to humans, there is still much debate concerning the theoretical generality of Lévy walk optimisation. Here, using a new and robust simulation environment, we investigate in the most detailed study to date (24×10(6) simulations) the foraging and search efficiencies of 2-D Lévy walks with a range of exponents, target resource distributions and several competing models. We find strong and comprehensive support for the predictions of the Lévy flight foraging hypothesis and in particular for the optimality of inverse square distributions of move step-lengths across a much broader range of resource densities and distributions than previously realised. Further support for the evolutionary advantage of Lévy walk movement patterns is provided by an investigation into the 'feast and famine' effect, with Lévy foragers in heterogeneous environments experiencing fewer long 'famines' than other types of searchers. Therefore overall, optimal Lévy foraging results in more predictable resources in unpredictable environments.

Scaling laws of ambush predator 'waiting' behaviour are tuned to a common ecology.

The decisions animals make about how long to wait between activities can determine the success of diverse behaviours such as foraging, group formation or risk avoidance. Remarkably, for diverse animal species, including humans, spontaneous patterns of waiting times show random 'burstiness' that appears scale-invariant across a broad set of scales. However, a general theory linking this phenomenon across the animal kingdom currently lacks an ecological basis. Here, we demonstrate from tracking the activities of 15 sympatric predator species (cephalopods, sharks, skates and teleosts) under natural and controlled conditions that bursty waiting times are an intrinsic spontaneous behaviour well approximated by heavy-tailed (power-law) models over data ranges up to four orders of magnitude. Scaling exponents quantifying ratios of frequent short to rare very long waits are species-specific, being determined by traits such as foraging mode (active versus ambush predation), body size and prey preference. A stochastic-deterministic decision model reproduced the empirical waiting time scaling and species-specific exponents, indicating that apparently complex scaling can emerge from simple decisions. Results indicate temporal power-law scaling is a behavioural 'rule of thumb' that is tuned to species' ecological traits, implying a common pattern may have naturally evolved that optimizes move-wait decisions in less predictable natural environments.