Defining Sex-Specific Habitat Suitability for a Northern Gulf of Mexico Shark Assemblage

Understanding the factors that influence species’ distributions is crucial for implementing effective management and conservation practices, yet difficult for highly vagile species like sharks. Many shark species demonstrate either spatial and/or temporal sexual segregation, further confounding accurate quantification of habitat suitability. Given the importance of understanding spatiotemporal patterns in the distribution of coastal shark assemblages, we sought to quantify sex-specific abiotic factors that influence seasonal variation in a coastal shark assemblage using data from a long-term fisheries-independent bottom longline program in the northern Gulf of Mexico. Catch data (individuals/100 hooks/hour) were coupled with a suite of potentially predictive variables: surface and bottom values for temperature (°C) and salinity (psu), sea surface height (m), three-dimensional surface and bottom current velocity (u, v, w, in m/s), bottom dissolved oxygen (mg/l), depth (m), substrate grain size (mm), daylength (min), and distance from shore (km). Data were analyzed using boosted regression trees (BRT) to describe the relationships between catch data and environmental factors potentially influencing sex-specific species distribution and abundance. Between May 2006 and November 2018, we conducted 1,226 bottom longline sets and caught 13,742 individuals encompassing 67 species. The majority of the animals captured (74%) were elasmobranchs, primarily sharks. Two species from each of the following three categories were selected for further analyses: small coastal sharks (Atlantic sharpnose shark Rhizoprionodon terraenovae and blacknose shark Carcharhinus acronotus), large coastal sharks (blacktip shark C. limbatus and sandbar shark C. plumbeus), and shelf-associated sharks (smoothhound sharks, Mustelus spp. and scalloped hammerhead Sphyrna lewini). Depth and distance from shore were the strongest predictors of distribution and relative abundance, followed by longitude and bottom salinity; other factors (e.g., temperature, daylength, substrate grain size) were less predictive. For the six species examined, predictive factors were often the same for males and females, although the range of preferred values varied. Surprisingly, the importance of these predictors varied little across seasons. Collectively, our findings demonstrate that sexual segregation is the norm for sharks in the north-central Gulf of Mexico. Long-term fishery-independent monitoring to further quantify these sex-based differences in habitat use should be prioritized, particularly in light of impending climate change.


