Population dynamics of the skipjack tuna (Katsuwonus pelamis) of the Eastern Pacific Ocean

ENGLISH: Catches of skipjack tuna supporting major fisheries in parts of the western, central and eastern Pacific Ocean have increased in recent years; thus, it is important to examine the dynamics of the fishery to determine man's effect on the abundance of the stocks. A general linear hypothesis model was developed to standardize fishing effort to a single vessel size and gear type. Standardized effort was then used to compute an index of abundance which accounts for seasonal variability in the fishing area. The indices of abundance were highly variable from year to year in both the northern and southern areas of the fishery but indicated a generally higher abundance in the south. Data from 438 fish tagged and recovered in the eastern Pacific Ocean were used to compute growth curves. A least-squares technique was used to estimate the parameters of the von Bertalanffy growth function. Two estimates of the parameters were made by analyzing the same data in different ways. For the first set of estimates, K= 0.819 on an annual instantaneous basis and L= 729 mm; for the second, K = 0.431 and L=881. These compared well with estimates derived using the Chapman-Richards growth function, which includes the von Bertalanffy function as a special case. It was concluded that the latter function provided an adequate empirical fit to the skipjack data since the more complicated function did not significantly improve the fit. Tagging data from three cruises involving 8852 releases and 1777 returns were used to compute mortality rates during the time the fish were in the fishery. Two models were used in the analyses. The best estimates of the catchability coefficient (q) in the north and south were 8.4 X 10- 4 and 5.0 X 10- 5 respectively. The other loss rate (X), which included losses due to emigration, natural mortality and mortality due to carrying a tag, was 0.14 on an annual instantaneous basis for both areas. To detect the possible effect of fishing on abundance and total yield, the relation between abundance and effort and between total catch and effort was examined. It was found that at levels of intensity observed in the fishery, fishing does not appear to have had any measurable effect on the stocks. It was concluded therefore that the total catch could probably be increased by substantially increasing total effort beyond the present level, and that the fluctuations in abundance are fishery-independent. The estimates of growth, mortality and fishing effort were used to compute yield-per-recruitment isopleths for skipjack in both the northern and southern areas. For a size at first entry of about 425 mm, the yield per recruitment was calculated at 3 pounds in the north and 1.5 pounds in the south. In both areas it would be possible to increase the yield per recruitment by increasing fishing effort. It was not possible to assess potential production of the skipjack stocks fished in the eastern Pacific, except to note that the fishery had not affected their abundance and that they were certainly under-exploited. It was concluded that the northern and southern stocks could support increased harvests, especially the latter. SPANISH: Las capturas de atún barrilete que sostienen las pesquerías principales de la parte occidental, central y oriental del Océano Pacífico han aumentado en los últimos años; así que es importante examinar la dinámica de la pesquería para determinar el efecto que pueda tener sobre la abundancia de los stocks. Se desarrolló un modelo hipotético, lineal para standardizar el esfuerzo de pesca a un solo tamaño de barco y tipo de arte. Luego se usó el esfuerzo standardizado para computar un índice de la abundancia que pueda dar razón de la variabilidad estacional en el área de pesca. Los índices de la abundancia variaron mucho de un año a otro tanto en el área septentrional como en el área meridional de la pesquería, pero indicaron una abundancia generalmente superior en el sur. Se emplearon los datos de 438 peces marcados y recuperados en el Océano Pacífico oriental para computar las curvas de crecimiento. Una técnica de mínimos cuadrados fue usada para estimar los parámetros de la función de crecimiento de van Bertalanffy. Se hicieron dos estimativos de los parámetros mediante el análisis de los mismos datos, de diferente manera. Para el primer juego de estimativos, K=0.819 sobre una base anual instantánea y L∞=729 mm; para el segundo, K=0.431 y L∞=881. Estos se correlacionaron bien con los estimativos obtenidos usando la función de crecimiento de Chapman-Richards, que incluye la de von Bertalanffy como un caso especial. Se decidió que la última función proveía un ajuste empírico, adecuado a los datos del barrilete, ya que la función más complicada no mejoró significativamente el ajuste. Los datos de marcación de tres cruceros incluyendo 8852 liberaciones y 1777 retornos, fueron usados para computar las tasas de mortalidad durante el tiempo en que los peces estuvieron en la pesquería. Se usaron dos modelos en los análisis. Los mejores estimativos del coeficiente de capturabilidad (q) en el norte y en el sur fueron 8.4 X 10-4 y 5.0 X 10-5 , respectivamente. La otra tasa de pérdida (X), la cual incluyó pérdidas debidas a la emigración, mortalidad natural y mortalidad debida a llevar una marca, fue 0.14 sobre una base anual instantánea para las dos áreas. Con el fin de descubrir el efecto que posiblemente pueda tener la pesca sobre la abundancia y el rendimiento total, se examinó la relación entre la abundancia y el esfuerzo y entre la captura total y el esfuerzo. Se encontró que a los niveles de la intensidad observada en la pesquería, la pesca no parece haber tenido ningún efecto perceptible en los stocks. Por lo tanto se decidió que mediante un aumento substancial del esfuerzo total, más allá del nivel actual, la captura total probablemente podría aumentarse, y que las fluctuaciones de la abundancia son independientes de la pesquería. Los estimativos del crecimiento, mortalidad y esfuerzo de pesca fueron usados para computar las isopletas del rendimiento por recluta del barrilete, tanto en las áreas del norte como del sur. Para una talla de primera entrada de unos 425 mm, el rendimiento por recluta fue calculado en 3 libras en el norte y 1.5 libras en el sur. En ambas áreas sería posible aumentar el rendimiento por recluta mediante un aumento del esfuerzo de pesca. No fue posible determinar la producción potencial de los stocks del barrilete pescado en el Pacífico oriental, excepto para observar que la pesquería no ha afectado su abundancia y que ciertamente se encuentran subexplotados. Se concluyó que los stocks norte y sur pueden soportar un aumento en el rendimiento, especialmente este último. (PDF contains 274 pages.)

