An Index of Abundance for Coastal Species of Juvenile Sharks from the Northeast Gulf of Mexico

A fishery-independent assessment of juvenile coastal shark populations in U.S. waters of the northeast Gulf of Mexico was conducted using two methods: gillnets and longlines. Surveys were conducted monthly during April–October in two fixed sampling areas from 1996 to 1998. The Atlantic sharpnose shark, Rhizoprionodon terraenovae, and the blacktip shark, Carcharhinus limbatus, were the most common species captured with either longlines or gillnets. An additional 14 shark species were captured, and juvenile indices of abundance were developed for 8 species with gillnets and 6 species of sharks with longlines. Trends in catch-per-unit-effort were found to vary depending on species. Length-frequency information revealed that the majority of sharks captured were juveniles. Given the direct relationship between stock and recruitment for sharks, continued monitoring of juvenile abundance will aid in determining the strength of the parental stock size and for predicting future population strength.

Influence of soak time and fish accumulation on catches of reef fishes in a multispecies trap survey

Catch rates from fishery-independent surveys often are assumed to vary in proportion to the actual abundance of a population, but this approach assumes that the catchability coefficient (q) is constant. When fish accumulate in a gear, the rate at which the gear catches fish can decline, and, as a result, catch asymptotes and q declines with longer fishing times. We used data from long-term trap surveys (1990–2011) in the southeastern U.S. Atlantic to determine whether traps saturated for 8 reef fish species because of the amount of time traps soaked or the level of fish accumulation (the total number of individuals of all fish species caught in a trap). We used a delta-generalized-additive model to relate the catch of each species to a variety of predictor variables to determine how catch was influenced by soak time and fish accumulation after accounting for variability in catch due to the other predictor variables in the model. We found evidence of trap saturation for all 8 reef fish species examined. Traps became saturated for most species across the range of soak times examined, but trap saturation occurred for 3 fish species because of fish accumulation levels in the trap. Our results indicate that, to infer relative abundance levels from catch data, future studies should standardize catch or catch rates with nonlinear regression models that incorporate soak time, fish accumulation, and any other predictor variable that may ultimately influence catch. Determination of the exact mechanisms that cause trap saturation is a critical need for accurate stock assessment, and our results indicate that these mechanisms may vary considerably among species.

Habitat and diet overlap of 4 piscivorous fishes: variation on the inner continental shelf off New Jersey

Piscivorous fishes, many of which are economically valuable, play an important role in marine ecosystems and have the potential to affect fish and invertebrate populations at lower trophic levels. Therefore, a quantitative understanding of the foraging ecology of piscivores is needed for ecosystem-based fishery management plans to be successful. Abundance and stomach contents of seasonally co-occurring piscivores were examined to determine overlap in resource use for Summer Flounder (Paralichthys dentatus; 206–670 mm total length [TL]), Weakfish (Cynoscion regalis; 80–565 mm TL), Bluefish (Pomatomus saltatrix; 55–732 mm fork length [FL]), and Striped Bass (Morone saxatilis; 422–920 mm FL). We collected samples from monthly, fishery-independent trawl surveys conducted on the inner continental shelf (5–27 m) off New Jersey from June to October 2005. Fish abundances and overlaps in diet and habitat varied over this study period. A wide range of fish and invertebrate prey was consumed by each species. Diet composition (determined from 1997 stomachs with identifiable contents) varied with ontogeny (size) and indicated limited overlap between most of the species size classes examined. Although many prey categories were shared by the piscivores examined, different temporal and spatial patterns in habitat use seemed to alleviate potential competition for prey. Nevertheless, the degree of overlap in both fish distributions and diets increased severalfold in the fall as species left estuaries and migrated across and along the study area. Therefore, the transitional period of fall migration, when fish densities are higher than at other times of the year, may be critical for unraveling resource overlap for these seasonally migrant predators.

Multiscale analysis of factors that affect the distribution of sharks throughout the northern Gulf of Mexico

Identification of the spatial scale at which marine communities are organized is critical to proper management, yet this is particularly difficult to determine for highly migratory species like sharks. We used shark catch data collected during 2006–09 from fishery-independent bottom-longline surveys, as well as biotic and abiotic explanatory data to identify the factors that affect the distribution of coastal sharks at 2 spatial scales in the northern Gulf of Mexico. Centered principal component analyses (PCAs) were used to visualize the patterns that characterize shark distributions at small (Alabama and Mississippi coast) and large (northern Gulf of Mexico) spatial scales. Environmental data on temperature, salinity, dissolved oxygen (DO), depth, fish and crustacean biomass, and chlorophyll-a (chl-a) concentration were analyzed with normed PCAs at both spatial scales. The relationships between values of shark catch per unit of effort (CPUE) and environmental factors were then analyzed at each scale with co-inertia analysis (COIA). Results from COIA indicated that the degree of agreement between the structure of the environmental and shark data sets was relatively higher at the small spatial scale than at the large one. CPUE of Blacktip Shark (Carcharhinus limbatus) was related positively with crustacean biomass at both spatial scales. Similarly, CPUE of Atlantic Sharpnose Shark (Rhizoprionodon terraenovae) was related positively with chl-a concentration and negatively with DO at both spatial scales. Conversely, distribution of Blacknose Shark (C. acronotus) displayed a contrasting relationship with depth at the 2 scales considered. Our results indicate that the factors influencing the distribution of sharks in the northern Gulf of Mexico are species specific but generally transcend the spatial boundaries used in our analyses.