INTRODUCTION
Understanding species distributions is central to realizing their ecological roles, fundamental for developing effective conservation and management measures, and increasingly necessary in the face of imminent climate change. Examining the mechanisms underpinning species distributions can help identify critical habitats and aid in the development of habitat suitability models. For some species, habitat associations shift ontogenetically and may be sex-specific; for these species, collecting data to accurately characterize spatially and temporally shifting distributions presents a challenge, particularly given changing environmental conditions. Sharks are highly mobile species whose distributions can be affected by a number of abiotic and environmental variables. Abiotic factors influencing shark distributions are complex, interactive, species-specific, and geographically variable (Schlaff et al., 2014). Perhaps the most common factor influencing the distribution of sharks is temperature. For example, the timing of massive seasonal migrations of blacktip shark (Carcharhinus limbatus) along the southeast coast of Florida is best predicted by water temperature, where large aggregations are only seen when water temperatures are below 25 • C (Kajiura and Tellman, 2016). Similarly, residency for migratory adult lemon sharks (Negaprion brevirostris) in southeast Florida coincides with water temperatures less than 24 • C (Kessel et al., 2014). While the influence of temperature, as well as salinity (e.g., Ubeda et al., 2009), dissolved oxygen (e.g., Heithaus et al., 2009), and interactions thereof (e.g., Drymon et al., 2014) are well known, other abiotic factors such as daylength (Grubbs et al., 2005) and distance from shore (Haig et al., 2018) may be important, and warrant further consideration.
Characterizing shark distributions is further complicated by tendencies toward sex-specific segregation (Springer, 1967;Sims, 2005;Wearmouth and Sims, 2008). A common type of sexspecific segregation is geographical segregation (Backus et al., 1956), defined as differences in habitat selection between sexes. This type of segregation can vary temporally; for example, females of many shark species move into discrete areas specifically for parturition (e.g., Feldheim et al., 2002). For other species, femaleonly aggregations have been recorded in shallow lagoons where water temperatures are 1-2 • C warmer than surrounding areas (Economakis and Lobel, 1998). According to the thermal-niche hypothesis (Sims, 2005), thermal advantages, including shorter gestation and/or larger size at birth, may be conferred to pregnant females during these periods (Jirik and Lowe, 2012). In some shark fisheries, sexual segregation may lead to disproportionate removals of males or females. This has been shown for male shortfin mako (Isurus oxyrinchus) in the South Pacific Ocean (Mucientes et al., 2009) and female spiny dogfish (Squalus acanthias) in the northwest Atlantic Ocean (Rago and Sosebee, 2010). For these species, understanding the factors determining sex-specific distributions is useful for reducing harvest pressure on the over-exploited portion of the population (Dell'Apa et al., 2017;Haugen et al., 2017).
The northern Gulf of Mexico (GoM) supports a diverse assemblage of coastal sharks, and while the environmental factors underpinning the distribution within this assemblage are welldescribed (e.g., Drymon et al., 2013), a detailed understanding of seasonal and sex-specific drivers within this assemblage is lacking. For example, the interactive effects of temperature, salinity, and proximity to tidal inlets influence habitat use for juvenile bull (Carcharhinus leucas), blacktip, and bonnethead (Sphyrna tiburo) sharks in coastal Texas (Froeschke et al., 2010;Plumlee et al., 2018). Similarly, temperature, depth, and salinity predict habitat use for juvenile blacktip, bonnethead, finetooth (C. isodon), spinner (C. brevipinna), scalloped hammerhead (S. lewini), and Atlantic sharpnose (Rhizoprionodon terraenovae) sharks in bay systems in north Florida (Ward-Paige et al., 2014). While these studies characterize estuarine habitat use, many coastal sharks in the GoM demonstrate across-shelf movements from coastal estuaries to offshore waters. Quantifying inshore and offshore habitat use and how that may differ between males and females, is therefore a critical data gap. Given the importance of understanding spatiotemporal patterns in the distribution of coastal shark assemblages, as well as the influential mechanisms, the goals of the current study are to (1) measure seasonal variation in the northern GoM coastal shark assemblage, specifically movements into deeper water, and (2) quantify sexspecific abiotic factors that influence seasonal variation in the coastal shark assemblage.

Bottom Longline Survey
Catch data were collected as part of inshore (Drymon et al., 2010) and offshore (Powers et al., 2018) fisheriesindependent bottom longline surveys conducted off the coasts of Louisiana, Mississippi and Alabama (Figure 1). For both surveys, bottom longline locations were selected using a stratifiedrandom sampling design. Bottom longline sampling followed standardized methods described in Drymon et al. (2013). The main line consisted of 1.85 km (1 nmi) of 4 mm monofilament (545 kg test) that was set with 100 gangions. Gangions consisted of a longline snap and a 15/0 circle hook baited with Atlantic mackerel (Scomber scombrus). Each gangion was made of 3.66 m of 3 mm monofilament (320 kg test). All sets were soaked for 1 h. Once the bottom longline was retrieved, all sharks that could be boated safely were removed from the main line, unhooked, and identified to species following Castro (2010). For each individual, length (precaudal, fork, and stretch in cm), weight (in kg), sex and maturity stage (when possible) were recorded. Maturity in males was assessed following Clark and Von Schmidt (1965). Sharks were tagged either on the first dorsal fin with a plastic Rototag (Premier 1) or just below the first dorsal fin with a metal dart tag. The tag type and location depended both on species and size of shark at capture (Kohler and Turner, 2001). Catch data were converted to catch per unit effort (CPUE), expressed as individuals 100 hooks −1 hour −1 .

Environmental Data
Multiple environmental variables were examined to explain the distribution of species captured on the bottom longline during the course of this study. Temperature and salinity were explored because they have been shown to commonly influence shark distributions (Schlaff et al., 2014), particularly for coastal species (e.g., blacktip shark). Sea surface height and current velocities were examined because they are oceanographic predictors that may aid in predicting habitat suitability for species that occur primarily offshore (e.g., scalloped hammerhead). Since the species assemblage in this region included both coastal and offshore species, the variables depth and distance from shore were also considered as potential predictors. Given the distribution of sampling effort across multiple seasons (i.e., variable photoperiods), daylength was also examined as a potential predictive variable following Grubbs et al. (2005). Finally, dissolved oxygen and substrate grain size were explored because these variables are potentially relevant to this study area given episodic, large-scale hypoxia (Rabalais et al., 2002) and transitions in sediment type (Balsam and Beeson, 2003).
Data for these environmental variables were obtained from multiple sources. Surface and bottom temperature ( • C), surface and bottom salinity (psu), sea surface height (m), and threedimensional surface and bottom current velocity (u, v, w in m s −1 ) were obtained from the Hybrid Coordinate Ocean Model (HYCOM) data server (4 km resolution). Bottom dissolved oxygen (mg l −1 ) was obtained from the National Oceanic and Atmospheric Administration (NOAA) 1 , and interpolated across ∼100 to 250 survey stations (number varied by year). Depth (m) and substrate grain size (mm) were obtained from the United States Geological Survey (USGS) 2 , 0.33 arc seconds, ∼10 m resolution). Daylength (min) was calculated in R (R Core Team, 2013) using code by SD 3 ; distance from shore (km) was calculated in QGIS (Quantum GIS Development Team, 2019).

