Maximum likelihood estimation of mortality and growth with individual variability from multiple length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock and derive the underlying length distribution of the population and the catch when there is individual variability in the von Bertalanffy growth parameter L∞. The model is flexible enough to accommodate 1) any recruitment pattern as a function of both time and length, 2) length-specific selectivity, and 3) varying fishing effort over time. The maximum likelihood method gives consistent estimates, provided the underlying distribution for individual variation in growth is correctly specified. Simulation results indicate that our method is reasonably robust to violations in the assumptions. The method is applied to tiger prawn data (Penaeus semisulcatus) to obtain estimates of natural and fishing mortality.

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.

Year-class formation in Pacific herring (Clupea pallasi) estimated from spawning-date distributions of juveniles in San Francisco Bay, California

Inter and intra-annual variation in year-class strength was analyzed for San Francisco Bay Pacific herring (Clupea pallasi) by using otoliths of juveniles. Juvenile herring were collected from March through June in 1999 and 2000 and otoliths from subsamples of these collections were aged by daily otolith increment analysis. The composition of the year classes in 1999 and 2000 were determined by back-calculating the birth date distribution for surviving juvenile herring. In 2000, 729% more juveniles were captured than in 1999, even though an estimated 12% fewer eggs were spawned in 2000. Spawning-date distributions show that survival for the 2000 year class was exceptionally good for a short (approximately 1 month) period of spawning, resulting in a large abundance of juvenile recruits. Analysis of age at size shows that growth rate increased significantly as the spawning season progressed both in 1999 and 2000. However, only in 2000 were the bulk of surviving juveniles a product of the fast growth period. In the two years examined, year-class strength was not predicted by the estimated number of eggs spawned, but rather appeared to depend on survival of eggs or larvae (or both) through the juvenile stage. Fast growth through the larval stage may have little effect on year-class strength if mortality during the egg stage is high and few larvae are available.

Modeling red sea urchin (Strongylocentrotus franciscanus) growth using six growth functions

The growth of red sea urchins (Strongylocentrotus franciscanus) was modeled by using tag-recapture data from northern California. Red sea urchins (n=211) ranging in test diameter from 7 to 131 mm were examined for changes in size over one year. We used the function Jt+1 = Jt + f(Jt) to model growth, in which Jt is the jaw size (mm) at tagging, and Jt+1 is the jaw size one year later. The function f(Jt), represents one of six deterministic models: logistic dose response, Gaussian, Tanaka, Ricker, Richards, and von Bertalanffy with 3, 3, 3, 2, 3, and 2 minimization parameters, respectively. We found that three measures of goodness of fi t ranked the models similarly, in the order given. The results from these six models indicate that red sea urchins are slow growing animals (mean of 7.2 ±1.3 years to enter the fishery). We show that poor model selection or data from a limited range of urchin sizes (or both) produces erroneous growth parameter estimates and years-to-fishery estimates. Individual variation in growth dominated spatial variation at shallow and deep sites (F=0.246, n=199, P=0.62). We summarize the six models using a composite growth curve of jaw size, J, as a function of time, t: J = A(B – e–Ct) + Dt, in which each model is distinguished by the constants A, B, C, and D. We suggest that this composite model has the flexibility of the other six models and could be broadly applied. Given the robustness of our results regarding the number of years to enter the fishery, this information could be incorporated into future fishery management plans for red sea urchins in northern California.

Comparative study of seasonal changes in some blood parameters and genetic polymorphisms of serum antibodies and antigens on red blood cells surface in immature Caspian salmon

