Horizontal and Vertical Movement Patterns and Habitat Use of Juvenile Porbeagles (Lamna nasus) in the Western North Atlantic

The porbeagle (Lamna nasus) is a large, highly migratory endothermic shark broadly distributed in the higher latitudes of the Atlantic, South Pacific, and Indian Oceans. In the North Atlantic, the porbeagle has a long history of fisheries exploitation and current assessments indicate that this stock is severely overfished. Although much is known of the life history of this species, there is little fisheries-independent information about habitat preferences and ecology. To examine migratory routes, vertical behavior, and environmental associations in the western North Atlantic, we deployed pop-up satellite archival transmitting tags on 20 porbeagles in late November, 2006. The sharks, ten males and ten females ranging from 128 to 154 cm fork length, were tagged and released from a commercial longline fishing vessel on the northwestern edge of Georges Bank, about 150 km east of Cape Cod, MA. The tags were programmed to release in March (n = 7), July (n = 7), and November (n = 6) of 2007, and 17 (85%) successfully reported. Based on known and derived geopositions, the porbeagles exhibited broad seasonally-dependent horizontal and vertical movements ranging from minimum linear distances of 937 to 3,310 km and from the surface to 1,300 m, respectively. All of the sharks remained in the western North Atlantic from the Gulf of Maine, the Scotian Shelf, on George's Bank, and in the deep, oceanic waters off the continental shelf along the edge of, and within, the Gulf Stream. In general, the population appears to be shelf-oriented during the summer and early fall with more expansive offshore radiation in the winter and spring. Although sharks moved through temperatures ranging from 2 to 26°C, the bulk of their time (97%) was spent in 6-20°C. In the summer months, most of the sharks were associated with the continental shelf moving between the surface and the bottom and remaining < 200 m deep. In the late fall and winter months, the porbeagles moved into pelagic habitat and exhibited two behavioral patterns linked with the thermal features of the Gulf Stream: “non-divers” (n = 7) largely remained at epipelagic depths and “divers” (n = 10) made frequent dives into and remained at mesopelagic depths (200–1000 m). These data demonstrate that juvenile porbeagles are physiologically capable of exploiting the cool temperate waters of the western North Atlantic as well as the mesopelagic depths of the Gulf Stream, possibly allowing exploitation of prey not available to other predators.

The porbeagle (Lamna nasus) is a large, highly migratory endothermic shark broadly distributed in the higher latitudes of the Atlantic, South Pacific, and Indian Oceans. In the North Atlantic, the porbeagle has a long history of fisheries exploitation and current assessments indicate that this stock is severely overfished. Although much is known of the life history of this species, there is little fisheries-independent information about habitat preferences and ecology. To examine migratory routes, vertical behavior, and environmental associations in the western North Atlantic, we deployed pop-up satellite archival transmitting tags on 20 porbeagles in late November, 2006. The sharks, ten males and ten females ranging from 128 to 154 cm fork length, were tagged and released from a commercial longline fishing vessel on the northwestern edge of Georges Bank, about 150 km east of Cape Cod, MA. The tags were programmed to release in March (n = 7), July (n = 7), and November (n = 6) of 2007, and 17 (85%) successfully reported. Based on known and derived geopositions, the porbeagles exhibited broad seasonally-dependent horizontal and vertical movements ranging from minimum linear distances of 937 to 3,310 km and from the surface to 1,300 m, respectively. All of the sharks remained in the western North Atlantic from the Gulf of Maine, the Scotian Shelf, on George's Bank, and in the deep, oceanic waters off the continental shelf along the edge of, and within, the Gulf Stream. In general, the population appears to be shelf-oriented during the summer and early fall with more expansive offshore radiation in the winter and spring. Although sharks moved through temperatures ranging from 2 to 26 • C, the bulk of their time (97%) was spent in 6-20 • C. In the summer months, most of the sharks were associated with the continental shelf moving between the surface and the bottom and remaining < 200 m deep. In the late fall and winter months, the porbeagles moved into pelagic habitat and exhibited two behavioral patterns linked with the thermal features of the Gulf Stream: "non-divers" (n = 7) largely remained at epipelagic depths and "divers" (n = 10) made frequent dives into and remained at

