Movement, Behavior, and Habitat Use of a Marine Apex Predator, the Scalloped Hammerhead

Conservation and management efforts of marine apex predators are more reliable when information on movement and habitat use patterns are known. The scalloped hammerhead (Sphyrna lewini) was the first shark species to be protected under the U.S. Endangered Species Act and has life history characteristics that make this species particularly at risk for local depletion. Consequently, the goal of this study was to better understand the movement dynamics of this species in the Gulf of Mexico (GOM) where discards through the longline fishery can be substantial. A total of 33 scalloped hammerheads were tagged with fin mounted satellite tags and tracked for an average of 146 days (ranging from 5 to 479 days) to examine horizontal movements and quantify space use. Scalloped hammerheads showed a wide range of movements throughout the GOM continental shelf with limited long-distance dispersal and females displayed a shelf-edge association relative to more mid-shelf use by males. A generalized additive model was developed to identify habitat suitability for scalloped hammerheads in the GOM, while state-space modeling was used to examine movement behaviors. Model results highlighted the use of continental shelf waters with high occurrence at close proximities to both artificial and hard-bottom habitat combined with low chlorophyll a concentrations (~0–4 mg m-3) and moderate salinities (33–35.5). Habitat suitability for scalloped hammerheads was predicted to be high on the mid to outer continental shelf inside the 200 m isobaths and state-space models support area-restricted behavior was most common relative to transient behavior. Findings from this study provide important information on movement of this species in the GOM and highlight their restricted use of continental shelf habitat and resident behavior that will need to be incorporated in future stock assessments and extinction risk analyses.

Conservation and management efforts of marine apex predators are more reliable when information on movement and habitat use patterns are known. The scalloped hammerhead (Sphyrna lewini) was the first shark species to be protected under the U.S. Endangered Species Act and has life history characteristics that make this species particularly at risk for local depletion. Consequently, the goal of this study was to better understand the movement dynamics of this species in the Gulf of Mexico (GOM) where discards through the longline fishery can be substantial. A total of 33 scalloped hammerheads were tagged with fin mounted satellite tags and tracked for an average of 146 days (ranging from 5 to 479 days) to examine horizontal movements and quantify space use. Scalloped hammerheads showed a wide range of movements throughout the GOM continental shelf with limited long-distance dispersal and females displayed a shelf-edge association relative to more mid-shelf use by males. A generalized additive model was developed to identify habitat suitability for scalloped hammerheads in the GOM, while state-space modeling was used to examine movement behaviors. Model results highlighted the use of continental shelf waters with high occurrence at close proximities to both artificial and hard-bottom habitat combined with low chlorophyll a concentrations (∼0-4 mg m −3 ) and moderate salinities (33-35.5). Habitat suitability for scalloped hammerheads was predicted to be high on the mid to outer continental shelf inside the 200 m isobath and state-space model results suggest area-restricted behavior was most common relative to transient behavior. Findings from this study provide important information on movement of this species in the GOM and highlight their restricted use of continental shelf habitat and resident behavior that will need to be incorporated in future stock assessments and extinction risk analyses.