Seasonal sampling from 40 immature Caspian salmon were performed in summer, autumn, winter and spring. The maximum ranges of RBC counts, Hct, Hb, WBC count and clotting times were observed in spring, summer, spring, spring and winter, respectively. The minimum amounts of these factors were counted in summer, winter, winter, winter and winter, respectively. Blood Samples were taken from healthy smolt, immature and adult Caspian salmon in spawning time. Hematological determinations and biochemical serum analysis were performed in 101 fish in the three samples. The ranges of hematological values for sample mean were counted. Red blood cell counts were 866600 mm3 and 1259400 mm3 in smolt and adult respectively. Hematocrit was 48.39% in smolt and 44.29% in adult. Hemoglobin was 8.85 gr/dl in smolt and 10.91 gr/dl in adult. White blood cell count was 8781.58 mm3 in smolt and 5217.55 mm3 in adult and mean were differential of WBC, Lymphocyte 90.57%in smolt and73.22% in adult. Neutrophil was 5.12% in smolt and 16.92% in adult, Monocyte were 1.27% in smolt and 4.24% in adult, Clotting time was 282.34 Seconds in smolt and 291.47 seconds in adult MCV, MCH and MCHC also meagered in smolt and adult. Biochemical parameter in immature and mature Caspian salmon meagered .Glucose concentration was 2.97 mmol.l- in immature and 1.99 mmol.l- in mature .Cholesterol concentration was 4.26 mmol.l- in immature and 7.06 mmol.l- in mature. Triglyceride amount was 2.35 mmol.l- in immature and 2.47 mmol.l- in mature and Calcium was 2.47 in immature and 2.61 mmol.l- in mature. An in situ study was made on erythrocytic isoantigens and hetero-antigen and their corresponding iso-and hetero-antibodies of sera by means of hemoagglutination tests on the blood sample, of 450 immature and 50 mature Caspian salmon. The absence of erythrocyte iso-antigens and hetero-antigen and their corresponding iso-and hetero-antibodies were shown by the experimental. It could be indicated an intra-specific variation and differences in species for kelardasht hatchery.

Variation in size of yellowfin tuna (Thunnus albacares) within individual purse-seine sets

ENGLISH:Length-frequency samples of yellowfin tuna from 276 individual purse-seine sets were examined. Evidence of schooling by size is presented. Yellowfin schooled with skipjack are smaller and more homogeneous in length than are yellowfin from pure schools. Yellowfin in schools associated with porpoise appear to be more variable in size than yellowfin from other types of schools. No relationship was found between the tonnage of yellowfin in a school and the mean length of the yellowfin. Despite the tendency to school by size, the size variation within individual schools was judged to be enough to complicate greatly any program of regulation aimed at maximizing the yield-per-recruit through increasing the minimum size of yellowfin at first capture. SPANISH: Fueron examinadas las muestras frecuencia-longitud de atún aleta amarilla, de 276 lances individuales de redes de cerco. Se presenta la evidencia de agrupación por tamaños. Los atunes aleta amarilla agrupados con barrilete, son más pequeños y más homogéneos en longitud, que los atunes aleta amarilla de cardúmenes puros. El atún aleta amarilla en cardúmenes asociados con delfines parece ser más variable en tamaño, que el atún aleta amarilla proveniente de otros tipos de cardúmenes. No se encontró relac¡'ón entre el tonelaje del atún aleta amarilla en un cardumen y la longitud media de esta especie. A pesar de la tendencia a agruparse por tamaño, se juzgó, que la variación de tamaño en cardúmenes individuales, sería suficiente para complicar grandemente cualquier programa de reglamentación, dirigido a obtener el máximo del rendimiento por recluta a través del incremento del tamaño mínimo del atún aleta amarilla en la primera captura.

Genetic variation survey of intra-specific (Sepia pharaonis) in the Persian Gulf and the Oman Sea using PCR-RFLP

This research was conducted to identify Cuttlefish population (Sepia pharaonis) in The Persian Gulf and the Oman Sea using PCR-RFLP. Specimens were collected from )0 different stations. Bottom trawling method was used for sampling from different zones of the Persian Gulf and the Oman Sea, and finally specimens from S. Pharaonis were collected at each station . DNA was extracted by phenol—Coloroform method. One pair primer was designed based on 1As rRNA gene nucleotide sequences. The results obtained from 1 As rRNA gene RFLP, which was reproduced by PCR technique, were analyzed and utilized for study of diversity of the Cuttlefish population. PCR product with o pair base in length achieved for all specimens, which was subjected to enzymatic digestion by A restriction action enzymes: Alu I-Taq I-Mnl I-Rsa I-Hind III-Dra I-vu II and Hae II DNA bands patterns in all specimens digested by those enzymen showed similarity with no any polymorphism. From this result, it can be concluded that there is not any possibility to isolate different populations in the studied Cuttlefish species under exploitation of rRNA gene.