Biscayne Bay commercial pink shrimp, Farfantepenaeus duorarum, fisheries, 1986-2005

The Biscayne Bay bait (1986–2005) and food (1989–2005) fisheries for pink shrimp were examined using dealer-reported individual vessel-trip landings data, separated by waterbody code to represent only catches from Biscayne Bay. Annual landings varied little during the 1980’s and early 1990’s, and landings of the bait shrimp fishery exceeded those of the food shrimp fishery. The number of trips and landings in both fisheries increased from the late 1990’s through 2002 and food shrimp landings exceeded landings of bait shrimp; landings in both fisheries decreased sharply in 2003. Landings in both fisheries increased in 2004 and 2005, but the increase in food shrimp landings was stronger. Annual catch per trip was much lower in the bait fishery than the food fishery. Each fishery exploited shrimp of a different size. The bait fishery targeted shrimp less than 19 mm carapace length (CL), whereas the food fishery caught shrimp greater than 19 mm CL. We compared monthly bait shrimp catch per unit of effort (CPUE) from the fishery to an estimate of shrimp density from a fishery-independent sampling effort over a 3-yr period and found a strong statistical relationship with the density estimate lagged by 3 mo. The relationship supported the use of bait shrimp fishery CPUE as an index of abundance in upcoming assessments of the effect of a massive water-management-based ecosystem restoration project on pink shrimp in Biscayne Bay. Project implementation will affect freshwater inflows to the bay and salinity patterns. An abundance index with a lengthy pre-implementation history that can be carried into the operational phase of the restoration project will be invaluable in assessing project effects and protecting an important fishery resource of Biscayne Bay. The bait shrimp fishery can provide a continuing index of shrimp abundance from late 1986 forward.

Analyzing and mapping fish assemblages off central California, USA

This study describes fish assemblages and their spatial patterns off the coast of California from Point Arena to Point Sal, by combining the results of the multivariate analyses of several fisheries datasets with a geographic information system. In order to provide comprehensive spatial coverage for the areas of inshore, continental shelf, and continental slope, three fisheries datasets were analyzed: 1) Inshore: the California Department of Fish and Game dataset of fishery-dependent commercial passenger fishing vessel trips that targeted rockfish; 2) Continental Shelf: the National Marine Fisheries Service (NMFS) fishery-independent bottom trawls; and 3) Continental Slope: the NMFS fishery-independent bottom trawls on the continental slope. One-hundred seven species were analyzed. These species represented those captured in at least 5% of the fishing trips or trawls in at least one of the three data sets. We analyzed each of the three datasets separately, and the three sets of results were combined to define 28 species assemblages and 23 site groups. A species assemblage consisted of species caught together, whereas a site group consisted of fishing trips or trawl locations that tended to have the same species assemblages. At the scale of these datasets, 97% of all site groups were significantly segregated by depth.

A sudden collapse in distribution of Pacific sardine (Sardinops sagax) off southwestern Australia enables an objective re-assessment of biomass estimates

Fishery-independent estimates of spawning biomass (BSP) of the Pacific sardine (Sardinops sagax) on the south and lower west coasts of Western Australia (WA) were obtained periodically between 1991 and 1999 by using the daily egg production method (DEPM). Ichthyoplankton data collected during these surveys, specifically the presence or absence of S. sagax eggs, were used to investigate trends in the spawning area of S. sagax within each of four regions. The expectation was that trends in BSP and spawning area were positively related. With the DEPM model, estimates of BSP will change proportionally with spawning area if all other variables remain constant. The proportion of positive stations (PPS), i.e., stations with nonzero egg counts — an objective estimator of spawning area — was high for all south coast regions during the early 1990s (a period when the estimated BSP was also high) and then decreased after the mid-1990s. There was a decrease in PPS from the mid-1990s to 1999. The particularly low estimates in 1999 followed a severe epidemic mass mortality of S. sagax throughout their range across southern Australia. Deviations from the expected relationship between BSP and PPS were used to identify uncertainty around estimates of BSP. Because estimation of spawning area is subject to less sampling bias than estimation of BSP, the deviation in the relation between the two provides an objective basis for adjusting some estimates of the latter. Such an approach is particularly useful for fisheries management purposes when sampling problems are suspected to be present. The analysis of PPS undertaken from the same set of samples from which the DEPM estimate is derived will help provide information for stock assessments and for the management of purse-seine fisheries.

