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.)

A comparison of seafloor habitats and associated benthic fauna in areas open and closed to bottom trawling along the central California Continental Shelf

Executive Summary: A number of studies have shown that mobile, bottom-contact fishing gear (such as otter trawls) can alter seafloor habitats and associated biota. Considerably less is known about the recovery of these resources following such disturbances, though this information is critical for successful management. In part, this paucity of information can be attributed to the lack of access to adequate control sites – areas of the seafloor that are closed to fishing activity. Recent closures along the coast of central California provide an excellent opportunity to track the recovery of historically trawled areas and to compare recovery rates to adjacent areas that continue to be trawled. In June 2006 we initiated a multi-year study of the recovery of seafloor microhabitats and associated benthic fauna inside and outside two new Essential Fish Habitat (EFH) closures within the Cordell Bank and Gulf of the Farallones National Marine Sanctuaries. Study sites inside the EFH closure at Cordell Bank were located in historically active areas of fishing effort, which had not been trawled since 2003. Sites outside the EFH closure in the Gulf of Farallones were located in an area that continues to be actively trawled. All sites were located in unconsolidated sands at equivalent water depths. Video and still photographic data collected via a remotely operated vehicle (ROV) were used to quantify the abundance, richness, and diversity of microhabitats and epifaunal macro-invertebrates at recovering and actively trawled sites, while bottom grabs and conductivity/temperature/depth (CTD) casts were used to quantify infaunal diversity and to characterize local environmental conditions. Analysis of still photos found differences in common seafloor microhabitats between the recovering and actively trawled areas, while analysis of videographic data indicated that biogenic mound and biogenic depression microhabitats were significantly less abundant at trawled sites. Each of these features provides structure with which demersal fishes, across a wide range of size classes, have been observed to associate. Epifaunal macro-invertebrates were sparsely distributed and occurred in low numbers in both treatments. However, their total abundance was significantly different between treatments, which was attributable to lower densities at trawled sites. In addition, the dominant taxa were different between the two sites. Patchily-distributed buried brittle stars dominated the recovering site, and sea whips (Halipteris cf. willemoesi) were most numerous at the trawled site though they occurred in only five of ten transects. Numerical classification (cluster analysis) of the infaunal samples also revealed a clear difference between benthic assemblages in the recovering vs. trawled areas due to differences in the relative abundances of component species. There were no major differences in infaunal species richness, H′ diversity, or J′ evenness between recovering vs. trawled site groups. However, total infaunal abundance showed a significant difference attributable to much lower densities at trawled sites. This pattern was driven largely by the small oweniid polychaete Myriochele gracilis, which was the most abundant species in the overall study region though significantly less abundant at trawled sites. Other taxa that were significantly less abundant at trawled sites included the polychaete M. olgae and the polychaete family Terebellidae. In contrast, the thyasirid bivalve Axinopsida serricata and the polychaetes Spiophanes spp. (mostly S. duplex), Prionospio spp., and Scoloplos armiger all had significantly to near significantly higher abundances at trawled sites. As a result of such contrasting species patterns, there also was a significant difference in the overall dominance structure of infaunal assemblages between the two treatments. It is suggested that the observed biological patterns were the result of trawling impacts and varying levels of recovery due to the difference in trawling status between the two areas. The EFH closure was established in June 2006, within a month of when sampling was conducted for the present study, however, the stations within this closure area are at sites that actually have experienced little trawling since 2003, based on National Marine Fishery Service trawl records. Thus, the three-year period would be sufficient time for some post-trawling changes to have occurred. Other results from this study (e.g., similarly moderate numbers of infaunal species in both areas that are lower than values recorded elsewhere in comparable habitats along the California continental shelf) also indicate that recovery within the closure area is not yet complete. Additional sampling is needed to evaluate subsequent recovery trends and persistence of effects. Furthermore, to date, the study has been limited to unconsolidated substrates. Ultimately, the goal of this project is to characterize the recovery trajectories of a wide spectrum of seafloor habitats and communities and to link that recovery to the dynamics of exploited marine fishes. (PDF has 48 pages.)

Benthic Macrofauna of the New York Bight, 1979-89

