Geographical distributions of effort and catches of tunas by purse-seine vessels in the Eastern Pacific Ocean during 1965-1998

ENGLISH: This report presents fine-scale spatial summaries of annual catch and effort information compiled by the IATTC staff. These data summaries are presented in a graphical format, and display only information collected from purse-seine vessels fishing in the eastern Pacific Ocean (EPO) during 1965-1998. Data collected from baitboat and longline vessels fishing in the EPO are not considered in this report. Equivalent data from Japanese longline vessels fishing in the EPO are presented by Uosaki and Bayliff(1999) for 1988-1992 and by' references cited in that publication for 1956-1987. SPANISH: Este informe presenta resúmenes espaciales a escala fina de información sobre captura y esfuerzo anuales compilada por el personal de la CIAT. Se presentan los resúmenes en formato gráfico, e incluyen solamente información tomada de barcos cerqueros pescando en el Océano Pacifico oriental (OPO) durante 1965-1998. No se consideran en el informe datos provenientes de barcos de camada y palangreros que pescan en el OPO. Se presentan datos equivalentes de barcos palangreros japoneses pescando en el OPO en Uosaki y Bayliff (1999) para 1988-1992 yen las referencias citadas en dicha publicación para 1956-1987. (PDF contains 104 pages.)

Technical workshop on calculating Nmin for the dolphin stocks of the eastern Pacific Ocean

ENGLISH: This report based on the minutes of a technical workshop carried out under the auspices of the Agreement on the International Dolphin Conservation Program, which took place in La Jolla, California, USA, on August 2-5, 2005. It is reproduced as an IATTC Special Report to make it more widely available to the general public. Some minor changes in formatting have been made, but nothing of scientific importance has been deleted from or added to the report. SPANISH: El presente informe se basa en el acta de una reunión técnica que se celebró en La Jolla, California (EE.UU.) del 2 al 5 de agosto de 2005, bajo los auspicios del Acuerdo sobre el Programa Internacional para la Conservación de los Delfines. Se reproduce como Informe Especial de la CIAT para difundirlo más ampliamente al público general. Se han cambiado unos detalles del formato, pero no se ha añadido ni sustraido nada de importancia científica.

An Evaluation of the area stratification used for sampling tunas in the eastern Pacific Ocean and implications for estimating total annual catches

The Inter-American Tropical Tuna Commission (IATTC) staff has been sampling the size distributions of tunas in the eastern Pacific Ocean (EPO) since 1954, and the species composition of the catches since 2000. The IATTC staff use the data from the species composition samples, in conjunction with observer and/or logbook data, and unloading data from the canneries to estimate the total annual catches of yellowfin (Thunnus albacares), skipjack (Katsuwonus pelamis), and bigeye (Thunnus obesus) tunas. These sample data are collected based on a stratified sampling design. I propose an update of the stratification of the EPO into more homogenous areas in order to reduce the variance in the estimates of the total annual catches and incorporate the geographical shifts resulting from the expansion of the floating-object fishery during the 1990s. The sampling model used by the IATTC is a stratified two-stage (cluster) random sampling design with first stage units varying (unequal) in size. The strata are month, area, and set type. Wells, the first cluster stage, are selected to be sampled only if all of the fish were caught in the same month, same area, and same set type. Fish, the second cluster stage, are sampled for lengths, and independently, for species composition of the catch. The EPO is divided into 13 sampling areas, which were defined in 1968, based on the catch distributions of yellowfin and skipjack tunas. This area stratification does not reflect the multi-species, multi-set-type fishery of today. In order to define more homogenous areas, I used agglomerative cluster analysis to look for groupings of the size data and the catch and effort data for 2000–2006. I plotted the results from both datasets against the IATTC Sampling Areas, and then created new areas. I also used the results of the cluster analysis to update the substitution scheme for strata with catch, but no sample. I then calculated the total annual catch (and variance) by species by stratifying the data into new Proposed Sampling Areas and compared the results to those reported by the IATTC. Results showed that re-stratifying the areas produced smaller variances of the catch estimates for some species in some years, but the results were not significant.

A-SCALA: an age-structured statistical catch-at-length analysis for assessing tuna stocks in the eastern tropical Pacific Ocean