Spawning cycles and habitats for ballyhoo (Hemiramphus brasiliensis) and balao (H. balao) in south Florida

Two halfbeak species, ballyhoo (Hemiramphus brasiliensis) and balao (H. balao), are harvested as bait in south Florida waters, and recent changes in fishing effort and regulations prompted this investigation of the overlap of halfbeak fishing grounds and spawning grounds. Halfbeaks were sampled aboard commercial fishing vessels, and during fishery-independent trips, to determine spatial and temporal spawning patterns of both species. Cyclic patterns of gonadosomatic indices (GSIs) indicated that both species spawned during spring and summer months. Histological analysis demonstrated that specific stages of oocyte development can be predicted from GSI values; for example, female ballyhoo with GSIs >6.0 had hydrated oocytes that were 2.0−3.5 mm diameter. Diel changes in oocyte diameters and histological criteria demonstrated that final oocyte maturation occurred over a 30- to 36-hour period and that ballyhoo spawned at dusk. Hydration of oocytes began in the morning, and ovulation occurred at sunset of that same day; therefore females with hydrated oocytes were ready to spawn within hours. We compared maps of all locations where fish were collected to maps of locations where spawning females (i.e. females with GSIs >6.0) were collected to determine the degree of overlap of halfbeak fishing and spawning grounds. We also used geographic information system (GIS) data to describe the depth and bottom type of halfbeak spawning grounds. Ballyhoo spawned all along the coral reef tract of the Atlantic Ocean, inshore of the reef tract, and in association with bank habitats within Florida Bay. In the Atlantic Ocean, balao spawned along the reef tract and in deeper, more offshore waters than did ballyhoo; balao were not found inshore of the coral reef tract or in Florida Bay. Both halfbeak species, considered together, spawned throughout the fishing grounds of south Florida.

Horizontal and vertical movements of juvenile bluefin tuna (Thunnus thynnus) in relation to oceanographic conditions of the western North Atlantic, determined with ultrasonic telemetry

We employed ultrasonic transmitters to follow (for up to 48 h) the horizontal and vertical movements of five juvenile (6.8–18.7 kg estimated body mass) bluefin tuna (Thunnus thynnus) in the western North Atlantic (off the eastern shore of Virginia). Our objective was to document the fishes’ behavior and distribution in relation to oceanographic conditions and thus begin to address issues that currently limit population assessments based on aerial surveys. Estimation of the trends in adult and juvenile Atlantic bluefin tuna abundance by aerial surveys, and other fishery-independent measures, is considered a priority. Juvenile bluefin tuna spent the majority of their time over the continental shelf in relatively shallow water (generally less then 40 m deep). Fish used the entire water column in spite of relatively steep vertical thermal gradients (≈24°C at the surface and ≈12°C at 40 m depth), but spent the majority of their time (≈90%) above 15 m and in water warmer then 20°C. Mean swimming speeds ranged from 2.8 to 3.3 knots, and total distance covered from 152 to 289 km (82–156 nmi). Because fish generally remained within relatively con-fined areas, net displacement was only 7.7–52.7 km (4.1–28.4 nmi). Horizontal movements were not correlated with sea surface temperature. We propose that it is unlikely that juvenile bluefin tuna in this area can detect minor horizontal temperature gradients (generally less then 0.5°C/km) because of the steep vertical temperature gradients (up to ≈0.6°C/m) they experience during their regular vertical movements. In contrast, water clarity did appear to influence behavior because the fish remained in the intermediate water mass between the turbid and phytoplankton-rich plume exiting Chesapeake Bay (and similar coastal waters) and the clear oligotrophic water east of the continental shelf.

Severe decline in abundance of the red porgy (Pagrus pagrus) population off the southeastern United States

