Oceanic Diel Vertical Movement Patterns of Blue Sharks Vary With Water Temperature and Productivity to Change Vulnerability to Fishing

fmars-08-688076 July 5, 2021 Time: 19:23 # 1 ORIGINAL RESEARCH published: 09 July 2021 doi: 10.3389/fmars.2021.688076 Edited by: Mourier Johann, Institut de Recherche Pour le Développement (IRD), France Reviewed by: Daniel M. Coffey, Texas A&M University–Corpus Christi, United States Yuuki Watanabe, National Institute of Polar Research, Japan *Correspondence: Nuno Queiroz [email protected] Specialty section: This article was submitted to Marine Megafauna, a section of the journal Frontiers in Marine Science Received: 30 March 2021 Accepted: 18 June 2021 Published: 09 July 2021 Citation: Vedor M, Mucientes G, Hernández-Chan S, Rosa R, Humphries N, Sims DW and Queiroz N (2021) Oceanic Diel Vertical Movement Patterns of Blue Sharks Vary With Water Temperature and Productivity to Change Vulnerability to Fishing. Front. Mar. Sci. 8:688076. doi: 10.3389/fmars.2021.688076 Oceanic Diel Vertical Movement Patterns of Blue Sharks Vary With Water Temperature and Productivity to Change Vulnerability to Fishing Marisa Vedor1,2, Gonzalo Mucientes1,3, Sofia Hernández-Chan1, Rui Rosa2, Nick Humphries4, David W. Sims4,5,6 and Nuno Queiroz1,4* 1Centro de Investigação em Biodiversidade e Recursos Genéticos/Rede de Investigação em Biodiversidade e Biologia Evolutiva (CIBIO/InBIO), Universidade do Porto, Vairão, Portugal, 2Laboratório Marítimo da Guia, Marine and Environmental Sciences Centre (MARE), Faculdade de Ciências, Universidade de Lisboa, Cascais, Portugal, 3Instituto de Investigaciones Marinas, Consejo Superior de Investigaciones Científicas (IIM-CSIC), Vigo, Spain, 4Marine Biological Association of the United Kingdom, Plymouth, United Kingdom, 5Ocean and Earth Science, National Oceanography Centre Southampton, University of Southampton, Southampton, United Kingdom, 6Centre for Biological Sciences, University of Southampton, Southampton, United Kingdom In the pelagic environment diel vertical movements (DVM) are widespread across taxa, from zooplankton ascending from day-time depths into surface layers at night to avoid visual predators, to apex predators following prey movements to maximise foraging opportunities. The drivers of DVM in large predators such as pelagic sharks have only recently begun to be investigated in detail with the advent of sophisticated archival tags and high-resolution oceanographic datasets. In this study, we satellite tagged adult [>180 cm fork length, (FL)] blue sharks (Prionace glauca) in the North Atlantic Ocean to examine behavioural changes in response to the encountered environment, and therefore, to determine potential risks of capture using pelagic longline fisheries data. Although blue sharks recurrently use surface waters, cyclic diel behaviours were observed, with >95% of night-time spent above 250 m depth and variable day-time depth use. Hence, three different diel behaviours were identified during the tracking period: (i) regular normal DVM (nDVM) (dawn descent – dusk ascent, with over 90% of nighttime spent above 250 m, and between 5 and 50% of the day below this threshold); (ii) surface-oriented behaviour (occupation of surface waters both day and night), and (iii) deep depth-oriented nDVM [dawn descent – dusk ascent, with the majority (>50%) of daytime spent at depth]. Importantly, diel behaviours generally occurred in different ocean regions with nDVM frequently observed in high latitudes, associated with cold, highly productive waters (e.g., North Atlantic Current/Labrador Current convergence zone, West African upwelling area), while depth-oriented nDVM was observed in warm, oligotrophic areas. Thus, day-time occupation of shallow waters significantly increased with lower water temperature at depth (100 m), and with increasing concentration (and decreasing depth) of the chlorophyll a maximum. During nights of full moon blue sharks spent significantly more time in the depth range of longline hooks, while fishing effort and catches were also higher. We demonstrate that increased occupancy of Frontiers in Marine Science | www.frontiersin.org 1July 2021 | Volume 8 | Article 688076 fmars-08-688076 July 5, 2021 Time: 19:23 # 2 Vedor et al. Vertical Movement Patterns of Blue Sharks surface layers driven by highly productive, cold waters and greater lunar illumination lead to higher capture risk. Understanding habitat-specific vulnerability to fishing in a commercially important pelagic shark species is essential for improving management and conservation measures. Keywords: diel vertical movement, longline fisheries, fishing risk, environment, lunar phase INTRODUCTION Apex pelagic predators, given their extensive movements and wide distributions, play an important role shaping the ecological structure and habitat use of oceanic communities (Sims, 2003; Campana, 2016;Boerder et al., 2019). Yet, over the past decades, populations of commercially important species have been declining due to overexploitation (Baum et al., 2003; Pacoureau et al., 2021). This has been partially fuelled by a lack of knowledge on the movement and habitat use patterns, and where and when fish interact with fishing vessels, which in turn hinder stock assessment models (Mejuto and García-Cortés, 2005;Queiroz et al., 2019). While original stock assessments assumed fish populations were evenly distributed at an oceanbasin scale, recent approaches that include a spatial structure have shown a largely improved estimation performance (Punt, 2019). However, most current stock assessments still account for catchability as a measure of spatial distribution (Punt, 2019), raising the need for movement-based stock assessments as opposed to fisheries-based (Baum et al., 2003;Dunn et al., 2019). Advances in satellite telemetry and the increasing availability of remotely sensed environmental data have, however, allowed the coupling of recorded behavioural data with aspects of the encountered environment (Whitford and Klimley, 2019; Francisco et al., 2020). This has permitted studies that describe habitat preferences of migratory marine vertebrates even in remote areas, such as the open ocean (Humphries et al., 2010;Guzman et al., 2018;Sequeira et al., 2018). In addition, understanding how environmental changes affect the movements and behaviour of pelagic fish has also enabled the identification of areas and periods of higher vulnerability to fisheries (Song et al., 2009;Queiroz et al., 2016, 2019;Hays et al., 2019). Taken together, this spatial and behavioural information is therefore essential to the development of more effective conservation measures (Hammerschlag et al., 2016;Sequeira et al., 2019; Hindell et al., 2020). Besides wide-ranging distributions, pelagic predators such as sharks and tunas (Thunnus spp.) are also efficient divers, exploring a multitude of vertical habitats from the surface to the bathyal (Schaefer and Fuller, 2002;Wilson et al., 2005; Howey-Jordan et al., 2013). A particularly distinctive cyclic behaviour has been described for a wide range of pelagic species, from zooplankton to apex predators – diel vertical movements (DVM) (e.g., Weng et al., 2009;Coffey et al., 2017; Hafker et al., 2017). DVM in zooplankton is characterised by occupation of surface waters at nighttime and, as a phototaxic response, diving into deeper waters during daytime where light levels are lower (Lampert, 1989). These cyclic movements have been primarily associated with zooplankton avoidance of visual predators, which are then mimicked throughout the food chain to optimise prey encounter success (Hays, 2003). However, drivers of DVM for pelagic predators have generally been associated with thermoregulation and/or foraging (Pade et al., 2009;Campana et al., 2011;Queiroz et al., 2012;Andrzejaczek et al., 2019). For example, Atlantic bluefin tunas in the Gulf Stream reduced time spent in shallow warm waters and exhibited deeper DVMs possibly as a behavioural mechanism to maintain optimal body temperature (Teo et al., 2006). However, DVM in ectothermic fish may also represent a behavioural mechanism that balances the higher energy costs associated with night-time foraging activity in shallow, warm waters, with reduced activity in deep cold waters during daytime that lowers metabolic rate (Sims et al., 2006). DVM associated with predator avoidance is particularly important for species at lower trophic levels, where day-time occupation of aphotic layers and night-time shoaling in the darkness decreases the predation success of visual predators, consequently increasing the chances of prey survival (Wirsing et al., 2011). However, for pelagic top predators, DVM has been mostly associated with foraging behaviour (Carey et al., 1990;Andrzejaczek et al., 2019). By following prey diel rhythms, predators increase the spatial and temporal overlap with prey distribution, increasing the success of prey encounter (Sims et al., 2005). Several species of marine megafauna are diel vertical migrators, and some of the most important components of their diet include mesopelagic fish and diel vertically migrating species, supporting the notion that DVM increases foraging success (Sims et al., 2005;Pade et al., 2009;Duffy et al., 2017;Braun et al., 2019a). Shifts in diel behavioural patterns have been linked to variations in the thermal profile of the water column and to the distribution and timing of peak productivity. For instance, both blue and salmon sharks have shown an increased surface occupancy in highly productive frontal regions, associated with cold surface temperatures and sharp water-column stratification (Queiroz et al., 2012;Coffey et al., 2017). Besides diel differences in vertical space-use, activity rates also follow a diel cycle, being generally higher at night than during the day (Andrews et al., 2009;Papastamatiou et al., 2018). For example, increased nighttime activity of sharks was linked to high food abundance (Sims et al., 2006). Hence, this increased nocturnal activity has been associated with foraging on prey aggregations occurring in warm surface waters during this period (Sims et al., 2006; Papastamatiou et al., 2018). Although changes in vertical behaviour are strongly linked to the physical and biological environment, they can also vary in response to the illuminated lunar fraction. Distributions of plankton on nights of full moon are generally deeper than on nights of new moon which manifests as a phototaxic response for predator avoidance (Tarling, 1999;Last et al., 2016). This Frontiers in Marine Science | www.frontiersin.org 2July 2021 | Volume 8 | Article 688076 fmars-08-688076 July 5, 2021 Time: 19:23 # 3 Vedor et al. Vertical Movement Patterns of Blue Sharks behaviour is driven by a similar light-avoidance mechanism to the one observed for DVM, but here deepening occurs at night in response to moonlight. This behaviour cascades through trophic levels, from mesopelagic fish (Olivar et al., 2017) up to top predators (e.g., tunas and swordfish Xiphias gladius) where deeper movements on nights of full moon have been hypothesised to improve foraging success within prey aggregated at deeper depths (Musyl et al., 2003;Dewar et al., 2011;Abascal et al., 2015). Furthermore, behavioural variations linked to the lunar phase affect the capture risk of pelagic predators by commercial fisheries. For example, catches of commercial species such as yellowfin and blackfin tuna, and blue sharks were higher in nights of full moon (Curran, 2014;Orbesen et al., 2017), while catch rates for swordfish were highest during both new and full moon phases (Poisson et al., 2010;Lerner et al., 2013;Orbesen et al., 2017). The vertical behaviour of blue sharks has been the subject of a large number of studies (e.g., Carey et al., 1990;Campana et al., 2011;Queiroz et al., 2012; Braun et al., 2019b), but knowledge of the fine-scale behaviour of oceanic adult blue sharks and their vulnerability to fishing as a consequence is scarce. Blue sharks are one of the most exploited species by pelagic longline fisheries, with populations declining by ∼40% since the 1970s (Baum et al., 2003;Pacoureau et al., 2021) largely as a result of high mortality risk from fishing (Queiroz et al., 2016, 2019). Despite recently established fishing quotas for blue sharks in some regions, the uncertainty in current stock assessments make understanding shark finescale behaviour a priority for effective conservation management (Hammerschlag et al., 2016;Robinson et al., 2017;Boerder et al., 2019). To address these knowledge gaps, we used pop-off satellitelinked archival transmitters (PSATs) attached to adult blue sharks in the North Atlantic Ocean, to investigate potential coupling between open-ocean shark DVM with environmental fields and to identify both the spatial and temporal risks posed by longline fishing in the region. MATERIALS AND METHODS Shark Tagging A total of 22 adult blue sharks (Prionace glauca) were tagged in two general oceanic locations (in the mid-Atlantic and the northwestern Atlantic regions) between June 2010 and August 2011. Sharks were captured on commercial baited surface longlines and brought alongside the vessel in the gear-hauling phase and tagged. Pop-off satellite-linked archival transmitter tags (PSATs; Mk-10 model, Wildlife Computers, Redmond, WA, United States) were rigged with a monofilament tether covered with silicone tubing and looped through a small hole made in the base of the first dorsal fin. Depth, external temperature, and light-level parameters were archived at 1 s intervals and stored as summary data over set intervals of 6 h (00:00, 06:00, 12:00, and 18:00). For each period, time-at-depth histograms (TAD; aggregated in eight depth bins, 50, 100, 150, 200, 250, 400, 600, >600 m), minimum and maximum depth and temperature, as well as profiles of water temperature at depth were available. All shark tagging procedures undertaken in this study were approved by institutional ethical review committees and completed by licenced, trained, and experienced personnel. Track Processing The movement of PSAT-tagged sharks was estimated using either satellite relayed data from each tag or from archival data after the tags were physically recovered. Positions of each shark between attachment and tag pop-up were reconstructed using software provided by the manufacturer (WC-GPE, global position estimator programme suite), where daily maximal rate-of-change in light intensity was used to estimate local time of midnight or midday for longitude calculations, and day-length estimation for determining latitude. Anomalous longitude estimates resulting from dive-induced shifts in the estimated timings of dawn and dusk from light curves were automatically discarded from the dataset using software provided by the manufacturer (WC-GPE); latitude estimates were subsequently iterated for the previously obtained longitudes. An integrated state-space model [unscented Kalman filter – UKFSST (Lam et al., 2008); using spatially complete NOAA Optimum Interpolation Quarter Degree Daily SST Analysis data] was then applied to correct the raw geolocation estimates and obtain the most probable track. A regular time-series of locations was then estimated using a continuous-time correlated random walk Kalman filter, CTCRW (Jonsen et al., 2005) performed in R [crawl package (Johnson et al., 2008)]. Subsequently, the CTCRW state-space model was applied to each individual track, producing a single position estimate per day. Argos positions were parameterised with the K error model parameters for longitude and latitude implemented in the crawl package (Jonsen et al., 2005). To obtain unbiased estimates of shark space use, gaps between consecutive dates in the raw tracking data were interpolated to one position per day. However, any tracks with gaps exceeding 20 days were split into segments prior to interpolation, thus avoiding the inclusion of unrepresentative location estimates (Queiroz et al., 2016, 2019). Diel Diving Behaviour To detect cyclical patterns in the behaviour of the sharks, maximum depths for each 6 h interval were analysed with a Lomb-Scargle (LS) periodogram (Lomb, 1976;Scargle, 1982), using a detection range between 10 and 30 h in R package lomb (Ruf, 2010;Azzurro et al., 2012). The maximum percentage of total data variance fitted by the corresponding periodicity was chosen as the peak value (Campbell et al., 2010; Azzurro et al., 2012). To analyse diel differences in the behaviour of blue sharks, the time between 12:00 GMT and 18:00 GMT (6 h period) was considered daytime, while night-time was considered between 00:00 GMT and 06:00 GMT (6 h period). Time bins which encompassed sunset and sunrise times (between 06:00– 12:00 GMT and 18:00–00:00 GMT, respectively; Supplementary Figure 1) were excluded from the analyses (Dewar et al., 2011; Abecassis et al., 2012). Following previous studies on the diel behaviour of blue sharks, three classes were defined based on individual TAD and maximum depth for each daytime and night-time period over 5-day periods (Campana et al., 2011; Queiroz et al., 2012). A 5-day period was classified as (i) Frontiers in Marine Science | www.frontiersin.org 3July 2021 | Volume 8 | Article 688076 fmars-08-688076 July 5, 2021 Time: 19:23 # 4 Vedor et al. Vertical Movement Patterns of Blue Sharks depth-oriented nDVM if a shark spent more than 90% of the night-time above 250 m and 50% of the daytime below this depth; (ii) regular nDVM was characterised by over 90% of the night above 250 m, and between 5 and 50% of the day below the threshold. Finally, (iii) surface-oriented behaviour was characterised by over 90% of both day and night-time above 250 m. This depth threshold was used because it corresponds to the depth above which >95% of night-time occupancy for tracked blue sharks occurs (Supplementary Figure 1). Lastly, periods without clear day and night TAD patterns were classified as “other” and excluded from further analyses (corresponding to 4.0% of the total analysed periods). Subsequently, timeweighted depth and temperature were calculated for each 6 h dayand night-time bins using the function weighted.mean in R software, using the middle point of the bins defined for TAD data and the maximum depth/temperature. Pairwise Wilcoxon Rank Sum tests were used to statistically compare day and night-time weighted values, within each diel class. To further explore the behavioural thermoregulation hypothesis, a Pearson correlation was used to test the relationship between timeweighted temperatures during the day and the consecutive night. Moreover, individual dives were identified using the diveMove R package for high-resolution archival dive data available for three sharks (S3, S11, and S15), following Queiroz et al. (2017). For each dive, we assessed its duration and respective timeweighted temperatures. A Spearman correlation was used to compare the average time-weighted temperatures between dayand night-time dives. Monthly mean modelled environmental data (0.25◦×0.25◦) were acquired from Copernicus Marine Environment Monitoring Service (CMEMS) global ocean biogeochemistry non-assimilative hindcast (PISCES 1998–2011) for chlorophyll a(mg m−3) and global ocean physics reanalysis (GLORYS2V3 1993–2012) for temperature (◦C) products for the North Atlantic Ocean. Environmental data was extracted from the surface to 1,750 m depth, and to account for the spatial error around real individual geolocations, data was averaged for 1.25◦in latitude and 0.75◦in longitude (using a 5 ×3 grid cell) around each position. The relationship between day-time occupation of shallow waters (TAD above 250 m depth) and environmental variables, (i) sea surface temperature (SST), (ii) temperature at depth 100 m, (iii) maximum concentration of chlorophyll ain the water column, and (iv) depth of the chlorophyll amaximum, was investigated using a general linear mixed model (GLMM) with binomial distribution and cauchit link function. Data exploration techniques were used to identify potential outliers and assess collinearity among independent variables. Hence, chlorophyll aconcentrations above 1 mg m−3were considered outliers and discarded from further analysis (corresponding to 2.56% of data); in addition, the depth of the chlorophyll amaximum showed a high collinearity with concentration – variance inflation factor (VIF) higher than three (Zuur et al., 2009) – and was removed from the model selection process. Individual sharks were considered an independent sampling unit and were included as random effects; however, individuals with less than 10 data points were discarded from the analyses [resulting in the exclusion of S12; (Zuur et al., 2009)]. General mixed models were constructed by backward selection of individual terms to allow for testing of biologically meaningful interactions. The selected terms were: (i) sea surface temperature (SST), (ii) temperature at depth 100 m, and (iii) maximum concentration of chlorophyll ain the water column. The model was trained with 75% of data (n= 456) and the optimal selection was based on Akaike Information Criterion (AIC). At each stage of the selection process, fitted models were compared to the null model. Consequently, the fitted model with the lowest AIC was chosen as an optimal structure, with which a GLMM was built using MASS R package. Normal quantile–quantile plots of deviance residuals were assessed for normality of residuals and fit, while homoscedasticity, model misspecification, and residual spatial autocorrelation were evaluated by inspecting plots of response residuals against fitted values and candidate explanatory variables, respectively. Spatial and temporal residual autocorrelation was further assessed by including respective covariate structures and comparing model performance. This resulted in the inclusion of an autocorrelation structure of order 1corAR1 (nmle R package) in the model. Lunar Phase and Hook Encounter Risk Lunar phase (visible fraction of the disc, a continuous variable where 0 is dark moon and 1 corresponds to full moon) was extracted for each estimated shark geolocation. To account for the spatial error around individual geolocations, visible fraction of the disc data was averaged for 1.25◦in latitude and 0.75◦in longitude around each shark position. Data was obtained using the function moonAngle in the oce R package. To analyse shark night-time depth occupation in relation to the moon phase, the illuminated lunar fraction was divided into four 0.25 bins representative of the four lunar phases: (i) new moon, (ii) first and (iii) third quarter, and (iv) full moon. Drifting pelagic longline fishing effort data at 0.01◦×0.01◦ grid resolution were obtained from Global Fishing Watch (GFW) and subsequently gridded using the same spatial resolution as the environmental data (0.25 ×0.25◦grid cells) and between the years 2012 to 2016. The number of fishing days in each month within each grid cell was summed and then averaged across years (Queiroz et al., 2019). Geo-referenced catch data for blue shark were also obtained from logbooks of Spanish commercial longliners between 2013 and 2017. Similarly, monthly catch per unit of effort (CPUE) was summed within each 0.25 ×0.25◦ grid cells and averaged across years. GFW and CPUE analyses also considered the spatial error around individual geolocations. Despite the lack of fishing data for the same period sharks were tracked (2010–2011), fishing patterns and the associated capture risk are consistent across years (Kroodsma et al., 2018; Queiroz et al., 2019). Hence, by averaging fishing effort and CPUE across multiple years, interannual variation was considered when calculating the mean vulnerability of tracked sharks. Since longlines are deployed at night at depths between 100 and 300 m (Domingo et al., 2016), the TAD between 100 and 250 m was also calculated (TADhook) for each shark location. This percentage of time at hook depth was initially compared between nights of new and full moon; GFW fishing effort and CPUE were also analysed Frontiers in Marine Science | www.frontiersin.org 4July 2021 | Volume 8 | Article 688076 fmars-08-688076 July 5, 2021 Time: 19:23 # 5 Vedor et al. Vertical Movement Patterns of Blue Sharks FIGURE 1 | Individual geolocations for 15 blue sharks tracked in the North Atlantic Ocean between 2010 and 2011. in relation to the lunar phase using an independent-sample Wilcoxon rank sum test. RESULTS Between 2010 and 2011, a total of 15 blue sharks were successfully tracked in the North Atlantic for a total of 1,325 cumulative days (Figure 1 and Table 1). Sharks tagged in the mid-Atlantic either remained in the same general area for the tracking duration (between 90 and 120 days; S1, S4, S8, and S14) or moved southward (S2, S3, S7, and S15; these sharks were tracked between 78 and 120 days). Of these, S2, S7, and S15 moved south/southeast into the Cape Verde islands area, while S3 moved southwest into oligotrophic waters (Figure 1). The spatial distribution pattern for sharks tagged in the West Atlantic was more complex. Of the seven individuals tagged, four blue sharks (S5, S10, S11, and S12) remained in the overall area associated with the Gulf Stream (the North Atlantic Current/Labrador Current convergence zone) generally moving north along the edges of the frontal area, for periods ranging from seven to 89 days (Figure 1). Two sharks (S6 and S9) moved south into warmer waters (for 56 and 120 days, respectively) and one shark (S13) moved west into the shelf area off the American east coast, south of Nova Scotia. This shark was tracked for 180 days and was the only to have moved into cold surface waters (Supplementary Figure 2). Diel Diving Behaviour Periodograms of individual maximum dive depth for each 6 h bin showed a significant cyclic behaviour for 12 (out of 15) sharks, of which 11 individuals peaked significantly within the 23–25 h range, which was considered to represent a diel rhythm (Shepard et al., 2006;Campbell et al., 2010). The periodogram of S3 peaked significantly at 12 h, while the analysis was non-significant for sharks S8, S11, and S12 (Table 1). Regular nDVM was the most observed behaviour (55.7% of the time; Figure 2A), followed by surface-oriented (22.8% of the time; Figure 2C) and finally, depth-oriented nDVM (17.5% of the time; Figure 2B). Overall, the time-weighted depth of blue sharks was significantly deeper during daytime and shallower during nighttime (pairwise Wilcoxon signed rank test, V= 27243, α= 0.05, p<0.001, n= 236), with an average day-time depth of 175 ±99 m and an average night-time depth of 59 ±46 m. The same pattern was observed for each diel behaviour; while performing regular nDVM, the average day-time depth of tracked blue sharks was 146 ±95 m and average night-time was 45 ±38 m (median: 147 m and 33 m, respectively), whereas in depth-oriented nDVM behaviour, average day-time depth was 246 ±74 m and night-time 92 ±42 m (median: 244 m and 95 m, respectively). For surface-oriented behaviour, a shallower average depth was observed, with an average depth of 101 ±49 m during the day and 19 ±12 m during the night (median: 116 m and 21 m, respectively). Hence, the median day and night-time depth for each diel behaviour was different (Kruskal– Wallis rank sum test, daytime: chi-squared = 77.276, α= 0.05, p<0.001, n= 236; night-time: chi-squared = 73.78, α= 0.05, p<0.001, n= 236). Moreover, during the day, time-weighted water temperature was significantly colder (19.34 ±3.16◦C) than at night (23.26 ±2.99◦C; pairwise Wilcoxon signed rank test, V= 29306, α= 0.05, p<0.001, n= 259). This general pattern was consistent across diel behaviours; during nDVM, average day-time temperature was 20.83 ±3.10◦C and average night-time temperature was 23.39 ±3.33◦C (paired Wilcoxon signed rank exact test, V= 3125, α= 0.05, p<0.001, n= 88); while sharks were performing depth-oriented nDVM, the average temperature was 17.53 ±2.66◦C and 23.02 ±2.50◦C for dayand night-time, respectively (paired Wilcoxon signed rank exact test, V= 2065, α= 0.05, p<0.001, n= 64). Finally, for surface-oriented behaviour, average day-time temperature was 18.85 ±2.42◦C and 23.28 ±3.11◦C at night (paired Wilcoxon signed rank exact test, V= 611, α= 0.05, p<0.001, n= 35). However, time weighted temperatures during the day were positively correlated with the time-weighted temperatures experienced in consecutive nights (Pearson correlation = 0.49, t= 8.66, α= 0.05, p<0.001, n= 244). Frontiers in Marine Science | www.frontiersin.org 5July 2021 | Volume 8 | Article 688076 fmars-08-688076 July 5, 2021 Time: 19:23 # 6 Vedor et al. Vertical Movement Patterns of Blue Sharks TABLE 1 | Summary data for satellite-tagged blue sharks. Shark ID Fork length (cm) Sex Tag type Location tagged Tagging date Days-at-liberty Pop-up date %DVM Periodogram Shark 1 260 M PAT-Mk10 Mid-Atlantic 21/08/2011 90 30/12/2011 93.75 24 Shark 2 250 M PAT-Mk10 Mid-Atlantic 22/08/2011 120 20/12/2011 94.74 24.02 Shark 3* 240 M PAT-Mk10 Mid-Atlantic 26/08/2011 78 24/12/2011 68.75 12 Shark 4 200 M PAT-Mk10 Mid-Atlantic 26/08/2011 120 24/12/2011 95.83 24.02 Shark 5 185 M PAT-Mk10 West Atlantic 26/06/2010 89 23/09/2010 100 23.99 Shark 6 192 F PAT-Mk10 West Atlantic 27/06/2010 56 25/09/2010 100 23.94 Shark 7 240 F PAT-Mk10 Mid-Atlantic 21/08/2011 90 19/11/2011 100 24.02 Shark 8 240 F PAT-Mk10 Mid-Atlantic 26/08/2011 120 25/12/2011 78.57 ns Shark 9 260 F PAT-Mk10 West Atlantic 29/06/2010 120 27/10/2010 70.83 24.01 Shark 10 240 F PAT-Mk10 West Atlantic 30/06/2010 27 27/12/2010 100 24.09 Shark 11* 200 M PAT-Mk10 West Atlantic 25/06/2010 14 23/10/2010 100 ns Shark 12 210 M PAT-Mk10 West Atlantic 25/06/2010 7 23/10/2010 100 ns Shark 13 235 M PAT-Mk10 West Atlantic 25/06/2010 180 22/12/2010 58.82 23.99 Shark 14 220 F PAT-Mk10 Mid-Atlantic 27/08/2011 120 25/12/2011 95.45 24.06 Shark 15* 220 F PAT-Mk10 Mid-Atlantic 28/08/2011 94 30/11/2011 100 23.97 DNR 1 165 F PAT-Mk10 West Atlantic 25/06/2010 120 23/10/2010 – – DNR 2 170 F PAT-Mk10 West Atlantic 30/06/2010 Failed – – – DNR 3 185 M PAT-Mk10 West Atlantic 25/06/2010 Failed – – – DNR 4 230 F PAT-Mk10 Mid-Atlantic 28/08/2011 180 24/02/2012 – – DNR 5 220 F PAT-Mk10 Mid-Atlantic 29/08/2011 180 25/02/2012 – – DNR 6 240 M PAT-Mk10 Mid-Atlantic 27/08/2011 180 23/02/2012 – – DNR 7 245 M PAT-Mk10 Mid-Atlantic 27/08/2011 180 23/02/2012 – – F, female; M, male; Ns, non-significant periodogram. Starred individuals indicate archival data was retrieved from the tag. Similarly, for a subset of archival data, average time-weighted temperatures during day-time dives were positively correlated with night-time dives (Spearman correlation = 0.18, S = 848744, α= 0.05, p<0.05, n= 184). Environmental-Linked Behavioural Changes The vertical behaviour of the sharks changed throughout the tracking period in relation to both encountered temperature and chlorophyll aconcentration in the water column, with regular nDVM (and surface-oriented) behaviour generally associated with colder, more productive waters (Figure 3). In contrast, depth-oriented diel behaviour was commonly observed in regions with less productive, warmer surface waters, but, strikingly, with deep chlorophyll amaxima layers (Figure 3). Similar patterns were evident from high-resolution archival data (Figure 4). Regular nDVM was generally characterised by shallower day-time dives, colder surface waters (23.55 ±3.76◦C) and higher chlorophyll aconcentrations near the surface (Figures 4A,C), while depth-oriented nDVM was performed in regions of warmer (27.35 ±0.80◦C), oligotrophic surface waters, and lower, deeper maximum concentrations of chlorophyll a (Figures 4B,D). Differences in habitat use were also clear between the two diel behaviours over 24 h periods (Figure 5). During both normal diel patterns, sharks showed greater nighttime occupation of waters above ca. 25◦C (Figures 5A,B), however, during nighttime, a consistent occupation in lower temperatures (<17◦C) was evident during depth-oriented nDVM (Figure 5D). Similarly, blue sharks maximised time at medium to high (>4.5 ×10−2mg m−3) chlorophyll a concentrations during nighttime (Figures 5C,D), although time spent at these concentrations was greater during regular nDVM (Figure 5C). Importantly, in depth-oriented nDVM, a peak in the amount of time spent in layers of low chlorophyll a concentrations (ca.<1.4 ×10−2mg m−3) was also observed during daytime (Figure 5D). The GLMM analysis revealed that day-time occupancy in shallow water (i.e., TAD above 250 m) significantly increased with decreasing (i) water temperature at depth (100 m), and (ii) with increasing concentration of chlorophyll a(Table 2, Figure 6, and Supplementary Figure 3). The depth of the maximum concentration of chlorophyll awas inversely correlated with the concentration (Spearman’s rank correlation: rho = −0.92, S= 30395248, p-value <0.001, n= 456); hence, the observed increased shark day-time occupancy in shallow water was also linked to the shoaling depth of the chlorophyll amaximum. Therefore, regular nDVM was predominantly observed associated with the Gulf Stream and the North Atlantic Current-Labrador Current convergence zone (NAC-LCCZ), but also in mid-Atlantic regions and in the proximity of archipelagos, such as the Azores and Cape Verde (Figure 7A). Surface-oriented behaviour was almost exclusively observed off the American east coast, in the vicinity of the Gulf Stream and in the midAtlantic (Supplementary Figure 4), whereas depth-oriented nDVM was generally observed during southward movements in mid-Atlantic regions (Figure 7B). Interestingly, almost no spatial overlap was observed between the two nDVM behaviours. Frontiers in Marine Science | www.frontiersin.org 6July 2021 | Volume 8 | Article 688076 fmars-08-688076 July 5, 2021 Time: 19:23 # 7 Vedor et al. Vertical Movement Patterns of Blue Sharks FIGURE 2 | Diel behaviour of blue sharks. Left panel shows the percentage of time-at-depth for day- (white) and night-time (grey) for (A) nDVM (day-time occupancy above 250 m: 66.0%), (B) depth-oriented nDVM (32.7%), and (C) surface-oriented (100%) diel diving behaviours; Right panels show the correspondent 5-day archival depth use. Horizontal red lines denote the depth threshold considered to define each diel class (250 m depth). This resulted in a latitudinal pattern, with surface-oriented and regular nDVM commonly observed in high latitudes of colder, productive waters (Figure 8). Regular nDVM was also frequent in low latitudes but only associated with cold, productive waters of frontal regions (in the western Africa upwelling region) (Figures 7A,8). The depth-oriented nDVM more frequently observed in lower latitudes of warmer, deep chlorophyll a maxima and oligotrophic waters (Figure 8). Lunar Phase and Overlap With Surface Longlining Hooks During nights of new moon, occupation of surface waters (less than 100 m; average: 81.9% ±22.4) was higher than in nights of full moon (70.6% ±25.5; Wilcoxon rank sum test, W= 15862, α= 0.05, p<0.001, n= 424). However, TADhook was significantly greater on nights of full moon when compared to nights of new moon (15.55% ±22.28 and 9.37% ±18.06, respectively; Wilcoxon rank sum test, W= 25475, α= 0.05, p<0.01, n= 424). In fact, along shark geolocations, both fishing effort (890 ±611 fishing days) as well as CPUE were also highest on nights of full moon (278 ±144 kg grid cell−1set−1) when compared to nights of new moon (687 ±443 days and 253 ±192 kg grid cell−1set−1, respectively; Wilcoxon rank sum test, fishing effort: W= 25475, α= 0.05, p<0.01, n= 284; CPUE: W= 5777, α= 0.05, p<0.01, n= 194). DISCUSSION Blue sharks performed extensive movements covering a large area of the North Atlantic, with results showing that diel Frontiers in Marine Science | www.frontiersin.org 7July 2021 | Volume 8 | Article 688076 fmars-08-688076 July 5, 2021 Time: 19:23 # 8 Vedor et al. Vertical Movement Patterns of Blue Sharks FIGURE 3 | Maximum daily dive depths for two individual blue sharks, shark #7 on the left (A,C) and #15 on the right (B,D). Dive series are overlaid on water temperature (A,B) and chlorophyll aconcentration (C,D). Coloured bar above represent the identified diel behaviour. behaviours shifted in response to environmental changes. Specifically, sharks (i) spent more time at depth in less productive, warmer, and deep chlorophyll amaximum regions, and (ii) performed behavioural patterns linked with surface occupancy in colder, more productive waters. The latter was also associated with a higher spatio-temporal overlap with hook depths and thus, higher fishing risk in these areas, particularly on nights of full moon. Horizontal Movements Tracked sharks performed typical, extensive movements in the open ocean of the North Atlantic and also into shelf waters. Such movements are generally thought to be part of long-distance, seasonal migrations linked with foraging and/or reproduction (Queiroz et al., 2005;Nakano and Stevens, 2008;Stevens, 2010). During summer, blue sharks typically migrate north following the warming of surface waters at higher latitudes (Queiroz et al., 2010), exemplified by known movements into productive waters of the northern Gulf Stream and North Atlantic Current (Queiroz et al., 2019). Movements into shelf waters also occur, which allows sharks to forage on shelf-associated seasonal schools of small pelagic fish and tuna (Henderson et al., 2001; Southward et al., 2004). Hence, the summer diet of blue sharks in northern latitudes has been described to shift from being typically comprised of cephalopods to a preference for teleosts (Stevens, 1973). Besides performing extensive horizontal movements, blue sharks we tracked also displayed high residency in specific regions, broadly associated with mesoscale oceanic features, such as thermal fronts and upwelling regions. For example, individuals tracked in the West Atlantic showed greater occupation of the Gulf Stream and the NAC-LCCZ, both frontal regions of high primary productivity and forage accumulation, where megafauna is known to aggregate (Campana et al., 2011;Scales et al., 2014;Braun et al., 2019b). Similarly, sharks tagged in the mid-Atlantic, near the Azores archipelago showed a longer residence within this area, which has been previously described as a preferred wintering ground for this species (Vandeperre et al., 2014). Three sharks displayed southward movements during the winter into the tropical Atlantic likely following warm SST isotherms toward southern waters (Queiroz et al., 2012). Female blue sharks S7 and S14 moved into the Eastern Tropical Atlantic (ETA), a highly productive area associated with the strong West African Frontiers in Marine Science | www.frontiersin.org 8July 2021 | Volume 8 | Article 688076 fmars-08-688076 July 5, 2021 Time: 19:23 # 9 Vedor et al. Vertical Movement Patterns of Blue Sharks FIGURE 4 | Archival 5-day period depth series for two individual blue sharks, during nDVM for shark #15 on the left (A,C) and depth-oriented nDVM for shark #3 on the right (B,D). Shark depth is overlaid with water temperature (A,B) and chlorophyll aconcentration (C,D). Coloured bar above represent the identified diel behaviour. upwelling regime (Karstensen et al., 2008;Stramma et al., 2010; Gilly et al., 2013). Diel Vertical Movements The majority of tracked sharks in this study displayed a cyclical behaviour, with 11 out of 15 individuals showing a diel rhythm [peak within 23–25 h range; (Campbell et al., 2010)]. The vertical movements for one individual were cyclic at 12 h (semi-diel rhythm; S3) indicating circatidal movements (Shepard et al., 2006). However, no diel behaviour was detected (i.e., the periodogram was non-significant) for three individuals. These were either tracked for short periods of time (seven to 14 days, for S12 and S11, respectively) or were frequently recorded diving to extreme depths (S8). In this last case, the irregularity of deep dives along the tracking period could have reduced the chance of detecting cyclical patterns in the maximum depth, and thus a significant diel pattern. Nonetheless, using Lomb-Scargle periodograms we were able to identify cyclical rhythms for the majority of the tracked sharks, likely because it is robust to the presence of non-regular and zero-inflated data, both of which are usually considered to be typical disadvantages of satellite relayed telemetry data (Ruf, 2010;Schaffeld et al., 2016). Interestingly, this study observed for the first time, to our knowledge, semi-diel movements for pelagic sharks in the open ocean. Instead, these have been often described for sharks associated with coastal and shelf areas, where the cyclical tidal transport of zooplankton influences the behaviour of predators in higher trophic levels (Shepard et al., 2006;RodríguezCabello et al., 2016). Nonetheless, shark S3 was tracked moving exclusively in the open ocean, with frequent excursions into the bathyal (maximum dive depth of 1,704 m), possibly associated with deep-water foraging on organisms influenced by semidiel patterns in deep-water currents. Indeed, semi-diel vertical migrations have also been observed in open-ocean deep-sea communities at over 1,000 m depth (Aguzzi et al., 2010, 2018), where tidal effects propagate throughout the water column from changes in strength and direction of deep-water currents (Uiblein et al., 2002;Trenkel et al., 2004;Lorance and Trenkel, 2006;Aguzzi et al., 2010) and by variations in temperature and salinity (Ratsimandresy et al., 2017). Furthermore, bathymetric features, such as deep-sea banks or ridges, create a displacement in the water column structure and on deep-water currents (Genin, 2004;Cotté and Simard, 2005), potentially favouring plankton advection from the bathyal to shallower depths, thereby acting as an upwelling area highly attractive for pelagic Frontiers in Marine Science | www.frontiersin.org 9July 2021 | Volume 8 | Article 688076 fmars-08-688076 July 5, 2021 Time: 19:23 # 16 Vedor et al. Vertical Movement Patterns of Blue Sharks marine megafauna conservation. Front. Mar. Sci. 6:639. doi: 10.3389/fmars. 2019.00639 Sequeira, A. M. M., Rodriguez, J. P., Eguiluz, V. M., Harcourt, R., Hindell, M., Sims, D. W., et al. (2018). Convergence of marine megafauna movement patterns in coastal and open oceans. Proc. Natl. Acad. Sci. U.S.A. 115, 3072–3077. doi: 10.1073/pnas.1716137115 Shepard, E. L. C., Ahmed, M. Z., Southall, E. J., Witt, M. J., Metcalfe, J. D., and Sims, D. W. (2006). Diel and tidal rhythms in diving behaviour of pelagic sharks identified by signal processing of archival tagging data. Mar. Ecol. Prog. Ser. 328, 205–213. doi: 10.3354/meps328205 Sims, D. W. (2003). Tractable models for testing theories about natural strategies: foraging behaviour and habitat selection of free-ranging sharks. J. Fish Biol. 63, 53–73. doi: 10.1111/j.1095-8649.2003.00207.x Sims, D. W., Southall, E. J., Tarling, G. A., and Metcalfe, J. D. (2005). Habitat-specific normal and reverse diel vertical migration in the planktonfeeding basking shark. J. Anim. Ecol. 74, 755–761. doi: 10.1111/j.13652656.2005.00971.x Sims, D. W., Wearmouth, V. J., Southall, E. J., Hill, J. M., Moore, P., Rawlinson, K., et al. (2006). Hunt warm, rest cool: bioenergetic strategy underlying diel vertical migration of a benthic shark. J. Anim. Ecol. 75, 176–190. doi: 10.1111/j. 1365-2656.2005.01033.x Song, L., Zhou, J., Zhou, Y., Nishida, T., Jiang, W., and Wang, J. (2009). Environmental preferences of bigeye tuna, Thunnus obesus, in the Indian Ocean: an application to a longline fishery. Environ. Biol. Fish. 85, 153–171. doi: 10.1007/s10641-009-9474-7 Southward, A. J., Langmead, O., Hardman-Mountford, N. J., Aiken, J., Boalch, G. T., Dando, P. R., et al. (2004). Long-term biological and environmental researches in the western English Channel. Adv. Mar. Biol 47, 1–105. Stevens, J. D. (1973). Stomach contents of the blue shark (Prionace glauca L.) off south-west England. J. Mar. Biol. Assoc. U.K. 53, 357–361. Stevens, J. D. (2010). “Epipelagic oceanic elasmobranchs,” in Sharks and Their Relatives II: Biodiversity, Adaptive Physiology and Conservation, eds J. C. Carrier, J. A. Musick, and M. R. Heithaus (Boca Raton, FL: CRC Press), 3–35. Stramma, L., Schmidtko, S., Levin, L. A., and Johnson, G. C. (2010). Ocean oxygen minima expansions and their biological impacts. Deep Sea Res. I Oceanogr. Res. Pap. 57, 587–595. doi: 10.1016/j.dsr.2010.01.005 Tarling, G. (1999). The effect of lunar eclipse on the vertical migration behaviour of Meganyctiphanes norvegica (Crustacea: Euphausiacea) in the Ligurian Sea. J. Plankton Res. 21, 1475–1488. doi: 10.1093/plankt/21.8.1475 Teo, S. L. H., Boustany, A., Dewar, H., Stokesbury, M. J. W., Weng, K. C., Beemer, S., et al. (2006). Annual migrations, diving behavior, and thermal biology of Atlantic bluefin tuna, Thunnus thynnus, on their Gulf of Mexico breeding grounds. Mar. Biol. 151, 1–18. doi: 10.1007/s00227-006-0447-5 Thums, M., Meekan, M., Stevens, J., Wilson, S., and Polovina, J. (2013). Evidence for behavioural thermoregulation by the world’s largest fish. J. R. Soc. Interface 10:20120477. doi: 10.1098/rsif.2012.0477 Trenkel, V. M., Lorance, P., and Mahévas, S. (2004). Do visual transects provide true population density estimates for deepwater fish? ICES J. Mar. Sci. 61, 1050–1056. Uiblein, F., Lorance, P., and Latrouite, D. (2002). Variation in locomotion behaviour in northern cutthroat eel (Synaphobranchus kaupi) on the Bay of Biscay continental slope. Deep Sea Res. I Oceanogr. Res. Pap. 49, 1689–1703. doi: 10.1016/s0967-0637(02)00065-1 Vandeperre, F., Aires-da-Silva, A., Fontes, J., Santos, M., Serrao Santos, R., and Afonso, P. (2014). Movements of blue sharks (Prionace glauca) across their life history. PLoS One 9:e103538. doi: 10.1371/journal.pone.0103538 Weng, K. C., Stokesbury, M. J., Boustany, A. M., Seitz, A. C., Teo, S. L., Miller, S. K., et al. (2009). Habitat and behaviour of yellowfin tuna Thunnus albacares in the Gulf of Mexico determined using pop-up satellite archival tags. J. Fish Biol. 74, 1434–1449. doi: 10.1111/j.1095-8649.2009.02209.x Whitford, M., and Klimley, A. P. (2019). An overview of behavioral, physiological, and environmental sensors used in animal biotelemetry and biologging studies. Anim. Biotelem. 7:26. doi: 10.1186/s40317-019-0189-z Wilson, S. G., Polovina, J. J., Stewart, B. S., and Meekan, M. G. (2005). Movements of whale sharks (Rhincodon typus) tagged at Ningaloo Reef, Western Australia. Mar. Biol. 148, 1157–1166. doi: 10.1007/s00227-005-0153-8 Wirsing, A. J., Heithaus, M. R., and Dill, L. (2011). Predator-induced modifications to diving behavior vary with foraging mode. Oikos 120, 1005–1012. doi: 10. 1111/j.1600-0706.2010.18844.x Zhang, X., and Dam, H. G. (1997). Downward export of carbon by diel migrant mesozooplankton in the central equatorial Pacific. Deep Sea Res. II Top. Stud. Oceanogr. 44, 2191–2202. doi: 10.1016/s0967-0645(97)00060-x Zuur, A. F., Ieno, E. N., Walker, N., Saveliev, A. A., and Smith, G. M. (2009). Mixed Effects Models and Extensions in Ecology with R, Statistics for Biology and Health. New York, NY: Springer-Verlag. Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. Copyright © 2021 Vedor, Mucientes, Hernández-Chan, Rosa, Humphries, Sims and Queiroz. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms. Frontiers in Marine Science | www.frontiersin.org 16 July 2021 | Volume 8 | Article 688076