Length-weight Relationships of Dolphinfish, Coryphaena hippurus, and Wahoo, Acanthocybium solandri: Seasonal Effects of Spawning and Possible Migration in the Central North Pacific

Weight-on-length (W-L) relationships for 2,482 dolphinfish, Coryphaena hippurus, and 1,161 wahoo, Acanthocybium solandri, were examined. Data on fork length, whole (round) weight, and sex were collected for dolphinfish at the Honolulu fish auction from March 1988 through November 1989. Unsexed weight and length data for wahoo were collected at the auction from July 1988 through November 1989. We also used sex specific weight and length data of 171 wahoo collected during 1977–1985 research cruises for analysis. Coefficients of W-L regressions were significantly different between the sexes for dolphinfish. Coefficients did not significantly differ between the sexes for wahoo based on research cruise data. In a general linear model evaluating month as a categorical factor, month was significant for female dolphinfish, male dolphinfish, and wahoo with sexes pooled. W-L and length-on-weight (L-W) relationships were fitted by nonlinear regression for all dolphinfish, female dolphinfish, male dolphinfish, and all wahoo sexes pooled. W-L relationships for monthly samples of female dolphinfish, male dolphinfish, and all wahoo with sexes pooled were also fitted by nonlinear regression. Predicted mean weight at length for wahoo was highest at the beginning of the spawning season in June and lowest after the spawning season in September. Maximum and minimum predicted mean weight at length for both sexes of dolphinfish did not correspond with the peak spawning period (March–May). Plausible migration models in conjunction with reproductive behavior were examined to explain the variability in monthly predicted mean weight at length for dolphinfish.

A Comparison of Two Methods for Enhancing the Recovery of Seagrasses into Propellor Scars: Mechanical Injection of a Nutrient and Growth Hormone Solution vs. Defecation by Roosting Seabirds: Final Report.

Based on the recovery rates for Thalassia testudinum measured in this study for scars of these excavation depths and assuming a linear recovery horizon, we estimate that it would take ~ 6.9 years (95% CI. = 5.4 to 9.6 years) for T. testudinum to return to the same density as recorded for the adjacent undisturbed population. The application of water soluble fertilizers and plant growth hormones by mechanical injection into the sediments adjacent to ten propellor scars at Lignumvitae State Botanical Site did not significantly increase the recovery rate of Thalassia testudinum or Halodule wrightii. An alternative method of fertilization and restoration of propellor scars was also tested by a using a method of “compressed succession” where Halodule wrightii is substituted for T. testudinum in the initial stages of restoration. Bird roosting stakes were placed among H.wrightii bare root plantings in prop scars to facilitate the defecation of nitrogen and phosphorus enriched feces. In contrast to the fertilizer injection method, the bird stakes produced extremely high recovery rates of transplanted H. wrightii. We conclude that use of a fertilizer/hormone injection machine in the manner described here is not a feasible means of enhancing T. testudinum recovery in propellor scars on soft bottom carbonate sediments. Existing techniques such as the bird stake approach provide a reliable, and inexpensive alternative method that should be considered for application to restoration of seagrasses in these environments. Document contains 40 pages)

Ecological condition of coastal ocean waters along the U.S. Western Continental Shelf: 2003

