Seasonal Movements and Habitat Use of Juvenile Smooth Hammerhead Sharks in the Western North Atlantic Ocean and Significance for Management

Upper trophic level predators dramatically impacted by fisheries include the large-bodied hammerhead sharks, which have become species of conservation concern worldwide. Implementing spatial management for conservation of hammerhead populations requires knowledge of temporal distribution patterns and habitat use, identification of essential habitat for protection, and quantification of interactions with human activities. There is little such information for the smooth hammerhead shark, Sphyrna zygaena. We used fin-mounted satellite tags to examine the movements and habitat use of juvenile smooth hammerheads, a demographic segment particularly threatened by exploitation. Six sharks were tagged off the US mid-Atlantic and tracked for 49–441 days (mean 187 ± 136 days). Sharks consistently showed area-restricted movements within a summer core area in waters of the New York Bight and a winter core area off Cape Hatteras, North Carolina, with directed movements between them in autumn. There was high overlap of shark winter core area use and the Mid-Atlantic Shark Area (MASA) – a 7 month per year, bottom-longline fishery closure – indicating that this area closure offers seasonal reduction in fishing pressure for this species. Based on timing of shark movements and the MASA closure, protection for juvenile smooth hammerheads may be increased by beginning the closure period 1 month earlier than currently scheduled. Generalized additive mixed models revealed that area-restricted movements of sharks in their summer and winter core areas coincided with high primary productivity, and elevated sea surface temperature. Consistency in use of summer and winter core areas suggests that the coastal waters of the New York Bight and Cape Hatteras, North Carolina could be considered for Essential Fish Habitat designation for this species. This study reveals the first high resolution movements and habitat use for smooth hammerheads in the western North Atlantic to inform management planning for this population.


