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.

Environmental drivers of biogeography and community structure in a Mid-Atlantic estuary.

Estuaries include some of the most productive yet anthropogenically impacted marine ecosystems on the planet, and provide critical habitat to many ecologically and economically important marine species. In order to elucidate ecological function in estuaries, we must understand what factors drive community dynamics. Delaware Bay is the third largest estuary in the United States and hosts over 200 species of migrant and resident fishes and invertebrates. The Delaware Division of Fish and Wildlife has conducted two long-term trawl surveys at monthly intervals in Delaware Bay since 1966. The two surveys collect data on environmental conditions, species composition, and number of fishes and macroinvertebrates across different size classes and life histories. Using a suite of multivariate approaches including hierarchical cluster analysis, canonical correlation analysis, and permutational multivariate analysis of variance, we characterized the fish and macroinvertebrate community in Delaware Bay and found that community composition and environmental conditions varied across spatial and seasonal scales. We identified four distinct biogeographic regions, based on environmental conditions and community composition, which were consistent across surveys. We found that the community was driven primarily by gradients in temperature and salinity and that abundant, frequently occurring species in the Bay have well-defined environmental associations. Our work represents the first attempt to use an existing historical survey to better understand how environmental parameters influence diversity and distribution of macrofauna within Delaware Bay, providing insight into how abiotic variables, influenced by climate, may impact the Delaware Bay ecosystem and similar estuarine ecosystems worldwide.

Evidence for dynamic resource partitioning between two sympatric reef shark species within the British Indian Ocean Territory.

Stable-isotope analyses (δ13 C, δ15 N and δ34 S) of multiple tissues (fin, muscle, red blood cells and plasma), revealed ontogenetic shifts in resource use by grey reef sharks Carcharhinus amblyrhynchos and resource partitioning with silvertip sharks Carcharhinus albimarginatus within the British Indian Ocean Territory marine protected area (MPA). Resource partitioning varied temporally, with C. albimarginatus feeding on more pelagic prey during October to January, potentially attributable to an influx of pelagic prey from outside the MPA at that time. Reef sharks may therefore be affected by processes outside an MPA, even if the sharks do not leave the MPA.

Stable isotope analysis reveals ontogenetic feeding shifts in Pacific blue marlin (Makaira nigricans) off eastern Taiwan.

To gain a better understanding of the trophic ecology of Pacific blue marlin Makaira nigricans off eastern Taiwan, nitrogen and carbon stable isotopes (δ15 N and δ13 C) and Bayesian mixing models were used to explore trophic dynamics and potential ontogenetic feeding shifts across M. nigricans of different size classes. Makaira nigricans samples from east of Taiwan (n = 213) and Palau (n = 37), as well as their prey (n = 70), were collected during 2012 and 2013. Results indicated increases in δ15 N with size, with values of larger size classes (> 200 cm eye-to-fork length; LEF ) significantly higher than those < 200 cm LEF . Values of δ13 C were negatively correlated with size. Makaira nigricans > 200 cm LEF had the highest estimated trophic position (4.44) and also exhibited ontogenetic changes in trophic position. Large M. nigricans fed more on dolphinfish Coryphaena hippurus and hairtail Trichiurus lepturus, while smaller M. nigricans consumed smaller forage fish (e.g., moonfish Mene maculata) and cephalopods. These changes may relate to greater swimming speeds and vertical habitat use in larger M. nigricans, allowing capture and consumption of larger prey items at higher trophic positions. The high trophic level of M. nigricans east of Taiwan confirms its important role as an apex predator in marine food webs and how ecological role changes with size.

Climate mediates hypoxic stress on fish diversity and nursery function at the land-sea interface.