English: We describe an age-structured statistical catch-at-length analysis (A-SCALA) based on the MULTIFAN-CL model of Fournier et al. (1998). The analysis is applied independently to both the yellowfin and the bigeye tuna populations of the eastern Pacific Ocean (EPO). We model the populations from 1975 to 1999, based on quarterly time steps. Only a single stock for each species is assumed for each analysis, but multiple fisheries that are spatially separate are modeled to allow for spatial differences in catchability and selectivity. The analysis allows for error in the effort-fishing mortality relationship, temporal trends in catchability, temporal variation in recruitment, relationships between the environment and recruitment and between the environment and catchability, and differences in selectivity and catchability among fisheries. The model is fit to total catch data and proportional catch-at-length data conditioned on effort. The A-SCALA method is a statistical approach, and therefore recognizes that the data collected from the fishery do not perfectly represent the population. Also, there is uncertainty in our knowledge about the dynamics of the system and uncertainty about how the observed data relate to the real population. The use of likelihood functions allow us to model the uncertainty in the data collected from the population, and the inclusion of estimable process error allows us to model the uncertainties in the dynamics of the system. The statistical approach allows for the calculation of confidence intervals and the testing of hypotheses. We use a Bayesian version of the maximum likelihood framework that includes distributional constraints on temporal variation in recruitment, the effort-fishing mortality relationship, and catchability. Curvature penalties for selectivity parameters and penalties on extreme fishing mortality rates are also included in the objective function. The mode of the joint posterior distribution is used as an estimate of the model parameters. Confidence intervals are calculated using the normal approximation method. It should be noted that the estimation method includes constraints and priors and therefore the confidence intervals are different from traditionally calculated confidence intervals. Management reference points are calculated, and forward projections are carried out to provide advice for making management decisions for the yellowfin and bigeye populations. Spanish: Describimos un análisis estadístico de captura a talla estructurado por edad, A-SCALA (del inglés age-structured statistical catch-at-length analysis), basado en el modelo MULTIFAN- CL de Fournier et al. (1998). Se aplica el análisis independientemente a las poblaciones de atunes aleta amarilla y patudo del Océano Pacífico oriental (OPO). Modelamos las poblaciones de 1975 a 1999, en pasos trimestrales. Se supone solamente una sola población para cada especie para cada análisis, pero se modelan pesquerías múltiples espacialmente separadas para tomar en cuenta diferencias espaciales en la capturabilidad y selectividad. El análisis toma en cuenta error en la relación esfuerzo-mortalidad por pesca, tendencias temporales en la capturabilidad, variación temporal en el reclutamiento, relaciones entre el medio ambiente y el reclutamiento y entre el medio ambiente y la capturabilidad, y diferencias en selectividad y capturabilidad entre pesquerías. Se ajusta el modelo a datos de captura total y a datos de captura a talla proporcional condicionados sobre esfuerzo. El método A-SCALA es un enfoque estadístico, y reconoce por lo tanto que los datos obtenidos de la pesca no representan la población perfectamente. Además, hay incertidumbre en nuestros conocimientos de la dinámica del sistema e incertidumbre sobre la relación entre los datos observados y la población real. El uso de funciones de verosimilitud nos permite modelar la incertidumbre en los datos obtenidos de la población, y la inclusión de un error de proceso estimable nos permite modelar las incertidumbres en la dinámica del sistema. El enfoque estadístico permite calcular intervalos de confianza y comprobar hipótesis. Usamos una versión bayesiana del marco de verosimilitud máxima que incluye constreñimientos distribucionales sobre la variación temporal en el reclutamiento, la relación esfuerzo-mortalidad por pesca, y la capturabilidad. Se incluyen también en la función objetivo penalidades por curvatura para los parámetros de selectividad y penalidades por tasas extremas de mortalidad por pesca. Se usa la moda de la distribución posterior conjunta como estimación de los parámetros del modelo. Se calculan los intervalos de confianza usando el método de aproximación normal. Cabe destacar que el método de estimación incluye constreñimientos y distribuciones previas y por lo tanto los intervalos de confianza son diferentes de los intervalos de confianza calculados de forma tradicional. Se calculan puntos de referencia para el ordenamiento, y se realizan proyecciones a futuro para asesorar la toma de decisiones para el ordenamiento de las poblaciones de aleta amarilla y patudo.

Growth of Skipjack tuna (Katsuwonus pelamis) in the eastern Pacific Ocean as estimated from tagging data