INTRODUCTION
The rapid expansion of elasmobranch fisheries and trade globally are principal drivers of population decline for many shark species (Dulvy et al., 2014). Some species, such as the large-bodied hammerhead sharks (great -Sphyrna mokarran, scalloped -S. lewini, and smooth -S. zygaena hammerheads), are especially vulnerable to fishing pressure because of their slow rates of population growth (Cortés et al., 2010) and high at-vessel and post-release mortality due to elevated stress response to capture (Morgan and Carlson, 2010;Eddy et al., 2016;Gallagher and Klimley, 2018). In addition, hammerhead sharks are taken in large numbers because of the high demand for their superiorquality fins (large size and high ceratotrichia count) in the global shark fin trade (Abercrombie et al., 2005;Clarke et al., 2006a,b;Cardeñosa et al., 2018). While population declines of the largebodied hammerhead shark species complex is thought to be largely driven by declines of scalloped hammerheads (Jiao et al., 2011), low catch rates in various parts of the world for all species suggest a significant historical decline in the abundance of all large-bodied hammerhead sharks (Baum et al., 2003;Baum and Blanchard, 2010;Ferretti et al., 2010).
Fishery exploitation of smooth hammerhead sharks via targeting or bycatch has been identified as the major threat to this species, particularly for juveniles (Casper et al., 2009;Cortés et al., 2010;Miller, 2016). The conservation of this species is an international priority, e.g., Vulnerable listing on the International Union for Conservation of Nature (IUCN) Red List (Casper et al., 2009); Appendix II listing on the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES); Appendix II listing on the Convention on Migratory Species of Wild Animals (CMS). Furthermore, an ecological risk assessment of sharks caught in Atlantic pelagic longline fisheries highlighted smooth hammerheads as a species in urgent need of biological data necessary for stock assessment (Cortés et al., 2010). However, conservation-relevant data on many aspects of the biology of smooth hammerheads are extremely limited, including information relating to their movements, seasonal distributions and habitat use (Miller, 2016;Gallagher and Klimley, 2018).
Achieving sustainable populations of fishery exploited species is critically dependent on the recruitment of immature individuals. Thus, identification and conservation of essential habitat for juvenile and sub-adult sharks is of paramount importance, requiring an understanding of this key demographic segment's spatiotemporal patterns of occurrence and associated oceanic environmental drivers (Kinney and Simpfendorfer, 2009;Schlaff et al., 2014). Equipped with adequate information on how environmental parameters influence seasonal movements, spatial management measures such as temporal closures of targeted areas can be enacted to promote recovery of overfished stocks. For example, based on understanding of temporal and spatial habitat use of dusky sharks (Carcharhinus obscurus), the Mid-Atlantic Shark Area (MASA) -a region closed to bottom longline fishing for 7 months each year -was established in 2005 as a means of reducing fishing mortality and enhancing recovery of this species (NMFS, 2009).
Smooth hammerheads have a circumglobal distribution in coastal and oceanic waters and occupy a wider latitudinal range than other sphyrnids (Compagno, 1984). Catch records from a variety of locations suggest that juveniles and sub-adults [<265 cm total length (TL)] are more common in inshore waters over coastal shelves, with larger individuals (>265 cm TL) found more frequently offshore (Diemer et al., 2011;Clarke et al., 2015;Francis, 2016;Deacy et al., 2020). This species is capable of long distance movements (e.g., 6,610 km over 150 days; Santos and Coelho, 2018), but also shows high levels of resident behavior within restricted areas (at least 80 days; Diemer et al., 2011).
To date, just two studies have reported on the finer-scale movements of smooth hammerheads via satellite telemetry. Horizontal and vertical movements of juveniles in the temperate western South Pacific appear to vary seasonally (Francis, 2016); in contrast, juvenile and adult smooth hammerheads in the tropical eastern Atlantic demonstrated relatively stable temporal diving behavior (Santos and Coelho, 2018). These studies illustrate the possibility that movements and behavior of smooth hammerheads may be influenced by thermal heterogeneity of their environment. Very little information on relationships between other environmental variables (primary productivity, oceanic fronts, water depth) and movements of smooth hammerheads exists (Couto et al., 2018), although these variables have been shown to greatly influence movement and habitat selection of highly mobile sharks (Block et al., 2011;Queiroz et al., 2016;Vaudo et al., 2017). Understanding interactions between patterns of habitat use and environmental variables would contribute toward construction of habitat models and an improved ability to predict the distribution of smooth hammerheads under climate change scenarios, as well as reveal potential interactions with human activities throughout their range.
Given minimal information on the spatial ecology of smooth hammerheads in general and conservation concerns for this species, our goal was to quantify habitat use and horizontal movements of juvenile smooth hammerheads in the western North Atlantic Ocean via satellite tag telemetry. Only very coarse scale information exists on movements of smooth hammerheads in this region, obtained from the recapture of just seven individuals out of 269 (0.02%) tagged with conventional identification tags over 52 years (Kohler and Turner, 2019). Our study objectives were to: (1) determine seasonal movements and distribution patterns; (2) identify core areas of habitat use; (3) evaluate the potential of the MASA seasonal closure for providing protection from fishing pressure, and (4) investigate relationships between movement behavior and environmental conditions, for juvenile smooth hammerheads.

Capture and Tagging
Between 22 July 2016 and 9 September 2017 six female, juvenile smooth hammerhead sharks were caught via rod and reel off the coast of Ocean City, Maryland United States (38.1 • N, 74.5 • W). Sharks were brought on board the fishing vessel where a saltwater hose was inserted into the mouth to irrigate the gills and then the sharks were measured, sexed, and fitted with a satellite-linked radio tag (SPOT-196 tag; Wildlife Computers, Redmond, WA, United States) on the dorsal fin. These tags directly communicate with the Argos tracking system 1 when the shark's dorsal fin breaks the sea surface exposing the tag to air, providing an estimated position (latitude and longitude) and an associated location class. Location class is determined by the number of transmissions received and the number of Argos satellites receiving transmissions, and categorized from most to least accurate as 3, 2, 1, 0, A and B. Estimated errors (1 SD) for each location class are LC 3: < 250 m, LC 2: 250-500 m; LC 1: 500-1500 m, and LC 0: > 1500 m; there is no spatial estimate of accuracy for LC A and B (CLS, 2016). The two tags deployed in 2016 were programmed to transmit for 1 h every other hour; the four tags deployed in 2017 were programmed to transmit for 1 h every 4 h to try to obtain longer duration tracks.