INTRODUCTION
Information on movement and habitat use of marine apex predators is needed to better understand population structure and to implement spatially explicit management strategies. Apex predators provide an important role in ecosystem stability since they are top consumers in the food web (Myers and Worm, 2003;Scheffer et al., 2005). Consequently, proper management of marine predators requires a fundamental understanding of the environmental processes driving their movements and behavior (Schlaff et al., 2014). Identifying important habitat and the factors responsible for movement are inherently difficult due to the mobility of large marine predators as they often traverse multiple jurisdictional boundaries, ecosystems, and habitats (Chin et al., 2017). Moreover, patterns of habitat use and residency are influenced by dynamic oceanographic conditions (e.g., shifting eddies or currents) and distribution and movement of prey resources (Brill et al., 1999;Queiroz et al., 2016). Physicochemical properties are increasingly used to delineate habitat boundaries and predict migratory trajectories (Block et al., 2001;Teo et al., 2007;Sequeira et al., 2012), and habitat utilization patterns have been described for several marine species, including large sharks (Block et al., 2011;Hammerschlag et al., 2011). Tagging methods (e.g., acoustic telemetry and archival tagging) have been used to discern habitat use and associated data provide fisheries managers with a more complete picture of the environmental and oceanographic conditions that influence movement and can be used in ecosystem-based fisheries management. Unfortunately, due to cost and logistical constraints, basic data on habitat use and movement are limited for many shark species, including many important stocks in the North Atlantic Ocean and associated basins, such as the Gulf of Mexico (GOM).
The scalloped hammerhead (Sphyrna lewini) is a circumglobal shark species that lives in warm temperate and tropical seas (Compagno, 1984). This species is highly mobile, capable of extensive migrations, and is generally found across a mosaic of ecosystems including bays and estuaries, continental shelves, and offshore pelagic waters (Compagno, 1984;Spaet et al., 2017). Current concerns over the population status of scalloped hammerheads include temporal declines in abundance (Baum and Blanchard, 2010;Powers et al., 2013) and high atvessel mortality rates (Morgan and Burgess, 2007). Scalloped hammerheads are primarily caught in targeted commercial shark fisheries, as well as bycatch in longline and coastal gillnet fisheries with mortality estimates sometimes exceeding 90% of individuals landed (Morgan and Burgess, 2007). Scalloped hammerheads are also captured in land-based shark fisheries in the GOM (Ajemian et al., 2016), and exhibit increased physiological sensitivity to capture stress compared to other pelagic sharks (Gallagher et al., 2014). This species has a high risk of overexploitation due to a late age of maturity and relatively slow growth (Cortes, 2000;Piercy et al., 2007). In addition, schooling behavior by this species makes them particularly vulnerable to capture in relatively large numbers (Hayes et al., 2009). Scalloped hammerhead fins are highly valued in the shark fin trade (Abercrombie et al., 2005) with estimates ranging between 1.3 and 2.7 million scalloped and smooth hammerheads (S. zygaena) used in the Hong Kong fin trade annually (Clark et al., 2006;Camhi et al., 2009). Following a petition to list scalloped hammerheads on the Endangered Species List in 2011, an extinction risk analysis (ERA) team provided a status review and assessed the risk of extinction (Miller et al., 2014), resulting in this species becoming the first shark to be protected under the U.S. Endangered Species Act for populations in the eastern Atlantic and Pacific Ocean. The risk of extinction for scalloped hammerheads in the GOM and northwest (NW) Atlantic Ocean was classified as low to moderate, but this population has suffered an 83% decline since 1981 due to commercial and recreational fishing (Hayes et al., 2009). As a result, the ERA suggested that the low levels of abundance combined with susceptible life history characteristics make this species particularly at risk for local depletion.
Population structure and connectivity studies on scalloped hammerheads have found variable patterns. Evidence suggests substantial genetic heterogeneity for a large marine predator on a global scale  and multiple subpopulations and cryptic speciation in the western North Atlantic Ocean (Quattro et al., 2006;Chapman et al., 2009). Tagging studies demonstrate that scalloped hammerheads have the capacity to travel distances exceeding 1,000 km (Kohler and Turner, 2001;Bessudo et al., 2011;Spaet et al., 2017). However, several studies have found that scalloped hammerheads have relatively small home ranges and display aggregative behaviors around bathymetric features, such as seamounts, potentially due to enhanced prey availability, geomagnetic associations from gradients in the seafloor, or simply short durations of tagged individuals (Klimley, 1993;Hearn et al., 2010). Other studies have observed ontogenetic differences in habitat use; for example, Duncan and Holland (2006) found juvenile scalloped hammerheads remained within a bay system for up to their first year of life before moving to offshore coastal areas. Consequently, a better understanding of shark movement patterns in relation to environmental conditions within the GOM is critical for proper conservation and successful management of this species.
The GOM is a highly productive marginal sea with unique oceanographic conditions influencing the distribution and abundance of sharks as well as other large marine predators (Richards et al., 1989;McKinney et al., 2012). Significant nutrient loading from the Mississippi River coupled with oceanic conditions (due to the proximity of slope waters to shore) results in optimal environments for many forage species and consequently large predators (Chesney et al., 2000). In addition, the basin is dominated by a unique physical feature, the loop current (LC), which influences the productivity and physicochemical parameters within the region (Biggs, 1992). The LC consists of warm surface waters from the Caribbean Sea entering the GOM through the Yucatan Strait, and flows clockwise around a surface "dome" of water before exiting through the Straits of Florida (Oey et al., 2005). Portions of the LC often spin off, forming cold and warm core eddies that drift into the NW GOM (Oey et al., 2005). As a result, this physical feature, combined with the Mississippi River plume, essentially divides the GOM into two unique oceanographic regions that influence habitat use patterns of pelagic fishes (Richards et al., 1989;Kraus and Rooker, 2007;Teo et al., 2007). Undoubtedly, differences in habitat use between the two regions is critical data that needs to be incorporated into habitat standardization procedures and it is therefore important to investigate the east-west barrier that the LC may present to mixing within the GOM, which may lead to differences in stock structure of sharks such as the scalloped hammerhead.
The objectives of this study were to characterize movement dynamics of scalloped hammerheads throughout the GOM. We were specifically interested in quantifying horizontal movements (i.e., space use), characterizing habitat associations in relation to environmental parameters, and developing predictive models to identify habitat suitability for scalloped hammerheads in the GOM based on their movement behaviors (i.e., area restricted, transiting).