The benthic macrofauna of the New York Bight has been monitored extensively, primarily to determine trends over space and time in biological effects of waste inputs. In the present study, from 44 to 48 stations were sampled each summer from 1980-1985. Data from other Bight benthic studies are included to· extend the temporal coverage from 1979 to 1989. Numbers of species and amphipods per sample, taken as relatively sensitive indicators of environmental stress, showed consistent spatial patterns. Lowest values were found in the Christiaensen Basin and other inshore areas, and numbers increased toward the outermost shelf and Hudson Shelf Valley stations. There were statistically significant decreases in species and amphipods at most stations from 1980 to 1985. (Preliminary data from a more recent study suggest numbers of species increased again between 1986 and 1989.) Cluster analysis of 1980-85 data indicated several distinct assemblages-sewage sludge dumpsite, sludge accumulation area, inner Shelf Valley, outer Shelf Valley, outer shelf-with little change over time. The "enriched" and "highly altered" assemblages in the Basin appear similar to those reported since sampling began there in 1968. No consistently defaunated areas have been found in any sampling programs over the past 20 years. On a gross level, therefore, recent faunal responses to any environmental changes are not evident, but the more sensitive measures used, i.e. numbers of species and amphipods, do indicate widespread recent effects. Causes of the faunal changes are not obvious; some possibilities, including increasing effects of sewage sludge or other waste inputs, natural factors, and sampling artifacts, are discussed. (PDF file contains 54 pages.)

Assessing the sustainability of brackish water aquaculture systems in the Phillipines. Working Paper 1: A typology of brackish-water pond aquaculture systems in the Philippines

Aquaculture in the Philippines is a long-standing activity but has witnessed relatively recent, rapid, technical change with the introduction of hatchery technology and commercial feed-mills changing the production possibilities for a fishpond operator. We are confronted with a diversity of aquaculture practices in the coastal areas of the Philippines, with new technologies being incorporated into more traditional systems. As a first step to understanding the sector, we therefore present a typology of farming systems with the motivation of generating domains (farm “types”) over which we can compare performance on a number of indicators. Our typology, restricted to brackish-water pond systems, is constructed using multivariate methods (principal components analysis, cluster analysis). Eight variables are used relating to the management of the farm across all the major factors of production. A stratified net sample of 136 observations provides the data for the analysis, from a farm-level survey carried out between January and June 2003 in the two main brackish-water production regions in the Philippines. We define five farm types from this analysis. In later work we will show how the use of this typology can be used for comparative study of economic, social and ecological performance at the farm-level. [PDF contains 42 pages]

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.

Mississippi River sustenance of brown shrimp (Farfantepenaeus aztecus) in Louisiana coastal waters

Brown shrimp (Farfantepenaeus aztecus) are abundant along the Louisiana coast, a coastline that is heavily influenced by one of the world’s largest rivers, the Mississippi River. Stable carbon, nitrogen, and sulfur (CNS) isotopes of shrimp and their proventriculus (stomach) contents were assayed to trace riverine support of estuarine-dependent brown shrimp. Extensive inshore and of fshore collections were made in the Louisiana coastal zone during 1999–2006 to document shrimp movement patterns across the bay and shelf region. Results showed an unexpectedly strong role for nursery areas in the river delta in supporting the offshore fishery, with about 46% of immigrants to offshore regions arriving from riverine marshes. Strong river influences also were evident offshore, where cluster analysis of combined CNS isotope data showed three regional station groups related to river inputs. Two nearer-river mid-shelf station groups showed isotope values indicating river fertilization and productivity responses in the benthic shrimp food web, and a deeper offshore station group to the south and west showed much less river inf luence. At several mid-shelf stations where hypoxia is common, shrimp were anomalously 15N depleted versus their diets, and this d15N difference or mismatch may be useful in monitoring shrimp movement responses to hypoxia.

Ichthyoplankton community in the Columbia River plume off Oregon: effects of fluctuating oceanographic conditions

