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.

Characterization of the U.S. Gulf of Mexico and South Atlantic penaeid and rock shrimp fisheries based on observer data

In July 2007, a mandatory Federal observer program was implemented to characterize the U.S. Gulf of Mexico penaeid shrimp (Farfantepenaeus aztecus, F. duorarum, and Litopenaeus setiferus) fishery. In June 2008, the program expanded to include the South Atlantic penaeid and rock shrimp, Sicyonia spp., fisheries. Data collected from 10,206 tows during 5,197 sea days of observations were analyzed by geographical area and target species. The majority of tows (~70%) sampled were off the coasts of Texas and Louisiana. Based on total hours towed, the highest concentrated effort occurred off South Texas and southwestern Florida. Gear information, such as net characteristics, bycatch reduction devices, and turtle excluder devices were fairly consistent among areas and target species. By species categories, finfish comprised the majority (≥57%) of the catch composition in the Gulf of Mexico and South Atlantic penaeid shrimp fisheries, while in the South Atlantic rock shrimp fishery the largest component (41%) was rock shrimp. Bycatch to shrimp ratios were lower than reported in previous studies for the Gulf of Mexico penaeid shrimp fishery. These decreased ratios may be attributed to several factors, notably decreased shrimp effort and higher shrimp catch per unit of effort (CPUE) in recent years. CPUE density surface plots for several species of interest illustrated spatial differences in distribution. Hot Spot Analyses for shrimp (penaeid and rock) and bycatch species identified areas with significant clustering of high or low CPUE values. Spatial and temporal distribution of protected species interactions were documented.

Drought frequency analysis for California from observed, synthetic, and proxy data

Drought frequency analysis can be performed with statistical techniques developed for determining recurrence intervals for extreme precipitation and flood events (Linsley et al 1992). The drought analysis method discussed in this paper uses the log-Pearson Type III distribution, which has been widely used in flood frequency research. Some of the difficulties encountered when using this distribution for drought analysis are investigated.

Year-to-year changes of vertical temperature distribution in the California Current region: 1954 to 1986

Studies by Enfield and Allen (1980), McLain et al (1985), and others have shown that anomalously warm years in the northern coastal California Current correspond to El Niño conditions in the equatorial Pacific Ocean. Ocean model studies suggest a mechanical link between the northern coastal California Current and the equatorial ocean through long waves that propagate cyclonically along the ocean boundary (McCreary 1976; Clarke 1983; Shriver et al 1991). However, distinct observational evidence of such an oceanic connection is not extensive. Much of the supposed El Niño variation in temperature and sea level data from the coastal California Current region can be associated with the effects of anomalously intense north Pacific atmospheric cyclogenesis, which is frequently augmented during El Niño years (Wallace and Gutzler 1981; Simpson 1983; Emery and Hamilton 1984). This study uses time series of ocean temperature data to distinguish between locally forced effects, initiated by north Pacific atmospheric changes, and remotely forced effects, initiated by equatorial Pacific atmospheric changes related to El Niño events.

Distribution and potential bioeffects of atrazine in coastal waters

Estuaries provide critical nursery habitat for many commercially and recreationally important fish and shellfish species. These productive, diverse ecosystems are particularly vulnerable to pollution because they serve as repositories for non–point-source contaminants from upland sources, such as pesticide runoff. Atrazine, among the most widely used pesticides in the United States, has also been one of the most extensively studied. There has not, however, been a specific assessment of atrazine in marine and estuarine ecosystems. This document characterizes the presence and transformation of atrazine in coastal waters, and the effects of atrazine on marine organisms. Review of marine and estuarine monitoring data indicate that atrazine is chronically present in U.S. coastal waters at relatively low concentrations. The concentrations detected have typically been below acute biological effects levels, and below the U.S. EPA proposed water quality criteria for atrazine. While direct risk of atrazine impacts are low, uncertainty remains regarding the effects of long-term low levels of atrazine in mixture with other contaminants. It is recommended that best management practices, such as the use of vegetative buffers and public education about pesticide use, be encouraged in the coastal zone to minimize runoff of atrazine into marine and estuarine waters.

The distribution of larval fishes of the Charleston Gyre Region off the southeastern United States in winter shaped by mesoscale, cyclonic eddies