Habitat Modeling
Boosted regression trees (BRTs) were used to describe the relationships between sex-specific CPUE from the bottom longline, and environmental variables potentially influencing species' distribution and abundance. Specifically, BRTs were fit for males and females of each species for each of three seasons (meteorological spring, summer, and autumn) resulting in six unique combinations per species. Boosted regression trees are a machine learning technique capable of fitting complex, non-linear relationships that offer predictive advantages over generalized linear or additive models (GLMs and GAMs). In particular, BRTs are relatively insensitive to multicollinearity (Abeare, 2009;Dedman et al., 2015Dedman et al., , 2017, especially in situations where outputs from the BRT are used to predict within the range of sampled data (Dormann et al., 2013). A thorough overview of BRTs is given in Elith et al. (2008).
Preliminary analyses indicated that catch data were zeroinflated (a high proportion of zero values). Zero-inflated catch data are characteristic for many sharks that naturally occur in relatively low abundance (i.e., true zeros, Martin et al., 2005). A detailed treatment of modeling zero-inflated data is given by Maunder and Punt (2004). To properly accommodate the preponderance of zeros, a 2-step (i.e., delta or hurdle) process was chosen to model catch data. Presence/absence probability was modeled using a BRT with a binary distribution and continuous, non-zero (i.e., abundance) probability was modeled using a BRT with a Gaussian distribution. Since the catch data were also long-tailed (some instances of anomalously high catch), nonzero data were natural-log-transformed. Predictions were reverse log-transformed; the final model is a product of the binary and Gaussian BRTs (Lo et al., 1992).
Sixteen variables were considered for the BRT models ( Table 1) and all models were fit using the package gbm.auto (Dedman et al., 2017) in R. Learning rate (lr), bag fraction (bf) and tree contribution (tc) are parameters that are used in concert to achieve minimum predictive error (Elith et al., 2008). These were optimized using gbm.auto for each species/season/sex model run.

Model Performance and Interpretation
The BRT modeling approach automatically partitioned the data into training and testing sets, with the ratio of these dictated by the bag fraction. Ten-fold cross-validation (CV) was then performed, with the members of the training/testing sets randomized each time. Performance metrics captured were training/testing correlation, CV deviance (and standard error (SE)) and correlation (and SE), as well as Area Under Receiver-Operator-Curve (AUC) and its CV and CV SE for the binary models (Parisien and Moritz, 2009). Final fitted functions from the binary BRT were visualized using marginal effect plots to show the effect of a variable on the response after accounting for the average effects of other model variables (Elith et al., 2008).

Habitat Suitability
The distribution of suitable habitat was predicted using the BRTs described above. Environmental data for model predictions were obtained in the same manner as above, except that HYCOM data were extracted for one representative date per season (monthly groupings per season, i.e., DJF/MAM/JJA/SON) at their native resolution of 4 km. Representative dates for environmental data were selected by ranking the absolute value of the differences of all sites' values for all variables against the mean for those variables, then selecting a single date within each season that most closely matched those values. The BRT models later generated predictive CPUE values for each 2 km × 2 km cell, which were then mapped in QGIS (Quantum GIS Development Team, 2019) using the heatmap setting. The heatmap setting produces a color point weighted by the predicted abundance generated from the BRT. The radius of each point was set at the default of 10 mm. Using gbm.auto, we also calculated the coefficient of variation for the predicted abundance values at each 2 km × 2 km cell to represent model error. Given the large differences in predicted abundance between species, the predicted abundance values and the coefficients of variation are not evenly scaled across all seasons, sexes, and species.

Bottom Longline Survey
Between May 2006 and November 2018, we conducted 1,226 bottom longline sets and caught 13,742 individuals across 67 species: 43 species of teleost, 23 species of elasmobranch, and 1 invertebrate (Supplementary Table S1). Survey effort was relatively well distributed across the three seasons sampled: spring (n = 460), summer (n = 405) and autumn (n = 361) (Figure 1). The majority of the individuals captured (74%) were elasmobranchs, primarily sharks. Atlantic sharpnose shark (n = 4,815) was the most common species, followed by blacknose shark (n = 1,532). Together, these two small coastal species represented nearly half of all individuals sampled (46%). The two most abundant large coastal sharks were blacktip shark (n = 1,207) and sandbar shark (C. plumbeus, n = 569). The two most abundant shelf-associated sharks were scalloped hammerhead (n = 209) and smoothhound sharks (Mustelus spp., n = 176). Collectively, individuals of both sexes for these six species occurred during all seasons sampled and were therefore selected for further analyses. Given the difficulty in distinguishing between three morphologically conserved species of Mustelus commonly found in the GoM (smooth dogfish M. canis, Florida smoothhound M. norrisi, and Gulf smoothhound M. sinusmexicanus, Giresi et al., 2015), NOAA Fisheries currently manages all individuals in the genus Mustelus in the GoM as a complex (SEDAR, 2015); therefore, all Mustelus will hereafter be aggregately referred to as "smoothhound sharks."