Study Area and Tagging
Tagging activities occurred throughout the northern GOM including the northeast (NE) off Alabama and Florida, northcentral (NC) off Louisiana, NW (NW) off Texas, United States. Tagging procedures occurred through hook and line fishing from sportfishing vessels or through fishery independent bottom longline (BLL) sampling. For hook and line fishing, once caught using either live or dead bait, the animal was immediately brought alongside the vessel and retained in the water for tagging procedures while some individuals were brought onboard with a cradle during the longline sampling. Scalloped hammerheads captured and tagged during fishery-independent BLL sampling followed methods outlined in Drymon et al. (2010). Briefly, a BLL consisting of 1.85 km of 4-mm monofilament (545-kg test) mainline was set with 100 gangions. Gangions were 3.66 m (3mm monofilament) in length and consisted of a longline snap and a 15/0 circle hook (Mustad model 39960D) baited with Atlantic mackerel (Scomber scombrus). All bottom longline BLLs were set for 1 h.
During tagging, two to four small holes were drilled into the distal portion of the leading edge of the dorsal fin for placement of a smart position or temperature (SPOT) transmitter tag (Wildlife Computers, Seattle, WA, United States). Tags were coated in black anti-fouling paint (Micron 66, Interlux R , NJ, United States) to prevent growth of fouling organisms, and sharks tagged in Texas were equipped with copper sensors and an extrastiff antenna. Transmitters were programmed with a maximum allotment up to 250 transmissions per day. Shark tracks were filtered to remove poor location class (Z) positions. For days in which multiple positions were recorded for an individual shark, the highest quality location class position for each day was retained. The SPOT tag was positioned on the dorsal fin to facilitate near real-time tracking via Argos satellites during surface interval periods. Prior to release, the fork length (FL, cm) was measured and sex of each shark were recorded.

Space Use
Activity space use and site fidelity estimates were calculated in the rhr package (Signer and Balkenhol, 2015) in R (R Development Core Team, 2011). Kernel density estimation (KDE) was used to quantify the activity space of individual sharks (n = 26; sharks with fewer than 20 total detections were excluded due to limited data) using reference bandwidth estimation. We defined areas of core use as those constrained by 50 % KDE isopleths (TinHan et al., 2018). Two movement metrics [mean squared distance (MSD) from center of activity and linearity index (LI)] were calculated for individual animal tracks and for 100 random trajectories. MSD from center of activity was calculated as a metric of dispersion of positions around the centroid of the track. LI was calculated as the linear distance between the first and last position of an individual track divided by the sum of distances between each position (total distance traveled), where values approaching 1 indicate linear paths and values approaching 0 indicate tortuous paths. Animal tracks significantly more constrained (MSD) or less linear (LI) (p < 0.05) than the distributions of MSD and LI values from randomly generated trajectories were considered indicative of site fidelity (Spencer et al., 1990). The relationship between core area size and shark size (FL), and core area size and sex were tested with a linear regression and Kruskal-Wallis one-way ANOVA, respectively.