Executive Summary: The western National Coastal Assessment (NCA-West) program of EPA, in conjunction with the NOAA National Ocean Service (NOS), conducted an assessment of the status of ecological condition of soft sediment habitats and overlying waters along the western U.S. continental shelf, between the target depths of 30 and 120 m, during June 2003. NCA-West and NOAA/NOS partnered with the West Coast states (Washington (WA), Oregon (OR), and California (CA)), and the Southern California Coastal Water Research Project (SCCWRP) Bight ’03 program to conduct the survey. A total of 257 stations were sampled from Cape Flattery, WA to the Mexican border using standard methods and indicators applied in previous coastal NCA projects. A key study feature was the incorporation of a stratified-random sampling design with stations stratified by state and National Marine Sanctuary (NMS) status. Each of the three states was represented by at least 50 random stations. There also were a total of 84 random stations located within NOAA’s five NMSs along the West Coast including the Olympic Coast NMS (OCNMS), Cordell Bank NMS (CBNMS), Gulf of Farallones NMS (GFNMS), Monterey Bay NMS (MBNMS), and Channel Islands NMS (CINMS). Collection of flatfish via hook-and-line for fish-tissue contaminant analysis was successful at 50 EMAP/NCA-West stations. Through a collaboration developed with the FRAM Division of the Northwest Fisheries Science Center, fish from an additional 63 stations in the same region and depth range were also analyzed for fish-tissue contaminants. Bottom depth throughout the region ranged from 28 m to 125 m for most stations. Two slightly deeper stations from the Southern California Bight (SCB) (131, 134 m) were included in the data set. About 44% of the survey area had sediments composed of sands (< 20% silt-clay), about 47% was composed of intermediate muddy sands (20-80% silt-clay), and about 9% was composed of muds (> 80% silt-clay). The majority of the survey area (97%) had relatively low percent total organic carbon (TOC) levels of < 2%, while a small portion (< 1%) had high TOC levels (> 5%), in a range potentially harmful to benthic fauna. Salinity of surface waters for 92% of the survey area were > 31 psu, with most stations < 31 psu associated with the Columbia River plume. Bottom salinities ranged only between 31.6 and 34.4 psu. There was virtually no difference in mean bottom salinities among states or between NMS and non-NMS stations. Temperatures of surface water (range 8.5 -19.9 °C) and bottom water (range 5.8 -14.7 °C) averaged several degrees higher in CA in comparison to WA and OR. The Δσt index of watercolumn stratification indicated that about 31% of the survey area had strong vertical stratification of the water column. The index was greatest for waters off WA and lowest for CA waters. Only about 2.6 % of the survey area had surface dissolved oxygen (DO) concentrations ≤ 4.8 mg/L, and there were no values below the lower threshold (2.3 mg/L) considered harmful to the survival and growth of marine animals. Surface DO concentrations were higher in WA and OR waters than in CA, and higher in the OC NMS than in the CA sanctuaries. An estimated 94.3% of the area had bottom-water DO concentrations ≤ 4.8 mg/L and 6.6% had concentrations ≤ 2.3 mg/L. The high prevalence of DO from 2.3 to 4.8 mg/L (85% of survey area) is believed to be associated with the upwelling of naturally low DO water across the West Coast shelf. Mean TSS and transmissivity in surface waters (excluding OR due to sample problems) were slightly higher and lower, respectively, for stations in WA than for those in CA. There was little difference in mean TSS or transmissivity between NMS and non-NMS locations. Mean transmissivity in bottom waters, though higher in comparison to surface waters, showed little difference among geographic regions or between NMS and non-NMS locations. Concentrations of nitrate + nitrite, ammonium, total dissolved inorganic nitrogen (DIN) and orthophosphate (P) in surface waters tended to be highest in CA compared to WA and OR, and higher in the CA NMS stations compared to CA non-sanctuary stations. Measurements of silicate in surface waters were limited to WA and CA (exclusive of the SCB) and showed that concentrations were similar between the two states and approximately twice as high in CA sanctuaries compared to OCNMS or nonsanctuary locations in either state. The elevated nutrient concentrations observed at CA NMS stations are consistent with the presence of strong upwelling at these sites at the time of sampling. Approximately 93% of the area had DIN/P values ≤ 16, indicative of nitrogen limitation. Mean DIN/P ratios were similar among the three states, although the mean for the OCNMS was less than half that of the CA sanctuaries or nonsanctuary locations. Concentrations of chlorophyll a in surface waters ranged from 0 to 28 μg L-1, with 50% of the area having values < 3.9 μg L-1 and 10% having values > 14.5 μg L-1. The mean concentration of chlorophyll a for CA was less than half that of WA and OR locations, and concentrations were lowest in non-sanctuary sites in CA and highest at the OCNMS. Shelf sediments throughout the survey area were relatively uncontaminated with the exception of a group of stations within the SCB. Overall, about 99% of the total survey area was rated in good condition (<5 chemicals measured above corresponding effect range low (ERL) concentrations). Only the pesticides 4,4′-DDE and total DDT exceeded corresponding effect range-median (ERM) values, all at stations in CA near Los Angeles. Ten other contaminants including seven metals (As, Cd, Cr, Cu, Hg, Ag, Zn), 2-methylnaphthalene, low molecular weight PAHs, and total PCBs exceeded corresponding ERLs. The most prevalent in terms of area were chromium (31%), arsenic (8%), 2-methylnaphthalene (6%), cadmium (5%), and mercury (4%). The chromium contamination may be related to natural background sources common to the region. The 2-methylnaphthalene exceedances were conspicuously grouped around the CINMS. The mercury exceedances were all at non-sanctuary sites in CA, particularly in the Los Angeles area. Concentrations of cadmium in fish tissues exceeded the lower end of EPA’s non-cancer, human-health-risk range at nine of 50 EMAP/NCA-West and nine of 60 FRAM groundfish-survey stations, including a total of seven NMS stations in CA and two in the OCNMS. The human-health guidelines for all other contaminants were only exceeded for total PCBs at one station located in WA near the mouth of the Columbia River. Benthic species richness was relatively high in these offshore assemblages, ranging from 19 to 190 taxa per 0.1-m2 grab and averaging 79 taxa/grab. The high species richness was reflected over large areas of the shelf and was nearly three times greater than levels observed in estuarine samples along the West Coast (e.g NCA-West estuarine mean of 26 taxa/grab). Mean species richness was highest off CA (94 taxa/grab) and lower in OR and WA (55 and 56 taxa/grab, respectively). Mean species richness was very similar between sanctuary vs. non-sanctuary stations for both the CA and OR/WA regions. Mean diversity index H′ was highest in CA (5.36) and lowest in WA (4.27). There were no major differences in mean H′ between sanctuary vs. nonsanctuary stations for both the CA and OR/WA regions. A total of 1,482 taxa (1,108 to species) and 99,135 individuals were identified region-wide. Polychaetes, crustaceans and molluscs were the dominant taxa, both by percent abundance (59%, 17%, 12% respectively) and percent species (44%, 25%, 17%, respectively). There were no major differences in the percent composition of benthic communities among states or between NMSs and corresponding non-sanctuary sites. Densities averaged 3,788 m-2, about 30% of the average density for West Coast estuaries. Mean density of benthic fauna in the present offshore survey, averaged by state, was highest in CA (4,351 m-2) and lowest in OR (2,310 m-2). Mean densities were slightly higher at NMS stations vs. non-sanctuary stations for both the CA and OR/WA regions. The 10 most abundant taxa were the polychaetes Mediomastus spp., Magelona longicornis, Spiophanes berkeleyorum, Spiophanes bombyx, Spiophanes duplex, and Prionospio jubata; the bivalve Axinopsida serricata, the ophiuroid Amphiodia urtica, the decapod Pinnixa occidentalis, and the ostracod Euphilomedes carcharodonta. Mediomastus spp. and A. serricata were the two most abundant taxa overall. Although many of these taxa have broad geographic distributions throughout the region, the same species were not ranked among the 10 most abundant taxa consistently across states. The closest similarities among states were between OR and WA. At least half of the 10 most abundant taxa in NMSs were also dominant in corresponding nonsanctuary waters. Many of the abundant benthic species have wide latitudinal distributions along the West Coast shelf, with some species ranging from southern CA into the Gulf of Alaska or even the Aleutians. Of the 39 taxa on the list of 50 most abundant taxa that could be identified to species level, 85% have been reported at least once from estuaries of CA, OR, or WA exclusive of Puget Sound. Such broad latitudinal and estuarine distributions are suggestive of wide habitat tolerances. Thirteen (1.2%) of the 1,108 identified species are nonindigenous, with another 121 species classified as cryptogenic (of uncertain origin), and 208 species unclassified with respect to potential invasiveness. Despite uncertainties of classification, the number and densities of nonindigenous species appear to be much lower on the shelf than in the estuarine ecosystems of the Pacific Coast. Spionid polychaetes and the ampharetid polychaete Anobothrus gracilis were a major component of the nonindigenous species collected on the shelf. NOAA’s five NMSs along the West Coast of the U.S. appeared to be in good ecological condition, based on the measured indicators, with no evidence of major anthropogenic impacts or unusual environmental qualities compared to nearby nonsanctuary waters. Benthic communities in sanctuaries resembled those in corresponding non-sanctuary waters, with similarly high levels of species richness and diversity and low incidence of nonindigenous species. Most oceanographic features were also similar between sanctuary and non-sanctuary locations. Exceptions (e.g., higher concentrations of some nutrients in sanctuaries along the CA coast) appeared to be attributable to natural upwelling events in the area at the time of sampling. In addition, sediments within the sanctuaries were relatively uncontaminated, with none of the samples having any measured chemical in excess of ERM values. The ERL value for chromium was exceeded in sediments at the OCNMS, but at a much lower percentage of stations (four of 30) compared to WA and OR non-sanctuary areas (31 of 70 stations). ERL values were exceeded for arsenic, cadmium, chromium, 2- methylnaphthalene, low molecular weight PAHs, total DDT, and 4,4′-DDE at multiple sites within the CINMS. However, cases where total DDT, 4,4′-DDE, and chromium exceeded the ERL values were notably less prevalent at CINMS than in non-sanctuary waters of CA. In contrast, 2-methylnaphthalene above the ERL was much more prevalent in sediments at the CINMS compared to non-sanctuary waters off the coast of CA. While there are natural background sources of PAHs from oil seeps throughout the SCB, this does not explain the higher incidence of 2-methylnaphthalene contamination around CINMS. Two stations in CINMS also had levels of TOC (> 5%) potentially harmful to benthic fauna, though none of these sites exhibited symptoms of impaired benthic condition. This study showed no major evidence of extensive biological impacts linked to measured stressors. There were only two stations, both in CA, where low numbers of benthic species, diversity, or total faunal abundance co-occurred with high sediment contamination or low DO in bottom water. Such general lack of concordance suggests that these offshore waters are currently in good condition, with the lower-end values of the various biological attributes representing parts of a normal reference range controlled by natural factors. Results of multiple linear regression, performed using full model procedures to test for effects of combined abiotic environmental factors, suggested that latitude and depth had significant influences on benthic variables regionwide. Latitude had a significant inverse influence on all three of the above benthic variables, i.e. with values increasing as latitude decreased (p< 0.01), while depth had a significant direct influence on diversity (p < 0.001) and inverse effect on density (p <0.01). None of these variables varied significantly in relation to sediment % fines (at p< 0.1), although in general there was a tendency for muddier sediments (higher % fines) to have lower species richness and diversity and higher densities than coarser sediments. Alternatively, it is possible that for some of these sites the lower values of benthic variables reflect symptoms of disturbance induced by other unmeasured stressors. The indicators in this study included measures of stressors (e.g., chemical contaminants, eutrophication) that are often associated with adverse biological impacts in shallower estuarine and inland ecosystems. However, there may be other sources of humaninduced stress in these offshore systems (e.g., bottom trawling) that pose greater risks to ambient living resources and which have not been captured. Future monitoring efforts in these offshore areas should include indicators of such alternative sources of disturbance. (137pp.) (PDF contains 167 pages)