Serial, cyclonic, mesoscale eddies arise just north of the Charleston Bump, a topographical rise on the continental slope and Blake Plateau, and characterize the U.S. outer shelf and upper slope in the region of the Charleston Gyre. This region was transected during the winters of 2000, 2001, and 2002, and hydrographic data and larval fishes were collected. The hydrodynamics of the cyclonic eddies of the Charleston Gyre shape the distribution of larval fishes by mixing larvae from the outer continental shelf and the Gulf Stream and entraining them into the eddy circulation at the peripheral margins, the wrap-around filaments. Over all years and transects (those that intercepted eddies and those that did not), chlorophyll a concentrations, zooplankton displacement volumes, and larval fish concentrations were positively correlated. Chlorophyll a concentrations were highest in filaments that wrapped around eddies, and zooplankton displacement volumes were highest in the continental shelf–Gulf Stream–frontal mix. Overall, the concentration of all larval fishes declined from inshore to offshore with highest concentrations occurring over the outer shelf. Collections produced larvae from 91 fish families representing continental shelf and oceanic species. The larvae of shelf-spawned fishes—Atlantic Menhaden Brevoortia tyrannus, Round Herring Etrumeus teres, Spot Leiostomus xanthurus, and Atlantic Croaker Micropogonias undulatus—were most concentrated over the outer shelf and in the continental shelf–Gulf Stream–frontal mix. The larvae of ocean-spawned fishes—lanternfishes, bristlemouths, and lightfishes—were more evenly dispersed in low concentrations across the outer shelf and upper slope, the highest typically in the Gulf Stream and Sargasso Sea, except for lightfishes that were highest in the continental shelf–Gulf Stream–frontal mix. Detrended correspondence analysis rendered groups of larval fishes that corresponded with a gradient between the continental shelf and Gulf Stream and Sargasso Sea. Eddies propagate northeastward with a residence time on the outer shelf and upper slope of ∼1 month, the same duration as the larval period of most fishes. The pelagic habitat afforded by eddies and fronts of the Charleston Gyre region can be exploited as nursery areas for feeding and growth of larval fishes within the southeastern Atlantic continental shelf ecosystem of the U.S. Eddies, and the nursery habitat they provide, translocate larvae northeastward.

Integration of data for nowcasting of harmful algal blooms

Harmful algal blooms (HABs) are a significant and potentially expanding problem around the world. Resource management and public health protection require sufficient information to reduce the impacts of HABs by response strategies and through warnings and advisories. To be effective, these programs can best be served by an integration of improved detection methods with both evolving monitoring systems and new communications capabilities. Data sets are typically collected from a variety of sources, these can be considered as several types: point data, such as water samples; transects, such as from shipboard continuous sampling; and synoptic, such as from satellite imagery. Generation of a field of the HAB distribution requires all of these sampling approaches. This means that the data sets need to be interpreted and analyzed with each other to create the field or distribution of the HAB. The HAB field is also a necessary input into models that forecast blooms. Several systems have developed strategies that demonstrate these approaches. These range from data sets collected at key sites, such as swimming beaches, to automated collection systems, to integration of interpreted satellite data. Improved data collection, particularly in speed and cost, will be one of the advances of the next few years. Methods to improve creation of the HAB field from the variety of data types will be necessary for routine nowcasting and forecasting of HABs.

An ecological characterization of the marine resources of Vieques, Puerto Rico. Part I: Historical data synthesis

From the 1940s until 2003, portions of the island of Vieques, a municipality within the Commonwealth of Puerto Rico, were used by the US Navy as a base and training facility, resulting in development and zoning history that differ in comparison to other Caribbean islands. The majority of former Navy lands are now under the jurisdiction of the Department of the Interior’s Fish and Wildlife Service as a National Wildlife Refuge, while a smaller percentage of land was transferred to the Vieques municipality and the Puerto Rico Conservation Trust. An analysis of the distribution and status of the marine resources is timely in light of the recent land transfer, increases in development and tourism, and potential changes in marine zoning around the island. To meet this need, NOAA’s Biogeography Branch, in cooperation with the Office of Response and Restoration and other local and regional partners, conducted Part I of an ecological characterization to integrate historical data and research into a synthesis report. The overall objective of this report is to provide resource managers and residents a comprehensive characterization of the marine resources of Vieques to support research, monitoring, and management. For example, knowledge of the spatial distribution of physical features, habitats, and biological communities is necessary to make an informed decision of the establishment and placement of a marine protected area (MPA). The report is divided into chapters based on the physical environment (e.g., climate, geology, bathymetry), habitat types (e.g., reefs and hardbottom, seagrasses, mangroves) and major faunal groups (e.g. fish, turtles, birds). Each section includes five subsections: an overview, description of the relevant literature, methods of analysis, information on the distribution, status and trends of the particular resource, and a discussion of ecological linkages with other components of the Vieques marine ecosystem and surrounding environment. The physical environment of Vieques is similar to other islands within the Greater Antilles chain, with some distinctions. The warm, tropical climate of Vieques, mediated by the northeasterly trade winds, is characterized by a dry season (December-April) and a rainy season (May-November), the latter of which is characterized by the occasional passage of tropical cyclones. Compared to mainland Puerto Rico, Vieques is characterized by lower elevation, less annual precipitation, and higher average temperatures. The amount of annual precipitation also varies spatially within Vieques, with the western portion of the island receiving higher amounts of rainfall than further east. While the North Equatorial Current dominates the circulation pattern in the Greater Antilles region, small scale current patterns specific to Vieques are not as well characterized. These physical processes are important factors mitigating the distribution and composition of marine benthic habitats around Vieques. In general, the topography of Vieques is characterized by rolling hills. Mt. Pirata, the tallest point at 301 m, is located near the southwest coast. In the absence of island wide sedimentation measurements, information on land cover, slope, precipitation, and soil type were used to estimate relative erosion potential and sediment delivery for each watershed. While slope and precipitation amount are the primary driving factors controlling runoff, land use practices such as urban development, military activity, road construction, and agriculture can increase the delivery of pollution and sediments to coastal waters. Due to the recent land transfer, increased development and tourism is expected, which may result in changes in the input of sediments to the coastal environment.