Coastal ecosystems provide numerous important ecological services, including maintenance of biodiversity and nursery grounds for many fish species of ecological and economic importance. However, human population growth has led to increased pollution, ocean warming, hypoxia, and habitat alteration that threaten ecosystem services. In this study, we used long-term datasets of fish abundance, water quality, and climatic factors to assess the threat of hypoxia and the regulating effects of climate on fish diversity and nursery conditions in Elkhorn Slough, a highly eutrophic estuary in central California (United States), which also serves as a biodiversity hot spot and critical nursery grounds for offshore fisheries in a broader region. We found that hypoxic conditions had strong negative effects on extent of suitable fish habitat, fish species richness, and abundance of the two most common flatfish species, English sole (Parophrys vetulus) and speckled sanddab (Citharichthys stigmaeus). The estuary serves as an important nursery ground for English sole, making this species vulnerable to anthropogenic threats. We determined that estuarine hypoxia was associated with significant declines in English sole nursery habitat, with cascading effects on recruitment to the offshore adult population and fishery, indicating that human land use activities can indirectly affect offshore fisheries. Estuarine hypoxic conditions varied spatially and temporally and were alleviated by strengthening of El Niño conditions through indirect pathways, a consistent result in most estuaries across the northeast Pacific. These results demonstrate that changes to coastal land use and climate can fundamentally alter the diversity and functioning of coastal nurseries and their adjacent ocean ecosystems.

Assessing niche width of endothermic fish from genes to ecosystem.

Endothermy in vertebrates has been postulated to confer physiological and ecological advantages. In endothermic fish, niche expansion into cooler waters is correlated with specific physiological traits and is hypothesized to lead to greater foraging success and increased fitness. Using the seasonal co-occurrence of three tuna species in the eastern Pacific Ocean as a model system, we used cardiac gene expression data (as a proxy for thermal tolerance to low temperatures), archival tag data, and diet analyses to examine the vertical niche expansion hypothesis for endothermy in situ. Yellowfin, albacore, and Pacific bluefin tuna (PBFT) in the California Current system used more surface, mesopelagic, and deep waters, respectively. Expression of cardiac genes for calcium cycling increased in PBFT and coincided with broader vertical and thermal niche utilization. However, the PBFT diet was less diverse and focused on energy-rich forage fishes but did not show the greatest energy gains. Ecosystem-based management strategies for tunas should thus consider species-specific differences in physiology and foraging specialization.

Diet Composition and Trophic Ecology of Northeast Pacific Ocean Sharks.

Although there is a general perception of sharks as large pelagic, apex predators, most sharks are smaller, meso- and upper-trophic level predators that are associated with the seafloor. Among 73 shark species documented in the eastern North Pacific (ENP), less than half reach maximum lengths >200cm, and 78% occur in demersal or benthic regions of the continental shelf or slope. Most small (≤200cm) species (e.g., houndsharks) and demersal, nearshore juveniles of larger species (e.g., requiem sharks) consume small teleosts and decapod crustaceans, whereas large species in pelagic coastal and oceanic environments feed on large teleosts and squids. Several large, pelagic apex predator species occur in the ENP, but the largest species (i.e., Basking Shark, Whale Shark) consume zooplankton or small nekton. Size-based dietary variability is substantial for many species, and segregation of juvenile and adult foraging habitats also is common (e.g., Horn Shark, Shortfin Mako). Temporal dietary differences are most pronounced for temperate, nearshore species with wide size ranges, and least pronounced for smaller species in extreme latitudes and deep-water regions. Sympatric sharks often occupy various trophic positions, with resource overlap differing by space and time and some sharks serving as prey to other species. Most coastal species remain in the same general region over time and feed opportunistically on variable prey inputs (e.g., season migrations, spawning, or recruitment events), whereas pelagic, oceanic species actively seek hot spots of prey abundance that are spatiotemporally variable. The influence of sharks on ecosystem structure and regulation has been downplayed compared to that of large teleosts species with higher per capita consumption rates (e.g., tunas, billfishes). However, sharks also exert indirect influences on prey populations by causing behavioural changes that may result in restricted ranges and reduced fitness. Except for food web modelling efforts in Alaskan waters, the trophic impacts of sharks are poorly incorporated into current ecosystem approaches to fisheries management in the NEP.

Stable isotope analysis of vertebrae reveals ontogenetic changes in habitat in an endothermic pelagic shark.