Characterization of the benthos, marine debris and bottom fish at Gray’s Reef National Marine Sanctuary

Executive Summary: Baseline characterization of resources is an essential part of marine protected area (MPA) management and is critical to inform adaptive management. Gray’s Reef National Marine Sanctuary (GRNMS) currently lacks adequate characterization of several key resources as identified in the 2006 Final Management Plan. The objectives of this characterization were to fulfill this need by characterizing the bottom fish, benthic features, marine debris, and the relationships among them for the different bottom types within the sanctuary: ledges, sparse live bottom, rippled sand, and flat sand. Particular attention was given to characterizing the different ledge types, their fish communities, and the marine debris associated with them given the importance of this bottom type to the sanctuary. The characterization has been divided into four sections. Section 1 provides a brief overview of the project, its relevance to sanctuary needs, methods of site selection, and general field procedures. Section 2 provides the survey methods, results, discussion, and recommendations for monitoring specific to the benthic characterization. Section 3 describes the characterization of marine debris. Section 4 is specific to the characterization of bottom fish. Field surveys were conducted during August 2004, May 2005, and August 2005. A total of 179 surveys were completed over ledge bottom (n=92), sparse live bottom (n=51), flat sand (n=20), and rippled sand (n=16). There were three components to each field survey: fish counting, benthic assessment, and quantification of marine debris. All components occurred within a 25 x 4 m belt transect. Two divers performed the transect at each survey site. One diver was responsible for identification of fish species, size, and abundance using a visual survey. The second diver was responsible for characterization of benthic features using five randomly placed 1 m2 quadrats, measuring ledge height and other benthic structures, and quantifying marine debris within the entire transect. GRNMS is composed of four main bottom types: flat sand, rippled sand, sparsely colonized live bottom, and densely colonized live bottom (ledges). Independent evaluation of the thematic accuracy of the GRNMS benthic map produced by Kendall et al. (2005) revealed high overall accuracy (93%). Most discrepancies between map and diver classification occurred during August 2004 and likely can be attributed to several factors, including actual map or diver errors, and changes in the bottom type due to physical forces. The four bottom types have distinct physical and biological characteristics. Flat and rippled sand bottom types were composed primarily of sand substrate and secondarily shell rubble. Flat sand and rippled sand bottom types were characterized by low percent cover (0-2%) of benthic organisms at all sites. Although the sand bottom types were largely devoid of epifauna, numerous burrows indicate the presence of infaunal organisms. Sparse live bottom and ledges were colonized by macroalgae and numerous invertebrates, including coral, gorgonians, sponges, and “other” benthic species (such as tunicates, anemones, and bryozoans). Ledges and sparse live bottom were similar in terms of diversity (H’) given the level of classification used here. However, percent cover of benthic species, with the exception of gorgonians, was significantly greater on ledge than on sparse live bottom. Percent biotic cover at sparse live bottom ranged from 0.7-26.3%, but was greater than 10% at only 7 out of 51 sites. Colonization on sparse live bottom is likely inhibited by shifting sands, as most sites were covered in a layer of sediment up to several centimeters thick. On ledge bottom type, percent cover ranged from 0.42-100%, with the highest percent cover at ledges in the central and south-central region of GRNMS. Biotic cover on ledges is influenced by local ledge characteristics. Cluster analysis of ledge dimensions (total height, undercut height, undercut width) resulted in three main categories of ledges, which were classified as short, medium, and tall. Median total percent cover was 97.6%, 75.1%, and 17.7% on tall, medium, and short ledges, respectively. Total percent cover and cover of macroalgae, sponges, and other organisms was significantly lower on short ledges compared to medium and tall ledges, but did not vary significantly between medium and tall ledges. Like sparse live bottom, short ledges may be susceptible to burial by sand, however the results indicate that ledge height may only be important to a certain threshold. There are likely other factors not considered here that also influence spatial distribution and community structure (e.g., small scale complexity, ocean currents, differential settlement patterns, and biological interactions). GRNMS is a popular site for recreational fishing and boating, and there has been increased concern about the accumulation of debris in the sanctuary and potential effects on sanctuary resources. Understanding the types, abundance, and distribution of debris is essential to improving debris removal and education efforts. Approximately two-thirds of all observed debris items found during the field surveys were fishing gear, and about half of the fishing related debris was monofilament fishing line. Other fishing related debris included leaders and spear gun parts, and non-gear debris included cans, bottles, and rope. The spatial distribution of debris was concentrated in the center of the sanctuary and was most frequently associated with ledges rather than at other bottom types. Several factors may contribute to this observation. Ledges are often targeted by fishermen due to the association of recreationally important fish species with this bottom type. In addition, ledges are structurally complex and are often densely colonized by biota, providing numerous places for debris to become stuck or entangled. Analysis of observed boat locations indicated that higher boat activity, which is an indication of fishing, occurs in the center of the sanctuary. On ledges, the presence and abundance of debris was significantly related to observed boat density and physiographic features including ledge height, ledge area, and percent cover. While it is likely that most fishing related debris originates from boats inside the sanctuary, preliminary investigation of ocean current data indicate that currents may influence the distribution and local retention of more mobile items. Fish communities at GRNMS are closely linked to benthic habitats. A list of species encountered, probability of occurrence, abundance, and biomass by habitat is provided. Species richness, diversity, composition, abundance, and biomass of fish all showed striking differences depending on bottom type with ledges showing the highest values of nearly all metrics. Species membership was distinctly separated by bottom type as well, although very short, sparsely colonized ledges often had a similar community composition to that of sparse live bottom. Analysis of fish communities at ledges alone indicated that species richness and total abundance of fish were positively related to total percent cover of sessile invertebrates and ledge height. Either ledge attribute was sufficient to result in high abundance or species richness of fish. Fish diversity (H`) was negatively correlated with undercut height due to schools of fish species that utilize ledge undercuts such as Pareques species. Concurrent analysis of ledge types and fish communities indicated that there are five distinct combinations of ledge type and species assemblage. These include, 1) short ledges with little or no undercut that lacked many of the undercut associated species except Urophycis earlii ; 2) tall, heavily colonized, deeply undercut ledges typically with Archosargus probatocephalus, Mycteroperca sp., and Pareques sp.; 3) tall, heavily colonized but less undercut with high occurrence of Lagodon rhomboides and Balistes capriscus; 4) short, heavily colonized ledges typically with Centropristis ocyurus, Halichoeres caudalis, and Stenotomus sp.; and 5) tall, heavily colonized, less undercut typically with Archosargus probatocephalus, Caranx crysos and Seriola sp.. Higher levels of boating activity and presumably fishing pressure did not appear to influence species composition or abundance at the community level although individual species appeared affected. These results indicate that merely knowing the basic characteristics of a ledge such as total height, undercut width, and percent cover of sessile invertebrates would allow good prediction of not only species richness and abundance of fish but also which particular fish species assemblages are likely to occur there. Comparisons with prior studies indicate some major changes in the fish community at GRNMS over the last two decades although the causes of the changes are unknown. Species of interest to recreational fishermen including Centropristis striata, Mycteroperca microlepis, and Mycteroperca phenax were examined in relation to bottom features, areas of assumed high versus low fishing pressure, and spatial dispersion. Both Mycteroperca species were found more frequently when undercut height of ledges was taller. They often were found together in small mixed species groups at ledges in the north central and southwest central regions of the sanctuary. Both had lower mode size and proportion of fish above the fishery size limit in heavily fished areas of the sanctuary (i.e. high boat density) despite the presence of better habitat in that region. Black sea bass, C. striata, occurred at 98% of the ledges surveyed and appeared to be evenly distributed throughout the sanctuary. Abundance was best explained by a positive relationship with percent cover of sessile biota but was also negatively related to presence of either Mycteroperca species. This may be due to predation by the Mycteroperca species or avoidance of sites where they are present by C. striata. Suggestions for monitoring bottom features, marine debris, and bottom fish at GRNMS are provided at the end of each chapter. The present assessment has established quantitative baseline characteristics of many of the key resources and use issues at GRNMS. The methods can be used as a model for future assessments to track the trajectory of GRNMS resources. Belt transects are ideally suited to providing efficient and quantitative assessment of bottom features, debris, and fish at GRNMS. The limited visibility, sensitivity of sessile biota, and linear nature of ledge habitats greatly diminish the utility of other sampling techniques. Ledges should receive the bulk of future characterization effort due to their importance to the sanctuary and high variability in physical structure, benthic composition, and fish assemblages. (PDF contains 107 pages.)

Age determination methods for northwest Atlantic species

The successful application of techniques to enhance detection of age marks in biological specimens is of vital importance in fisheries research. This manual documents age determination techniques used by staff at the Woods Hole Laboratory, National Marine Fisheries Service. General information on procedures for preparing anatomical structures is described, together with criteria used to interpret growth patterns and assign ages. Annotated photographs of age structures are provided to illustrate criteria. Detailed procedures are given for the following species: Atlantic herring (Clupea harengus), haddock (Melanogrammus aeglefinus), Atlantic cod (Gadus morhua), pollock (Pollachius virens), silver hake (Merluccius bilinearis), red hake (Urophycis chuss), black sea bass (Centropristis striata), weakfish (Cynoscion regalis), Atlantic mackerel (Scomber scombrus), butterfish (Peprilus triacanthus), redfish (Sebastes fasciatus), summer flounder (Paralichthys dentatus), winter flounder (Pseudopleuronectes americanus), witch flounder (Glyptocephalus cynoglossus), American plaice (Hippoglossoides platessoides), yellowtail flounder (Limanda ferruginea), surf clam (Spisula solidissima), and ocean quahog (Arctica islandica). (PDF file contains 142 pages.)

The application of methods of ecosystem analysis to the evaluation of water quality

Guided by experience and the theoretical development of hydrobiology, it can be considered that the main aim of water quality control should be the establishment of the rates of the self-purification process of water bodies which are capable of maintaining communities in a state of dynamic balance without changing the integrity of the ecosystem. Hence, general approaches in the elaboration of methods for hydrobiological control are based on the following principles: a. the balance of matter and energy in water bodies; b. the integrity of the ecosystem structure and of its separate components at all levels. Ecosystem analysis makes possible a revelation of the whole totality of factors which determine the anthropogenic evolution of a water body. This is necessary for the study of long-term changes in water bodies. The principles of ecosystem analysis of water bodies, together with the creation of their mathematical models, are important because, in future, with the transition of water demanding production into closed cycles of water supply, changes in water bodies will arise in the main through the influence of 'diffuse' pollution (from the atmosphere, with utilisation in transport etc.).

Assessing trends in the density of Atlantic croaker (Micropogonias undulatus): a comparison of passive acoustic and trawl methods

Using data collected simultaneously from a trawl and a hydrophone, we found that temporal and spatial trends in densities of juvenile Atlantic croaker (Micropogonias undulatus) in the Neuse River estuary in North Carolina can be identified by monitoring their sound production. Multivariate analysis of covariance (MA NCOVA) revealed that catch per unit of effort (CPUE) of Atlantic croaker had a significant relationship with the dependent variables of sound level and peak frequency of Atlantic croaker calls. Tests of between-subject correspondence failed to detect relationships between CPUE and either of the call parameters, but statistical power was low. Williamson’s index of spatial overlap indicated that call detection rate (expressed by a 0–3 calling index) was correlated in time and space with Atlantic croaker CPUE. The correspondence between acoustic parameters and trawl catch rates varied by month and by habitat. In general, the calling index had a higher degree of overlap with this species’ density than did the received sound level of their calls. Classification and regression tree analysis identified calling index as the strongest correlate of CPUE. Passive acoustics has the potential to be an inexpensive means of identifying spatial and temporal trends in abundance for soniferous fish species.

Restricted Access vs. Open Access Methods of Management: Toward More Effective Regulation of Fishing Effort

This paper gives an overview of the economic rationale for limited entry as a method of fishery management and discusses general advantages and disadvantages of license limitation and catch rights as the two primary methods of restricting access to marine fisheries. Traditional open-access methods of regulation (e.g., gear restrictions, size limits, trip limits, quotas, and closures) can be temporarily effective in protecting fish populations, but they generally fail to provide lasting biological or economic benefits to fishermen because they do not restrict access to the fishery. The general result of regulation with unrestricted access to a fishery is additional and more costly and complex regulations as competition increases for dwindling fishery resources. Regulation that restricts access to a fishery in conjunction with selected traditional methods of regulation would encourage efficient resource usage and minimize the need for future regulatory adjustments, provided that enforcement and monitoring costs are not too great. In theory, catch rights are superior to license limitation as a means of restricting access to a fishery.

A general framework for integrating environmental time series into stock assessment models: model description, simulation testing, and example

We present a method to integrate environmental time series into stock assessment models and to test the significance of correlations between population processes and the environmental time series. Parameters that relate the environmental time series to population processes are included in the stock assessment model, and likelihood ratio tests are used to determine if the parameters improve the fit to the data significantly. Two approaches are considered to integrate the environmental relationship. In the environmental model, the population dynamics process (e.g. recruitment) is proportional to the environmental variable, whereas in the environmental model with process error it is proportional to the environmental variable, but the model allows an additional temporal variation (process error) constrained by a log-normal distribution. The methods are tested by using simulation analysis and compared to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. In the traditional method, the estimates of recruitment were provided by a model that allowed the recruitment only to have a temporal variation constrained by a log-normal distribution. We illustrate the methods by applying them to test the statistical significance of the correlation between sea-surface temperature (SST) and recruitment to the snapper (Pagrus auratus) stock in the Hauraki Gulf–Bay of Plenty, New Zealand. Simulation analyses indicated that the integrated approach with additional process error is superior to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. The results suggest that, for the snapper stock, recruitment is positively correlated with SST at the time of spawning.

Determining the role of linear and nonlinear interactions in records of climate change: possible solar influences in annual to century-scale records

Time series analysis methods have traditionally helped in identifying the role of various forcing mechanisms in influencing climate change. A challenge to understanding decadal and century-scale climate change has been that the linkages between climate changes and potential forcing mechanisms such as solar variability are often uncertain. However, most studies have focused on the role of climate forcing and climate response within a strictly linear framework. Nonlinear time series analysis procedures provide the opportunity to analyze the role of climate forcing and climate responses between different time scales of climate change. An example is provided by the possible nonlinear response of paleo-ENSO-scale climate changes as identified from coral records to forcing by the solar cycle at longer time scales.

A preliminary account and review of the simple methods for determining the operational parameters of fishing gear, underwater, with notes on its application

An account and review of the simple methods for determining the operational parameters of fishing gear, underwater, such a tilt of otter boards (outwards or inwards, forwards or afterwards), vertical height of net, its horizontal spread, angle of divergence at bosom, spread between wing tips, angle of inclination of danlenos, butterfly, slope of legs and sweep-line has been given. The relationship of distance between the otter boards spread and the vertical height of net has been obtained as generally linear. The possibilities of regulating the vertical height of net (dependent variate) and spread of otter boards (independent variate) for increasing the fishing efficiency has been discussed. The angle of attack of oval shaped otter boards used during the operations still remain undetermined, however, it has been explained how the best angle of attack for increasing the efficiency of gear can be obtained by regulating the ratio of depth to warp for a given net. The inadequacy of the mere indices of catch per hour of trawling in comparing the relative efficiency of trawls in gear research studies has been indicated. The importance of estimating the operational parameters and its application to commercial fisheries depending upon the distribution pattern of fish and in gear research has been discussed. The efficiency of the jelly bottle method has been compared statistically with the observations made on the trawl gear underwater with instruments.

On the use of combined compact scheme for numerical solution of the two-layer oceanic shallow water equations

Many types of oceanic physical phenomena have a wide range in both space and time. In general, simplified models, such as shallow water model, are used to describe these oceanic motions. The shallow water equations are widely applied in various oceanic and atmospheric extents. By using the two-layer shallow water equations, the stratification effects can be considered too. In this research, the sixth-order combined compact method is investigated and numerically implemented as a high-order method to solve the two-layer shallow water equations. The second-order centered, fourth-order compact and sixth-order super compact finite difference methods are also used to spatial differencing of the equations. The first part of the present work is devoted to accuracy assessment of the sixth-order super compact finite difference method (SCFDM) and the sixth-order combined compact finite difference method (CCFDM) for spatial differencing of the linearized two-layer shallow water equations on the Arakawa's A-E and Randall's Z numerical grids. Two general discrete dispersion relations on different numerical grids, for inertia-gravity and Rossby waves, are derived. These general relations can be used for evaluation of the performance of any desired numerical scheme. For both inertia-gravity and Rossby waves, minimum error generally occurs on Z grid using either the sixth-order SCFDM or CCFDM methods. For the Randall's Z grid, the sixth-order CCFDM exhibits a substantial improvement , for the frequency of the barotropic and baroclinic modes of the linear inertia-gravity waves of the two layer shallow water model, over the sixth-order SCFDM. For the Rossby waves, the sixth-order SCFDM shows improvement, for the barotropic and baroclinic modes, over the sixth-order CCFDM method except on Arakawa's C grid. In the second part of the present work, the sixth-order CCFDM method is used to solve the one-layer and two-layer shallow water equations in their nonlinear form. In one-layer model with periodic boundaries, the performance of the methods for mass conservation is compared. The results show high accuracy of the sixth-order CCFDM method to simulate a complex flow field. Furthermore, to evaluate the performance of the method in a non-periodic domain the sixth-order CCFDM is applied to spatial differencing of vorticity-divergence-mass representation of one-layer shallow water equations to solve a wind-driven current problem with no-slip boundary conditions. The results show good agreement with published works. Finally, the performance of different schemes for spatial differencing of two-layer shallow water equations on Z grid with periodic boundaries is investigated. Results illustrate the high accuracy of combined compact method.