Report on Marine Mammal Unusual Mortality Event UMESE0501Sp: Multispecies Mass Stranding of Pilot Whales (Globicephala macrorhynchus), Minke Whale (Balaenoptera acutorostrata), and Dwarf Sperm Whales (Kogia sima) in North Carolina on 15-16 January 2005

On 15-16 January 2005, three offshore species of cetaceans (33 short-finned pilot whales, Globicephala macrorhynchus, one minke whale, Balaenoptera acutorostrata, and two dwarf sperm whales, Kogia sima) stranded alive on the beaches of North Carolina. The pilot whales stranded near Oregon Inlet, the minke whale in northern North Carolina, and the dwarf sperm whales near Cape Hatteras. Live strandings of three species in one weekend was unique in North Carolina and qualified as an Unusual Mortality Event. Gross necropsies were conducted on 16-17 January 2005 on 27 pilot whales, two dwarf sperm whales, and the minke whale. Samples were collected for clinical pathology, parasitology, gross pathology, histopathology, microbiology and serology. There was variation in the number of animals sampled for each collection type, however, due to carcasses washing off the beach or degradation in carcass condition during the course of the response. Comprehensive histologic examination was conducted on 16 pilot whales, both dwarf sperm whales, and the minke whale. Limited organ or only head tissue suites were obtained from nine pilot whales. Histologic examination of tissues began in February 2005 and concluded in December 2005 when final sampling was concluded. Neither the pilot whales nor dwarf sperm whales were emaciated although none had recently ingested prey in their stomachs. The minke whale was emaciated; it was likely a dependent calf that became separated from the female. Most serum biochemistry abnormalities appear to have resulted from the stranding and indicated deteriorating condition from being on land for an extended period. Three pilot whales had clinical evidence of pre-existing systemic inflammation, which was supported by histopathologic findings. Although gross and histologic lesions involving all organ systems were noted, consistent lesions were not observed across species. Verminous pterygoid sinusitis and healed fishery interactions were seen in pilot whales but neither of these changes were causes of debilitation or death. In three pilot whales and one dwarf sperm whale there was evidence of clinically significant disease in postcranial tissues which led to chronic debilitation. Cardiovascular disease was present in one pilot whale and one dwarf sperm whale; musculoskeletal disease and intra-abdominal granulomas were present in two pilot whales. These lesions were possible, but not definitive, causal factors in the stranding. Remaining lesions were incidental or post-stranding. The minke whale and three of five tested pilot whales had positive morbillivirus titers (≥1:8 with one at >1:256), but there was no histologic evidence of active viral infection. Parasites (nematodes, cestodes, and trematodes) were collected from 26 pilot whales and two dwarf sperm whales. Sites of collection included stomach, nasal/pterygoid, peribullar sinuses, blubber, and abdominal cavity. Parasite species, locations and loads were within normal limits for free-ranging cetaceans and were not considered causative for the stranding event. Gas emboli lesions which were considered consistent with or diagnostic of sonarassociated strandings of beaked whales or small cetaceans were not found in the whales stranded as part of UMESE0501Sp. Twenty-five heads were examined with nine specific anatomic locations of interest: extramandibular fat, intramandibular fat, auditory meatus, peribullar acoustic fat, peribullar soft tissue, peribullar sinus, pterygoid sinus, melon, and brain. The common finding in all examined heads was verminous pterygoid sinusitis. Intramandibular adipose tissue reddening, typically adjacent to the vascular plexus, was observed in some individuals and could represent localized hemorrhage resulting from vascular rete rupture, hypostatic congestion, or erythrocyte rupture during the freeze/thaw cycle. One cetacean had peracute to acute subdural hemorrhage that likely occurred from thrashing on the beach post-stranding, although its occurrence prior to stranding cannot be excluded. Information provided to NMFS by the U.S. Navy indicated routine tactical mid-frequency sonar operations from individual surface vessels over relatively short durations and small spatial scales within the area and time period investigated. No marine mammals were detected by marine mammal observers on operational vessels; standard operating procedure for surface naval vessels operating mid-frequency sonar is the use of trained visual lookouts using high-powered binoculars. Sound propagation modeling using information provided to NMFS indicated that acoustic conditions in the vicinity likely depended heavily on position of the receivers (e.g., range, bearing, depth) relative to that of the sources. Absent explicit information on the location of animals meant that it was not possible to estimate received acoustic exposures from active sonar transmissions. Nonetheless, the event was associated in time and space with naval activity using mid-frequency active sonar. It also had a number of features in common (e.g., the “atypical” distribution of strandings involving multiple offshore species, all stranding alive, and without evidence of common infectious or other disease process) with other sonar-related cetacean mass stranding events. Given that this event was the only stranding of offshore species to occur within a 2-3 day period in the region on record (i.e., a very rare event), and given the occurrence of the event simultaneously in time and space with a naval exercise using active sonar, the association between the naval sonar activity and the location and timing of the event could be a causal rather than a coincidental relationship. However, evidence supporting a definitive association is lacking, and, in particular, there are differences in operational/environmental characteristics between this event and previous events where sonar has apparently played a role in marine mammal strandings. This does not preclude behavorial avoidance of noise exposure. No harmful algal blooms were present along the Atlantic coast south of the Chesapeake Bay during the months prior to the event. Environmental conditions, including strong winds, changes in upwelling- to downwelling-favorable conditions, and gently sloping bathymetry, were consistent with conditions which have been correlated with other mass strandings. In summary, we did not find commonality in gross and histologic lesions that would indicate a single cause for this stranding event. Three pilot whales and one dwarf sperm whale had debilitating conditions identified that could have contributed to stranding, one pilot whale had a debilitating condition (subdural hemorrhage) that could have been present prior to or resulting from stranding. While the pilot and dwarf sperm whale strandings may have had a common cause, the minke whale stranding was probably just coincidental. On the basis of examination of physical evidence in the affected whales, however, we cannot definitively conclude that there was or was not a causal link between anthropogenic sonar activity or environmental conditions (or a combination of these factors) and the strandings. Overall, the cause of UMESE0501Sp in North Carolina is not and likely will not be definitively known. (PDF contains 240 pages)