Movements, Distribution, and Behavior
Because Argos positions of sharks varied in temporal frequency and spatial accuracy, we obtained standardized positions (hereafter "positions") at 12 h intervals that were comparable between individuals and over time by processing Argos locations using a behavioral switching state-space model (SSM) within a Bayesian framework developed by Jonsen et al. (2005). Since parameter estimation is improved when conducted jointly across multiple individual datasets (Jonsen, 2016), we produced most probable tracks using a hierarchical joint estimation model (hSSM) that produced temporally regular positional estimates based on the Argos location class, mean turning angle, and autocorrelation in speed and direction. Previous research has shown that the accuracy of the hSSM parameter estimates declines in response to outlier locations (from poor quality satellite positions) and long gaps in detection data (Bailey et al., 2008); therefore, prior to fitting hSSMs, each track was filtered using the argosfilter package (Freitas et al., 2008) in R Core Team (2014) with parameters listed in Vaudo et al. (2017). To reduce spurious results associated with long detection gaps, tracks were broken into multiple segments when gaps between Argos locations were >10 days. Resulting segments <20 days in duration were excluded from the hSSM (Block et al., 2011). Given that 84.7% of gaps between positions in our tracks were <12 h (Supplementary Figure S1), we used a time step of 12 h in the hSSM to produce two positions per day for each shark.
The hSSM model was fit by running two Markov Chain Monte Carlo (MCMC) chains in parallel for a total of 60,000 samples, with the first 50,000 being discarded as burn-in, and the remaining 10,000 samples thinned by retaining every 10th sample to reduce autocorrelation (n = 1,000 per chain). Each MCMC iteration provides not only a most probable track but also assigns each estimated location to one of two possible behavior modes (resident and transient). The final estimated track is the average of all 2,000 MCMC samples and the final output for each behavioral state represents the proportion of samples for a given position classified as resident (MCMC diagnostics given in Supplementary Figure S5). When the proportion is high (resident) or low (transient) the classification can be confidently assessed. Consequently, following Breed et al. (2009), we classified proportions ≤0.3 as transient, ≥0.7 as resident, and proportions of 0.3-0.7 as uncertain. The hSSM was fit using the bsam package (Jonsen et al., 2015) in R.
Using the hSSM positions, a seasonal utilization distribution (UD) was calculated for all sharks pooled across the meteorological seasons (summer: June-August, autumn: September-November, winter: December-February, and spring: March-May) using the adehabitat package in R (Calenge, 2006). The UD estimate was calculated following methods described in Vaudo et al. (2017).
To investigate vertical diel behavior in the absence of transmitted depth data, we used successful Argos transmissions as a proxy for surfacing behavior since locations are only obtained when sharks are at the surface (Doyle et al., 2015). Using the Argos Satellite Pass Prediction tool 2 , satellite pass data was obtained for all six available satellites from June 2017-September 2018. Because each satellite can simultaneously detect all transmitters within an approximately 5,000 km diameter circle below it (CLS, 2016), satellite pass data was obtained for 40 • N and 74 • W, which encompassed all shark positions received. Because the number of satellites passing overhead varies by hour of the day (in effect increasing the amount of listening effort when more satellites are present; Supplementary Figure S2), surfacing behavior was determined by summing the number of Argos locations obtained per shark per hour (Eastern Standard Time), and dividing by the cumulative amount of time that all satellites were overhead during each hour (in general, each satellite takes roughly 10 min to pass over a stationary object). The resulting value represents a standardized number of Argos locations per hour of satellite coverage (hereafter termed "surfacing index"), providing information on temporal patterns of surfacing, regardless of the number of satellites overhead. The surfacing index (square root transformed) was compared among hours of the day using a linear mixed effects (LME) model as surfacing index ∼ hour + ID, where surfacing index was the response variable, hour of day was the explanatory variable and shark ID was a random factor using the lmer function in the lme4 R package (Bates et al., 2014). Tests of multiple comparisons were obtained using the glht function in the multcomp package (Hothorn et al., 2008). This analysis of diel surfacing behavior was limited to Sharks 3-6 because satellite pass data is only retained by the Argos system for 1 year and this analysis was not undertaken until 2018; thus, satellite pass data could only be obtained for the four sharks tagged in 2017.