Characterizing a cyanobacterial bloom in western Lake Erie using satellite imagery and meteorological data

The distribution and intensity of a bloom of the toxic cyanobacterium, Microcystis aeruginosa, in western Lake Erie was characterized using a combination of satellite ocean-color imagery, field data, and meteorological observations. The bloom was first identified by satellite on 14 August 2008 and persisted for more than 2 months. The distribution and intensity of the bloom was estimated using a satellite algorithm that is sensitive to near-surface concentrations of M. aeruginosa. Increases in both area and intensity were most pronounced for wind stress less than 0.05 Pa. Area increased while intensity did not change for wind stresses of 0.05–0.1 Pa, and both decreased for wind stress greater than 0.1 Pa. The recovery in intensity at the surface after strong wind events indicated that high wind stress mixed the bloom through the water column and that it returned to the surface once mixing stopped. This interaction is consistent with the understanding of the buoyancy of these blooms. Cloud cover (reduced light) may have a weak influence on intensity during calm conditions. While water temperature remained greater than 15°C, the bloom intensified if there were calm conditions. For water temperature less than 15°C, the bloom subsided under similar conditions. As a result, wind stress needs to be considered when interpreting satellite imagery of these blooms.

Distribution and biology of prowfish (Zaprora silenus) in the northeast Pacific

The prowfish (Zaprora silenus) is an infrequent component of bottom trawl catches collected on stock assessment surveys. Based on presence or absence in over 40,000 trawl catches taken throughout Alaskan waters southward to southern California, prowfish are most frequently encountered in the Gulf of Alaska and the Aleutian Islands at the edge of the continental shelf. Based on data from two trawl surveys, relative abundance indicated by catch per swept area reaches a maximum between 100 m and 200 m depth and is much higher in the Aleutian Islands than in the Gulf of Alaska. Females weigh 3.7% more than males of the same length. Weight-length functions are W (g) = 0.0164 L2.92 (males) and W = 0.0170 L2.92 (females). Length at age does not differ between sexes and is described by L = 89.3(1 – e–0.181(t+0.554)), where L is total length in cm and t is age in years. Females reached 50% maturity at a length of 57.0 cm and an age of 5.1 years. Prowfish diet is almost entirely composed of gelatinous zooplankton, primarily scyphozoa and salps.

Distribution and abundance of humpback whales (Megaptera novaeangliae) and other marine mammals off the northern Washington coast

We examined the summer distribution of marine mammals off the northern Washington coast based on six ship transect surveys conducted between 1995 and 2002, primarily from the NOAA ship McArthur. Additionally, small boat surveys were conducted in the same region between 1989 and 2002 to gather photographic identification data on humpback whales (Megaptera novaeangliae) and killer whales (Orcinus orca) to examine movements and population structure. In the six years of ship survey effort, 706 sightings of 15 marine mammal species were made. Humpback whales were the most common large cetacean species and were seen every year and a total of 232 sightings of 402 animals were recorded during ship surveys. Highest numbers were observed in 2002, when there were 79 sightings of 139 whales. Line-transect estimates for humpback whales indicated that about 100 humpback whales inhabited these waters each year between 1995 and 2000; in 2002, however, the estimate was 562 (CV= 0.21) whales. A total of 191 unique individuals were identified photographically and mark recapture estimates also indicated that the number of animals increased from under 100 to over 200 from 1995 to 2002. There was only limited interchange of humpback whales between this area and feeding areas off Oregon and California. Killer whales were also seen on every ship survey and represented all known ecotypes of the Pacific Northwest, including southern and northern residents, transients, and offshore-type killer whales. Dall’s porpoise (Phocoenoides dalli) were the most frequently sighted small cetacean; abundance was estimated at 181−291 individuals, except for 2002 when we observed dramatically higher numbers (876, CV= 0.30). Northern fur seals (Callorhinus ursinus) and elephant seals (Mirounga angustirostris) were the most common pinnipeds observed. There were clear habitat differences related to distance offshore and water depth for different species.

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.