English: Data obtained from tagging experiments initiated during 1953-1958 and 1969-1981 for skipjack tuna from the coastal eastern Pacific Ocean (EPO) are reanalyzed, using the Schnute generalized growth model. The objective is to provide information that can be used to generate a growth transition matrix for use in a length-structured population dynamics model. The analysis includes statistical approaches to include individual variability in growth as a function of length at release and time at liberty, measurement error, and transcription error. The tagging data are divided into northern and southern regions, and the results suggest that growth rates differ between the two regions. The Schnute model provides a significantly better fit to the data than the von Bertalanffy model, a sub-model of the Schnute model, for the northern region, but not for the southern region. Individual variation in growth is best described as a function of time at liberty and as a function of growth increment for the northern and southern regions, respectively. Measurement error is a significant part of the total variation, but the results suggest that there is no bias caused by the measurement error. Additional information, particularly for small and large fish, is needed to produce an adequate growth transition matrix that can be used in a length-structured population dynamics model for skipjack tuna in the EPO. Spanish: Los datos obtenidos de los experimentos de marcado iniciados durante los períodos de 1953- 1958 y de 1969-1981 para el atún barrilete en las costas del Océano Pacífico Oriental (OPO) fueron analizados nuevamente, utilizando el modelo de crecimiento generalizado de Schnute. El objetivo es brindar información que sea útil para producir una matriz sobre la tran-sición de crecimiento que pueda utilizarse en un modelo de dinámica poblacional estructurado por talla. El análisis usa enfoques estadísticos para poder incluir la variabilidad individual del crecimiento como función de la talla de liberación y tiempo en libertad, el error de medición, y el error de transcripción. Los datos de marcado son divididos en regiones norte y sur, y los resultados sugieren que las tasas de crecimiento en las dos regiones son diferentes. En la región norte, pero no en la región sur, el modelo de Schnute se ajusta significativamente mejor a los datos que el modelo von Bertalanffy, un sub-modelo del modelo de Schnute. La mejor descripción de la variación individual en el crecimiento es como una función del tiempo en libertad y como una función del incremento de crecimiento para las regiones norte y sur, respectivamente. El error de medición es una parte significativa de la variación total, pero los resultados sugieren que no existe un sesgo causado por el error de medición. Se necesita información adicional, particularmente para peces pequeños y grandes, para poder producir una matriz de transición de crecimiento adecuada que pueda utilizarse en el modelo de dinámica poblacional estructurado por tallas para el atún barrilete en el OPO.

A model of the pelagic ecossystem in the eastern tropical Pacific Ocean