Variation in the isotopic signatures of juvenile gray snapper (Lutjanus griseus) from five southern Florida region

Measurements of 18O/16O and 13C/12C ratios in the carbonate of juvenile gray snapper (Lutjanus griseus) sagittal otoliths collected during 2001–2005 from different southern Florida regions indicated significant variations in the ratios between Florida Bay and surrounding areas. Annual differences in isotopic composition were also observed. Classification accuracy of individual otoliths to a region averaged 80% (63% to 96%), thereby enabling the probability of assigning an unknown individual to the appropriate juvenile nursery habitat. Identification of isotopic signatures in the otoliths of gray snapper from Florida Bay and adjacent ecosystems may be important for distinguishing specific portions of the bay that are crucial nursery grounds for juveniles. Separation of gray snapper between geographic regions and nursery sites is possible and has the potential to establish a link between adult gray snapper present on offshore reefs and larvae and juveniles at nursery habitats in Florida Bay or adjacent areas.

Temporal variation in growth of yellowfin tuna (Thunnus albacares) larvae in the Panama Bight, 1990

Tuna larvae (at flexion, postflexion, and transformation stages) were collected by dip net and light traps at night in the northwestern Panama Bight during the season of reduced upwelling (June−September) of 1990, 1991, 1992, and 1997. The larvae were identified as yellowfin tuna (Thunnus albacares) by mtDNA analysis. Ichthyoplankton data from bongo and Tucker trawl tows were used to examine the potential prey abundance in relation to the mean size-at-age and growth rates of the yellowfin tuna larvae and their otoliths. The most rapid growth rates occurred during June 1990 when plankton volumes were at their highest levels. The lowest plankton volumes coincided with the lowest growth rates and mean sizes-at-age during the August−September 1991 period. High densities of larval fish were prevalent in the ichthyoplankton tows during the 1991 period; therefore intra- and interspecific competition for limited food resources may have been the cause of slower growth (density-dependent growth) in yellowfin tuna larvae The highest mean seasurface temperature and the lowest mean wind stress occurred during an El Niño-Southern Oscillation (ENSO) event during the 1997 period. There appeared to be no clear association between these environmental factors and larval growth rates, but the higher temperatures may have caused an increase in the short-term growth of otoliths in relat