Environmental Variables
Water depth (m) and sea surface temperature (SST; • C) values were obtained using the NOAA ETOPO1 Global Relief Model (one arc-minute resolution) and the Multiscale Ultra-high Resolution (MUR) SST dataset 3 (0.01 • resolution), respectively, using the "xtractomatic" package in R (Mendelssohn, 2017). SST gradient (a proxy for temperature fronts) was calculated as the maximum difference in SST across a moving window of a 15 × 15 grid cell matrix (totaling ∼0.15 • area covered) using the raster package (Hijmans et al., 2017) in R. Finally, using the rerddap (Chamberlain et al., 2019) and rerddapXtracto 4 packages in R, we obtained 8-day composite primary productivity (PP) (mg C/m 2 /day; 0.0125 • resolution) data from the National Aeronautics and Space Administration's (NASA) Aqua satellite with its Moderate Resolution Imaging Spectroradiometer sensor (MODIS-Aqua).
Generalized additive mixed models (GAMMs) were used to determine the best environmental predictors of smooth hammerhead shark resident behavior. Prior to inclusion in the global model, univariate models were constructed with potential environmental predictors standardized by their mean and standard deviation, and tested using a likelihood ratio test.
Significant predictor variables were then tested for collinearity using a Pearson's rank correlation matrix (Zuur et al., 2009) and all non-spatial combinations were <0.7 (Supplementary Figure S3). The proportion of behavioral states categorized as resident for each position by the hSSM was used as the response variable and all predictor variables were included in the global model. The model was run using a Gaussian response distribution and identity link. The importance of various combinations of autocorrelation structures was tested while holding other variables constant. Similarly, to determine the best random effects structure aimed at accounting for any temporal effect (e.g., increasing temporal gaps between positions since tagging or season) or individual effect imposed on the sharks' behavior, we considered shark ID, season and days at liberty as possible random effects. The performance of the final model output was assessed using the C index, where values closer to 1 indicate better performance, and the corresponding Somers' Dxy rank correlation, which is a measure of ordinal association between the response and predictor variables (Lea et al., 2018).