English: Recent calls for a more holistic approach to fisheries management have motivated development of trophic mass-balance models of ecosystems that underlie fisheries production. We developed a model hypothesis of the pelagic ecosystem in the eastern tropical Pacific Ocean (ETP) to gain insight into the relationships among the various species in the system and to explore the ecological implications of alternative methods of harvesting tunas. We represented the biomasses of and fluxes between the principal elements in the ecosystem with Ecopath, and examined the ecosystem's dynamic, time-series behavior with Ecosim. We parameterized the model for 38 species or groups of species, and described the sources, justifications, assumptions, and revisions of our estimates of the various parameters, diet relations, fisheries landings, and fisheries discards in the model. We conducted sensitivity analyses with an intermediate version of the model, for both the Ecopath mass-balance and the dynamic trajectories predicted by Ecosim. The analysis showed that changes in the basic parameters for two components at middle trophic levels, Cephalopods and Auxis spp., exert the greatest influence on the system. When the Cephalopod Q/B and Auxis spp. P/B were altered from their initial values and the model was rebalanced, the trends of the biomass trajectories predicted by Ecosim were not sensitive, but the scaling was sensitive for several components. We described the review process the model was subjected to, which included reviews by the IATTC Purse-seine Bycatch Working Group and by a working group supported by the National Center for Ecological Analysis and Synthesis. We fitted the model to historical time series of catches per unit of effort and mortality rates for yellowfin and bigeye tunas in simulations that incorporated historical fishing effort and a climate driver to represent the effect of El Niño-Southern Oscillation-scale variation on the system. The model was designed to evaluate the possible ecological implications of fishing for tunas in various ways. We recognize that a model cannot possibly represent all the complexity of a pelagic ocean ecosystem, but we believe that the ETP model provides insight into the structure and function of the pelagic ETP. Spanish: Llamamientos recientes hacia un enfoque más holístico al ordenamiento de la pesca han motivado el desarrollo de modelos tróficos de balance de masas de los ecosistemas que sostienen la producción pesquera. Desarrollamos una hipótesis modelo del ecosistema pelágico en el Océano Pacífico oriental tropical (POT) con miras a mejorar los conocimientos de las relaciones entre las distintas especies en el sistema y explorar las implicaciones ecológicas de métodos alternativos de capturar atunes. Con Ecopath representamos las biomasas de los elementos principales en el ecosistema, y los flujos entre los mismos, y con Ecosim examinamos el comportamiento dinámico del ecosistema con el tiempo. Parametrizamos el modelo para 38 especies o grupos de especies (denominados “componentes” del modelo), y describimos las fuentes, justificaciones, supuestos, y revisiones de nuestras estimaciones de los distintos parámetros, relaciones basadas en dieta, capturas retenidas de las pesquerías, y descartes de las mismas en el modelo. Realizamos análisis de sensibilidad con una versión intermedia del modelo, para el balance de masas de Ecopath y las trayectorias dinámicas predichas por Ecosim también. El análisis demostró que cambios en los parámetros básicos para dos componentes en niveles tróficos medianos, Cefalópodos y Auxis spp., ejercieron la mayor influencia sobre el sistema. Cuando se alteraron el Q/B de los Cefalópodos y el P/B de los Auxis spp. de sus valores iniciales y se balanceó el modelo de nuevo, las tendencias de las trayectorias de la biomasa predichas por Ecosim no fueron sensibles, pero la escala fue sensible para varios componentes. Describimos el proceso de revisión al que fue sujeto el modelo, inclusive revisiones por el Grupo de Trabajo sobre Captura Incidental de la CIAT y un grupo de trabajo apoyado por el Centro Nacional para Síntesis y Análisis Ecológicos. Ajustamos el modelo a series de tiempo históricas de capturas por unidad de esfuerzo y tasas de mortalidad de atunes aleta amarilla y patudo en simulaciones que incorporaron esfuerzo de pesca histórico e impulsos climáticos para representar el efecto de variaciones a escala de El Niño-Oscilación del Sur sobre el sistema. El modelo fue diseñado para evaluar las posibles implicaciones ecológicas de la pesca atunera de varias formas. Reconocemos la imposibilidad de que el modelo represente toda la complejidad de un ecosistema oceánico pelágico, pero creemos que el modelo del POT mejora los conocimientos de la estructura y función del POT pelágico.

A review of the Japanese longline fishery for tunas and billfishes in the eastern Pacific Ocean, 1993-1997

English: This report reviews the Japanese longline fishery in the eastern Pacific Ocean during the 1993-1997 period, extending the studies for the 1956-1992 period made by other investigators. The spatial and temporal distributions of fishing effort, catch, apparent abundance, sexual maturity, and size composition are examined for the principal species of tunas and billfishes taken by that fishery. Some information on the catches of sharks by the Japanese longline fishery is given. The interactions between the surface and longline fisheries are discussed. Spanish: En este informe se presenta un análisis de la actividad pesquera de buques palangreros japoneses en el Océano Pacífico oriental durante el período de 1993-1997, extendiendo los estudios del período de 1956-1992 realizados por otros investigadores. Se examinan las distribuciones espacial y temporal del esfuerzo de pesca, la captura, la abundancia aparente, la madurez sexual, y la composición por talla de las principales especies de atunes y picudos capturadas por dicha pesquería. Se presenta cierta información sobre las capturas de tiburones por la pesquería palangrera japonesa. Se describen las interacciones entre las pesquerías de superficie y palangrera.

Bycatches of sharks in the tuna purse-seine fishery of the eastern Pacific Ocean reported by observers of the Inter-American Tropical Tuna Commission, 1993-2004

Information on bycatches of sharks collected by observers of the Inter-American Tropical Tuna Commission (IATTC) between 1993 and 2004 is presented in this data report. This report contains two sections. The first section summarizes information used by the staff of the IATTC to review and revise IATTC observers’ at-sea species identifications of Carcharhinus falciformis, C. limbatus, and C. longimanus. The revisions were based on 1) data collected on species-specific diagnostic characteristics as part of a special sampling program conducted between March 2000, and March 2001 and 2) a review of observers’ archival field notes for the 1993-2004 period. The second section summarizes the shark bycatches reported by IATTC observers between 1993 and 2004, incorporating the revisions of observers’ at-sea identifications. The IATTC-observed shark bycatch data are summarized as tables with annual tallies of observed bycatches and maps of the spatial distributions of the average bycatches per set and size compositions of the bycatches.