Model Performance
Model performance was assessed for the two small coastal sharks, the two large coastal sharks, and the two shelf-associated sharks for each sex/season combination (n = 36 model configurations, Table 2). Training data AUC scores ranged from very good (0.84 for Atlantic sharpnose shark, both sexes for all seasons) to excellent (0.99 for smoothhound sharks, both sexes for all seasons) according to criteria defined in Lane et al. (2009). Crossvalidated AUC scores were only slightly lower than training data scores, indicating that model overfitting was negligible (Hijmans and Elith, 2013). For most species/sex/season combinations, distributions were modeled using Gaussian BRTs; for female scalloped hammerheads and male smoothhound sharks, binary BRTs were used, as they had relatively fewer data available.

Small Coastal Sharks
Atlantic sharpnose sharks of both sexes were encountered across a range of similar sizes (Figure 2A) and were captured across the extent of the study area ( Figure 3A). Females were more common than males (1.1:1), which differed from the expected sex ratio (x 2 = 12.95, p < 0.01). Final models for Atlantic sharpnose shark season/sex combinations had training data AUC scores of 0.84, and cross validation AUC scores ranging from 0.78 (±0.01) to 0.79 (±0.01) ( Table 2). For females, distance from shore (47%), depth (13%), and longitude (10%) were the three most influential predictors of distribution, and this was consistent across seasons ( Table 2). Female Atlantic sharpnose sharks were consistently distributed greater than 30 km offshore, in depths ranging from 20 to 80 m, particularly west of 88 • W longitude ( Figure 4A). These trends were consistent across seasons. Wide expanses of suitable habitat were predicted for female Atlantic sharpnose shark, predominantly west of 88 • W. Subtle seasonal shifts in predicted habitat suitability were seen, primarily in deeper waters ( Figure 5A). Coefficients of variation were highest in the deepest portion of the study region (Supplementary Figure S1A). The factors influencing male Atlantic sharpnose sharks were similar to those influencing females, yet the magnitude of these effects was different. The two predictors with the strongest influence for males were longitude and distance from shore, with little seasonal variation ( Table 2). Like females, an east-west gradient was apparent for males, yet was even more pronounced (Figure 4A), such that relatively large areas east of 88 • W longitude were predicted to be unsuitable for male Atlantic sharpnose sharks ( Figure 5A). Coefficients of variation were highest for males in the western region of the study area, yet comparatively lower for males than females (Supplementary Figure S1A).
Similar to Atlantic sharpnose shark, blacknose sharks were captured across a wide range of sizes ( Figure 2B), primarily in nearshore areas ( Figure 3B). Females were encountered less frequently than males (0.77:1), which differed from the expected sex ratio (x 2 = 26.21, p < 0.01). Final models for blacknose shark season/sex combinations had training data AUC scores ranging from 0.85 to 0.87, and cross validation AUC scores ranging from 0.74 (±0.01) to 0.82 (±0.01) ( Table 2). For females, depth was the most influential predictor of abundance (29% across seasons), followed by distance from shore and daylength ( Table 2). Large areas of suitable habitat were identified for female blacknose shark (Figure 5B), and coefficients of variation were highest in the fall (Supplementary Figure S1B).
Male blacknose sharks were encountered in depths shallower than female blacknose sharks ( Figure 4B). Predicted suitable habitat for male blacknose sharks was less widely dispersed across the survey region but showed more seasonality relative to females ( Figure 5B). The primary factors influencing male blacknose sharks were the same as those influencing females (depth and distance from shore), yet also included salinity ( Table 2), with male blacknose sharks encountered primarily  in waters between 30 and 33 psu ( Figure 4B). Consequently, predicted suitable habitat for male blacknose sharks more closely followed the contours of the Chandeleur Islands, Louisiana and the barrier islands south of Mississippi and Alabama ( Figure 5B). Coefficients of variation were relatively low for male blacknose sharks (Supplementary Figure S1B).