Environmental Data
A suite of environmental variables were used to relate to scalloped hammerhead occurrence including chlorophyll a concentration (chl a), salinity, sea surface height anomaly (SSHA), sea surface temperature (SST), bathymetric slope, depth, distance to artificial structure, and distance to hard-bottom habitat. Oceanographic data (chl a, salinity, SSHA, SST) were extracted over the duration of the study period (19 September 2012 to 19 July 2016) using the Marine Geospatial Ecology Toolbox (v. 0.8a68, Roberts et al., 2010) in ArcMap 10.4.1 (ESRI). Mean 8-day and monthly chl a data were obtained from the Moderate Resolution Imagine Spectroradiometer (MODIS) at a 4 km resolution. For grid cells obscured by cloud cover in the 8-day data, monthly mean values of chl a were used. Daily mean SST (4 km resolution) and salinity (7.8 km resolution) were obtained from the Hybrid Coordinate Ocean Model (HYCOM) data server. SSHA (24 km resolution) data were obtained from the Aviso altimetry dataset, and the Copernicus Marine Environment Monitoring Service (CMEMS).
Bathymetry data (0.003 km resolution) for the GOM were downloaded from the Scripps Institution of Oceanography 1 (Becker et al., 2009). Kriging interpolation methods were used to create a continuous raster surface to estimate benthic sediment composition (% rock, gravel, sand, mud) using sediment data from the United States Geological Survey (USGS 2 ) usSeabed data set (Buczkowski et al., 2006). Hard-bottom habitat was then classified as areas comprised of greater than 25% rock substrate. Locations of standing oil and gas platforms were obtained from the Bureau of Ocean Energy Management (BOEM), while the locations of known artificial reefs were obtained from data sets compiled by both BOEM and National Oceanic and Atmospheric Administration (NOAA) (Supplementary Figure S1). For the purpose of our models, these two data sets (oil and gas platforms as well as artificial reefs) were combined and collectively referred to as artificial habitat. Distance to the nearest artificial and hardbottom habitat was calculated as the shortest in-water distance from a sampling location using the Cost Distance tool in ArcGIS.

Habitat Modeling
A generalized additive modeling (GAM) framework similar to Aarts et al. (2008) was used to examine species-habitat relationships and to predict the distribution of suitable habitat in the northern GOM. Daily positions of sharks represented a sample of telemetry positions visited by each shark (presence). Each presence was complemented by a sample position representing a location that the shark did not visit (absence) randomly drawn from a surface of accessibility based on the elapsed time between positions and the estimated maximum sustained swimming speed of scalloped hammerheads (∼1 body length s −1 ; Lowe, 1996;Ryan et al., 2015). This process ensured that each absence was in a location that could be realistically reached by scalloped hammerheads (based on the previous position). In addition, because each absence is matched temporally to a telemetry position, we can assume that the shark did not visit that position in space at that time (Aarts et al., 2008). A binomial GAM with individual included as a random factor was then used to evaluate the effects of the eight environmental variables (chl a, salinity, SSHA, SST, depth, bathymetric slope, distance to artificial structure, and distance to hard-bottom habitat) on the presence of scalloped hammerheads. Binomial GAMs were fit in the mgcv package (Wood, 2006) in R, with a logit link function and cubic regression splines restricted to a maximum of four degrees of freedom in order to avoid overfitting (Ciannelli et al., 2008). A manual stepwise backward selection procedure based on minimizing Akaike's information criterion (AIC; Akaike, 1974) was used to select variables explaining the presence of scalloped hammerheads. Non-significant variables (p-value < 0.05) were removed from the final model if variable exclusion resulted in a AIC < 2 (Burnham and Anderson, 2002). Prior to backward selection, all variables were tested for collinearity using variance inflation factor (VIF), and variables with VIF > 2 were removed from analysis. There was no collinearity detected among the variables selected. The relative importance of each predictor variable in the final model was assessed by examining the difference in AIC ( AIC) and percent deviance explained ( DE) when each variable was removed from the final model (Rooker et al., 2012).
The best-fit model was then used to assess the spatial distribution of suitable habitat for scalloped hammerheads within the northern GOM. Oceanographic and habitat data for model predictions were obtained in the same manner as previously described for telemetry data, with the exception that oceanographic data (chl a, salinity, SSHA, and SST) were based on cumulative means across summer months over a single year (01 July 2016 to 30 September 2016). Environmental data were extracted to prediction grids (resolution = 0.05 • ) and used to predict the probability of occurrence of scalloped hammerheads at each grid point using the predict.gam function in R. For visualization, grid points were kriged in ArcGIS to create a continuous raster surface and smoothed using bilinear interpolation. For the purposes of this study, habitat suitability was then classified based on probability of occurrence as: high (>0.66), moderate (0.33-0.66), and low (<0.33).