Postrelease survival, vertical and horizontal movements, and thermal habitats of five species of pelagic sharks in the central Pacific Ocean

From 2001 to 2006, 71 pop-up satellite archival tags (PSATs) were deployed on five species of pelagic shark (blue shark [Prionace glauca]; shortfin mako [Isurus oxyrinchus]; silky shark [Carcharhinus falciformis]; oceanic whitetip shark [C. longimanus]; and bigeye thresher [Alopias superciliosus]) in the central Pacific Ocean to determine species-specific movement patterns and survival rates after release from longline fishing gear. Only a single postrelease mortality could be unequivocally documented: a male blue shark which succumbed seven days after release. Meta-analysis of published reports and the current study (n=78 reporting PSATs) indicated that the summary effect of postrelease mortality for blue sharks was 15% (95% CI, 8.5–25.1%) and suggested that catch-and-release in longline fisheries can be a viable management tool to protect parental biomass in shark populations. Pelagic sharks displayed species-specific depth and temperature ranges, although with significant individual temporal and spatial variability in vertical movement patterns, which were also punctuated by stochastic events (e.g., El Niño-Southern Oscillation). Pelagic species can be separated into three broad groups based on daytime temperature preferences by using the unweighted pair-group method with arithmetic averaging clustering on a Kolmogorov-Smirnov Dmax distance matrix: 1) epipelagic species (silky and oceanic whitetip sharks), which spent >95% of their time at temperatures within 2°C of sea surface temperature; 2) mesopelagic-I species (blue sharks and shortfin makos, which spent 95% of their time at temperatures from 9.7° to 26.9°C and from 9.4° to 25.0°C, respectively; and 3) mesopelagic-II species (bigeye threshers), which spent 95% of their time at temperatures from 6.7° to 21.2°C. Distinct thermal niche partitioning based on body size and latitude was also evident within epipelagic species.

Accuracy of sex determination for northeastern Pacific Ocean thornyheads (Sebastolobus altivelis and S. alascanus)

Determining the sex of thornyheads (Sebastolobus alascanus and S. altivelis) can be difficult under field conditions. We assessed our ability to correctly assign sex in the field by comparing results from field observations to results obtained in the laboratory through both macroscopic and microscopic examination of gonads. Sex of longspine thornyheads was more difficult to determine than that of shortspine thornyheads and correct determination of sex was signif icantly related to size. By restricting the minimum size of thornyheads to 18 cm for macroscopic determination of sex we reduced the number of fish with misidentified sex by approximately 65%.

A modeling approach to identify optimal habitat and suitable fishing grounds for neon flying squid (Ommastrephes bartramii) in the Northwest Pacific Ocean

We developed a habitat suitability index (HSI) model to understand and identify the optimal habitat and potential fishing grounds for neon f lying squid (Ommastrephes bartramii) in the Northwest Pacific Ocean. Remote sensing data, including sea surface temperature, sea surface salinity, sea surface height, and chlorophyll-a concentrations, as well as fishery data from Chinese mainland squid f leets in the main fishing ground (150–165°E longitude) from August to October, from 1999 to 2004, were used. The HSI model was validated by using fishery data from 2005. The arithmetic mean modeling with three of the environmental variables—sea surface temperature, sea surface height anomaly, and chlorophyll- a concentrations—was defined as the most parsimonious HSI model. In 2005, monthly HSI values >0.6 coincided with productive fishing grounds and high fishing effort from August to October. This result implies that the model can reliably predict potential f ishing grounds for O. bartramii. Because spatially explicit fisheries and environmental data are becoming readily available, it is feasible to develop a dynamic, near real-time habitat model for improving the process of identifying potential fishing areas for and optimal habitats of neon flying squid.

Reproductive biology of blue marlin (Makaira nigricans) in the western Pacific Ocean