Large Coastal Sharks
Blacktip sharks were the most common large coastal shark, with females and males encountered across a wide range of sizes ( Figure 2C), particularly in shallow waters ( Figure 3C). Females were more common than males (1.36:1), which differed from the expected sex ratio (x 2 = 26.51, p < 0.01). Final models for blacktip shark season/sex combinations had training data AUC scores ranging from 0.87 to 0.91 and cross validation AUC scores ranging from 0.79 (±0.02) to 0.87 (±0.01) ( Table 2). For females, the three most influential predictors of abundance were bottom salinity (36%), distance from shore (16%), and longitude (13%); there was no seasonal variability in the strength of these predictors ( Table 2). Female blacktip sharks preferred inshore regions with bottom salinity values between 26 and 33 psu, particularly west of 88 • W longitude (Figure 6A), preferences which did not change seasonally. Predicted suitable habitat for female blacktip sharks was concentrated in the Mississippi Sound and west of the Chandeleur Islands ( Figure 7A). Coefficients of variation were relatively high for females in the nearshore, western region of the study area in the spring, but not summer or autumn (Supplementary Figure S2A).
Male blacktip shark distributions were influenced by daylength (17%), bottom salinity (13%), and surface temperature (11%) ( Table 2); similar to females, no seasonal variation in the strength of these predictors was evident. Males were most abundant during periods when daylength exceeded 750 min, in waters between 26 and 33 psu, and at sea surface temperatures between 20 and 24 • C ( Figure 6A). Suitable habitat was predicted across the continental shelf for male blacktip sharks, particularly in the offshore waters south of Mississippi ( Figure 7A). In contrast to females, coefficients of variation were highest for males during autumn and comparatively lower in spring and summer (Supplementary Figure S2A).
Sandbar sharks were the second most abundant large coastal shark. Males and females were encountered across a relatively narrow range of sizes, encompassing primarily sub-adult and adult individuals (Baremore and Hale, 2012; Figure 2D). Individuals were captured offshore in waters greater than 20 m, yet were practically absent from nearshore areas ( Figure 3D). Females were encountered more often than males (1.38:1), which differed from the expected sex ratio (x 2 = 12.18, p < 0.01). Final models for sandbar shark season/sex combinations had training data AUC scores of 0.94 and cross validation AUC scores ranging from 0.88 (±0.02) to 0.89 (±0.02) ( Table 2). The two most influential predictors of abundance for female sandbar sharks were distance from shore (41%) and bottom salinity (10%), with little to no seasonal variation ( Table 2). Females were only encountered greater than 20 km from shore, where bottom salinity was greater than 34 psu ( Figure 6B). Consequently, areas with suitable habitat for female sandbar shark were predominately located past the 20 m isobath ( Figure 7B). In general, coefficients of variation increased with depth (Supplementary Figure S2B).
Male sandbar shark distributions were influenced by distance from shore (∼30%), longitude (∼17%), and substrate grain size (11%) ( Table 2). The relative influence of these predictors varied little across seasons. Male sandbar sharks were encountered farther offshore than females (20-45 km), almost exclusively east of 88 • W longitude in areas with increasing substrate grain size ( Figure 6B). Areas where suitable habitat was predicted for male sandbar sharks were restricted to the southeastern region of our sampling region (Figure 7B), and coefficients of variation were relatively small (Supplementary Figure S2B).

Shelf-Associated Sharks
Scalloped hammerheads were the most abundant shelf-associated species; males were encountered across a wider range of sizes than females (Figure 2E). Males and females were encountered broadly across the continental shelf ( Figure 3E). Females were much less common than males (0.36:1), which differed from the expected sex ratio (x 2 = 42.63, p < 0.01). Final models for scalloped hammerhead season/sex combinations had training data AUC scores ranging from 0.92 to 0.95 and cross validation AUC scores ranging from 0.83 (±0.02) to 0.87 (±0.03) ( Table 2). Female scalloped hammerhead abundance was influenced by distance from shore (28-29%), depth (∼23%), and bottom salinity (7-8%) ( Table 2), with little seasonal variation. Female scalloped hammerheads were encountered 75-85 km offshore, at depths between 50 and 100 m ( Figure 8A). Suitable habitat for female scalloped hammerheads was restricted to a small core area directly south of Mobile Bay, Alabama ( Figure 9A). Coefficients of variation were highest for females during summer (Supplementary Figure S3A).
Like females, male scalloped hammerhead abundance was influenced by distance from shore (∼17%), as well as bottom salinity (11-12%) and bottom velocity (11%) ( Table 2), with marginal seasonal variation. Male scalloped hammerheads were encountered closer to shore (15-75 km offshore) relative to females, at depths between 25 and 100 m ( Figure 8A). Consequently, greater amounts of suitable habitat were predicted for male scalloped hammerheads compared to females (Figure 9A), and coefficients of variation associated with these predictions were comparatively higher, particularly in spring (Supplementary Figure S3A).
Smoothhound sharks were the second most abundant shelf-associated shark and were encountered across a range of sizes, primarily encompassing sub-adults and adults (SEDAR 39, Figure 2F). Smoothhound sharks were encountered on the continental shelf, but across a relatively narrow depth range ( Figure 3F). Females were markedly more common than males (5.42:1), which differed from the expected sex ratio (x 2 = 79.19, p < 0.01). Final models for smoothhound shark season/sex combinations had training data AUC scores of 0.99, and cross validation AUC scores ranging from 0.93 (±0.03) to 0.97 (±0.01) ( Table 2). Female smoothhound shark abundance was influenced by depth (37-44%), distance from shore (30-33%), and daylength (7-10%) ( Table 2). The relative strength of these predictors varied slightly between seasons. Female smoothhound sharks were encountered at depths between 50 and 100 m, between 75 and 85 km offshore ( Figure 8B). Predicted habitat suitable for female smoothhound sharks was similar to predicted suitable habitat for female scalloped hammerheads ( Figure 9B). Coefficients of variation were highest for females in the deepest portion of the study area, particularly in spring and summer (Supplementary Figure S3B).
Male smoothhound shark abundance was influenced by distance from shore (35-37%), depth (17-18%), and substrate grain size (8-9%) ( Table 2). Male smoothhound sharks were encountered 75-85 km offshore, at depths ranging from 50 to 100 m (Figure 8B), similar to what was shown for females, however, larger areas of suitable habitat were predicted for males FIGURE 5 | Predicted relative abundance from BRT models for female (red) and male (blue) (A) Atlantic sharpnose shark and (B) blacknose shark during spring, summer and autumn. Darker shades indicate higher predicted abundance, and lighter shades indicate lower predicted abundance.
relative to females (Figure 9B). Coefficients of variation were similar between males and females (Supplementary Figure S3B).