Ontogenetic changes in habitat are driven by shifting life-history requirements and play an important role in population dynamics. However, large portions of the life history of many pelagic species are still poorly understood or unknown. We used a novel combination of stable isotope analysis of vertebral annuli, Bayesian mixing models, isoscapes and electronic tag data to reconstruct ontogenetic patterns of habitat and resource use in a pelagic apex predator, the salmon shark (Lamna ditropis). Results identified the North Pacific Transition Zone as the major nursery area for salmon sharks and revealed an ontogenetic shift around the age of maturity from oceanic to increased use of neritic habitats. The nursery habitat may reflect trade-offs between prey availability, predation pressure and thermal constraints on juvenile endothermic sharks. The ontogenetic shift in habitat coincided with a reduction of isotopic niche, possibly reflecting specialization upon particular prey or habitats. Using tagging data to inform Bayesian isotopic mixing models revealed that adult sharks primarily use neritic habitats of Alaska yet receive a trophic subsidy from oceanic habitats. Integrating the multiple methods used here provides a powerful approach to retrospectively study the ecology and life history of migratory species throughout their ontogeny.

Reconstructing transoceanic migration patterns of Pacific bluefin tuna using a chemical tracer toolbox.

Large pelagic predators play important roles in oceanic ecosystems, and may migrate vast distances to utilize resources in different marine ecoregions. Understanding movement patterns of migratory marine animals is critical for effective management, but often challenging, due to the cryptic habitat of pelagic migrators and the difficulty of assessing past movements. Chemical tracers can partially circumvent these challenges by reconstructing recent migration patterns. Pacific bluefin tuna (Thunnus orientalis; PBFT) inhabit the western and eastern Pacific Ocean, and are in steep decline due to overfishing. Understanding age-specific eastward transpacific migration patterns can improve management practices, but these migratory dynamics remain largely unquantified. Here, we combine a Fukushima-derived radiotracer (134Cs) with bulk tissue and amino acid stable isotope analyses of PBFT to distinguish recent migrants from residents of the eastern Pacific Ocean. The proportion of recent migrants to residents decreased in older year classes, though the proportion of older PBFT that recently migrated across the Pacific was greater than previous estimates. This novel toolbox of biogeochemical tracers can be applied to any species that crosses the North Pacific Ocean.

Novel CTD tag establishes shark fins as ocean observing platforms.

Animal-borne tags are effective instruments for collecting ocean data and can be used to fill spatial gaps in the observing network. We deployed the first conductivity, temperature, and depth (CTD) satellite tags on the dorsal fin of salmon sharks (Lamna ditropis) to demonstrate the potential of sharks to monitor essential ocean variables and oceanographic features in the Gulf of Alaska. Over 1360 km and 36 days in the summer of 2015, the salmon shark collected 56 geolocated, temperature-salinity profiles. The shark swam through a plume of anomalously salty water that originated from the "Blob" and encountered several mesoscale eddies, whose subsurface properties were altered by the marine heatwave. We demonstrate that salmon sharks have the potential to serve as submesoscale-resolving oceanographic platforms and substantially increase the spatial coverage of observations in the Gulf of Alaska.

Environmental stress reduces shark residency to coral reefs.

Coral reef ecosystems are highly threatened and can be extremely sensitive to the effects of climate change. Multiple shark species rely on coral reefs as important habitat and, as such, play a number of significant ecological roles in these ecosystems. How environmental stress impacts routine, site-attached reef shark behavior, remains relatively unexplored. Here, we combine 8 years of acoustic tracking data (2013-2020) from grey reef sharks resident to the remote coral reefs of the Chagos Archipelago in the Central Indian Ocean, with a satellite-based index of coral reef environmental stress exposure. We show that on average across the region, increased stress on the reefs significantly reduces grey reef shark residency, promoting more diffuse space use and increasing time away from shallow forereefs. Importantly, this impact has a lagged effect for up to 16 months. This may have important physiological and conservation consequences for reef sharks, as well as broader implications for reef ecosystem functioning. As climate change is predicted to increase environmental stress on coral reef ecosystems, understanding how site-attached predators respond to stress will be crucial for forecasting the functional significance of altering predator behavior and the potential impacts on conservation for both reef sharks and coral reefs themselves.

Integrating multiple chemical tracers to elucidate the diet and habitat of Cookiecutter Sharks.