Ichthyoplankton samples were collected at approximately 2-week intervals, primarily during spring and summer 1999−2004, from two stations located 20 and 30 km from shore near the Columbia River, Oregon. Northern anchovy (Engraulis mordax) was the most abundant species collected, and was the primary species associated with summer upwelling conditions, but it showed significant interannual and seasonal fluctuations in abundance and occurrence. Other abundant taxa included sanddabs (Citharichthys spp.), English sole (Parophrys vetulus), and blacksmelts (Bathylagidae). Two-way cluster analysis revealed strong species associations based primarily on season (before or after the spring transition date). Ichthyoplankton abundances were compared to biological and environmental data, and egg and larvae abundances were found to be most correlated with sea surface temperature. The Pacific Decadal Oscillation changed sign (from negative to positive) in late 2002 and indicated overall warmer conditions in the North Pacific Ocean. Climate change is expected to alter ocean upwelling, temperatures, and Columbia River flows, and consequently fish eggs and larvae distributions and survival. Long-term research is needed to identify how ichthyoplankton and fish recruitment are affected by regional and largescale oceanographic proces

Diets of and trophic relationships among dominant marine nekton within the northern California Current ecosystem

In this study we analyzed the diets of 26 nekton species collected from two years (2000 and 2002) off Oregon and northern California to describe dominant nekton trophic groups of the northern California Current (NCC) pelagic ecosystem. We also examined interannual variation in the diets of three nekton species. Cluster analysis of predator diets resulted in nekton trophic groups based on the consumption of copepods, euphausiids, brachyuran larvae, larval juvenile fishes, and adult nekton. However, many fish within trophic groups consumed prey from multiple trophic levels—euphausiids being the most widely consumed. Comparison of diets between years showed that most variation occurred with changes in the contribution of euphausiids and brachyuran larvae to nekton diets. The importance of euphausiids and other crustacean prey to nekton indicates that omnivory is an important characteristic of the NCC food web; however it may change during periods of lower or higher upwelling and ecosystem production.

An objective classification of climatic regions in the Pacific and Indian oceans

We have applied a number of objective statistical techniques to define homogeneous climatic regions for the Pacific Ocean, using COADS (Woodruff et al 1987) monthly sea surface temperature (SST) for 1950-1989 as the key variable. The basic data comprised all global 4°x4° latitude/longitude boxes with enough data available to yield reliable long-term means of monthly mean SST. An R-mode principal components analysis of these data, following a technique first used by Stidd (1967), yields information about harmonics of the annual cycles of SST. We used the spatial coefficients (one for each 4-degree box and eigenvector) as input to a K-means cluster analysis to classify the gridbox SST data into 34 global regions, in which 20 comprise the Pacific and Indian oceans. Seasonal time series were then produced for each of these regions. For comparison purposes, the variance spectrum of each regional anomaly time series was calculated. Most of the significant spectral peaks occur near the biennial (2.1-2.2 years) and ENSO (~3-6 years) time scales in the tropical regions. Decadal scale fluctuations are important in the mid-latitude ocean regions.

The Hudson-Raritan Estuary as a crossroads for distribution of blue (Callinectes sapidus), lady (Ovalipes ocellatus), and Atlantic rock (Cancer irroratus) crabs

Blue (Callinectes sapidus)(Portunidae),lady (Ovalipes ocellatus)(Portunidae), and Atlantic rock (Cancer irroratus) (Cancridae) crabs inhabit estuaries on the northeast United States coast for parts or all of their life cycles. Their distributions overlap or cross during certain seasons. During a 1991–1994 monthly otter trawl survey in the Hudson-Raritan Estuary between New York and New Jersey, blue and lady crabs were collected in warmer months and Atlantic rock crabs in colder months. Sex ratios, male:female, of mature crabs were 1:2.0 for blue crabs, 1:3.1 for lady crabs, and 21.4:1 for Atlantic rock crabs. Crabs, 1286 in total, were subsampled for dietary analysis, and the dominant prey taxa for all crabs, by volume of foregut contents, were mollusks and crustaceans. The proportion of amphipods and shrimp in diets decreased as crab size increased. Trophic niche breadth was widest for blue crabs, narrower for lady crabs, and narrowest for Atlantic rock crabs. Trophic overlap was lowest between lady crabs and Atlantic rock crabs, mainly because of frequent consumption of the dwarf surfclam (Mulinia lateralis) by the former and the blue mussel (Mytilus edulis) by the latter. The result of cluster analysis showed that size class and location of capture of predators in the estuary were more influential on diet than the species or sex of the predators.

