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.

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.

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.

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.

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.

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.

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.

Environmental context explains Lévy and Brownian movement patterns of marine predators.

An optimal search theory, the so-called Lévy-flight foraging hypothesis, predicts that predators should adopt search strategies known as Lévy flights where prey is sparse and distributed unpredictably, but that Brownian movement is sufficiently efficient for locating abundant prey. Empirical studies have generated controversy because the accuracy of statistical methods that have been used to identify Lévy behaviour has recently been questioned. Consequently, whether foragers exhibit Lévy flights in the wild remains unclear. Crucially, moreover, it has not been tested whether observed movement patterns across natural landscapes having different expected resource distributions conform to the theory's central predictions. Here we use maximum-likelihood methods to test for Lévy patterns in relation to environmental gradients in the largest animal movement data set assembled for this purpose. Strong support was found for Lévy search patterns across 14 species of open-ocean predatory fish (sharks, tuna, billfish and ocean sunfish), with some individuals switching between Lévy and Brownian movement as they traversed different habitat types. We tested the spatial occurrence of these two principal patterns and found Lévy behaviour to be associated with less productive waters (sparser prey) and Brownian movements to be associated with productive shelf or convergence-front habitats (abundant prey). These results are consistent with the Lévy-flight foraging hypothesis, supporting the contention that organism search strategies naturally evolved in such a way that they exploit optimal Lévy patterns.

Foraging success of biological Lévy flights recorded in situ.

It is an open question how animals find food in dynamic natural environments where they possess little or no knowledge of where resources are located. Foraging theory predicts that in environments with sparsely distributed target resources, where forager knowledge about resources' locations is incomplete, Lévy flight movements optimize the success of random searches. However, the putative success of Lévy foraging has been demonstrated only in model simulations. Here, we use high-temporal-resolution Global Positioning System (GPS) tracking of wandering (Diomedea exulans) and black-browed albatrosses (Thalassarche melanophrys) with simultaneous recording of prey captures, to show that both species exhibit Lévy and Brownian movement patterns. We find that total prey masses captured by wandering albatrosses during Lévy movements exceed daily energy requirements by nearly fourfold, and approached yields by Brownian movements in other habitats. These results, together with our reanalysis of previously published albatross data, overturn the notion that albatrosses do not exhibit Lévy patterns during foraging, and demonstrate that Lévy flights of predators in dynamic natural environments present a beneficial alternative strategy to simple, spatially intensive behaviors. Our findings add support to the possibility that biological Lévy flight may have naturally evolved as a search strategy in response to sparse resources and scant information.

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.

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 marine predator search behaviour.

Many free-ranging predators have to make foraging decisions with little, if any, knowledge of present resource distribution and availability. The optimal search strategy they should use to maximize encounter rates with prey in heterogeneous natural environments remains a largely unresolved issue in ecology. Lévy walks are specialized random walks giving rise to fractal movement trajectories that may represent an optimal solution for searching complex landscapes. However, the adaptive significance of this putative strategy in response to natural prey distributions remains untested. Here we analyse over a million movement displacements recorded from animal-attached electronic tags to show that diverse marine predators-sharks, bony fishes, sea turtles and penguins-exhibit Lévy-walk-like behaviour close to a theoretical optimum. Prey density distributions also display Lévy-like fractal patterns, suggesting response movements by predators to prey distributions. Simulations show that predators have higher encounter rates when adopting Lévy-type foraging in natural-like prey fields compared with purely random landscapes. This is consistent with the hypothesis that observed search patterns are adapted to observed statistical patterns of the landscape. This may explain why Lévy-like behaviour seems to be widespread among diverse organisms, from microbes to humans, as a 'rule' that evolved in response to patchy resource distributions.

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.

High activity and Levy searches: jellyfish can search the water column like fish.

Over-fishing may lead to a decrease in fish abundance and a proliferation of jellyfish. Active movements and prey search might be thought to provide a competitive advantage for fish, but here we use data-loggers to show that the frequently occurring coastal jellyfish (Rhizostoma octopus) does not simply passively drift to encounter prey. Jellyfish (327 days of data from 25 jellyfish with depth collected every 1 min) showed very dynamic vertical movements, with their integrated vertical movement averaging 619.2 m d(-1), more than 60 times the water depth where they were tagged. The majority of movement patterns were best approximated by exponential models describing normal random walks. However, jellyfish also showed switching behaviour from exponential patterns to patterns best fitted by a truncated Lévy distribution with exponents (mean µ=1.96, range 1.2-2.9) close to the theoretical optimum for searching for sparse prey (µopt≈2.0). Complex movements in these 'simple' animals may help jellyfish to compete effectively with fish for plankton prey, which may enhance their ability to increase in dominance in perturbed ocean systems.

Lévy flight and Brownian search patterns of a free-ranging predator reflect different prey field characteristics.

1. Search processes play an important role in physical, chemical and biological systems. In animal foraging, the search strategy predators should use to search optimally for prey is an enduring question. Some models demonstrate that when prey is sparsely distributed, an optimal search pattern is a specialised random walk known as a Lévy flight, whereas when prey is abundant, simple Brownian motion is sufficiently efficient. These predictions form part of what has been termed the Lévy flight foraging hypothesis (LFF) which states that as Lévy flights optimise random searches, movements approximated by optimal Lévy flights may have naturally evolved in organisms to enhance encounters with targets (e.g. prey) when knowledge of their locations is incomplete. 2. Whether free-ranging predators exhibit the movement patterns predicted in the LFF hypothesis in response to known prey types and distributions, however, has not been determined. We tested this using vertical and horizontal movement data from electronic tagging of an apex predator, the great white shark Carcharodon carcharias, across widely differing habitats reflecting different prey types. 3. Individual white sharks exhibited movement patterns that predicted well the prey types expected under the LFF hypothesis. Shark movements were best approximated by Brownian motion when hunting near abundant, predictable sources of prey (e.g. seal colonies, fish aggregations), whereas movements approximating truncated Lévy flights were present when searching for sparsely distributed or potentially difficult-to-detect prey in oceanic or shelf environments, respectively. 4. That movement patterns approximated by truncated Lévy flights and Brownian behaviour were present in the predicted prey fields indicates search strategies adopted by white sharks appear to be the most efficient ones for encountering prey in the habitats where such patterns are observed. This suggests that C. carcharias appears capable of exhibiting search patterns that are approximated as optimal in response to encountered changes in prey type and abundance, and across diverse marine habitats, from the surf zone to the deep ocean. 5. Our results provide some support for the LFF hypothesis. However, it is possible that the observed Lévy patterns of white sharks may not arise from an adaptive behaviour but could be an emergent property arising from simple, straight-line movements between complex (e.g. fractal) distributions of prey. Experimental studies are needed in vertebrates to test for the presence of Lévy behaviour patterns in the absence of complex prey distributions.

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.