The Cookiecutter shark (Isistius brasiliensis) is an ectoparasitic, mesopelagic shark that is known for removing plugs of tissue from larger prey, including teleosts, chondrichthyans, cephalopods, and marine mammals. Although this species is widely distributed throughout the world's tropical and subtropical oceanic waters, like many deep-water species, it remains very poorly understood due to its mesopelagic distribution. We used a suite of biochemical tracers, including stable isotope analysis (SIA), fatty acid analysis (FAA), and environmental DNA (eDNA), to investigate the trophic ecology of this species in the Central Pacific around Hawaii. We found that large epipelagic prey constituted a relatively minor part of the overall diet. Surprisingly, small micronektonic and forage species (meso- and epipelagic) are the most important prey group for Cookiecutter sharks across the studied size range (17-43 cm total length), with larger mesopelagic species or species that exhibit diel vertical migration also being important prey. These results were consistent across all the tracer techniques employed. Our results indicate that Cookiecutter sharks play a unique role in pelagic food webs, feeding on prey ranging from the largest apex predators to small, low trophic level species, in particular those that overlap with the depth distribution of the sharks throughout the diel cycle. We also found evidence of a potential shift in diet and/or habitat with size and season. Environmental DNA metabarcoding revealed new prey items for Cookiecutter sharks while also demonstrating that eDNA can be used to identify recent prey in stomachs frozen for extended periods. Integrating across chemical tracers is a powerful tool for investigating the ecology of elusive and difficult to study species, such as meso- and bathypelagic chondrichthyans, and can increase the amount of information gained from small sample sizes. Better resolving the foraging ecology of these mesopelagic predators is critical for effective conservation and management of these taxa and ecosystems, which are intrinsically vulnerable to overfishing and exploitation.

Analysing detection gaps in acoustic telemetry data to infer differential movement patterns in fish.

A wide array of technologies are available for gaining insight into the movement of wild aquatic animals. Although acoustic telemetry can lack the fine-scale spatial resolution of some satellite tracking technologies, the substantially longer battery life can yield important long-term data on individual behavior and movement for low per-unit cost. Typically, however, receiver arrays are designed to maximize spatial coverage at the cost of positional accuracy leading to potentially longer detection gaps as individuals move out of range between monitored locations. This is particularly true when these technologies are deployed to monitor species in hard-to-access locations.Here, we develop a novel approach to analyzing acoustic telemetry data, using the timing and duration of gaps between animal detections to infer different behaviors. Using the durations between detections at the same and different receiver locations (i.e., detection gaps), we classify behaviors into "restricted" or potential wider "out-of-range" movements synonymous with longer distance dispersal. We apply this method to investigate spatial and temporal segregation of inferred movement patterns in two sympatric species of reef shark within a large, remote, marine protected area (MPA). Response variables were generated using network analysis, and drivers of these movements were identified using generalized linear mixed models and multimodel inference.Species, diel period, and season were significant predictors of "out-of-range" movements. Silvertip sharks were overall more likely to undertake "out-of-range" movements, compared with gray reef sharks, indicating spatial segregation, and corroborating previous stable isotope work between these two species. High individual variability in "out-of-range" movements in both species was also identified.We present a novel gap analysis of telemetry data to help infer differential movement and space use patterns where acoustic coverage is imperfect and other tracking methods are impractical at scale. In remote locations, inference may be the best available tool and this approach shows that acoustic telemetry gap analysis can be used for comparative studies in fish ecology, or combined with other research techniques to better understand functional mechanisms driving behavior.

Predicted hotspots of overlap between highly migratory fishes and industrial fishing fleets in the northeast Pacific.

Many species of sharks and some tunas are threatened by overexploitation, yet the degree of overlap between industrial fisheries and pelagic fishes remains poorly understood. Using satellite tracks from 933 industrial fishing vessels and predictive habitat models from 876 electronic tags deployed on seven shark and tuna species, we developed fishing effort maps across the northeast Pacific Ocean and assessed overlap with core habitats of pelagic fishes. Up to 35% of species' core habitats overlapped with fishing effort. We identified overlap hotspots along the North American shelf, the equatorial Pacific, and the subtropical gyre. Results indicate where species require international conservation efforts and effective management within national waters. Only five national fleets (Mexico, Taiwan, China, Japan, and the United States) account for >90% of overlap with core habitats of our focal sharks and tunas on the high seas. These results inform global negotiations to achieve sustainability on the high seas.

The political biogeography of migratory marine predators.