The reproductive biology of blue marlin (Makaira nigricans) was assessed from 1001 fish (ranging from 121 to 275 cm in eye-to-fork length; EFL) caught by Taiwanese offshore longliners in the western Pacific Ocean from September 2000 to December 2001 and from 843 gonad samples from these fish, The overall sex ratio of the catch was approximately 1:1 dur ing the sampling period, but blue marlin are sexually dimorphic; females are larger than males. Reproductive activity (assessed by histology), a gonadosomatic index, and the distribution of oocyte diameters, indicated that spawning occurred predominantly from May to September. The estimated sizes-at-maturity (EFL50) were 179.76 ±1.01 cm (mean ±standard error) for females and 130 ±1 cm EFL for males. Blue marlin are multiple spawners and oocytes develop asynchronously. The proportion of mature females with ovaries containing postovulatory follicles (0.41) and hydrated oocytes (0.34) indicated that the blue marlin spawned once every 2–3 days on average. Batch fecundity (BF) for 26 females with the most advanced oocytes (≥1000 μm), but without postovulatory follicles, ranged from 2.11 to 13.50 million eggs (6.94 ± 0.54 million eggs). The relationships between batch fecundity (BF, in millions of eggs) and EFL and round weight (RW, kg) were BF = 3.29 × 10 –12 EFL5.31 (r2 = 0.70) and BF = 1.59 × 10–3 RW 1.73 (r2= 0.67), respectively. The parameters estimated in this study are key information for stock assessments of blue marlin in the western Pacific Ocean and will contribute to the conservation and sustainable yield of

Abundance, condition, and diet of juvenile Pacific ocean perch (Sebastes alutus) in the Aleutian Islands

T he relative value of pelagic habitat for three size classes of juvenile Pacific ocean perch (Sebastes alutus) was investigated by comparing their abundance and condition in two areas of the Aleutian Islands. Diet, zooplankton biomass, and water column temperatures were examined as potential factors affecting observed differences. Juvenile Pacific ocean perch abundance and condition, and zooplankton biomass varied significantly between areas, whereas juvenile Pacific ocean perch diet varied only by size class. Observed differences in fish condition may have been due to the quantity or quality of pelagic prey items consumed. For the delineation of essential demersal fish habitat, important ecological features of the pelagic habitat must therefore be considered.

Prevalence of the commensal barnacle Xenobalanus globicipitis on cetacean species in the eastern tropical Pacific Ocean, and a review of global occurrence

Distribution and prevalence of the phoretic barnacle Xenobalanus on cetacean species are reported for 22 cetaceans in the eastern tropical Pacific Ocean (21 million km2). Four cetacean species are newly reported hosts for Xenobalanus: Bryde’s whale (Balaenoptera edeni), long-beaked common dolphin (Delphinus capensis), humpback whale (Megaptera novaeangliae), and spinner dolphin (Stenella longirostris). Sightings of Xenobalanus in pelagic waters are reported for the first time, and concentrations were located within three productive zones: near the Baja California peninsula, the Costa Rica Dome and waters extending west along the 10°N Thermocline Ridge, and near Peru and the Galapagos Archipelago. Greatest prevalence was observed on blue whales (Balaenoptera musculus) indicating that slow swim speeds are not necessary for effective barnacle settlement. Overall, prevalence and prevalence per sighting were generally lower than previously reported. The number of barnacles present on an individual whale was greatest for killer whales, indicating that Xenobalanus larvae may be patchily distributed. The broad geographic distribution and large number of cetacean hosts, indicate an extremely cosmopolitan distribution. A better understanding of the biology of Xenobalanus is needed before this species can be used as a biological tag.

An ecological analysis of rockfish (Sebastes spp.) assemblages in the North Pacific Ocean along broad-scale environmental gradients

Environmental variability affects the distributions of most marine fish species. In this analysis, assemblages of rockfish (Sebastes spp.) species were defined on the basis of similarities in their distributions along environmental gradients. Data from 14 bottom trawl surveys of the Gulf of Alaska and Aleutian Islands (n=6767) were used. Five distinct assemblages of rockfish were defined by geographical position, depth, and temperature. The 180-m and 275-m depth contours were major divisions between assemblages inhabiting the shelf, shelf break, and lower continental slope. Another noticeable division was between species centered in southeastern Alaska and those found in the northern Gulf of Alaska and Aleutian Islands. The use of environmental variables to define the species composition of assemblages is different from the use of traditional methods based on clustering and nonparametric statistics and as such, environmentally based analyses should result in predictable assemblages of species that are useful for ecosystem-based management.