State-Space Modeling
A first-difference correlated random walk switching (DCRWS) model was used to initially identify whether scalloped hammerheads exhibited multiple modes of movement behavior (area restricted, unclassified, transiting) over the course of a track. Bayesian DCRWS models permit estimation of states (i.e., biological parameters, including discrete behavioral modes) at regular intervals throughout potentially irregular time series data, and the probability of switching between behavioral modes at each interval, while making allowances for observational errors such as those inherent in satellite tag data. Models were run on individual sharks with at least 20 location estimates (n = 26) using JAGS 4.0.0 (Just Another Gibbs Sampler) via package bsam in R. Markov chain Monte Carlo (MCMC) techniques were used to produce mean estimates of location ( latitude , longitude ), and random walk index (turning angle and speed). Random walk index estimates were then used to estimate behavioral mode (b t ) where values range from 1 to 2. Values of b t approaching 1 indicate area-restricted movement, and values approaching 2 indicate transiting movement. Continuous estimates of b t were classified into discrete behavioral modes where values < 1.25 were considered as area-restricted movements and values > 1.75 were considered as transiting movements (Jonsen et al., 2007). Where switches among multiple behavioral modes could not be identified, a first difference correlated random walk (DCRW) model was fit to all shark tracks. To characterize the movements of tagged scalloped hammerheads, random walk index γ was estimated in the DCRW model, where 0 < γ < 1. Similar to behavioral mode classification, values of γ approaching 0 indicate uncorrelated turning angle and speed at each location (simple random walk), while values approaching 1 indicate more linear tracks, with autocorrelation in turning angle and speed (correlated random walk). DCRWS and DCRW models regularized samples to 12 h time intervals for each animal, and were fit using 45,000 MCMC samples, with a burn in of 5,000, and thinned by a factor of 4 (final MCMC length: 10,000 samples). The Raftery-Lewis diagnostic (Raftery and Lewis, 1992) was used to estimate minimum MCMC lengths and thinning factors to precisely estimate posterior sample quantiles, and Gelman-Rubin's diagnostic (Gelman and Rubin, 1992) was used to confirm model convergence.

RESULTS
A total of 33 scalloped hammerheads were captured and tagged throughout the northern GOM and consisted of 11 in the NE, 17 in the NC, and five individuals in the NW GOM (Table 1). Eighty five percent (n = 28) of tagged sharks were males and 15% (n = 5) were females. Mean size at the time of tagging was 159.4 cm FL (±30.6 cm SD) ranging from 102 to 220 cm FL. Tagging activities occurred over a 4-year period from 2012 to 2015 with a mean number of days at liberty of 146 days (±139.7 SD) ranging from 5 to 479 total days until the tag stopped communicating. Four tags communicated over one full year, five additional tags for over six months, and seven more tags for over four months with a total of 4,543 days of communication among the 33 total tagged sharks ( Table 1).

Space Use
Movement patterns of scalloped hammerheads tagged over the northern GOM continental shelf showed a wide range of movement tracks with limited long-distance dispersal (Figures 1A,B). Core use area (50% KDE) sizes for scalloped hammerheads were relatively small (9,159 ± 10,735 km 2 ; mean ± SD; Figure 1A and Table 2). MSDs to centroids and LIs indicated site fidelity in nearly half (MSD: 65%, LI: 50%) of tagged scalloped hammerheads ( Table 2). Core area size was not significantly related to sex (Kruskal-Wallis test; H = 0.852, p = 0.356) or shark size (F 1,24 = 0.138, p = 0.714). Though no statistically significant relationship was found between sex and horizontal distribution, male shark detections were located across the continental shelf primarily in depths less than 200 m, whereas detections of female sharks were concentrated at the 200 m isobaths and deeper ( Figure 1B).