DISCUSSION
A diverse assemblage of sharks, dominated by the Atlantic sharpnose shark, was identified across the extent of the area we surveyed. Our findings are remarkably similar to those of Murawski et al. (2018), who identified the continental shelf of the northern GoM as the most speciose region of the Gulf and documented high relative abundance of the Atlantic sharpnose shark in this region. Our study illustrates a strong degree of sexual segregation between female and male Atlantic sharpnose sharks, as first recorded by Springer (1967). This segregation is most pronounced during spring and autumn, when adult males are abundant in the inshore waters north of the barrier islands while females inhabit offshore areas. Parsons and Hoffmayer (2005) also noted this segregation, and suggested the possibility of a reproductively-motivated summer exodus of mature male Atlantic sharpnose sharks to deeper waters for mating. Our findings support this assertion; spatial and temporal overlap between adult male and female Atlantic sharpnose sharks was only predicted offshore during the summer. Collectively, these patterns support the notion that this fast-growing species likely does not depend on discrete, inshore nurseries (Knip et al., 2010;Heupel et al., 2018); rather, mating and parturition for this species may take place in deeper waters during the summer months, as suggested by Drymon et al. (2010).
Notable differences in population structure were apparent between Atlantic sharpnose shark and the other small coastal shark examined in detail, the blacknose shark. First, the composition of the blacknose shark population differed from that of Atlantic sharpnose shark. A smaller size range of blacknose FIGURE 6 | Line plots for the three most influential marginal effects from the binary BRT models for female (red) and male (blue) (A) blacktip shark and (B) sandbar shark during spring, summer, and autumn.
sharks was noted, and most individuals were identified as adults (Driggers et al., 2004); moreover, the sex ratio was biased in favor of males. Second, the factors best predicting distribution and abundance of blacknose sharks appear to act independently of sex, unlike the factors influencing Atlantic sharpnose shark distribution and abundance. Interestingly, males were predicted to exist in higher abundance shallower and closer to shore than females. Third, although the blacknose shark was the second most abundant species sampled in our study, its relative abundance was less than one third of the relative abundance of Atlantic sharpnose shark. This trend aligns with both fishery-dependent (Scott-Denton et al., 2011) and fishery-independent (Murawski et al., 2018) catch data; both of these datasets also suggest that blacknose sharks are even less abundant in areas outside of the northern GoM. The most recent assessment for blacknose shark was unable to determine the status of the GoM stock (SEDAR, 2011); consequently, commercial retention of blacknose sharks in the GoM is currently prohibited. Despite preliminary signs of shark recovery as indicated by increasing population trends for GoM Atlantic sharpnose shark, the relative abundance of GoM blacknose shark is declining (Peterson et al., 2017). Yet, of the six species we examined, blacknose shark possesses the largest FIGURE 7 | Predicted relative abundance from BRT models for female (red) and male (blue) (A) blacktip shark and (B) sandbar shark during spring, summer and autumn. Darker shades indicate higher predicted abundance, and lighter shades indicate lower predicted abundance.
predicted range of suitable habitat. We are therefore optimistic that ongoing efforts to rebuild this stock will not be impeded by a lack of suitable habitat.
Blacktip shark was the most abundant large coastal species observed during the survey. This species was captured across a wide range of sizes (neonate through mature individuals) and showed strong sexual segregation. Gravid females and neonates were encountered throughout the nearshore regions of the study area, in accordance with the classical shark population model proposed by Springer (1967). Despite the presence of gravid females inshore, temperature was not a significant predictor of female presence or relative abundance, suggesting that the thermal-niche hypothesis (Wearmouth and Sims, 2008) does not explain this case of sexual segregation. Instead, we speculate that gravid females may be preferentially occupying nearshore waters to provide predator refugia for neonate blacktip sharks following parturition. Surprisingly, the clear sexual segregation demonstrated by this species has not resulted in a sex-biased fishery, as has been shown for shortfin mako in the equatorial Pacific (Mucientes et al., 2009) and spiny dogfish and smooth dogfish in the northwest Atlantic (Dell'Apa et al., 2014, 2018Haugen et al., 2017). Rather, the most recent update to the stock assessment (SEDAR, 2018) indicates that GoM blacktip sharks represent an outcome of successful fisheries management. Effective quota allocation may contribute to this accomplishment. The commercial quota for the federal GoM blacktip shark fishery is disproportionately divided between the areas east and west of 88 • W longitude. The habitat suitability predictions from our BRTs are in striking agreement with this longitudinal quota delineation. Large areas of highly suitable blacktip shark habitat are predicted in the western GoM subregion (i.e., west of 88 • W longitude); in contrast, little suitable habitat is predicted in the eastern GoM sub-region (i.e., east of 88 • W longitude). Accordingly, the 2019 quota for western GoM blacktip sharks is approximately ten times greater than the available quota for eastern GoM blacktip sharks. Our BRT results support the assertion that GoM blacktip shark can support sustainable fisheries when strong science-based management strategies are applied (Simpfendorfer and Dulvy, 2017).
Sandbar sharks were encountered across a relatively narrow range of sizes; most individuals were adults approximately 2 m in length. This species has been protected in the GoM since 2010, and populations in the GoM seem to be increasing as a result (Peterson et al., 2017). Our findings show very clear and interesting associations between sandbar sharks and Alabama's Artificial Reef Zone (AARZ), where depredation (i.e., the partial or complete removal of a hooked species by a non-target species) by sandbar sharks is common in hook and line fisheries (Drymon et al., 2019). In Western Australia, depredation by sharks is higher in areas receiving higher fishing pressure, perhaps indicting learned behavior by the sharks (Mitchell et al., 2018), which may also be the case for sandbar sharks inhabiting the AARZ. Although limited, recapture data from four adult sandbar sharks tagged in the AARZ indicates strong site fidelity. Each of these FIGURE 9 | Predicted relative abundance from BRT models for female (red) and male (blue) (A) scalloped hammerhead and (B) smoothhound sharks during spring, summer and autumn. Darker shades indicate higher predicted abundance, and lighter shades indicate lower predicted abundance.
individuals was recaptured less than 30 km from the location where it was originally tagged. In particular, one individual was at large for nearly 5 years after release and was recaptured only 18 km from the tagging location (Drymon unpublished data). We postulate that the strong associations of sandbar sharks with, and their possible site fidelity to, the AARZ may be driven by a factor not examined in our BRT analysis, namely red snapper (Lutjanus campechanus). The AARZ represents a regional hotspot of red snapper abundance (Karnauskas et al., 2017) and experiences heavy fishing pressure during relatively short recreational fishing seasons. The recent compression of the federal red snapper recreational fishing season has resulted in a derby-style fishery, particularly in the AARZ, where the number of anglers per day has increased in accordance with shortened seasons (Powers and Anson, 2016). Consequently, the spatially and temporally limited recreational fishery for red snapper may explain the high predicted abundance of sandbar sharks in the AARZ. This and similar competition between anglers and sharks over a shared resource (in this case, red snapper) exemplifies an emerging challenge to be addressed as shark populations continue to recover (Carlson et al., 2019).
Male scalloped hammerheads were encountered across a wider range of sizes and were approximately three times more abundant than females. Sexual segregation has previously been described for scalloped hammerheads in the GoM (Branstetter, 1987), and our results support those findings. Similar dynamics have been documented for scalloped hammerheads in the Pacific Ocean, where sexual segregation may allow faster growth among females resulting from increased prey availability in offshore waters (Klimley, 1987). Our findings are consistent with this concept, termed the "forage-selection hypothesis" (Sims, 2005), although dietary studies would be required to confirm this as the mechanism driving sexual segregation in scalloped hammerheads in the northern GoM. Despite the sexual segregation shown for scalloped hammerheads in the current study, distance from shore was the strongest predictor of presence for both sexes, with little seasonal variation. Specifically, females' predicted suitable habitat is deeper than that of males. Habitat use determined from satellite-tagged scalloped hammerheads support these findings; females more frequently use shelf-edge (>200 m isobath) areas compared to shallower depths used by males (Wells et al., 2018). Regardless of the mechanism driving sexual segregation in scalloped hammerheads, additional efforts to understand habitat use by mature females are needed to reveal habitat areas of particular concern for this species, such as potential areas for parturition (Wells et al., 2018).
Individuals in the genus Mustelus are morphologically conserved, with three genetically distinct species known to occur across the GoM (Giresi et al., 2015). As a result, the smooth dogfish, the Gulf smoothhound and the Florida smoothhound are currently managed as a single, multi-species complex. Accordingly, our sex ratio and BRT analyses apply to all sharks in the genus Mustelus. However, our field identifications since publication of Giresi et al. (2015) demonstrate that nearly all smoothhound sharks represented in our survey are the Gulf smoothhound. Gulf smoothhound sharks are thought to be endemic to the GoM (Compagno et al., 2005) and limited data on their horizontal distribution indicates that they occur predominantly between the Florida panhandle and Louisiana (Murawski et al., 2018). In addition to this horizontal stratification, Gulf smoothhound shark also occupy specific depth ranges. Giresi et al. (2015) report the mean depth of capture for Gulf smoothhound to be 112 m, strikingly similar to the depth range predicted by the BRTs in the current study. Collectively, these habitat preferences place the Gulf smoothhound at risk; species with narrow habitat ranges (e.g., limited geographic range, depth range) are most vulnerable to the synergistic effects of climate change and overfishing (Cheung et al., 2018). In fact, Cheung et al. (2018) report that the two species at highest risk from the combined effects of overfishing and climate change are in the genus Mustelus (M. canis and M. lenticulatus). Therefore, while Gulf smoothhound sharks (and the smoothhound shark complex) are not currently overfished or experiencing overfishing (NOAA, 2019), the relatively restrictive habitat preferences identified for this species indicate that it should be monitored closely and managed in the future with a precautionary approach.
The shark assemblage detailed in this study straddles a unique zoogeographic province in the northern GoM. Portnoy and Gold (2012) suggest that multiple vicariant events formed a marine suture-zone in the vicinity of Mobile Bay, Alabama, which hosts at least fifteen pairs of sister taxa. It is thought that historical events created a physical barrier for sessile organisms that led to speciation of these organisms; however, these events also resulted in modern-day longitudinal differences that appear to be influencing the shark assemblage. One such difference is a gradient in sediment type. Specifically, the area south of Mobile Bay represents a transition from terrigenous sediments in the west to carbonate sediments in the east (Balsam and Beeson, 2003). This transition occurs at ∼88 • W longitude and stems from a variety of other factors that influence the distribution of sharks in this assemblage in different ways. For example, the terrigenous sediments that characterize habitats west of Mobile Bay stem from increases in primary production associated with nutrient-rich outflow from Mobile Bay and the Mississippi River. The resultant increases in secondary production and forage fish biomass provide an explanation for the abundance of suitable habitat for the Atlantic sharpnose shark in areas west of Mobile Bay (Drymon et al., 2012(Drymon et al., , 2013. In addition, increased freshwater output lowers salinity, which the current study suggests increases habitat suitability for blacktip sharks. Unrelated to these abiotic factors, dense concentrations of artificial habitat occur just east of 88 • W longitude in the AARZ, and may explain the high predicted concentrations for sandbar sharks, particularly males. Thus, our findings align with previous work indicating a regional zoogeographic break in the juvenile shark assemblage in this region (Bethea et al., 2015) and expand upon those findings by revealing trends in predicted habitat among adult sharks across deeper waters.
Our findings are consistent with previous studies (summarized in Wearmouth and Sims, 2008) and indicate that sexual segregation is the norm for sharks in the northern GoM. Similar to the "sexual line in the sea" demonstrated for shortfin mako in the South Pacific Ocean (Mucientes et al., 2009), the strong sexual segregation demonstrated in this study needs to be considered when forming state and federal management strategies, particularly for species like blacktip shark, which are currently managed under different quotas in the eastern and western GoM. The factors influencing relative abundance among the shark assemblage showed surprisingly little seasonal variation, suggesting that future efforts to characterize the dynamics of this shark community could take place any time of the year. These efforts would be aided by the application of acoustic and/or satellite telemetry, which would provide insights that complement traditional monitoring programs. For example, electronic tags would serve to verify habitat suitability predictions for animals outside of the range of the current study area (e.g., Lee et al., 2019). In addition, telemetry would allow investigations of long-term population dynamics, including partial migrations, which have been shown for larger, more vagile species examined in the current study (e.g., scalloped hammerhead, Hoyos-Padilla et al., 2014). Ultimately, the species-and sex-specific differences in habitat use by the shark assemblage in this region provide further evidence that long-term fishery-independent monitoring should be prioritized.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

AUTHOR CONTRIBUTIONS
JD and SP contributed to all aspects of the research including design, implementation, analysis, and writing. SD, JF, ES, and AJ contributed to data analysis and writing. AK and JM assisted with implementation of the research and writing.

FUNDING
Funding for this work was provided via grants from NOAA SEAMAP and the Alabama Department of Conservation and Natural Resources, Marine Resources Division through the National Fish and Wildlife Foundation, Gulf Environmental Benefits Fund.

ACKNOWLEDGMENTS
We would like to thank Dauphin Island Sea Lab captains Tom Guoba, Jonathan Wittmann, Rodney Collier, Willie Johnston, and Russell Wilson for keeping us safe on the vessels. Additional thanks to the numerous graduate students, technicians, and interns who assisted with field work, in particular Crystal Hightower and Trey Spearman. This manuscript was greatly improved by thoughtful comments from the two reviewers.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars. 2020.00035/full#supplementary-material FIGURE S1 | Coefficient of variation of predicted relative abundance from BRT models for female (red) and male (blue) (A) Atlantic sharpnose shark and (B) blacknose shark during spring, summer and autumn. Darker shades indicate higher variance, and lighter shades indicate lower variance.
FIGURE S2 | Coefficient of variation of predicted relative abundance from BRT models for female (red) and male (blue) (A) blacktip shark and (B) sandbar shark during spring, summer and autumn. Darker shades indicate higher variance, and lighter shades indicate lower variance.
FIGURE S3 | Coefficient of variation of predicted relative abundance from BRT models for female (red) and male (blue) (A) scalloped hammerhead and (B) smoothhound sharks during spring, summer and autumn. Darker shades indicate higher variance, and lighter shades indicate lower variance.