INTRODUCTION
The porbeagle, Lamna nasus, is a large, highly migratory pelagic shark broadly distributed in the temperate latitudes of the Atlantic, South Pacific, and South Indian Oceans (Compagno, 2002). Historically, the species has been subjected to both commercial and recreational fisheries throughout its range (Curtis et al., 2016;Haugen, 2020). In the North Atlantic, the porbeagle is currently managed as two stocks (i.e., east and west) and the most recent stock assessment indicates that both are heavily overfished (ICCAT, 2009;Curtis et al., 2016). As a result, the IUCN Red List of Threatened Species assessed the species as Vulnerable in 2018 (Rigby et al., 2018). The porbeagle is also currently listed on Annex I of the United Nations Convention on Law of the Sea, Appendix II of the Convention for the Conservation of Migratory Species (CMS, 2008), and, since 2014, Appendix II of the Convention on International Trade in Endangered Species of Wild Fauna and Flora (www.cites.org; Curtis et al., 2016). In 2015, the International Commission for the Conservation of Atlantic Tunas (ICCAT) passed a recommendation requiring contracting and cooperating parties to release unharmed porbeagles caught in association with ICCAT fisheries. In 2016, the porbeagle was proposed to be listed for protection under the U.S. Endangered Species Act, but an evaluation of all available information determined that the listing was not warranted due to a low risk of extinction in the foreseeable future (Curtis et al., 2016). Although major harvesters of porbeagle have either stopped (European Union, Canada) or capped landings (USA) in an effort to allow for the rebuilding of stocks (Campana et al., 2015;reviewed by Haugen, 2020), this species occurs in the highest bycatch risk zone in the North Atlantic (Queiroz et al., 2019) and mortality associated with bycatch may be high (Haugen, 2020). The most optimistic projections predict that stock recovery on both sides of the Atlantic will require decades (ICCAT, 2009). Given the current status of porbeagle stocks, the Standing Committee on Research and Statistics in ICCAT proffered a suite of research recommendations that included a better understanding of stock structure, vertical and migratory movements, and habitat associations with oceanographic features (ICCAT, 2009). Such information on the temporal and spatial distribution of porbeagles in relation to environmental/ecosystem features would not only facilitate a better understanding of catch trends, but also spatial management to reduce potential fishing mortality (ICCAT, 2009).
Although the life history of the porbeagle has been wellstudied Joyce et al., 2002;Natanson et al., 2019), much of what is known about stock structure, movements, habitat use, and spatial ecology has come from fisheries-dependent data (Campana and Joyce, 2004;Kohler and Turner, 2019). These studies (using conventional tags and commercial catch data) suggest that porbeagles exhibit seasonal shifts in abundance and are highly migratory in the western North Atlantic (WNA). However, relatively little is known about the three-dimensional movements and habitat use of these fish as well as the biotic and abiotic factors (e.g., ambient temperature, oceanic fronts, prey availability, and reproductive behaviors) that may influence their movements.
While fisheries-dependent data suggest that porbeagles do indeed have a wide thermal distribution (i.e., captured in waters with a sea surface temperature [SST] ranging from 2 to 23 • C), most catches near the continental shelves are associated with the thermocline and with oceanographic frontal zones in a relatively narrow temperature range between 5 and 10 • C (Campana and Joyce, 2004). Based on these data, Campana and Joyce (2004) suggested that the vertical distribution of the porbeagle in the WNA is relatively narrow (<200 m), although they did not rule out the possibility of vertical migrations to greater depths.
During the last two decades, the advent of electronic tagging technologies has allowed researchers to gain a more holistic, fisheries-independent understanding of the temporal and spatial movement patterns and stock structure of highly migratory fishes (reviewed by Costa et al., 2012). In addition, these new technologies allow for the correlation of movement patterns with oceanographic features (Braun et al., 2019) and, therefore, provide a more robust understanding of habitat use. The most commonly used technology to study the three-dimensional movement ecology in commercially important teleosts and elasmobranchs is the pop-up satellite archival transmitting (PSAT) tag (e.g., Block, 2005;Skomal et al., 2009Skomal et al., , 2017Campana et al., 2010;Galuardi et al., 2010;Musyl et al., 2011). To date, several studies have deployed PSAT tags to examine porbeagle movements in the Atlantic (Pade et al., 2009;Campana et al., 2010;Saunders et al., 2011;Biais et al., 2017). Work in the eastern North Atlantic (ENA) off the coast of the United Kingdom showed that porbeagles appear to exhibit deep (>300 m) and, in some cases, diel diving behavior and highly variable horizontal movements (Pade et al., 2009;Saunders et al., 2011). However, conclusions from those studies were based on relatively low sample sizes (total n = 7) and short track durations (22-122 days). In the WNA, the tracks of 21 porbeagles showed deep diving (up to 1,360 m) and broad horizontal movements as far as the Sargasso Sea (Campana et al., 2010). Although track durations spanned almost a year (60-348 days), the porbeagle sample in that study was dominated by large, sexually mature females.
Based on those data, Campana et al. (2010) hypothesized that the Sargasso Sea may be an important parturition area for porbeagles in the WNA. Most recently, Biais et al. (2017) deployed PSAT tags on nine subadult/adult porbeagles in the Bay of Biscay for up to a year, and showed broad migrations to the Arctic Circle, Madeira, and the mid-Atlantic Ridge, as well as site fidelity to the tagging location. However, the lack of juveniles from both sexes in these studies limits the extent to which these findings are applicable to all age classes.
The porbeagle and other members of the family Lamnidae are unique among sharks in having the capacity to retain metabolically produced heat to significantly elevate tissue (i.e., eyes, brain, viscera and swimming muscle) temperatures above ambient (i.e., regional endothermy, Carey and Teal, 1969a;Goldman, 1997;Bernal et al., 2001;Goldman et al., 2004;Sepulveda et al., 2004). Among lamnids, the porbeagle and its congener the salmon shark (Lamna ditropis) are able to maintain their eyes, brain, viscera, and swimming muscles up to 20 • C above ambient (Goldman et al., 2004;Bernal et al., 2005), resulting in a relatively stable elevated tissue operating temperature. This may allow lamnid sharks to make long duration or more frequent dives below the thermocline or to inhabit very cold, highly productive, subpolar waters (2-10 • C) for prolonged periods (i.e., numerous months) of time without compromising locomotor, sensory, and digestive function. In the porbeagle, regional endothermy may allow them to expand their thermal niche to forage in cooler, more productive waters at greater depths and/or higher latitudes (Block and Carey, 1985;Carey et al., 1985;Bernal et al., 2001;Dickson and Graham, 2004). While the vertical movement patterns of large, sexually mature females support a wide thermal (i.e., depth and latitude) distribution (Campana et al., 2010), the vertical movement patterns of juvenile porbeagles remain unknown.
Information on porbeagle movement ecology, habitat use, and stock structure as it relates to juveniles is warranted given relatively little is known about the three-dimensional movements of these fish. As this stock is highly overexploited, fishery managers need more empirical data for the development of enhanced conservation measures to facilitate the restoration of porbeagle populations in the Atlantic. In this study, we used PSAT technology to examine long-term, three-dimensional movements and to characterize habitat use by juvenile porbeagles in the WNA.