Status of the southeastern U.S. stock of red porgy (Pagrus pagrus) was estimated from fishery-dependent and fishery-independent data, 1972–97. Annual population numbers and fishing mortality rates at age were estimated from virtual population analysis (VPA) calibrated with fishery-independent data. For the VPA, a primary matrix of catch at age was based on age-length keys from fishery-independent samples; an alternate matrix was based on fishery-dependent keys. Additional estimates of stock status were obtained from a surplus-production model, also calibrated with fishery-independent indices of abundance. Results describe a dramatic increase in exploitation of this stock and concomitant decline in abundance. Estimated fully recruited fishing mortality rate (F) from the primary catch matrix increased from 0.10/yr in 1975 to 0.88/yr in 1997, and estimated static spawning potential ratio (SPR) declined from about 67% to about 18%. Estimated recruitment to age 1 declined from a peak of 3.0 million fish in 1973–74 to 94,000 fish in 1997, a decline of 96.9%. Estimated spawning-stock biomass declined from a peak of 3530 t in 1979 to 397 t in 1997, a decline of 88.8%. Results from the alternate catch matrix were similar. Retrospective patterns in the VPA suggest that the future estimates of this population decline will be severe, but may be less than present estimates. Long-term and marked declines in recruitment, spawning stock, and catch per unit of effort (both fishery-derived and fishery-independent)are consistent with severe overexploitation during a period of reduced recruitment. Although F prior to 1995 has generally been estimated at or below the current management criterion for overfishing (F equivalent to SPR=35%), the recent spawning-stock biomass is well below the biomass that could support maximum sustainable yield. Significant reductions in fishing mortality will be needed for rebuilding the southeastern U.S. stock.

Role of the wider Azores region as a nursery ground for North Atlantic blue shark (Prionace glauca)

Tese de Doutoramento, Ciências do Mar (Ecologia Marinha)

Changes in catch rates and length and age at maturity, but not growth, of an estuarine plotosid (cnidoglanis macrocephalus) after heavy fishing

The hypothesis that heavy fishing pressure has led to changes in the biological characteristics of the estuary cobbler (Cnidoglanis macrocephalus) was tested in a large seasonally open estuary in southwestern Australia, where this species completes its life cycle and is the most valuable commercial fish species. Comparisons were made between seasonal data collected for this plotosid (eeltail catfish) in Wilson Inlet during 2005-08 and those recorded with the same fishery-independent sampling regime during 1987-89. These comparisons show that the proportions of larger and older individuals and the catch rates in the more recent period were far lower, i.e., they constituted reductions of 40% for fish ≥430 mm total length, 62% for fish ≥4 years of age, and 80% for catch rate. In addition, total mortality and fishing-induced mortality estimates increased by factors of ~2 and 2.5, respectively. The indications that the abundance and proportion of older C. macrocephalus declined between the two periods are consistent with the perception of long-term commercial fishermen and their shift toward using a smaller maximum gill net mesh to target this species. The sustained heavy fishing pressure on C. macrocephalus between 1987-89 and 2005-08 was accompanied by a marked reduction in length and age at maturity of this species. The shift in probabilistic maturation reaction norms toward smaller fish in 2005-08 and the lack of a conspicuous change in growth between the two periods indicate that the maturity changes were related to fishery-induced evolution rather than to compensatory responses to reduced fish densities.

Aerial surveys to monitor bluefin tuna abundance and track efficiency of management measures

Conservation and management measures for exploited fish species rely on our ability to monitor variations in population abundance. In the case of the eastern stock of Atlantic bluefin tuna (ABFT), recent changes in management policies have strongly affected the reliability of fishery-dependent indicators due to drastic changes in fishing season/area, fisheries selectivity and strategy. However, fishery-independent indices of abundance are rare for large pelagic fish, and obtaining them is often costly and labor intensive. Here, we show that scientific aerial surveys are an appropriate tool for monitoring juvenile bluefin tuna abundance in the Mediterranean. We present an abundance index based on 62 aerial surveys conducted since 2000, using 2 statistical approaches to deal with the sampling strategy: line and strip transects. Both approaches showed a significant increase in juvenile ABFT abundance in recent years, resulting from the recovery plan established in 2007. Nonetheless, the estimates from the line transect method appear to be more robust and stable. This study provides essential information for fisheries management. Expanding the spatial coverage to other nursery grounds would further increase the reliability and representativeness of this index.

A generalized estimating equations approach for analysis of the impact of new technology on a trawl fishery

The article describes a generalized estimating equations approach that was used to investigate the impact of technology on vessel performance in a trawl fishery during 1988-96, while accounting for spatial and temporal correlations in the catch-effort data. Robust estimation of parameters in the presence of several levels of clustering depended more on the choice of cluster definition than on the choice of correlation structure within the cluster. Models with smaller cluster sizes produced stable results, while models with larger cluster sizes, that may have had complex within-cluster correlation structures and that had within-cluster covariates, produced estimates sensitive to the correlation structure. The preferred model arising from this dataset assumed that catches from a vessel were correlated in the same years and the same areas, but independent in different years and areas. The model that assumed catches from a vessel were correlated in all years and areas, equivalent to a random effects term for vessel, produced spurious results. This was an unexpected finding that highlighted the need to adopt a systematic strategy for modelling. The article proposes a modelling strategy of selecting the best cluster definition first, and the working correlation structure (within clusters) second. The article discusses the selection and interpretation of the model in the light of background knowledge of the data and utility of the model, and the potential for this modelling approach to apply in similar statistical situations.