Individual growth rates and movement of juvenile white shrimp (Litopenaeus setiferus) in a tidal marsh nursery

We measured growth and movements of individually marked free-ranging juvenile white shrimp (Litopenaeus setiferus) in tidal creek subsystems of the Duplin River, Sapelo Island, Georgia. Over a period of two years, 15,974 juvenile shrimp (40−80 mm TL) were marked internally with uniquely coded microwire tags and released in the shallow upper reaches of four salt marsh tidal creeks. Subsequent samples were taken every 3−6 days from channel segments arranged at 200-m intervals along transects extending from the upper to lower reach of each tidal creek. These collections included 201,384 juvenile shrimp, of which 184 were marked recaptures. Recaptured shrimp were at large an average of 3−4 weeks (range: 2−99 days) and were recovered a mean distance of <0.4 km from where they were initially marked. Mean residence times in the creek subsystems ranged from 15.2 to 25.5 days and were estimated from exponential decay functions describing the proportions of marked individuals recaptured with increasing days at large. Residence time was not significantly correlated with creek length (Pearson=−0.316, P=0.684 ), but there was suggestive evidence of positive associations with either intertidal (Pearson r=0.867, P=0.133) or subtidal (Pearson r=0.946, P=0.054) drainage area. Daily mean specific growth rates averaged 0.009 to 0.013 among creeks; mean absolute growth rates ranged from 0.56−0.84 mm/d, and were lower than those previously reported for juvenile penaeids in estuaries of the southeastern United States. Mean individual growth rates were not significantly different between years (t-test, P>0.30) but varied significantly during the season, tending to be greater in July than November. Growth rates were size-dependent, and temporal changes in size distributions rather than temporal variation in physical environmental factors may have accounted for seasonal differences in growth. Growth rates differed between creeks in 1999 (t-test, P<0.015), but not in 1998 (t-test, P>0.5). We suggest that spatial variation in landscape structure associated with access to intertidal resources may have accounted for this apparent interannual difference in growth response.

Allozyme variation in four hatchery populations of Thai pangas, Pangasius hypophthalmus in Bogra, Bangladesh

Genetic variation of four hatchery stocks of Thai pangas, Pangasius hypophthalmu [sic] of Bogra region, Bangladesh was studied from 1 January 2002 to 31 December 2003. Muscle samples were collected for allozyme analysis from four (Bhai-Bhai, Jahangir, Belal and Bhai-Bon) different hatchery populations. For allozyme electrophoresis, eight enzymes were used and 11 loci viz. Adh-1*, Est-1*, GJpdh-1*, Gpi-1*, Gpi-2*, Jdhp-1*, Ldh-1*, Ldh-2*, Mdh-1*, Mdh-2*and Pgm* were identified, of which three loci (Est-1*, Gpi-2*, G3pdh-1 and Pgm*) were polymorphic in all the four populations. The mean proportion of polymorphic loci per population and the mean proportion of heterozygous loci per individual was 36.36% and 13.33, respectively for all the population studied. The highest variability measured by the mean number of alleles per locus was 1.545 in Bhai-Bon hatchery population. Based on Nei's (1972) genetic distance, the dendrogram (UPGMA) shows that four populations have made two clusters by D-value (D=0.043). Bhai-Bhai and Jahangir hatchery populations have made cluster-I, and Belal Uddin and Bhai-Bon hatchery populations formed cluster-II. Among the four populations, BhaiBhai and Jahangir hatchery populations were differentiated from each other by the D-value of 0.013, and Belal Uddin and Bhai-Bon populations were differentiated from each other by the D-value of 0.002, which suggests that the four populations may be fallen into the local population or race.