MATERIALS AND METHODS
All methods associated with the capture and handling of sharks were conducted in accordance to standards from the Institutional Animal Care and Use Committee of the University of Massachusetts, Dartmouth (Protocol #05-07).

Tagging
During this study, 20 porbeagles were captured and tagged in the western North Atlantic Ocean east of Cape Cod, MA near Cultivator Shoal, George's Bank (41.844 • N, 68.095 • W) in November, 2006 aboard the F/V Eagle Eye II. Sharks were caught on pelagic longline fishing gear comprised of monofilament gangions (9.1 m; 180 kg test) with 18/0 circle hooks (Mustad Corp., Miami, FL, USA) baited with whole Atlantic mackerel (Scomber scombrus). Two longline sets of 210 and 315 hooks were allowed to soak for 3 and 5 h, respectively, before hauling. During each set, a temperature profile of the water column was obtained using a conductivity-temperature-depth meter (model SBE911, Sea-Bird Electronics, Bellevue, WA, USA).
All hooked sharks were brought alongside the vessel and inspected for physical trauma (e.g., hook damage); only healthy, actively swimming, mouth-hooked specimens that showed little or no bleeding were selected for this study. Each shark was guided by the monofilament leader into a mesh-lined cradle and quickly lifted to the deck for sex determination, accurate length (over the body fork length, FL) measurement, and PSAT tag attachment (n = 20, Model Mk10-PAT, Wildlife Computers Inc., Redmond, WA, USA). Each PSAT tag was tethered to an intramuscular nylon umbrella dart, which was implanted into the musculature and through the cartilaginous supports of the dorsal fin (Kohler and Turner, 2001). All PSAT tags were programmed to release if the shark showed no significant vertical movement (± 1 m) for 96 continuous hours.

Data Analyses
Tags were programmed to collect depth (minimum resolution: 0.5 m), ambient temperature (minimum resolution: 0.05 • C), and light intensity (measured as irradiance at 550 nm wavelength) every 10 s for the duration of each deployment. Three subsets of deployed tags were programmed to detach from the sharks after 120 days (n = 7; T 120 ; March, 2007), 240 days (n = 7; T 240 ; July, 2007), and 360 days (n = 6; T 360 ; November, 2007). Upon release, pre-processed binned data (summarized every 6 h [120 and 240 day deployments] or 12 h [360 day deployments]; proportion of time spent at depth [m]: <5, 10, 25, 50, 75, 100, 150, 200, 300, 400, 500, 600, 800, >800; proportion of time spent at ambient temperature: 2-26 • C at 2 • C intervals) were transmitted through the Argos satellite system over a period of 6-10 days. Each tag also transmitted a temperature-depth profile of the water column inhabited by the shark during each 6-or 12-h interval. This comprised minimum and maximum water temperatures in eight depth bins, ranging from the shallowest to the deepest swimming depths for that time interval. Finally, the tag processed the light level data onboard to correct for depth and subsequently estimated times of dawn, dusk, and midnight or midday.
Temperature profiles of the water column and each shark's possible association with a thermocline, were assessed by plotting depth vs. temperature for discrete periods of time when diving behavior indicated repeated use of a consistent depth range. Given our findings that porbeagles exhibited two distinct diving behaviors [i.e., epipelagic (0-200 m) and mesopelagic (200-1000 m)], we grouped them as "non-divers" and "divers, " respectively. This allowed the examination of depth associations, averaged across individuals within each behavior group, using contoured probability plots generated with the ffgrid package in MATLAB (version R2010a).
To estimate the horizontal movements of each shark, we then used the HMMoce package , which uses light-based position estimates (from the tag), sea surface temperature (SST), and depth-temperature profiles, to construct gridded likelihoods of occupancy using a state space hidden Markov model framework. Observation-based likelihoods were derived from remotely-sensed SST, light-based longitude, and depth-temperature profile data collected by the tags, using four separate likelihood calculations: (1) an SST likelihood was generated for tag-based SST values integrated according to an error term (+/-1%) and compared to remotely-sensed SST from the Multi-scale Ultra-high Resolution (MUR) sea surface temperature analysis (MUR-SST, 0.01 • resolution) (NASA JPL, 2015); (2) light-based longitudinal likelihood was derived using estimates of longitude from GPE2 software (Wildlife Computers, Inc., Redmond, WA, USA) that allowed visual checking of light curves; (3) depth-temperature profile summaries recorded by the tag were compared to daily reanalysis model depth-temperature products from the HYbrid Coordinate Ocean Model (Bleck, 2002;Chassignet et al., 2007; HYCOM, 0.08 • resolution) at the standard depth levels available in the model and individual likelihood surfaces for each depth level were then multiplied together for an overall profile likelihood at that time point; and (4) Ocean Heat Content (OHC) was obtained by integrating the heat content of the water column above the minimum daily tagmeasured temperature to the most shallow depth recorded by the PSAT and included in the HYCOM fields (Luo et al., 2015). Start and end locations were fixed in all model runs.
The resulting observation likelihoods, in all reasonable pairwise and triplicate combinations, were convolved with a diffusive movement kernel that allowed swim speeds up to 2 m/s (see Braun et al., 2018 for full details of the convolution, filtering, and smoothing components of the model). Model selection comparing different combinations of observation likelihoods used Akaike Information Criterion (AIC). Results from the final smoothing step of the selected model represent the posterior distribution of each state (position) over time. The spatial means of the daily posterior distributions were used to calculate a most probable track.
Residency distribution (RD) plots were generated using the posterior likelihood surfaces from the position estimates. RDs were determined by aggregating the daily posterior probabilities across all tagged sharks, seasonally. Gridded probabilities were calculated with a 0.1 • resolution and seasons were based on the lunar calendar: fall, September 20-December 19; winter, December 20-March 19; spring, March 20-June 19; summer, June 20-September 19.

Horizontal Movements
The total minimum straight-line distance traveled by each of the tagged sharks ranged from 937 to 1,903 km for T 120 , 1,579 to 2,857 km for T 240 , and 2,517 to 3,310 km for T 360 ( Table 1). The porbeagles largely remained in the temperate WNA yearround, occurring in the Gulf of Maine, the Scotian Shelf, on George's Bank, and in the deep, oceanic waters off the continental shelf along the edge of, and within, the Gulf Stream (Figure 1). Geoposition estimates of latitude ranged from 35 • N (off the Outer Banks, North Carolina, USA) to 45 • N; longitude ranged from shelf waters east to 55 • W (Figure 1). Residency distribution plots show that tagged porbeagles exhibited seasonal changes in distribution (Figure 2). During the fall, the sharks were broadly distributed from the Gulf of Maine to oceanic areas, but largely concentrated in deeper regions of the former. In the late fall and winter, porbeagles moved out of the Gulf and off the shelf into oceanic habitat along the northern edge of and within the Gulf Stream (Figure 2). In the spring, the porbeagles remained in oceanic habitat along the shelf edge prior to returning to the Gulf of Maine and Scotian Shelf in the summer (Figure 2).

Vertical Movements
For all sharks combined, ambient water temperature and swimming depth ranged from 2 to 24 • C and from the surface to 1,304 m, respectively. During the late fall and winter months, two general patterns of vertical behavior were observed: (1) "nondivers" (2 females and 5 males) largely remained at shallow depths and rarely penetrated depths >600 m (Figure 3); and (2) "divers" (8 females and 2 males) made frequent dives into and remained at mesopelagic depths (200-1000 m) for 38-109 days (average 71 ± 26 days from fall into winter), rarely moving into the upper epipelagic zone (Figure 4). Overall, these two distinct dive patterns did not appear to be influenced by the broad-scale horizontal (coastal to offshore) movements of these fish as there were no significant differences in total minimum straight-line distances traveled between divers (2,450 ± 659 km) and non-divers (2,288 ± 715 km; p > 0.05, Student's paired Ttest). When stratified by dive behavior, the estimated geopositions of the porbeagles during the late fall and winter indicate that non-divers were almost exclusively associated with the cold side of the northern edge of the Gulf Stream, while divers primarily occupied the warm side of this edge, within the Gulf Stream, or in warm water south of the Gulf Stream (Figure 5).
When analyzed together, all fish, regardless of season, spent 75% and 50% of their time at depths of <200 m and 75-200 m, respectively, and 97% of their time in temperatures of 6-20 • C. In general, porbeagles occupied similar depth and temperature distributions, but divers and non-divers diverged during the fall and winter months as they moved off the shelf (Figures 6, 7). In the fall, divers occupied a broad range of depths from 5 to 800 meters, spending 31% of their time >200 m (Figure 6), and 68% of the time in water temperatures ranging from 8 to  14 • C (Figure 7). In contrast, non-divers did not move frequently to mesopelagic depths in the fall spending only 11% of the time >200 m, but in a similar temperature range (79% at 8-14 • C) when compared to divers. During the winter months, divers exhibited a bimodal depth and temperature distribution comprised of a shallow (40% at 75-200 m), warm (17% at 18-20 • C) component and a deep (43% at 300-800 m), cool (22% at 12-14 • C) component (Figures 6, 7). However, non-divers remained at shallow depths in the winter (89% <200 m), yet occupied a similar temperature range to that of the divers (Figure 7). In addition, these fish were not associated with the warmer (18-20 • C) waters encountered by the divers during the same time period. In the spring, both groups returned to epipelagic depths along the edge of the continental shelf with 70% and 92% of the time between 75 and 200 m and in an ambient temperature range of 6-14 • C, respectively (Figures 3,  4, 6, 7). In the summer, all fish were <200 m (99% of the time) experiencing a wider range of ambient temperatures (from 4 to 18 • C) for 97% of the time, but showed a bimodal depth distribution with shallow (59% <25 m) and deeper (25% at 100-200 m) components (Figures 3, 4, 6, 7).

DISCUSSION
Based on fisheries-independent data derived from 17 PSAT tags, we show that juvenile porbeagles tagged in the western North Atlantic exhibited broad, seasonally-dependent horizontal and vertical movements, including movement into the Gulf Stream.
In the fall, juvenile porbeagles were associated with shelf waters in the Gulf of Maine and George's Bank (Figure 2). The sharks tagged in this study were captured on the northern side of Georges Bank in 90-140 m of water as part of 88 juveniles hooked on only two longline sets in this area-this was clearly indicative of high abundance in the area. As fall transitions to winter, juvenile porbeagles moved out of the Gulf of Maine, off the shelf, and into the oceanic waters of the WNA for the winter months, during which some of the sharks moved into the Gulf Stream and exhibited mesopelagic residence. As northeastern shelf waters warm during the spring months, these sharks transitioned back to the shelf where they spent the summer in the Gulf of Maine and on the Scotian Shelf moving between the surface and the bottom. Although some of these results are consistent with historical fisheries-dependent conventional tagging data (Kohler and Turner, 2019), this study provided evidence of more expansive offshore movement in the WNA. All the porbeagles tagged in this study remained in temperate waters north of 35 • N, which differs from other pelagic shark species that seasonally move into more southerly latitudes, including the blue shark These results are in sharp contrast to previous findings for juvenile porbeagles tagged with this technology in the WNA. Campana et al. (2010) concluded that males and juvenile sharks of both sexes remained primarily in cool temperate waters on the continental shelf, always north of latitude 38 • N, for periods of up to 348 days after tagging. However, their findings are likely linked to sample size and composition. Of the 21 sharks tagged in that study, 52% were adult females, two were adult males, and the remaining eight were juvenile females. Most of the tracks of the latter were relatively short (<76 days) and did not span the winter months during which our tagged juveniles moved  offshore. Campana et al. (2010) did show that adult females move into the Gulf Stream and as far as the Sargasso Sea. Our findings demonstrate that juveniles also move into and south of the Gulf Stream, but not as far south as adult females (∼20 • N). To date, three PSAT tagging studies have been conducted on porbeagles in the ENA: Pade et al. (2009;n = 4;UK), Saunders et al. (2011;n = 3;Ireland), and Biais et al. (2017;n = 9; Bay of Biscay). Although samples sizes were relatively low, collectively these studies show that porbeagles are associated with relatively shallow shelf waters during the spring and summer months and move to offshore, deeper waters in the fall and winter months. Most recently, Biais et al. (2017) found that porbeagles tagged in the Bay of Biscay moved extensively (up to 2,000 km) in the ENA as far as north as the Arctic Circle, west to the mid-Atlantic Ridge, and south to Madeira in late summer and returned to the Bay of Biscay in spring of the following year. However, like those tagged by Campana et al. (2010), these fish were primarily sub-adult and adult females, which appear to migrate more extensively than juveniles.

Diving Behavior and Water Temperature
As porbeagles migrated seasonally from shelf to offshore habitat, the juveniles tracked during this study also exhibited seasonal changes in their vertical and thermal distribution (Figures 6, 7). Not surprisingly, other pelagic fish have been shown to exhibit similar seasonal shifts in their depth distribution, which are associated with foraging behavior. For example, the movement of blue sharks from shallow to deep habitat is thought to be correlated with prey availability (Queiroz et al., 2010). Similarly, the movements of juvenile porbeagles that we observed are likely linked to foraging behavior mediated by temperature.
In the WNA, Joyce et al. (2002) examined the stomach contents of 1,022 porbeagles and concluded that this species is an opportunistic feeder on a diverse assemblage of pelagic, epipelagic, and benthic fish species, as well as cephalopods, depending on availability. In juveniles (<150 cm), they found that a variety of fishes and cephalopods constituted major portions of the diet. They also observed a strong seasonal diet shift linked to spatial changes in shark distribution, although they grouped the stomach samples coarsely into spring (January-June) and fall (July-December) "seasons" and, hence, do not coincide with true seasons. Regardless, there are some indications that the movements observed in our study are associated with feeding behavior. In stomachs sampled during the months of January to June, Joyce et al. (2002) found a high proportion of cephalopods and pelagic teleosts (e.g., lancetfish, Alepisaurus ferox; Atlantic herring, Clupea harengus). In our study, this roughly coincides with offshore foraging behavior at epipelagic and mesopelagic depths during the winter. With progression into the autumn months, juvenile porbeagles transitioned off the shelf (Figure 2), expanding their vertical behavior. During the late fall and winter, they were broadly distributed along the shelf edge and in offshore regions of the WNA, including the Gulf Stream. Those that we qualify as non-divers remained north of the Gulf Stream or along its northern edge and moved through epipelagic depths (<200 m), spending 84% of the time in a temperature range of 8-16 • C. Those juveniles that migrated further offshore into the warm waters of the Gulf Stream and beyond (i.e., divers, Figure 5) shifted from epipelagic to mesopelagic habitat and exhibited a bimodal temperature distribution reflective of warm, shallow water and cooler waters at depth. This suggests that, in a manner similar to blue sharks (Carey et al., 1990;Braun et al., 2019) and white sharks (Gaube et al., 2018), these juvenile sharks may potentially dive deep within the Gulf Stream and associated features (i.e., fronts, eddies) in order to exploit more abundant food resources, including pelagic teleosts and cephalopods.
The importance of cephalopods in the diets of large pelagic fishes is well documented in the world's oceans (Smale, 1996). Off the coast of New Zealand, cephalopods were found to be a major component of the porbeagle's diet (Horn et al., 2013). In the offshore waters of the Gulf Stream and Central North Atlantic, Logan et al. (2013) examined stomach samples from nine species of large pelagic tunas, swordfish, and sharks and found that ommastrephid squids were the most common prey group across predator species. Other important prey included octopod, histioteuthid, and architeuthid cephalopods, as well as mesopelagic fishes. Although porbeagle sharks were not sampled by that study, the predominance and abundance of cephalopod species in the region, and reflected in the diets of other pelagic predators, make them a probable component of the porbeagle's diet. During the late fall and winter months, some of the juvenile FIGURE 7 | Seasonal time-at-temperature histograms for porbeagle divers (solid bar; n = 10) and non-divers (hatched bar, n = 7) tagged in the western North Atlantic.
porbeagles tagged in this study moved into these offshore areas (Figure 5), which are dominated by cephalopod prey (Logan et al., 2013). Thus, the temporal and spatial distribution of these key prey species may be driving the seasonal differences in depth distribution observed in juvenile porbeagle sharks in this study.
During the spring months, when the porbeagles moved back to the shelf edge and George's Bank (Figure 2), they all inhabited epipelagic depths and a predominant temperature range of 6-14 • C. During the summer, the porbeagles were shelf-oriented (Figure 2), spending all of their time <200 m and 59% of the time in the top 25 m (Figure 6). Moving from the surface to the bottom through a well-stratified water column, these sharks spent 97% of the time in the temperature range of 4-18 • C (Figures 3A, 4A). In those fish sampled from July to December, Joyce et al. (2002) found a greater proportion of benthic teleosts (groundfish) and fewer pelagics in their stomachs. This reflects summer movement onto the Northeast Shelf Ecosystem, which includes the Gulf of Maine and Georges Bank, and is historically one of the most productive ecosystems in the world, supporting large numbers of fish, invertebrate, and cetacean species (Bigelow andSchroeder, 1953: Kenney andWinn, 1986;German, 1987;Sherman et al., 1988;Fogarty and Murawski, 1998).
The seasonal depth distribution that we observed in porbeagle sharks differs markedly from those previously described in the WNA based on fisheries-dependent catch data (Campana and Joyce, 2004) and PSAT tags (Campana et al., 2010), which suggested that juvenile sharks do not commonly dive >200 m. However, such behavior was observed in the large adult females tagged by Campana et al. (2010). In a pattern remarkably similar to what we observed in our "divers, " these mature females encountered the Gulf Stream between 22 December and 9 March, which resulted in an "almost instantaneous initiation of deep diving behavior (daily maximum depth of <248 m before entry and a mean of 845 m after entry)" (Campana et al., 2010). These authors suggested that the porbeagles were migrating beneath the main flow of the Gulf Stream to maximize their net swimming speed, to minimize their ambient temperature, and to ultimately reach the Sargasso Sea for parturition (Campana et al., 2010). Our findings that juvenile porbeagles also exhibit this behavior suggest that foraging, and not reproduction, is likely a key driver of this behavior. As noted by Pade et al. (2009), it is becoming increasingly apparent that the porbeagle exhibits considerable plasticity in horizontal and vertical habitat use and the variable behaviors likely reflect the need for different search strategies depending on habitat and prey types encountered.
Collectively, the results presented in this and previous studies indicate that while porbeagles have broad thermal tolerance, their vertical and horizontal movements are driven by an optimal temperature range. Although the juvenile porbeagles tagged in this study exhibited extensive seasonal shifts in habitat use, both vertically and horizontally, and occupied a broad temperature range of 2-24 • C, they spent 97% of the time at 6-20 • C. The lower limits reported herein are colder than previous findings using PSAT technology in the WNA (8.0 • C, Campana et al., 2010), ENA (9 • C, Saunders et al., 2011), and Southern Hemisphere (4.6 • C, Francis et al., 2015), but it should be noted that our tagged fish spent only 2% of the time in water <6 • C. The sharks in this study also moved into water that was warmer than the maximum reported in the ENA (18.5 • C, Pade et al., 2009), but similar to that reported by Campana et al. (2010;25.4 • C). This is associated with movement into the Gulf Stream, yet temperatures in excess of 20 • C constituted <1% of the observations. Campana et al. (2010) suggested that the deep diving behavior observed in adult females while in Gulf Stream and Sargasso Sea was driven by water temperature because the surface waters (>22 • C) were uninhabitable. Our findings also suggest that the surface waters (<200 m) of the Gulf Stream and temperatures in excess of 20 • C represent a thermal barrier that influence the vertical behavior of juvenile porbeagles.
It is important to consider the extent to which, if any, the depth and temperature distribution of porbeagles are influenced (or limited) by their unique thermal physiology. While the body temperature of most sharks closely matches ambient water temperature (i.e., ectothermy) due to loss of all metabolically produced heat via conductive and convective transfer to the surrounding water (Brill et al., 1994), a few sharks (i.e., lamnids and one alopiid) have evolved a suite of morphological adaptations that allow certain regions of their body (i.e., eye, brain, viscera, aerobic swimming muscles) to be maintained warmer relative to ambient temperature (i.e., regional endothermy; Teal, 1966, 1969a,b;Carey et al., 1971Carey et al., , 1985Block and Carey, 1985). Regional endothermy may allow these sharks to maintain a more thermally stable internal operating environment during their forays across thermal fronts (Bernal et al., 2001(Bernal et al., , 2018. Thus, the selective advantages of regional endothermy may have allowed these species, including porbeagles, to expand their thermal niche and make accessible the additional food resources of the cooler, more productive waters at both a greater depth and at higher latitudes (Block and Carey, 1985;Carey et al., 1985;Bernal et al., 2001;Dickson and Graham, 2004). Carey et al. (1985) compiled in vivo muscle and visceral temperatures and compared a suite of anatomical attributes associated with metabolic heat production and conservation (e.g., relative heart size, the amount and distribution of red muscle, the number of vessels in the lateral cutaneous rete) in lamnid sharks and ranked each species by its capacity to elevate its temperature. These endothermic rankings of Carey et al. (1985) are now being validated by PSAT technology, which is generating new horizontal and vertical data about habitat use by these species. For example, it has now been demonstrated that the Pacific salmon shark, ranked first, occupies subarctic winter habitat and demonstrates the greatest tolerance of cold waters, spending 68% of its time in waters cooler than 10 • C (Weng et al., 2005). The endothermic capacity of the porbeagle, ranked second by Carey et al. (1985) likely allows this species to remain in northern latitudes year-round in the North Atlantic, spending 46% of the time in waters <12 • C based on our findings. In sharp contrast, the white shark (ranked third by Carey et al., 1985) migrates seasonally into temperate latitudes of the WNA, including the Gulf of Maine and the Scotian Shelf, but returns south as water temperatures cool in the late fall; it spends 95% of its time at temperatures >12 • C (Skomal et al., 2017). Similarly, the shortfin mako (ranked fifth by Carey et al., 1985), migrates seasonally to northeastern latitudes (59 • N) in the WNA (Vaudo et al., 2016;Kohler and Turner, 2019), but habitat use is strongly associated with warmer features and water temperatures between 10 and 15 • C represent the lower limit for this species, thereby constraining them horizontally and vertically over the course of the year (Vaudo et al., 2016).

Management
Although juvenile, subadult, and adult porbeagles have been PSAT-tagged on both sides of the Atlantic, none, including the juveniles tagged in this study, have provided evidence of trans-Atlantic movement. The juvenile porbeagle sharks tagged in this study demonstrated an inshore-offshore movement pattern consistent with fidelity to the WNA. Similarly, Biais et al. (2017) found evidence for porbeagle site fidelity and consistent return migrations to and from the Bay of Biscay in the ENA. These porbeagles exhibited migrations as far west as the Mid-Atlantic Ridge, with subsequent movements back to the Bay of Biscay, thereby supporting the idea of a separate ENA stock (Biais et al., 2017). Collectively, these findings support the two stock hypotheses put forth by Pade et al. (2009) and Saunders et al. (2011), and are consistent with conventional tagging data (Kohler and Turner, 2019), in which there is only one documented instance of an immature female porbeagle moving 4,260 km during a trans-Atlantic crossing from Ireland to Canada (Cameron et al., 2018). Nonetheless, molecular analyses have found that porbeagles sampled in the WNA and ENA are not genetically distinct, suggesting some level of migration and mixing (Testerman, 2014). Most recently, Haugen (2020) conducted an interdisciplinary review of all available information (i.e., life history, genetics, and movement) to evaluate North Atlantic stock structure in this species and concluded that geographic variation in life-history traits and movement patterns indicate one phenotypic stock in the WNA and another in the ENA, despite genetic analyses indicating gene flow and genetic homogeneity between the two regions. Based on our findings, coupled with this information, management of the porbeagle should be based on two stocks.
In the North Atlantic, major harvesters of porbeagle have either ceased landing (European Union, Canada) or capped landings (USA) in an effort to allow for the rebuilding of stocks (Campana et al., 2015;reviewed by Haugen, 2020). However, porbeagle bycatch occurs in several fisheries including pelagic longlines, gillnets, and trawls (Curtis et al., 2016;Haugen, 2020). In this study, we found that juvenile porbeagles ranged from shelf waters to the Gulf Stream and are, therefore, susceptible to these fisheries. A recent analysis of the overlap of shark distribution and pelagic longline fishing effort concluded that porbeagle occur in the highest-risk zone in the North Atlantic with high potential for mortality as a result of incidental bycatch (Queiroz et al., 2019). Although the ICCAT (2009) assessment concluded that longline bycatch on the high seas is minor and does not pose a significant threat to the species, at-vessel and post-release mortality is thought to be high (average = 47%; reviewed by Curtis et al., 2016). In the current study, we captured 88 juvenile porbeagles on two pelagic longline sets (525 hooks total) with relatively short soak times (3-5 h) and 21 (24%) sharks were dead at haulback. This not only indicates that porbeagle catch rates can be high on pelagic longlines when these sharks are associated with the shelf, but could lead to significant at-vessel mortality.
In addition to pelagic longline gear, porbeagles are also taken by trawl and gillnet fisheries while in shelf waters. A recent analysis of observer data collected from US gillnet and trawl fisheries indicated that annual bycatch and discards of porbeagles can be on the order of several hundred metric tons in recent years (Haugen, 2020). Unfortunately, these estimates have not been incorporated into stock assessments for this species. Clearly, bycatch mortality might be more substantial than previously thought and all sources need to be taken into consideration in future stock assessments. Our findings show a relatively high risk for juveniles in multiple fisheries across the distribution of this species and highlight the need for continued protection for stock rebuilding.
In this study, we demonstrate that juvenile porbeagles are more broadly distributed in the WNA using near-coastal, shelf habitat during the warm months and migrating to pelagic habitat from the shelf edge to the Gulf Stream during the winter. When in the latter, these sharks exhibited mesopelagic foraging behavior likely driven by cephalopod prey abundance and water temperature. The residency of juvenile porbeagles in temperate latitudes throughout the year is likely facilitated by its endothermic capacity. Additional studies that link these movements with the foraging ecology and physiology of this species are warranted.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

ETHICS STATEMENT
All methods associated with the capture and handling of sharks were reviewed and approved by Institutional Animal Care and Use Committee of the University of Massachusetts, Dartmouth (Protocol #05-07).

AUTHOR CONTRIBUTIONS
GS and DB conceived the study. GS, HM, and LN performed the field work. GS, HM, BG, and CB conducted the data analyses. GS, HM, and DB wrote the initial draft and all authors provided edits and approval. All authors contributed to the article and approved the submitted version.