Habitat Modeling
The final model describing the influence of environmental variables on the presence of scalloped hammerheads in the northern GOM had an AIC of 3,153 and DE of 33.5%. Seven of the eight environmental variables were retained in the final model: chl a, salinity, SSHA, depth, bathymetric slope, distance to artificial structure, and distance to hard-bottom habitat. The most influential environmental variables according to AIC were depth ( AIC = 278.8), distance to artificial habitat ( AIC = 265.1), and chl a concentration ( AIC = 55.8; Supplementary Table S1). The same three variables were most influential according to DE. Response plots indicated that scalloped hammerhead presence increased at depths shallower than 1,500 m (Figure 2) in close proximity to both artificial and hard-bottom habitat. The presence of scalloped hammerheads was also greatest in waters with low chl a concentrations (∼0-4 mg m −3 ), moderate salinities (33-35.5), and neutral SSHA (−0.05-0 cm; Figure 2). Habitat suitability for scalloped hammerheads was predicted to be very high on the mid to outer continental shelf of the GOM inside the 200 m isobath (Figure 3).
Habitat suitability decreased in most nearshore areas (moderate to suitable) with the exception of the NC GOM offshore of Mobile Bay, Alabama to Cape San Blas, Florida, and south Texas, where highly suitable habitat spanned the continental shelf.

State-Space Modeling
First difference correlated random walk switching modeling identified behavioral switches in two scalloped hammerheads ( Fish ID 4 and 28, Figures 4, 5). Both sharks spent the majority of the time making area-restricted movements centralized in two patches separated by a distance of ∼350 km (Figures 5A,B). Both sharks exhibited a small number of switches from area-restricted movement to transiting movement when moving between these patches and across the Mississippi River plume. Another individual scalloped hammerhead (Fish ID 12, Figure 5C) made   directed movements (b t < 1.5) toward and along the shoreline of the Florida Panhandle for ∼7 days after tagging, before making less directionally persistent movements (b t > 1.5) over an area of 1,280 km 2 for ∼30 days. With the exception of two scalloped hammerheads (Fish ID 4 and 28), values of b t for tagged scalloped hammerheads did not meet criteria for classification into distinct behavioral modes (e.g., Fish ID 12). Non-switching (DCRW) models run for individual scalloped hammerheads did not converge for three sharks (Fish ID 2,22,27). Mean correlated random walk (γ) estimates from DCRW models on scalloped hammerheads indicated sharks typically followed uncorrelated random walk (highly tortuous, area-restricted) movements (γ = 0.08 ± 0.09 mean ± SD). No relationship was found between shark sex, size (FL), or random walk index.