RESULTS
The six juvenile smooth hammerheads TL (mean ± SD) 184.2 ± 18.5 cm were tracked for periods of 49-441 days and generated a total of 3,488 Argos locations. The number of Argos locations d −1 ranged from 0 to 21 (mean 3.1 ± 3.3). The mean time interval between Argos positions was 7.7 ± 33.8 h 4 https://github.com/rmendels/rerddapXtracto (median = 2.5 h). The number of days with Argos locations for each shark ranged from 46 to 263 days (mean 131 ± 72.5 days), resulting in a total of 786 days with locations out of 1121 days at liberty (mean 187 ± 136 days) ( Table 1). Among all sharks, this equates to being detected on 77.2 ± 0.1% of days at liberty. Once the Argos locations were filtered and standardized to a 12 h time interval using the hSSM, positions were removed for days lacking an Argos location. As a result, 1,531 positions remained, which served as the basis of subsequent analyses.
Most individuals displayed similar movements and habitat use throughout the course of the study (Figure 1A). In general, the hSSM indicated that sharks were resident in shallow water off southern Long Island, New York during the summer, with some southern movement to the waters off New Jersey, Delaware and Maryland in late summer. During autumn, directed southern movements through the mid-Atlantic region were common to all sharks, showing little affinity to any one region in the area, as indicated by the observation that 57% of positions were classified as transient during autumn. During winter and early spring, sharks displayed area restricted movements, primarily focused near the southeastern outer banks of Pamlico Sound, North Carolina ( Figure 1A and Supplementary Figure S4).
One shark (shark #5), tagged 17 June 2017 was tracked for 441 days (Figure 1b and Table 1). This shark was tagged near Ocean City, Maryland and remained there for almost 2 weeks after tagging, then moved north into the New York Bight in early July, where it remained until September when it began a 3 month journey south reaching the area off Cape Hatteras, North Carolina on 15 November. It remained in this area until 3 May 2018 and was not detected again until 9 July 2018 off New Jersey; the shark then moved to waters of the New York Bight until 1 September 2018. Another individual (Shark #1) moved south from the New York Bight similarly to other sharks, but continued moving south east Cape Hatteras eventually reached the coast of central Florida on 12 December 2016 at the time of the last detection 144 days after tagging (Figure 1a).
Seasonal utilization distributions showed similar seasonal movements as indicated by the hSSM (Figure 2). Core areas (50% UD) of the seasonal distributions were primarily centered in the New York Bight in the summer, expanded southward during the autumn as sharks moved south, and were concentrated off Cape Hatteras, North Carolina in the winter (Figure 2). The individual tracked for greater than 1 year moved north in late April/early May, in a similar manner to the northward movements of sharks shortly after tagging off Ocean City, Maryland (Figures 1b, 2d).   (Figure 1b), and transmissions from all other sharks stopped prior to exiting the MASA, so time spent within the MASA could not be assessed.
Because tags (sharks 3-6) deployed in 2017 were programmed to transmit just one out of every 4 h, diel vertical behavior is only described for the hours tags were set to transmit (0000-0100, 0400-0500, 0800-0900, 1200-1300, 1600-1700, and 2000-2100 h). Significant fixed effects for the 0400-0500 and 2000-2100 h blocks (0400-0500 LME Estimate = 0.25, SE = 0.07, t = 3.4, p = 0.004; 2000-2100 LME Estimate = 0.23, SE = 0.07, t = 3.1, p = 0.007) indicated that surfacing index varied over the 24 h diel period, and multiple comparisons revealed that sharks surfaced most frequently just before dawn, at midday, and just after dusk (Figure 4). The total number of Argos locations for all sharks pooled was greatest during the time interval 0400-0500 and 2000-2100 h (308 and 700 total Argos locations, respectively). The 2000-2100 h time interval coincided with the greatest number of satellite passes in the study area (all satellites combined = 794 passes, Supplementary Figure S2), resulting in 165.3 h of listening time, which yielded an overall surfacing index of 4.23 for all sharks combined. However, although the 0400-0500 h block only had 40.8 total h of satellite listening time, it had the highest surfacing index of 4.35 for all sharks combined (i.e., on average, there was roughly one position per shark per hour of satellite coverage just before and just after dawn and dusk, respectively; Figure 4).
After testing the importance of various combinations of autocorrelation structures while holding other variables constant, we found that the GAMM without an autocorrelation term was deemed more robust with better wAIC and AIC (wAIC = 0.31; Supplementary Table S1); thus, no autocorrelation structure was used in the final model. Additionally, wAIC and AIC revealed that treating shark ID and season as random effects resulted in the most parsimonious model, thus, days at liberty was not included (Supplementary Table S2). The final GAMM predicting residency behavior explained 34% of the sample variance (C index = 0.72, Dxy = 0.45, SD = 0.001, n = 1432). Mean SST, log of primary productivity and water depth were included in the best fit model, while SST gradient was removed given its lack of significance (p = 0.2) and improved model fit after removal ( AIC = 2.4, wAIC = 0.76). Model output indicated that most of the variation in the observed resident behavior was attributable to geographical location, followed by primary productivity concentration and depth ( Table 2). Probability of displaying resident behavior was highest at latitudes associated with the New York Bight (>40 • N), high primary productivity concentration [7.82 log(mg C/m 2 /day) ( Figure 5B)], and inshore neritic waters (<100 m; Figure 5C). In addition, SST of ∼18, 23 and >26 • C resulted in increased probability of sharks displaying resident behavior as these represented the temperatures experienced in core habitat areas ( Figure 5A and Table 2). FIGURE 3 | Locations of juvenile smooth hammerhead sharks (n = 4 for all groupings) within the mid-Atlantic Shark Area (MASA) during the closure period (January 1-July 1; orange points), and outside of the closure period (black points). Nearly all December locations (green points) fall within the MASA boundary but not within the closure period.

DISCUSSION
We provide the first detailed view of the movement dynamics of smooth hammerhead sharks in the western North Atlantic. Previous work has reported that smooth hammerhead sharks spend a large proportion of their time in surface or near-surface (<10 m) waters (Francis, 2016;Santos and Coelho, 2018), potentially making them good candidates for SPOT tags which only transmit data when exposed to air. Indeed, the sharks tracked here were detected on average 3.1 ± 3.3 times per day, and 84.7% of Argos locations occurred within 12 h of a previous location. The high frequency of satellite transmissions and Argos locations allowed for reconstruction of smooth hammerhead movements at a much higher resolution than has previously been described.
The sharks we tracked in the western North Atlantic Ocean displayed consistent seasonal movements between core areas of activity off Long Island, New York in summer and off Cape Hatteras, North Carolina in winter. Although seasonal movements of this species have been hypothesized previously based on surface sightings (Couto et al., 2018) and fisheries catch per unit effort data (Santos and Coelho, 2019) in the eastern North Atlantic, the telemetry results here provide a direct, fishery independent demonstration of this behavior by smooth hammerheads. Based on environmental characteristics of the core areas, sea surface temperature and productivity appear to be important drivers of their seasonal movement patterns, as has been demonstrated in other highly migratory marine megafauna (Weng et al., 2008;Block et al., 2011;Curtis et al., 2014;Kajiura and Tellman, 2016;Vaudo et al., 2017).
Seasonal movements and habitat use in other hammerhead species have been documented, but thus far suggest they are driven more so by foraging or reproduction, rather than dynamic oceanographic processes. For example, seasonal changes in abundance of scalloped hammerheads at offshore islands in the eastern tropical Pacific have been suggested as possibly related to movements for reproductive purposes and/or parturition, but currents and chlorophyll concentrations may also play a role in long term movements (Bessudo et al., 2011;Ketchum et al., 2014;Nalesso et al., 2019). Wells et al. (2018) found scalloped hammerhead movements in the northern Gulf of Mexico to be primarily driven by static bathymetric features rather than dynamic environmental variables and did not observe any seasonal patterns in shark movements. Furthermore, great hammerhead repeated seasonal presence and residency within the Bahamas is also believed to be related to reproduction or foraging, rather than climatic processes (Guttridge et al., 2017). However, sharks tracked here represent the juvenile to sub-adult size class of this species, and physiological tolerances to environmental conditions vary across ontogeny and may result in juveniles selecting different habitats than adults (Grubbs, 2010). Given that only one smooth hammerhead in our study was tracked for over a full year, it remains unclear how typical seasonal migratory behavior and environmental driven movement is in this species throughout its geographic and size range. Seasonal movement patterns of smooth hammerheads along the US East Coast was characterized by resident behavior during the summer and late winter/early spring. The timing of resident behavior coincided with increased levels of primary productivity, presumably tied to prey availability (Ware and Thomson, 2005;Priede and Miller, 2009). Stomach contents show that the major prey of smooth hammerheads is cephalopods (mainly ommastrephid squid) and small schooling fishes (Smale, 1991;Rogers et al., 2012;Bornatowski et al., 2014). The longfin squid Doryteuthis pealeii and shortfin squid Illex illecebrosus are the most common species of squid in the western North Atlantic from Georges Bank to Cape Hatteras, and both species undergo seasonal spawning migrations at northern and inshore locations in late spring/early summer and deeper, southern locations along the continental shelf edge in late autumn/early winter (Dawe et al., 2007;Jereb and Roper, 2010). Seasonal movements and aggregations in relation to high prey abundance has been reported in several species of sharks (Klimley et al., 1992;Heyman et al., 2001;Mourier et al., 2016), however, little information exists on smooth hammerhead diet in the study region to determine if they are taking advantage of this potential resource. Nevertheless, spawning and seasonal movements of these squid in the western North Atlantic spatially and temporally overlap with core areas used by smooth hammerheads tracked in this study.
The diel surfacing behavior patterns of smooth hammerheads tracked in our study may also be related to foraging. Highest surfacing indices were recorded shortly before dawn and after dusk, similar to the pattern observed in a juvenile smooth hammerhead (139 cm TL) tracked off the coast of New Zealand (Francis, 2016). Francis (2016) also reported diel differences in depth distribution of another juvenile smooth hammerhead tracked with a popup satellite transmitter, with a shallower distribution at night compared to daylight hours. In contrast, scalloped hammerhead sharks have been observed to remain in shallow waters during the day and dive at night presumably to forage (Klimley and Nelson, 1984;Hoffmayer et al., 2013), or show continuous deep diving behavior throughout the 24 h cycle (Spaet et al., 2017). Similarly, Santos and Coelho (2018) found that similarly sized smooth hammerheads to those in this study [T-test; T (6 . 5) = 1.1, p = 0.3] tracked using depth and temperature archival transmitters off the west coast of equatorial Africa occupied deeper, cooler water during the night compared to day. Our findings contrast somewhat with those of Santos and Coelho (2018); however, this difference may be an artifact of study location, where sharks tracked in Santos and Coelho (2018) were experiencing temperatures at depth several degrees warmer (26-27 • C) than SSTs observed here in the temperate western North Atlantic (mean 21.7 ± 2.2 • C). In addition, the surfacing index presented here is limited to when a sharks' dorsal fin breaks the surface and appropriate satellite coverage is overhead, so patterns observed here may not be fully representative of smooth hammerhead diel depth distribution in the western North Atlantic.
Decreasing population trends of smooth hammerheads have prompted conservation listings (e.g., IUCN, CITES, CMS) and calls for additional management. Because of high at-vessel (Coelho et al., 2012) and estimated post-release mortality of smooth hammerheads caught in fisheries (Braccini et al., 2012), reducing exposure to capture rather than relying on release after capture is a more effective management method for reduced fishing mortality. While acknowledging that our inferences are based on the four animals with long enough tracks, the consistent finding of the winter core area of activity largely falling within the boundaries of the MASA management zone during winter and spring, and high proportion of transmissions occurring within the MASA during the shark bottom longline fishery closure period (1 January-31 July), suggests the potential of the MASA for reducing fishing mortality of this species. Furthermore, as reported for sand tiger sharks (Carcharias taurus) (Teter et al., 2015), the smooth hammerhead spatial and temporal patterns of movement suggest that beginning the MASA closure on 1 December, rather than 1 January, would provide additional and extended protection from commercial fisheries for this species also.
Though there was some individual variability in movements of smooth hammerheads tracked in our study with a limited number of individuals, the high degree of spatial and temporal consistency demonstrated by the sharks in use of both summer and winter core areas as well as behaviors associated with foraging suggest that the coastal waters of the New York Bight and Cape Hatteras, North Carolina could be considered for designation of Essential Fish Habitat (EFH) for this species in the western North Atlantic, an important designation for protection consideration in U.S. fisheries management practices (NMFS, 2009) 5 . Seasonal movement between southern areas of increased presence in winter and northern areas of concentrated activity in summer have been reported for other species of sharks in the western North Atlantic, including sandbar sharks Carcharhinus plumbeus (Grubbs et al., 2007;McCandless et al., 2007;Conrath and Musick, 2008), dusky sharks Carcharhinus obscurus (Musick and Colvocoresses, 1986), sand tiger sharks Carcharias taurus (Teter et al., 2015) and white sharks Carcharodon carcharias (Curtis et al., 2018), and has led to delineation of nurseries and designation of EFH for several of these species (NMFS, 2009). Likely due to the lack of available data, there is currently no EFH in U.S. waters for smooth hammerhead sharks. 5 https://www.fisheries.noaa.gov/resource/map/essential-fish-habitat-mapper Successful management of populations is dependent on the survival of young individuals and recruitment to reproductive stock; therefore, understanding movement patterns, habitat use and EFH of juveniles is vital. In this study we have identified both winter and summer core areas of concentrated activity for juvenile smooth hammerheads in the western North Atlantic, as well as pathways traveled between those seasonal core areas. In addition, environmental conditions associated with resident behavior within these core areas and timing of directed movements between them enables improved ability to predict inter-and intra-annual distribution of smooth hammerheads, and how this may change over time with changing environmental conditions (e.g., increasing sea surface temperatures). These advances in understanding patterns of distribution and habitat use of juvenile smooth hammerheads in the western North Atlantic are directly applicable to effective management of this demographic component of their population. Future work should include studying the movement ecology of adult smooth hammerhead sharks of both sexes since their movements and habitat use patterns are likely to be different from those of juveniles.

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

ETHICS STATEMENT
The animal study was reviewed and approved by the Nova Southeastern University IACUC #DB1.

AUTHOR CONTRIBUTIONS
BW and MSS designed and implemented the research. MS led the fieldwork. RL led the analysis of the data with assistance from LS and JV. RL led the writing of the manuscript with assistance from all authors.

FUNDING
This work was supported by grants to MSS from the Guy Harvey Ocean Foundation (GHOF18-1), Save Our Seas Foundation (SOSF157), Shark Foundation/Hai-Stiftung (HF-19-2), the Levitetz Family Foundation , and a scholarship to RL from Fish Florida.