Abstract
Introduction:
Nursery areas for sandbar sharks (Carcharhinus plumbeus) are well-delineated along the Atlantic coast of the United States, with only a single nursery area identified in the Gulf of Mexico on the west coast of Florida.
Methods:
Fishery-independent surveys and fishery-dependent data were used to explore the frequency and seasonality of young-of-the-year (YOY) sandbar sharks off the coast of Texas.
Results:
Data presented in this study demonstrate YOY are caught off the Texas coast, suggesting a potential nursery in the region. Data collected by the Texas Parks and Wildlife Department estuarine gillnet surveys and recreational shore-based shark anglers documented the presence of YOY sandbar sharks in nearshore and estuarine waters of Texas.
Discussion:
While the increase in YOY individuals were detected in the decades long fishery-independent surveys (i.e., estuarine gillnet, SEAMAP coastal and NMFS offshore longlines), fishery-dependent data collected via a shore-based recreational shark tournament documented a substantial rise in YOY sandbar sharks over a much shorter time period (10 years), supporting the usefulness of data collected by citizen scientists. For shark species, especially those with rebuilding plans, such as sandbar sharks, identifying and conserving nursery areas is important, as the decreasing amount of suitable habitat could be a limiting factor for population recovery.
1 Introduction
Sandbar sharks (Carcharhinus plumbeus) are a large, highly migratory, coastal species that inhabits temperate, subtropical and tropical waters of the western North Atlantic Ocean, including the Gulf of Mexico (also known as the ‘Gulf of America’ in U.S. federal documents, hereafter Gulf; , ; ) and have been classified as globally endangered on the IUCN Red List of Threatened Species (Rigby et al., 2021). Sharks like sandbar sharks are highly susceptible to overexploitation due to their life history characteristics (e.g., slow growth and low reproductive output; Musick et al., 2000). Sandbar sharks were historically the primary targeted species of the U.S. commercial shark fishery, leading to overfishing as identified by the 2006 stock assessment (SEDAR, 2006). As a result, directed fisheries were closed in 2008 and remained closed today, except for a limited entry research fishery (NMFS, 2008; ). This closure significantly limited the availability of new and updated biological data (e.g., reproductive variables, litter size, and parturition) for the species. However, effective management strategies rely on updated and accurate data, especially for those species that are in rebuilding fishery management plans.
Periodic reassessment of reproduction, continued monitoring of the age structure, and habitat use of sandbar sharks will ensure that assessment scientists and managers have the best data possible. Size and age at 50% maturity were 151.6 cm FL and 12.1 years for males while females were 154.9 cm or 13.1 years (). Peak mating and parturition occur from April through July, with an average of 8 pups per litter (, Musick and Colvocoresses, 1988, ; ). The reproductive cycle is at least biennial, although there is evidence that some females have triennial cycles (Springer, 1960, , ; ; ).
Adult sandbar sharks are known to aggregate offshore and make large north–south migrations off the east coast of the United States, including the Gulf (; ; , Stunz, unpublished data). Neonate and juvenile sandbar sharks frequent coastal nursery areas (), with the largest and most well-known nursery area in the western North Atlantic Ocean for this species in Chesapeake Bay (Springer, 1960, Musick and Colvocoresses, 1988), with smaller nurseries reported in New York, Delaware, Virginia, South Carolina, and Florida (Springer, 1960; ; ; ; ; ). Juvenile sandbar sharks have been reported in Chesapeake Bay and eastern shore of Virginia in depths less than 10 m (; ), and in Bulls Bay, South Carolina (). Pratt and Merson (1996) have also reported Delaware Bay, New Jersey, as a major nursery area for neonate and juvenile sharks.
Although the nursery areas for this species in U.S. waters are largely considered to be off the east coast, a nursery area has been suggested in the western Gulf, based on a few females collected with full-term pups near the mouth of the Mississippi River and few juveniles captured off Texas and Louisiana (; Springer, 1960; , ). This nursery area was hypothesized to result from gravid females occasionally making their way into the Gulf but not be self-sustaining (Springer, 1960). documented the capture of neonates and juvenile sandbar sharks (n = 105) in Indian Pass and St. Andrew Sound, and because neonates (age <1) reside through their first summer in the parturition areas (Pratt and Merson, 1996), he reported these individuals were unlikely to have migrated from other areas. This description of a nursery area in the eastern Gulf suggests that nursery areas in the Gulf might be more extensive than previously reported (; ). Texas bays have previously been described as nursery habitat for several other shark species, including blacktip sharks (Carcharhinus limbatus) and bull sharks (Carcharhinus leucas), with the highest catches historically in Matagorda and San Antonio Bay (, ; , ).
The abundance estimates and diversity of shark populations have traditionally been described using fishery independent sampling (Xu et al., 2015; ). These standardized survey methods allow for more statistically robust assessments because their continuity and random stratified designs allow for more accurate assessments of population dynamics, climate change, management actions, and migratory dynamics (; ). One such fishery independent survey is the National Marine Fisheries Service bottom longline (NMFS BLL) survey, which samples shark species in offshore waters (>9 nautical miles (nm) or 16.67 km). A complementary survey is the Southeast Area Monitoring and Assessment Program (SEAMAP), which is a collaboration between federal/state/universities that focus on collecting bottom longline data in coastal waters (<9 nm; ). State sampling programs (e.g., Texas Parks and Wildlife Department [TPWD] gillnet and BLL) are used to monitor state waters and bays and estuaries (hereafter, estuaries). Despite their value to the assessment process, fishery-independent surveys are relatively expensive and may be restricted both temporally and geographically by funding (, ).
An increasing trend in fisheries science is the incorporation of recreational anglers to provide data that otherwise would not have been logistically and financially feasible (Williams et al., 2015; ). The Texas Shark Rodeo (TSR, texassharkrodeo.com) is an annual 10-month tournament, where anglers target sharks from the shore (<500 m from shore and <5 m depth) following best practices for catch-photo-release. Anglers participating in the TSR tag and submit a photograph of their catch along with biological measurements. Beginning in 2014, this tournament has recorded over 13,000 shark submissions, with anglers correctly identifying shark species 97.2% of the time, providing a unique long-term dataset for the Texas coast (). reported that sandbar sharks were the third most commonly caught species in the TSR. The purpose of this study was to compare fishery-independent (e.g., BLL and gillnet surveys) and dependent (e.g., TSR) datasets for sandbar shark seasonality and size composition, particularly for YOY individuals associated with the Texas coast.
2 Methods
This research was approved by the Texas A&M University-Corpus Christi Institutional Animal Care and Use Committee under protocols #08–15, #08–18, and #2023-007, and also by National Park Service permits PAIS-2010-SCI-0009, PAIS-2015-SCI-0001, and PAIS-2016-SCI-0018. Tournament participants captured and tagged sharks following the rules and regulations of the TPWD and the TSR.
2.1 NMFS Offshore BLL
Since 1995, the NOAA National Marine Fisheries Service (NMFS, ; , https://www.fisheries.noaa.gov/inport/item/28636) conducts fisheries-independent BLL surveys along the US continental shelf waters from Cape Hatteras, NC to Brownsville, TX, from late July to the beginning of October. The longline, which is 1,852 meters long, is deployed at depths ranging from 9 to 366 m. It consists of a 536 kg test monofilament mainline and 100 gangions, each with a 3.7 m, 3.0 mm diameter monofilament leader and a hook. Bait used was Atlantic Mackerel (Scomber scombrus). Weights are placed at the beginning, middle, and end of the line. Hook types varied over time, with J-hooks used from 1995 to 1998 and 15/0 circle hooks from 2001 to the present. In 1999 and 2000, both hook types were used. Given changes in hook type prior to 2001 and the subsequent standardization of longline gear and sampling procedures, this study includes only data collected from 2001 onward. Soak times are generally limited to one hour to minimize mortality rates of all captured organisms. Survey locations are randomly selected by stratified-random sampling with proportional allocation based on depth and continental shelf area. Catch data, including species identification, length, weight, and sex, as well as environmental data, such as temperature, salinity, and dissolved oxygen, are collected at each site ().
2.2 SEAMAP Coastal BLL
In 2008, recognizing the gap in the NMFS offshore BLL survey, which only samples at depths greater than 9 m, SEAMAP launched a complementary fishery-independent BLL survey in the coastal waters of the U.S. Gulf. SEAMAP is a collaborative program between federal, state, and university partners that focuses on collecting and sharing fisheries-independent, gear-standardized data. The sampling methodology used by SEAMAP mirrors that of the NMFS BLL survey ().
From 2008 to 2014, SEAMAP survey design and effort allocation parameters varied between states. In 2015, several changes were implemented: the depth range was modified to 3–10 m to include previously unsampled waters while minimizing overlap with the NMFS BLL survey’s range of 9–366 m, and sampling efforts were redistributed to create a more spatially balanced survey. Additionally, prior to 2015, SEAMAP surveys were conducted from March to October, but in 2015, sampling was reorganized into three seasonal periods: Spring (April–May), Summer (June–July), and Fall (August–September; ). For this study, all months were analyzed separately to avoid differences between seasonal sampling and allow comparable results between datasets.
2.3 TPWD Estuarine Gillnet
Gillnet data were obtained from the TPWD fishery-independent gillnet monitoring program. Since 1982, 45 gillnets were deployed in each major estuarine system during the spring (April–June) and fall (September–November) seasons. During each 10-week sampling period, a minimum of three gillnet deployments were conducted per week, to ensure sampling was distributed throughout the full sampling period (Zapp Sluis et al., 2025). Gill nets (182.9 m × 1.2 m) were comprised of four sections with varying mesh sizes: 76 mm, 102 mm, 127 mm, and 152 mm. These nets were deployed overnight, beginning approximately one hour before sunset and retrieved three to four hours after sunrise (average soak time ± SD = 13.7 ± 1.4 hours). Gillnets were set perpendicular to the shoreline, with the 76 mm mesh closest to the shore, at randomly selected locations (; ; ; Plumlee et al., 2018).
All organisms larger than 5 mm caught in gillnets were identified to the lowest taxonomic level (typically species), counted, and measured to the nearest millimeter for stretched total length (TL). At each sampling location, environmental data were recorded while nets soaked, including date, location, water temperature (°C), salinity (psu), dissolved oxygen (mg L−¹), and water depth (m). Environmental parameters such as temperature, salinity, and dissolved oxygen were measured approximately 0.15 m below the water surface (; Plumlee et al., 2018). Environmental conditions were only measured at deployment and retrieval of the gillnets, limiting the ability to assess how overnight variations during sampling influenced shark captures (). Thus, the average environmental conditions of each sampling event at deployment and retrieval were used for analyses.
All major estuaries in the TPWD dataset were included in analyses, with the exception of Cedar Lakes and East Matagorda Bay, which were excluded due to irregular sampling and low sample sizes over the sampling period (Plumlee et al., 2018). Moreover, data prior to 1983 have also been excluded, due to variation in sample sizes (). Data from 2024 was excluded due to sampling inconsistency in the fall. Thus, the temporal scope covered by this dataset is from 1983 to 2023, except for Sabine Lake where sampling began in 1986 (Zapp Sluis et al., 2025).
2.4 Shore-based TSR
Texas Shark Rodeo has been on-going since 2014, providing a long-term, fishery-dependent dataset. Each year, a mean of 325 teams (range: 180–560 teams/year) participate with teams consisting of up to 6 anglers. Although variations exist, the general strategy for shore-based shark fishing involves the use of large reels spooled with 800-1,000 m of 50-lb to 100-lb test of either monofilament or braided line with a top shot of approximately 100 m of monofilament of increased strength. A leader made of either wire or monofilament, a weight, and a hook (13/0 to 24/0) is attached to the top shot line and baited with typically sections of stingray (Rhinoptera spp. or Dasyatis spp.), crevalle jack (Caranx hippos), or striped mullet (Mugil cephalus), which is either surf casted or kayaked out 100–400 m offshore ().
Captured sharks were identified to species, measured, photographed with a tournament issued ruler, and then released. Date of capture, location, stretched total length (STL; measured from the tip of the snout to the tip of the stretched upper caudal lobe), fork length (FL), sex, species, and tag number if applicable, were recorded. Photos submitted to TSR were viewed for species confirmation (, ). Environmental data were not submitted as part of the tournament.
2.5 Data analyses
Statistical analyses were completed in R. Pups, a.k.a. young-of-the-year (YOY) were defined as less than 71 cm FL (age-1 length as used in SEDAR 54 [2017]). For NMFS offshore BLL surveys where FL was not measured but natural total length (TL) data were recorded, FL was estimated using the equation from :
For the SEAMAP coastal BLL and estuarine gillnet surveys, where only stretched total length was recorded, the NOAA calculator was used to determine the relationship between stretched total length and FL (Natanson et al., 2022).
For all fishery-independent datasets (i.e., offshore BLL, coastal BLL and estuarine gillnet), total length to fork length conversions were calculated to ensure consistency across all datasets. When FL was already measured, the recorded value was retained. Fork length was reported for all sharks in TSR, so no conversion was necessary.
Sandbar shark catch data were then standardized as catch-per-unit-effort (CPUE). For the shore-based TSR dataset, CPUE was calculated as:
Where catch is the number of sandbar sharks captured, and A is the total number of angler-days. An angler-day was defined as each unique combination of angler and date in which any shark was reported. In the case of estuarine gillnet data, CPUE was calculated as catch per net-hour (soak time):
Estuarine gillnet and shore-based TSR CPUE calculations were visualized along the Texas coast using ggplot2, ggspatial, and sf packages in R.
For the coastal and offshore BLL datasets, CPUE was calculated as:
Where c is the number of sharks captured, h is the number of hooks deployed, and t is the soak time in minutes. The multipliers 60 min and 100 hooks were used to standardize CPUE data as number of sharks caught per 100 hook hours (; Pickens et al., 2022). The fishery-independent datasets (i.e., offshore BLL, coastal BLL, and estuarine gillnet) included environmental data, therefore, generalized additive models (GAMs) with CPUE as response variable and environmental, spatial, and temporal predictors with a Tweedie distribution were applied. The Tweedie distribution was selected for the GAMs due to its flexibility in handling both zero-inflated data and overdispersion (), which are common features of ecological datasets. Model parametrization ranged from simpler to more complex GAMs incorporating smooth cubic regression spline functions of continuous predictors (year, depth and environmental variables), cyclic cubic splines for month predictors and thin-plate regression splines of spatial coordinates. In the case of NMFS BLL dataset, as sampling only occurred during the summer months, month was not included in the models. To account for potential effects associated with the sampling station, survey was included as a random effect in the BLL models and bay was included as a random effect in the gillnet models. Random effects were retained only when they improved model performance. Models included a gamma value of 1.4 to penalize overfitting and k was left as default to allow for model flexibility.
Multicollinearity and pairwise correlations among predictor variables were assessed prior to model fitting using scatter plot matrices and variance inflation factors (VIFs). All correlations were below 0.7 and VIFs were below 5, allowing retention of all variables in the models (O’Brien et al., 2025; Peterson et al., 2017), except for the depth variable in the case of the offshore NMFS BLL dataset, which was removed due to high correlation (0.83) with other environmental variables.
Model selection was based on Akaike’s Information Criterion (AIC; ; ). Manual backwards stepwise model fitting was performed to select explanatory variables influencing CPUE, based on minimizing the AIC. At each step of the backward selection process, the variable with the highest p-value (p > 0.05) was removed from the model. This process was repeated until further removal of variables led to an increase in AIC of >1%, at which point the variable was retained and the selection process was finished (; ; ; Sluis et al., 2021).
Model performance was evaluated using 5-fold cross-validation, assessing root mean square error (RMSE), mean absolute error (MAE), and the Pearson correlation between predicted and observed values (O’Brien et al., 2025). Temporal autocorrelation in model residuals was initially assessed through visual inspection of autocorrelation function (ACF) plots (; ; ). Formal testing of temporal autocorrelation was conducted using the testTemporalAutocorrelation() function from the DHARMa package in R (). In addition to temporal checks, DHARMa diagnostics were also used for broader model validation, with diagnostic plots evaluated to detect any systematic trends or overdispersion (; ; O’Brien et al., 2025). Spatial autocorrelation in residuals was assessed separately using a variogram of the residuals (Zuur et al., 2009) and Monte Carlo simulation of Moran’s I (; O’Brien et al., 2025; Potts and Rose, 2018; Rooker et al., 2012). Test results were considered acceptable when no significant autocorrelation was detected (p-value > 0.05). Finally, partial effects of individual covariates on CPUE were visualized using the ggplot2 package in R.
3 Results
The total number of sandbar sharks and the number of YOY individuals were recorded for each of the four datasets used in this study. A total of 1,333 individuals were sampled across the four datasets (offshore BLL, coastal BLL, estuarine gillnet, and shore-based TSR; Table 1) spanning the Texas coast. Based on length, 267 individuals were classified as YOY, with YOY catches concentrated between Corpus Christi Bay to north of Galveston Bay (Figure 1).
Table 1
| Year | All | <71 cm FL | ||||||
|---|---|---|---|---|---|---|---|---|
| Offshore BLL | Coastal BLL | Gillnet | Shore-based TSR | Offshore BLL | Coastal BLL | Gillnet | Shore-based TSR | |
| 1983 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1984 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1985 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1986 | NA | NA | 1 | NA | NA | NA | 1 | NA |
| 1987 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1988 | NA | NA | 1 | NA | NA | NA | 0 | NA |
| 1989 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1990 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1991 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1992 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1993 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1994 | NA | NA | 1 | NA | NA | NA | 0 | NA |
| 1995 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1996 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1997 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1998 | NA | NA | 0 | NA | NA | NA | 0 | NA |
| 1999 | NA | NA | 1 | NA | NA | NA | 1 | NA |
| 2000 | NA | NA | 3 | NA | NA | NA | 3 | NA |
| 2001 | 5 | NA | 0 | NA | 0 | NA | 0 | NA |
| 2002 | 11 | NA | 0 | NA | 0 | NA | 0 | NA |
| 2003 | 8 | NA | 1 | NA | 1 | NA | 1 | NA |
| 2004 | 10 | NA | 0 | NA | 0 | NA | 0 | NA |
| 2005 | 0 | NA | 0 | NA | 0 | NA | 0 | NA |
| 2006 | 5 | NA | 0 | NA | 0 | NA | 0 | NA |
| 2007 | 7 | NA | 2 | NA | 0 | NA | 2 | NA |
| 2008 | 1 | 0 | 0 | NA | 0 | NA | 0 | NA |
| 2009 | 13 | 0 | 4 | NA | 1 | 0 | 4 | NA |
| 2010 | 10 | 22 | 1 | NA | 0 | 0 | 1 | NA |
| 2011 | 0 | 0 | 0 | NA | 0 | 0 | 0 | NA |
| 2012 | 15 | 0 | 0 | NA | 0 | 0 | 0 | NA |
| 2013 | 18 | 1 | 0 | NA | 0 | 0 | 0 | NA |
| 2014 | 5 | 0 | 1 | 54 | 0 | 0 | 1 | 2 |
| 2015 | 8 | 0 | 0 | 102 | 0 | 0 | 0 | 74 |
| 2016 | 11 | 0 | 4 | 93 | 0 | 0 | 4 | 53 |
| 2017 | 18 | 0 | 4 | 45 | 2 | 0 | 4 | 30 |
| 2018 | 5 | 0 | 2 | 98 | 0 | 0 | 2 | 13 |
| 2019 | 12 | 0 | 1 | 88 | 0 | 0 | 1 | 16 |
| 2020 | 0 | 0 | 4 | 48 | 0 | 0 | 3 | 6 |
| 2021 | 4 | 5 | 12 | 105 | 0 | 0 | 12 | 9 |
| 2022 | 4 | 2 | 0 | 201 | 0 | 2 | 0 | 8 |
| 2023 | 9 | 0 | 0 | 130 | 0 | 0 | 0 | 9 |
| 2024 | 5 | 0 | 0 | 92 | 0 | 0 | 0 | 2 |
| TOTAL | 184 | 30 | 43 | 1056 | 4 | 2 | 40 | 222 |
Individual sandbar shark catches for fishery independent surveys (i.e., offshore BLL, coastal BLL, estuarine gillnet survey) and fishery dependent surveys (i.e., TSR) in and off Texas.
NA means the survey had not been implemented yet or were excluded from these analyses (see methods for details). Bold numbers are totals.
Figure 1
3.1 Offshore BLL
The offshore BLL surveys (since 2001) reported 184 sandbar sharks and suggested an increase in individuals overtime (Figure 2). To formally evaluate temporal patterns in CPUE, two candidate models were developed and compared based on statistical performance and residual diagnostics. Model Offshore A_1 (AIC = 484.23, deviance explained = 46.6%; Supplementary Material Table 1) showed no signs of overfitting, with training and cross-validation metrics remaining consistent (CV RMSE = 0.82, CV MAE = 0.37, CV R = 0.365; Supplementary Material Table 2), indicating reliable out-of-sample predictive performance. DHARMa residual diagnostics revealed no significant violations of model assumptions, although the outlier and dispersion tests returned significant results (Supplementary Material Figure 1). The dispersion test indicated mild under dispersion (φ = 0.96), meaning the observed variance in the response variable was slightly lower than expected under the model. This pattern can increase sensitivity to a small number of extreme observations (Zuur et al., 2009), consistent with the significant outlier test result. Given the characteristically right-skewed and zero-inflated nature of CPUE data for relatively rare size classes such as neonates, these deviations are not unexpected and are unlikely to meaningfully bias the directional trends inferred from the model. A Durbin-Watson test indicated the presence of temporal autocorrelation (p < 0.002), this was acknowledged as a limitation affecting the precision of standard errors rather than the validity of the model structure itself.
Figure 2
Model B_1, despite achieving a lower AIC (466.42) and higher deviance explained (67.6%), exhibited clear signs of overfitting, with a substantial divergence between training and cross-validation R values (0.868 vs. 0.375) (Supplementary Material Table 2). Additionally, DHARMa diagnostics for Model Offshore B_1 revealed a significant combined adjusted quantile test (Supplementary Material Figure 2), indicating misspecification of the residual distribution. Given that a lower AIC value is insufficient justification for model selection when accompanied by overfitting and residual misspecification, Model A_1 was retained as the preferred model for inference. The temporal autocorrelation identified in Model A_1 was accounted for in the interpretation of results, and temporal trends are discussed descriptively with appropriate caution regarding the precision of significance tests.
Model_Offshore_A_1 <- gam(CPUE ~.
s(YEAR, bs=“cr”) +.
s(TEMP, bs=“cr”) + s(SALINITY, bs=“cr”) +.
s(FLYLAST_OUTLON, FLYLAST_OUTLAT, bs=“tp”),
family = tw(link = “log”), gamma=1.4,
data = offshore_data_C_plumbeus, method = “REML”).
Adult sandbar shark CPUE increased from 2000 to 2015 and then gradually declined thereafter, supported by the GAM analysis, in which Year was a significant predictor of CPUE (p = 8.18 × 10−7; Figure 2B), though this result should be interpreted cautiously given temporal autocorrelation in the model residuals. Size‐composition patterns indicate that YOY sandbar sharks were rare offshore (Figure 3D), with only four YOY individuals recorded (mean: 56.85 cm FL; range: 54.5 - 59.6 cm FL), three females and one male. These occurred almost exclusively in September, coinciding with sampling times for Texas (Figure 4D). Larger juveniles and adults size structure was relatively consistent across months. The significant smooth term for location in the GAM model of the BLL data (s(Lat,Lon), p < 2e-16) confirmed that CPUE varied geographically offshore.
Figure 3
Figure 4

Fork length size composition of sandbar sharks for each data series: shore-based (A), estuaries (B), coastal (C) and offshore (D).
3.2 Coastal BLL
The coastal BLL (since 2008) recorded 30 sandbar sharks, of which two were YOY individuals (mean: 70.45 cm FL; range: 70.4 - 70.5 cm FL), one female and one male (Figure 3C). The YOY individuals were caught in May and July (Figure 4C). The GAM model for BLL coastal data did not converge due to low sample size and outliers, and therefore, it was not presented here.
3.3 Estuarine gillnet
The estuarine gillnet survey recorded a much higher proportion of YOY sandbar sharks, with 40 out of 43 (93%) individuals measuring <71 cm FL (mean: 57.66 cm FL; range: 40.1 - 68.0 cm FL; Figure 3B). Of these 40 individuals, 12 were females, 13 were males and no sex was recorded for 15 individuals. These catches first appeared in the 2000s and showed an increase after 2015 (Figure 4B). The observed increase in sandbar and YOY individuals was supported by the GAM analysis, in which Year was a significant predictor of CPUE (p = 0.0033; Figure 2A). The estuarine gillnet model selected through backward stepwise selection included smooth terms for year, month (cyclic cubic spline, k = 5), salinity, and depth, a thin-plate regression spline of spatial coordinates (longitude and latitude), and bay as a random effect (Model Gillnet C_1 in Supplementary Material Table 3). QQ plot and DHARMa residual diagnostics indicated no significant deviations from expected distributional assumptions, and the model was therefore retained (Supplementary Material Figure 3).
Model_Gillnet_C_1<- gam (CPUE ~.
s(YEAR, bs=“cr”) + s(MONTH, bs=“cc”, k=5) +.
s(SALINITY, bs=“cr”) +.
s(START_DEEP_DEPTH, bs=“cr”) +.
s(START_LONGITUDE, START_LATITUDE, bs=“tp”) +.
s(MAJOR_AREA_CODE, bs=“re”),
family = tw(link = “log”), gamma=1.4,
data = gillnet_data_C_plumbeus, method = “REML”).
Size‐composition patterns indicate that YOY sandbar sharks were more common in estuaries than adults, with a more distinct temporal pattern in size composition, with the largest number of sandbar sharks captured in September, indicating a strong seasonal peak late in the summer (Figure 4B). However, the smallest individuals observed in the gillnet data were captured earlier in the year, specifically in June, suggesting an early-summer presence of younger juveniles in estuaries. Relative abundance of YOY has been generally increasing annually (Figure 5). Spatially, YOY sandbar shark captures were concentrated along the central Texas coast, particularly in Matagorda Bay and Corpus Christi Bay, with only a single catch recorded near south Galveston Bay and another in Baffin Bay (Figure 1A). The significant smooth term for location in the GAM model of the estuaries (s(Lat, Lon), p = 0.011637) confirmed that YOY CPUE varied geographically along the coast.
Figure 5

Annual mean CPUE of young-of-the-year (YOY) sandbar sharks (<71 cm FL) from the TPWD gillnet survey conducted in Texas estuaries.
3.4 Shore-based TSR
From 2014 to 2024, 1,056 sandbar sharks were caught and reported through the shore-based TSR tournament. Of those, 849 were female and 207 were male. Of the total, 222 individuals classified as YOY (<71 cm FL; mean: 57.37 cm FL; range: 43.18 - 69.85 cm FL; Figure 3A) based on reported FL, with 131 of those females and 91 males (Table 1). An increased proportion of YOY was observed in summer months (June to September), coinciding with a decrease in the number of larger individuals starting in May (Figure 4A). Years 2015 to 2017 had the highest number of YOY sandbar sharks reported in the shore-based fishery despite the highest total number of all sandbar sharks reported in later years (2022-2023; Figure 6). Relative abundance (CPUE) peaked for YOY in 2016, declined until 2020, and has since been steady (Figure 6). Overall, CPUE of sandbar sharks has been increasing since the inception of TSR in 2014 (Figure 6). Highest CPUE was observed between Matagorda Bay and Galveston Bay systems with higher abundance also observed near Corpus Christi Bay. Lower CPUE was observed along the rest of the coast. For YOY, CPUE was also highest between Matagorda Bay and Galveston Bay systems; however, YOY CPUE was lower along the remaining coastline, especially south of Matagorda Bay despite significant fishing effort in that region (Figure 1B).
Figure 6

Mean CPUE trends for sandbar sharks caught in the shore-based TSR. Young-of-the-year (YOY) are individuals <71 cm FL.
4 Discussion
While most sandbar shark nurseries have been documented along the Atlantic coast (
In this study, sandbar shark CPUE has increased over time across three datasets (coastal, estuarine, and shore-based), suggesting an increasing trend in abundance in the western Gulf and supporting the effectiveness of the rebuilding plan. The offshore BLL predominantly captured adult sandbar sharks, reporting 184 sandbar sharks encountered in the 22-year survey period. The coastal BLL survey had even less encounters with a total of 30 individuals reported in the 16-year survey period. The estuarine gillnet survey had the longest time period of 41 years and reported 43 sandbar shark encounters, most of which were YOY. The catch has continued to increase in this survey starting in 2016. Data collected through the shore-based TSR, which has the shortest survey period of 10 years, sampled 1,056 sandbar sharks of all age classes, with about 21% being YOY. This sample size was ~5x more than the offshore BLL, ~35x more than the coastal BLL, and ~25x more than the estuarine gillnet data. These results may be interpreted in a few ways. Directly, as the sandbar shark population rebuilds, more individuals will be observed. However, the significant difference in catch numbers between the surveys may be explained spatially. The inshore surveys (i.e., estuarine gillnet and shore-based TSR) may be better suited to monitoring potential nursery habitat than the coastal and offshore BLL surveys (
Although not standardized, the intensive sampling conducted by citizen scientists during events like the shore-based TSR demonstrates that these programs can help monitor stocks of species that are routinely encountered (e.g.,
Seasonal trends were also evident in this study. Catches of YOY sandbar sharks were highest during the summer months, from June to September, which overlaps with the known pupping time for sandbar sharks in the northwestern Atlantic (
Relative abundance of sandbar sharks in both the shore-based TSR and TPWD estuarine gillnet surveys, despite their statewide coverage, were concentrated near the Matagorda Bay and Galveston Bay systems, suggesting that this might be the location of a potential nursery in Texas. For this area to be a nursery, three proposed criteria must be met as described by
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The requirement of ethical approval was waived by Texas A&M University-Corpus Christi Institutional Animal Care and Use Committee for the studies involving animals because data was collected under state and federal permits or through citizen submitted data. No protocol was deemed to be needed. The studies were conducted in accordance with the local legislation and institutional requirements.
Author contributions
KG-B: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. PDR-L: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. JS: Data curation, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review & editing. JD: Conceptualization, Data curation, Formal analysis, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. RW: Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Supervision, Visualization, Writing – original draft, Writing – review & editing. FM-A: Data curation, Investigation, Methodology, Resources, Validation, Writing – original draft, Writing – review & editing. WD: Data curation, Investigation, Methodology, Resources, Validation, Writing – original draft, Writing – review & editing. MS: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. Funding was provided by Coastal Conservation Association-Texas, the Harte Research Institute, and Texas Parks and Wildlife Department (CA-0005812).
Acknowledgments
We thank Texas Shark Rodeo, Sharkathon, and their participants for helping to collect this data, as we would not have been able to collect it without them. Thanks to the staff and students of the Center for Sportfish Science and Conservation for their help in processing data. We also thank the personnel from both TPWD and NOAA for their assistance with sample collection.
Conflict of interest
The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author JD declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Correction note
6 May 2026 This article has been corrected with minor changes. These changes do not impact the scientific content of the article.
26 May 2026 A correction has been made to this article. Details can be found at: 10.3389/fmars.2026.1879610.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2026.1814009/full#supplementary-material
References
1
AjemianM. J.JoseP. D.FroeschkeJ. T.WildhaberM. L.StunzG. W. (2016). Was everything bigger in Texas? Characterization and trends of a land-based recreational shark fishery. Mar. Coast. Fish.8, 553–566. doi: 10.1080/19425120.2016.1227404. PMID:
2
AkaikeH. (1973). Maximum likelihood identification of Gaussian autoregressive moving average models. Biometrika60, 255–265. doi: 10.1093/biomet/60.2.255
3
AltobelliA. N.SzedlmayerS. T. (2020). Migration and residency of sandbar, Atlantic sharpnose, bull, and nurse sharks in the Northern Gulf of Mexico. North Am. J. Fish. Manage.40, 1324–1343. doi: 10.1002/nafm.10501. PMID:
4
AraujoG.LabajaJ.SnowS.HuveneersC.PonzoA. (2020). Changes in diving behaviour and habitat use of provisioned whale sharks: implications for management. Sci. Rep.10, 16951. doi: 10.1038/s41598-020-73416-2. PMID:
5
BaremoreI. E.HaleL. F. (2012). Reproduction of the sandbar shark in the western North Atlantic Ocean and Gulf of Mexico. Mar. Coast. Fish.4, 560–572. doi: 10.1080/19425120.2012.700904. PMID:
6
BigelowH. B.Pérez FarfanteI.SchroederW. C. (1948). Fishes of the Western North Atlantic: 1. Lancelets. Cyclostomes. Sharks. Memoir Sears Foundation for Marine Research, I (New Haven: Sears Foundation for Marine Research, Yale University), xvii, 576.
7
BonarS. A.Contreras-BalderasS.IlesA. C. (2009). “ An introduction to standardized sampling,” in Standard methods for sampling north american freshwater fishes. Eds. BonarS. A.HubertW. A.WillisD. W. ( American Fisheries Society, Bethesda, MD).
8
BranstetterS. (1987). Age and growth validation of newborn sharks held in laboratory aquaria, with comments on the life history of the Atlantic sharpnose shark, Rhizoprionodon terraenovae. Copeia1987, 291–300. doi: 10.2307/1445764
9
BurnhamK. P.AndersonD. R. (2002). Model selection and multimodel inference: a practical information-theoretic approach (New York, NY: Springer New York).
10
CarlsonJ. K. (1999). Occurrence of neonate and juvenile sandbar sharks, Carcharhinus plumbeus, in the northeastern Gulf of Mexico. Fishery Bull.97, 387–391.
11
CarvalhoF. J.de SantanaD. G.SampaioM. V. (2020). Modeling overdispersion, autocorrelation, and zero-inflated count data via generalized additive models and bayesian statistics in an aphid population study. Neotrop. Entomol.49, 40–51. doi: 10.1007/s13744-019-00729-x. PMID:
12
CaseyJ. G.KohlerN. E. (1990). “ Long distance movements of Atlantic sharks from NMFS cooperative shark tagging program,” in Elasmobranchs as living resources: advances in the biology, ecology, systematics, and the status of the fisheries. Eds. PrattH. L. J.GruberS. H.TaniuchiT. ( U.S. Dept. Commer), 87–91. NOAA Tech. Rep. NMFS 90.
13
CaseyJ. G.NatansonL. J. (1992). Revised estimates of age and growth of the sandbar shark (Carcharhinus plumbeus) from the western North Atlantic. Can. J. Fish. Aquat.Sci.49, 1474–1477. doi: 10.1139/f92-162. PMID:
14
CastroJ. I. (1993). The shark nursery of Bulls Bay, South Carolina, with a review of the shark nurseries of the southeastern coast of the United States. Environ. Biol. Fishes38, 37–48. doi: 10.1007/978-94-017-3450-9_4. PMID:
15
ClarkE.vonSchmidtK. (1965). Sharks of the central gulf coast of florida. Bull. Mar. Sci.15, 13–83.
16
CliffG.DudleyS. F. J.DavisB. (1988). Sharks caught in the protective gill nets off Natal, South Africa: 1. the sandbar shark Carcharhinus plumbeus (Nardo). South. Afr. J. Mar. Sci.7, 255–265.
17
CornicM.RookerJ. R. (2018). Influence of oceanographic conditions on the distribution and abundance of blackfin tuna (Thunnus atlanticus) larvae in the Gulf of Mexico. Fish. Res.201, 1–10. doi: 10.1016/j.fishres.2017.12.015. PMID:
18
CurtisB.Brewster-GeiszK.WilliamsonA. (2025). Highly migratory species sandbar sharks management history (North Charleston, SC: SEDAR). SEDAR101-DW-01.
19
DanceM. A.RookerJ. R. (2019). Cross-shelf habitat shifts by red snapper (Lutjanus campechanus) in the Gulf of Mexico. PLoS One14, e0213506. doi: 10.1371/journal.pone.0213506. PMID:
20
DennisD.PlagányiÉ. (2015). Cost benefit of fishery-independent surveys: are they worth the money? Mar. Policy58, 108–115. doi: 10.1016/j.marpol.2015.04.016. PMID:
21
DriggersW. B. I.IngramG. W.GraceM. A.GledhillC. T.HenwoodT. A.HortonC. N.et al. (2008). Pupping areas and mortality rates of young tiger sharks Galeocerdo cuvier in the western North Atlantic Ocean. Aquat. Biol.2, 161–170. doi: 10.3354/ab00045. PMID:
22
DrymonJ. M.DedmanS.FroeschkeJ. T.SeubertE.JeffersonA. E.KroetzA. M.et al (2020). Defining sex-specific habitat suitability for a northern Gulf of Mexico shark assemblage. Frontiers in Marine Science.7, 35. doi: 10.3389/fmars.2020.00035
23
DrymonJ. M.PowersS. P.DindoJ.DzwonkowskiB.HenwoodT. A. (2010). Distributions of sharks across a continental shelf in the northern Gulf of Mexico. Mar. Coast. Fish.2, 440–450. doi: 10.1577/c09-061.1
24
EllisJ. K.MusickJ. A. (2007). Ontogenetic changes in the diet of the sandbar shark, Carcharhinus plumbeus, in lower Chesapeake Bay and Virginia (USA) coastal waters. Environ. Biol. Fishes80, 51–67. doi: 10.1007/s10641-006-9116-2. PMID:
25
FlowersK.BabcockE.PapastamatiouY.BondM.LambN.MirandaA.et al. (2022). Varying reef shark abundance trends inside a marine reserve: evidence of a Caribbean reef shark decline. Mar. Ecol. Prog. Ser.684, 97–107. doi: 10.3354/meps13954. PMID:
26
FontaineP.JensenC. C.MatichP.RookerJ. R.WellsR. D. (2024). Predicting habitat suitability for the co-occurrence of an estuarine mesopredator and two top predatory fishes. Front. Fish. Sci.2, 1443923. doi: 10.5194/oos2025-385. PMID:
27
FroeschkeJ.StunzG. W.WildhaberM. L. (2010). Environmental influences on the occurrence of coastal sharks in estuarine waters. Mar. Ecol. Prog. Ser.407, 279–292. doi: 10.3354/meps08546. PMID:
28
GibsonK. J.StreichM. K.ToppingT. S.StunzG. W. (2019). Utility of citizen science data: A case study in land-based shark fishing. PLoS One14, e0226782. doi: 10.1371/journal.pone.0226782. PMID:
29
Gibson BanksK.StreichM. K.DrymonJ. M.ScyphersS. B.MohanJ. A.WellsR. D.et al. (2023). Talk is cheap: Direct evidence of conservation-based changes in angler behavior. Conserv. Sci. Pract.5, e13001. doi: 10.1111/csp2.13001. PMID:
30
Gibson-BanksK.StreichM.StunzG.SmithJ. (2025). Partnering with recreational anglers to characterize the Texas shore-based shark fishery. Texas State Wildlife Grant Program TX T-266-R-1 /F23AF03057-00 Final Report44 pg.
31
GilchristR.DrinkwaterD. (2000). “ The use of the Tweedie distribution in statistical modelling,” in COMPSTAT: Proceedings in computational statistics 14th symposium held in Utrecht, the Netherlands ( Physica-Verlag HD, Heidelberg), 313–318.
32
GraceM.HenwoodT. (1997). Assessment of the distribution and abundance of coastal sharks in the U.S. Gulf of Mexico and Eastern Seaboard 1995 and 1996. Mar. Fish. Rev.59, 23–32.
33
GrubbsR. D.MusickJ. A.ConrathC. L.RomineJ. G. (2007). “ Long term movements, migration, and temporal delineation of a summer nursery for juvenile sandbar sharks in the Chesapeake Bay region,” in Shark nursery grounds of the Gulf of Mexico and the east coast waters of the United States. Eds. McCandlessC. T.KohlerN. E.PrattH. L. ( American Fisheries Society, Bethesda, Maryland), 87–107. Symposium 50.
34
HartigF.AbregoN.BushA.ChaseJ. M.Guillera-ArroitaG.LeiboldM. A.et al. (2024). Novel community data in ecology-properties and prospects. Trends Ecol. Evol.39, 280–293. doi: 10.1016/j.tree.2023.09.017. PMID:
35
HeistE. J.GoldJ. R. (1999). Microsatellite DNA variation in sandbar sharks (Carcharhinus plumbeus) from the Gulf of Mexico and mid-Atlantic Bight. Copeia1999 (1), 182–186. doi: 10.2307/1447399
36
HeistE. J.GravesJ. E.MusickJ. A. (1995). Population genetics of the sandbar shark (Carcharhinus plumbeus) in the Gulf of Mexico and Mid-Atlantic Bight. Copeia1995 (3), 555–562. doi: 10.2307/1446752. PMID:
37
HendonJ. M.HoffmayerE. R.PollackA. G.MareskaJ.Martinez-AndradeF.ResterJ.et al. (2025). Impacts of survey design on a Gulf of Mexico bottom longline survey and the transition to a unified, stratified-random design. Front. Mar. Sci.11. doi: 10.3389/fmars.2024.1426756. PMID:
38
HeupelM. R.CarlsonJ. K.SimpfendorferC. A. (2007). Shark nursery areas: concepts, definition, characterization and assumptions. Mar. Ecol. Prog. Ser.337, 287–297. doi: 10.3354/meps337287. PMID:
39
HowardR. A.CiannelliL.WakefieldW. W.HaltuchM. A. (2023). Comparing fishery-independent and fishery-dependent data for analysis of the distributions of Oregon shelf groundfishes. Fish. Res.258, 106553. doi: 10.1016/j.fishres.2022.106553. PMID:
40
HueterR. E.TyminskiJ. P. (2007). “ Species-specific distribution and habitat characteristics of shark nurseries in Gulf of Mexico waters off Peninsular Florida and Texas,” in Shark nursery grounds of the Gulf of Mexico and the east coast waters of the United States. Eds. McCandlessC. T.KohlerN. E.PrattH. L. ( American Fisheries Society, Bethesda, Maryland), 193–223. Symposium 50.
41
JacksonM. C.HuangL.XieQ.TiwariR. C. (2010). A modified version of Moran’s I. Int. J. Health Geographics9, 33. doi: 10.1186/1476-072x-9-33. PMID:
42
JonesL. M.GraceM. A. (2002). Shark nursery areas in the bay systems of Texas. In Shark nursery grounds of the Gulf of Mexico and the east coast waters of the United States: an overview, Eds. McCandlessC. T.PrattH. L.KohlerN. E. (Bethesda, Maryland: American Fisheries Society, Symposium 50), 209–220.
43
KingB.MorrisC.GreenJ.GregoryR.SnelgroveP.CoteD.et al. (2024). Acoustic telemetry and network analysis reveal seasonal spatial overlap between gadid species in a subarctic coastal marine protected area. Can. J. Fish. Aquat.Sci.81, 1547–1559. doi: 10.1139/cjfas-2023-0272. PMID:
44
KohlerN. E.CaseyJ. G.TurnerP. A. (1996). “ Length-length and length-weight relationships for 13 shark species from the Western North Atlantic,” in NOAA technical memorandum NMFS-NE-110.
45
KohlerN. E.CaseyJ. G.TurnerP. A. (1998). NMFS cooperative shark tagging program 1962–93: an atlas of shark tag and recapture data—National Marine Fisheries Service—statistical data included. Mar. Fish. Rev.60, 1–87.
46
MakwinjaR.MengistouS.KaundaE.AlemiewT.PhiriT. B.KosamuI. B. M.et al. (2021). Modeling of lake Malombe annual fish landings and catch per unit effort (CPUE). Forecasting3, 39–55. doi: 10.3390/forecast3010004. PMID:
47
Martinez-AndradeF. (2018). Trends in relative abundance and size of selected finfishes and shellfishes along the Texas coast: November 1975-2016 (Austin, Texas: Texas Parks and Wildlife Department, Coastal Fisheries Division), 140.
48
MatichP.AultJ. S.BoucekR. E.BryanD. R.GastrichK. R.HarveyC. L.et al. (2017a). Ecological niche partitioning within a large predator guild in a nutrient-limited estuary. Limnol. Oceanogr.62, 934–953. doi: 10.1002/lno.10477. PMID:
49
MatichP.BigelowC. L.ChambersB.DoddsJ. J.HebertJ. A.LemieuxA.et al. (2022). Delineation of blacktip shark (Carcharhinus limbatus) nursery habitats in the north‐western Gulf of Mexico. J. Fish Biol.101, 236–248. doi: 10.1111/jfb.15103. PMID:
50
MatichP.MohanJ. A.PlumleeJ. D.TinHanT.WellsR. D.FisherM. (2017b). Factors shaping the co-occurrence of two juvenile shark species along the Texas Gulf Coast. Mar. Biol.164, 141. doi: 10.1007/s00227-017-3173-2. PMID:
51
MatichP.MullinsL.PlumleeJ. D.FisherM. R.MareskaJ.DrymonJ. M. (2026). The dynamic nature of nearshore shark nurseries in the northern Gulf of Mexico. J. Biogeogr. 35, e70210. doi: 10.1111/jbi.70210. PMID: . In revision.
52
McAuleyR. B.SimpfendorferC. A.HyndesG. A.LenantonR. C. J. (2007). Distribution and reproductive biology of the sandbar shark, Carcharhinus plumbeus (Nardo), in western Australian waters. Mar. Freshw. Res.58, 116–126. doi: 10.1071/mf05234. PMID:
53
MedvedR. J.MarshallJ. A. (1981). Feeding behavior and biology of young sandbar sharks, Carcharhinus plumbeus (Pisces, Carcharhinidae), in Chincoteague Bay, Virginia. Fishery Bull.79, 441–447.
54
MedvedR. J.StillwellC. E.CaseyJ. J. (1985). Stomach contents of young sandbar sharks, Carcharhinus plumbeus, in Chincoteague Bay, Virginia. Fishery Bull.83, 395–402.
55
MersonR. R. (1998). Nursery grounds and maturation of the sandbar shark in the western North Atlantic (University of Rhode Island, Kingston). Doctoral dissertation.
56
MersonR. R.PrattH. L. (2001). Distribution, movements and growth of young sandbar sharks, Carcharhinus plumbeus, in the nursery grounds of Delaware Bay. Environ. Biol. Fishes61, 13–24. doi: 10.1023/a:1011017109776. PMID:
57
MillerT. J.SkalskiJ. R.IanelliJ. N. (2007). Optimizing a stratified sampling design when faced with multiple objectives. ICES J. Mar. Sci.64, 97–109. doi: 10.1093/icesjms/fsl013. PMID:
58
MurawskiS. A.PebblesE. B.GraciaA.TunnellJ. W.ArmenterosM. (2018). Comparative abundance, species composition, and demographics of continental shelf fish assemblages throughout the Gulf of Mexico. Mar. Coast. Fish.: Dyn. Manage. Ecosyst. Sci.10, 325–346. doi: 10.1002/mcf2.10033. PMID:
59
MusickJ. A.BranstetterS.ColvocoressesJ. A. (1993). “ Trend in shark abundance from 1974–1991 for the Chesapeake Bight region of the U.S. mid-Atlantic Coast,” in Conservation biology of elasmobranchs. Ed. BranstetterS. ( U.S. Dept. Commer. NOAA), 1–18. Tech. Rep. NMFS 115.
60
MusickJ. A.BurgessG.CaillietG.CamhiM.FordhamS. (2000). Management of sharks and their relatives (Elasmobranchii). Fisheries25, 9–13. doi: 10.1577/1548-8446(2000)025<0009:mosatr>2.0.co;2
61
MusickJ. A.ColvocoressesJ. A. (1988). “ Seasonal recruitment of subtropical sharks in Chesapeake Bight, USA,” in Intergovernmental Oceanographic Commission/Food and Agriculture Organization workshop on recruitment in tropical coastal demersal communities. Eds. Yáñez-ArancibiaA.PaulyD. ( UNESCO, Paris), 301–311.
62
NatansonL. J.McCandlessC. T.PasserottiM. S.BelcherC. N.BowlbyH.DriggersW. B.et al. (2022). Morphometric conversions for 33 shark species from the Western North Atlantic Ocean. Mar. Fish. Rev.84, 1–65. doi: 10.7755/MFR.84.3-4.1
63
NMFS (National Marine Fisheries Service) (2008). Final amendment 2 to the consolidated Atlantic highly migratory species fishery management plan (Silver Spring, MD: NOAA, Office of Sustainable Fisheries, Highly Migratory Species Management Division).
64
O’BrienK.CarlsonJ. K.CortésE.DriggersW. B.FrazierB. S.LatourR. J. (2025). Evaluation of a spatiotemporal index standardization method for coastal shark species; implications for future stock assessments. Front. Mar. Sci.12, 1621720. doi: 10.3389/fmars.2025.1621720
65
PetersonC. D.BelcherC. N.BetheaD. M.DriggersW. B.FrazierB. S.LatourR. J. (2017). Preliminary recovery of coastal sharks in the south-east United States. Fish Fish.18, 845–859. doi: 10.51952/9781529234985.ch011
66
PickensB. A.TaylorJ. C.CampbellM. D.DriggersW. B. (2022). Offshore snapper and shark distributions are predicted by prey and area of nearby estuarine environments in the Gulf of Mexico, USA. Mar. Ecol. Prog. Ser.682, 169–189. doi: 10.3354/meps13925. PMID:
67
PlumleeJ. D.DanceK. M.MatichP.MohanJ. A.RichardsT. M.TinHanT. C.et al. (2018). Community structure of elasmobranchs in estuaries along the northwest Gulf of Mexico. Estuar. Coast. Shelf Sci.204, 103–113. doi: 10.1016/j.ecss.2018.02.023. PMID:
68
PottsS. E.RoseK. A. (2018). Evaluation of GLM and GAM for estimating population indices from fishery independent surveys. Fish. Res.208, 167–178. doi: 10.1016/j.fishres.2018.07.016. PMID:
69
PrattH. L.MersonR. R. (1996). Delaware Bay sandbar shark nursery pilot study ( Narragansett Laboratory, National Marine Fisheries Service, NOAA), 23. Progress report.
70
RigbyC. L.DerrickD.DickenM.HarryA. V.PacoureauN.SimpfendorferC. (2021). “ Carcharhinus plumbeus,” in The IUCN red list of threatened species 2021: e.T3853A2874370. doi: 10.2305/IUCN.UK.2021-2.RLTS.T3853A2874370.en
71
RookerJ. R.SimmsJ. R.WellsR. D.HoltS. A.HoltG. J.GravesJ. E.et al. (2012). Distribution and habitat associations of billfish and swordfish larvae across mesoscale features in the Gulf of Mexico. PLoS One7, e34180. doi: 10.1371/journal.pone.0034180. PMID:
72
SEDAR (Southeast Data, Assessment, and Review) (2006). HMS sandbar shark stock assessment report (North Charleston, SC).
73
SEDAR (Southeast Data, Assessment, and Review) (2017). HMS sandbar shark stock assessment report (North Charleston, SC).
74
SluisM. Z.FujiwaraM.Martinez-AndradeF.WellsR. D. (2025). Spatiotemporal shifts and influence of environmental parameters on estuarine-dependent fishes in Texas bays. Estuar. Coast. Shelf Sci.312, 109034. doi: 10.1016/j.ecss.2024.109034. PMID:
75
SluisM. Z.JudkinsH.DanceM. A.VecchioneM.CornicM.SuttonT.et al. (2021). Taxonomic composition, abundance and habitat associations of squid paralarvae in the northern Gulf of Mexico. Deep Sea Res. Part. I: Oceanogr. Res. Papers174, 103572. doi: 10.1016/j.dsr.2021.103572. PMID:
76
SpringerS. (1960). Natural history of the sandbar shark Eulamia milberti. Fishery Bull.61, 1–38.
77
WilliamsS. M.HolmesB. J.PepperellJ. G. (2015). The novel application of non-lethal citizen science tissue sampling in recreational fisheries. PLoS One10. doi: 10.1371/journal.pone.0135743. PMID:
78
XuB.ZhangC.XueY.RenY.ChenY. (2015). Optimization of sampling effort for a fishery-independent survey with multiple goals. Environ. Monit. Assess.187, 252. doi: 10.1007/s10661-015-4483-9. PMID:
79
ZuurA. F.IenoE. N.WalkerN. J.SavelievA. A.SmithG. M. (2009). Mixed effects models and extensions in ecology with R (New York: Springer).
Summary
Keywords
citizen scientist, nursery, pups, sandbar sharks, Texas, young of the year
Citation
Gibson-Banks K, Dominguez Rein-Loring P, Smith JM, Drymon JM, Wells RJD, Martinez-Andrade F, Driggers III WB and Streich MK (2026) Evidence of a potential sandbar shark (Carcharhinus plumbeus) nursery in the Western Gulf of Mexico. Front. Mar. Sci. 13:1814009. doi: 10.3389/fmars.2026.1814009
Received
19 February 2026
Revised
31 March 2026
Accepted
01 April 2026
Published
04 May 2026
Corrected
26 May 2026
Volume
13 - 2026
Edited by
Adrian C. Gleiss, Murdoch University, Australia
Reviewed by
Eric Reyier, Herndon Solutions Group, LLC, United States
Jorge Manuel Morales-Saldaña, McGill University, Canada
Updates

Check for updates
Copyright
© 2026 Gibson-Banks, Dominguez Rein-Loring, Smith, Drymon, Wells, Martinez-Andrade, Driggers and Streich.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kesley Gibson-Banks, Kesley.Banks@tamucc.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.