DISCUSSION
Knowledge of movement patterns of marine apex predators is essential when attempting to understand connectivity within and across ocean basins for conservation and management efforts. In our study, scalloped hammerhead distribution was primarily restricted to the continental shelf, though track data revealed some cross-shelf movements.  determined population genetic structure of scalloped hammerheads was substantial for a large marine predator, and that restricted dispersal over deep open ocean habitat may be the primary driver, but found high connectivity and dispersal across continuous coastlines along continental shelves. Tagging studies have found similar overall patterns of limited long-distance movement for this species (Kohler and Turner, 2001); however, individuals are  capable of extensive movements (>1,000 km; Kohler and Turner, 2001;Bessudo et al., 2011). Scalloped hammerhead movements were monitored for just over a single year (479 d maximum, Table 1) making it difficult to assess long-term movement in this study. There was little indication of seasonal patterns in habitat use as demonstrated for other species of elasmobranchs in the northern GOM (Drymon et al., 2010;Bethea et al., 2015). Nevertheless, results show restricted movement across the continental shelf with limited evidence of open ocean habitat use beyond the shelf edge. Though only five females were tagged in this study, evidence suggests these individuals more frequently use shelf-edge (>200 m isobath) areas in contrast to shallower depths along the mid-shelf used by males. While this is in stark contrast to the presence of periodically high abundances of neonates in multiple coastal areas across the GOM (Hueter and Tyminski, 2007;Bethea et al., 2015), additional tagging studies that target mature females may better reveal inshore-offshore habitat connectivity as well as potential areas of parturition. Sex-specific movement patterns and segregation are relatively common in elasmobranchs (Sims, 2005), and sexual segregation has been reported for scalloped hammerheads in the GOM (Branstetter, 1987). Klimley (1987) found females in the Gulf of California moved further offshore at a smaller size than males, and provided evidence that this strategy may be to rapidly increase in size to reach sexual maturity due to enhanced prey availability. Size and sex segregated populations may be disproportionately affected by interactions with fisheries operating in specific habitats or regions. In the GOM, the majority of reported fishing mortality of scalloped hammerheads occurs in coastal-shelf BLL fisheries, with comparatively few scalloped hammerheads being taken in other sectors (e.g., pelagic longline; Miller et al., 2014). If scalloped hammerheads that most frequently interact with BLL sets in coastal waters are presumed to be males or sexually immature, GOM populations may not necessarily face the potential risk of mate limitation resulting from the loss of larger females.
Habitat preferences of scalloped hammerheads appeared to be primarily driven by bathymetric features such as depth and bottom type, rather than dynamic oceanographic processes (e.g., SST and SSHA). The preference for mid-shelf or shelf edge habitats (<200 m depth) is consistent with aforementioned fisheries dependent data that indicate this species is rarely captured in oceanic waters of the GOM (Miller et al., 2014). The positive association of tagged scalloped hammerheads with natural hard-bottom habitat in the current study is not surprising given that large schools are observed annually at mid-shelf and shelf edge banks (e.g., Flower Garden Banks) in the northern GOM (Childs, 2001). Bathymetric features such as seamounts and associated currents are thought to enhance prey availability to scalloped hammerheads (Klimley, 1993;Hearn et al., 2010;Ketchum et al., 2014a), and thus natural reefs and hard-bottom outcroppings may be important habitats or foraging areas for this species in the GOM. Still, scalloped hammerheads in the current study were more closely linked to artificial habitats, particularly (standing and toppled) oil and gas platforms. Oil and gas platforms are productive habitats (Claisse et al., 2014) that provide structure throughout the water column, and often hold higher fish densities compared to natural habitats (Streich et al., 2017). Given the high fish densities and number of such structures in the northern GOM (>2,000), it is likely that large predators such as scalloped hammerheads frequent oil and gas platforms to increase foraging opportunities. While habitat use in this study was primarily associated with bathymetric features, scalloped hammerheads were also negatively associated with high chlorophyll a, suggesting that this species may avoid nearshore areas of high nutrient input (e.g., Mississippi River plume, low salinity estuaries). This finding was supported by movement patterns from our switching models that suggest sharks avoided surface waters of the Mississippi River plume. Likewise, our habitat suitability maps indicated that habitat suitability was highest from the mid to outer shelf and over nearshore areas with relatively low freshwater inflow (NW Florida, south Texas). While other oceanographic variables contributed very little to our final models, it remains a possibility that oceanographic drivers may yet influence the movements of female, and larger, scalloped hammerheads, which are presumed to spend a greater proportion of time in off-shelf pelagic habitats with fewer bathymetric features.
The majority of tagged scalloped hammerheads displayed area-restricted movements, with limited evidence of behavioral switches to transiting movements. Several potential factors have been suggested to influence the site fidelity, home range size, and proportion of time a predator spends engaged in area-restricted searching versus transiting movements, including proximate cues such as prey abundance and inter-specific competition (Kittle et al., 2016), and ultimate drivers such as reproductive success (Patrick and Weimerskirch, 2017). The DCRWS modeling identified behavioral switches between transiting and arearestricted movements in two scalloped hammerheads. Though these switches did not appear to correspond to abiotic factors investigated, behavioral switches for both scalloped hammerheads were observed in proximity to the Mississippi River plume suggesting this feature may not act as a barrier for longitudinal movement of this species. The median detection interval for both sharks was less than 2 h, while movements across the Mississippi River plume were associated with gaps in detection greater than 5 days, similar to a lack of detections found around this plume for great hammerheads (Sphyrna mokarran; Drymon and Wells, 2017). Vertical stratification of the water column is common where the Mississippi River meets the GOM (Androulidakis et al., 2015), and prolonged intervals observed between SPOT transmissions in this area suggest that sharks avoided surface waters when traversing the plume, potentially in response to vertical gradients in environmental conditions (e.g., salinity, temperature, and dissolved oxygen) or prey availability. At the Mississippi River mouth, seasonal hypoxia begins forming in the spring and persists throughout the summer (Bianchi et al., 2010). Shark ID 4 exhibited transiting behavior across the "dead zone" regions during the summer when hypoxia is most severe, possibly suggesting avoidance of low oxygen water. Alternatively, as hypoxia can aggregate prey near the edges of the dead zone (Craig, 2012), it is possible the sharks are experiencing enhanced foraging on stressed prey in hypoxic environments, as demonstrated in teleost fish (Pihl et al., 1992;Long and Seitz, 2008). A satellite tagged scalloped hammerhead in the Gulf of California made deep dives exceeding 250 m into severely hypoxic water that was potentially related to pursuing deepwater squid prey (Jorgensen et al., 2009). Thus, if scalloped hammerheads are able to tolerate low oxygen conditions, they may experience an energetic tradeoff with enhanced prey consumption. Given the rate (∼3 km/h) at which sharks traversed the Mississippi River plume (when movements were classified as transiting behavior), it is likely that the DCRWS model did not identify true behavioral switches, but rather sharks simply continued to make area-restricted movements while tags were submerged and unable to transmit positional data. Alternatively, it is possible these sharks made transiting movements that were not identified by DCRWS models. Nevertheless, behavioral state cannot be accurately inferred for movements made over spatiotemporal scales smaller than those detectable by the tags in this study, due to positional error of tags (poor location classes), or tags being submerged.
Performance based on tag type can be species-specific and this study found fin-mounted SPOT tags appear to be a viable option to obtain daily position estimates of scalloped hammerheads, owing to what is likely a considerable amount of time spent at the surface. In contrast, a recent study by Drymon and Wells (2017) used a double tagging approach on great hammerheads in the GOM and found fin-mounted SPOT tags were a poor choice, despite known use of surface waters by this species (Queiroz et al., 2016). A similar study on juvenile smooth hammerheads (Sphyrna zygaena) also had limited success using fin-mounted SPOT tags in New Zealand (Francis, 2016). Several studies examining movement of scalloped hammerheads have used different tag technology including acoustic telemetry (Klimley, 1993;Bessudo et al., 2011;Ketchum et al., 2014a,b) and PAT tags (Jorgensen et al., 2009;Hoffmayer et al., 2013). The goal of this study was not to evaluate tag type; however, our results propose the utility of fin-mounted SPOT tags to provide information on habitat use for scalloped hammerheads in this region of the GOM.
Results of this study provide valuable movement information needed to properly conserve this species and to guide fishery management. Similar findings for the great hammerhead were found by Guttridge et al. (2017) highlighting the need for local conservation and management due to the seasonal residency and site fidelity of this species in the southeastern United States and Bahamas. A recent ERA (Miller et al., 2014) evaluated the six distinct population segments of scalloped hammerheads and determined that the NW Atlantic and GOM stock as well as the Central Pacific stock have declined to relatively low levels of abundance but were at low risk of current extinction. In contrast, the Central and Southwest Atlantic stock and the Indo-West Pacific stock were at moderate risk of extinction and the Eastern Atlantic and Eastern Pacific stocks were at high risks of extinction. Findings from this study provide important information on movement of this species in the GOM and highlight the restricted use of continental shelf habitat and resident behavior that will need to be incorporated in future stock assessments and extinction risk analyses. Areas where concentrated fishing effort is relatively high may have significant impacts on this species given the sex-specific differences in habitat use and deserves careful consideration for future conservation measures.

ETHICS STATEMENT
Sharks tagged from 2014 and beyond were approved by the Institutional Animal Care and Use Committee (IACUC) reference number 009413 at Texas A&M University. Sharks tagged prior to 2014 did not require approval in accordance with the Louisiana Department of Wildlife and Fisheries (LDWF) and National Marine Fisheries Service (NMFS).

AUTHOR CONTRIBUTIONS
All authors contributed to the design and implementation of the research. RW, TT, and MD contributed to the analysis of the results. RW led writing of the manuscript with assistance from all authors.

FUNDING
The open access publishing fees for this article have been covered by the Texas A&M University Open Access to Knowledge Fund (OAKFund), supported by the University Libraries and the Office of the Vice President for Research. and % deviance explained (% DE). The relative importance of each predictor variable to the final model is given as the difference in model AIC ( AIC) and % DE ( DE) if this variable is removed from the final model.