During their migrations, marine predators experience varying levels of protection and face many threats as they travel through multiple countries' jurisdictions and across ocean basins. Some populations are declining rapidly. Contributing to such declines is a failure of some international agreements to ensure effective cooperation by the stakeholders responsible for managing species throughout their ranges, including in the high seas, a global commons. Here we use biologging data from marine predators to provide quantitative measures with great potential to inform local, national and international management efforts in the Pacific Ocean. We synthesized a large tracking data set to show how the movements and migratory phenology of 1,648 individuals representing 14 species-from leatherback turtles to white sharks-relate to the geopolitical boundaries of the Pacific Ocean throughout species' annual cycles. Cumulatively, these species visited 86% of Pacific Ocean countries and some spent three-quarters of their annual cycles in the high seas. With our results, we offer answers to questions posed when designing international strategies for managing migratory species.

Influence of ontogeny and environmental exposure on mercury accumulation in muscle and liver of male Round Stingrays.

Mercury tissue distribution and its dynamics are poorly understood in elasmobranchs. Total mercury was measured in liver and muscle of male Round Stingrays (Urobatis halleri) from Seal Beach, California, an anthropogenically impacted site, and from the offshore island of Santa Catalina, a less impacted site. Stable isotope analysis was also performed on the muscle and red blood cells (RBCs) of a subset of rays over a range of age classes to investigate mercury accumulation with respect to trophic ecology. Mercury in both tissues was found to be significantly greater in adults than juveniles in mainland rays; however, liver mercury accumulation drastically increased after maturity and was significantly greater in mainland adult rays than Catalina rays. There were no patterns in δ15N or δ13C with size in muscle; however, there were indications of seasonal changes in RBC δ15N, suggesting short term shifts in diet or behavior is likely linked to reproductive status as is mercury accumulation.

Oceanographic drivers of the vertical distribution of a highly migratory, endothermic shark.

Salmon sharks Lamna ditropis are highly migratory, upper trophic level predators in North Pacific ecosystems. We analysed a multi-year satellite tag dataset to investigate the habitat use of female salmon sharks across their broad range in the eastern North Pacific (NEP) and identified key environmental factors that influence vertical distribution. Salmon sharks displayed remarkable plasticity in habitat use across disparate oceanographic regions in the NEP and increased utilization of deeper waters in offshore habitats. Diel shifts in vertical distribution and behaviour were consistently observed across their range and likely reflect shifts in their foraging ecology. Salmon sharks utilized a broad thermal niche and exhibited submergence behaviour, possibly for thermal refuge, when encountering sea surface temperatures outside their preferred temperature distribution. Moreover, the vertical distribution of salmon sharks indicates they were able to exploit low dissolved oxygen environments (<1-3 ml l-1), occasionally for extended periods of time in offshore habitats. However, salmon sharks generally reduced their use of deeper waters when encountering the combination of cold temperatures (<6 °C) and low dissolved oxygen concentrations (<1-3 ml l-1). Combining vertical distribution with high-resolution horizontal movements furthers our understanding of the ecological and environmental drivers of movement across short (diel) and long-term (migratory) scales.

Quantifying overlap between the Deepwater Horizon oil spill and predicted bluefin tuna spawning habitat in the Gulf of Mexico.

Atlantic bluefin tuna (Thunnus thynnus) are distributed throughout the North Atlantic and are both economically valuable and heavily exploited. The fishery is currently managed as two spawning populations, with the GOM population being severely depleted for over 20 years. In April-August of 2010, the Deepwater Horizon oil spill released approximately 4 million barrels of oil into the GOM, with severe ecosystem and economic impacts. Acute oil exposure results in mortality of bluefin eggs and larvae, while chronic effects on spawning adults are less well understood. Here we used 16 years of electronic tagging data for 66 bluefin tuna to identify spawning events, to quantify habitat preferences, and to predict habitat use and oil exposure within Gulf of Mexico spawning grounds. More than 54,000 km2 (5%) of predicted spawning habitat within the US EEZ was oiled during the week of peak oil dispersal, with potentially lethal effects on eggs and larvae. Although the oil spill overlapped with a relatively small portion of predicted spawning habitat, the cumulative impact from oil, ocean warming and bycatch mortality on GOM spawning grounds may result in significant effects for a population that shows little evidence of rebuilding.