The statistical properties of recreational catch rate data for some fish stocks off the northeast U.S. coast

Recreational fisheries in the waters off the northeast U.S. target a variety of pelagic and demersal fish species, and catch and effort data sampled from recreational fisheries are a critical component of the information used in resource evaluation and management. Standardized indices of stock abundance developed from recreational fishery catch rates are routinely used in stock assessments. The statistical properties of both simulated and empirical recreational fishery catch-rate data such as those collected by the National Marine Fisheries Service (NMFS) Marine Recreational Fishery Statistics Survey (MRFSS) are examined, and the potential effects of different assumptions about the error structure of the catch-rate frequency distributions in computing indices of stock abundance are evaluated. Recreational fishery catch distributions sampled by the MRFSS are highly contagious and overdispersed in relation to the normal distribution and are generally best characterized by the Poisson or negative binomial distributions. The modeling of both the simulated and empirical MRFSS catch rates indicates that one may draw erroneous conclusions about stock trends by assuming the wrong error distribution in procedures used to developed standardized indices of stock abundance. The results demonstrate the importance of considering not only the overall model fit and significance of classification effects, but also the possible effects of model misspecification, when determining the most appropriate model construction.

Changes over time in the spatial distribution of walleye pollock (Theragra chalcogramma) in the Gulf of Alaska, 1984-1996

Triennial bottom trawl survey data from 1984 to 1996 were used to evaluate changes in the summer distribution of walleye pollock in the western and central Gulf of Alaska. Differences between several age groups of pollock were evaluated. Distribution was examined in relation to several physical characteristics, including bottom depth and distance from land. Interspecies associations were also analyzed with the Bray-Curtis clustering technique to better understand community structure. Our results indicated that although the population numbers decreased, high concentrations of pollock remained in the same areas during 1984–96. However, there was an increase in the number of stations where low-density pollock concentrations of all ages were observed, which resulted in a decrease in mean population density of pollock within the GOA region. Patterns emerging from our data suggested an alternative to Mac-Call’s “basin hypothesis” which states that as population numbers decrease, there should be a contraction of the population range to optimal habitats. During 1984–96 there was a concurrent precipitous decline in Steller sea lions in the Gulf of Alaska. The results of our study suggest that decreases in the mean density of adult pollock, the main food in the Steller sea lion diet, combined with slight changes in the distribution of pollock (age-1 pollock in particular) in the mid-1980s, may have contributed to decreased foraging efficiency in Steller sea lions. Our results support the prevailing conceptual model for pollock ontogeny, although there is evidence that substantial spawning may also occur outside of Shelikof Strait.

Age, growth, mortality, and distribution of pinfish (Lagodon rhomboides) in Tampa Bay and adjacent Gulf of Mexico waters

Life-history dynamics of pinfish (Lagodon rhomboides) were examined from data derived from random station surveys conducted in Tampa Bay and adjacent Gulf of Mexico waters during 1993–97. In addition, patterns in spatial distribution and abundance in Gulf of Mexico waters were investigated. Ages determined from whole otoliths ranged from 0 to 7 years, and von Bertalanffy growth models for males and females were not significantly different. Von Bertalanffy growth model parameters were L∞=219.9 mm SL, k =0.33/yr, and t0 =–1.10 years for all fish combined. High gonadosomatic indices during October–December indicated that some spawning may occur in Tampa Bay. Estimated lengths at 50% maturity were 132 mm SL for males and 131 mm SL for females. Total instantaneous mortality rates derived from the Chapman-Robson estimator ranged from 0.88 to 1.08/yr, and natural mortality was estimated to be 0.78/yr. In Gulf of Mexico waters, pinfish catch rates declined with increasing depth, and most pinfish were caught in <17 m of water. Length distributions showed that pinfish segregate by size with increasing depth.