Population structure of king mackerel (Scomberomorus cavalla) around peninsular Florida, as revealed by microsatellite DNA

A total of 1006 king mackerel (Scomberomorus cavalla) representing 20 discrete samples collected between 1996 and 1998 along the east (Atlantic) and west (Gulf) coasts of Florida and the Florida Keys were assayed for allelic variation at seven nuclear-encoded microsatellites. No significant deviations from Hardy-Weinberg equilibrium expectations were found for six of the microsatellites, and genotypes at all microsatellites were independent. Allele distributions at each microsatellite were independent of sex and age of individuals. Homogeneity tests of spatial distributions of alleles at the microsatellites revealed two weakly divergent “genetic” subpopulations or stocks of king mackerel in Florida waters—one along the Atlantic coast and one along the Gulf coast. Homogeneity tests of allele distributions when samples were pooled along seasonal (temporal) boundaries, consistent with the temporal boundaries used currently for stock assessment and allocation of the king mackerel resource, were nonsignificant. The degree of genetic divergence between the two “genetic” stocks was small: on average, only 0.19% of the total genetic variance across all samples assayed occurred between the two regions. Cluster analysis, assignment tests, and spatial autocorrelation analysis did not generate patterns that were consistent with either geographic or spatial-temporal boundaries. King mackerel sampled from the Florida Keys could not be assigned unequivocally to either “genetic” stock. The genetic data were not consistent with current spatial-temporal boundaries employed in stock assessment and allocation of the king mackerel resource. The genetic differences between king mackerel in the Atlantic versus those in the Gulf most likely stem from reduced gene flow (migration) between the Atlantic and Gulf in relation to gene flow (migration) along the Atlantic and Gulf coasts of peninsular Florida. This difference is consistent with findings for other marine fishes where data indicate that the southern Florida peninsula serves (or has served) as a biogeographic boundary.

Population diversity studies of pharaoh cuttlefish (Sepia pharaonis) in the Persian Gulf and the Gulf of Oman

The pharaoh Cuttlefish (Sepia Pharaonis) is the most abundant species in the persian gulf and oman sea. The stock patterns of this species was studied conserning biological, morphological and electrophoretical aspects. In addition to measuring the biological patterns; 21 quantitative and qualitative factors were measured or counted. There tissues namely muscle, eye and liver were also used for further polyacrylamide electrophoresis analysis (SDS-Page method). The densitograms of protein bands of each tissue were prepared by Gel-Scanner and also the amount of bands' areas were estimated. The results of LSD test showed that tentacle length (TL), Tentacle club length (TCL) and TCL were indicative factors which they showed significant difference between male TL and female specimen of Bushehr and Balouchestan regions. Regarding to length and weight frequencies data the results indicated that males are always bigger than females and also, the cuttlefishes of the gulf of oman are bigger than persian gulf's samples. There were found a significant sample correlation (95%) between different quantitative parameters and the most correlation (0.963) was found to be between TL and TCL; whereas the less correlation (0.384) was observed between GL and LA3. However the gill length factor illustrated the less correlation with the other factors. The results of Cluster analysis for both sexes showed that the cuttlefishes of both studied regions belongs to seperate stocks. The electrophoretic experiments of proteins showed that the protein bands of the sample tissue of muscle and eye from both sampling areas revealed significant differences; whereas the tissue of liver wasn't recognized as an indicative tissue for population studies. Taking into consideration the findings of the present study including: (1) difference in spawning season, (2) results of dendrograms, (3) observed Significant differences in one-way analysis of variance (ANOVA) for morphometric measurements, (4) differences in body length and weight, (5) ecological variations of the persian gulf and oman sea, (6) as well as the results of electrophoresis of proteins have indicatend that: "The Pharaoh cuttlefishes of Bushehr and Balouchestan waters belong to two seperate stocks. Also, it is believed that each region has propably two different populations and more studies are needed to approve this result.