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.

An environmentally based growth model that uses finite difference calculus with maximum likelihood method: its application to the brackish water bivalve Corbicula japonica in Lake Abashiri, Japan

We present a growth analysis model that combines large amounts of environmental data with limited amounts of biological data and apply it to Corbicula japonica. The model uses the maximum-likelihood method with the Akaike information criterion, which provides an objective criterion for model selection. An adequate distribution for describing a single cohort is selected from available probability density functions, which are expressed by location and scale parameters. Daily relative increase rates of the location parameter are expressed by a multivariate logistic function with environmental factors for each day and categorical variables indicating animal ages as independent variables. Daily relative increase rates of the scale parameter are expressed by an equation describing the relationship with the daily relative increase rate of the location parameter. Corbicula japonica grows to a modal shell length of 0.7 mm during the first year in Lake Abashiri. Compared with the attain-able maximum size of about 30 mm, the growth of juveniles is extremely slow because their growth is less susceptible to environmental factors until the second winter. The extremely slow growth in Lake Abashiri could be a geographical genetic variation within C. japonica.

A model for assessing the likelihood of self-sustaining populations resulting from commercial production of triploid Suminoe oysters (Crassostrea ariakensis) in Chesapeake Bay

Culture of a non-native species, such as the Suminoe oyster (Crassostrea ariakensis), could offset the harvest of the declining native eastern oyster (Crassostrea virginica) fishery in Chesapeake Bay. Because of possible ecological impacts from introducing a fertile non-native species, introduction of sterile triploid oysters has been proposed. However, recent data show that a small percentage of triploid individuals progressively revert toward diploidy, introducing the possibility that Suminoe oysters might establish self-sustaining populations. To assess the risk of Suminoe oyster populations becoming established in Chesapeake Bay, a demographic population model was developed. Parameters modeled were salinity, stocking density, reversion rate, reproductive potential, natural and harvest-induced mortality, growth rates, and effects of various management strategies, including harvest strategies. The probability of a Suminoe oyster population becoming self-sustaining decreased in the model when oysters are grown at low salinity sites, certainty of harvest is high, mini-mum shell length-at-harvest is small, and stocking density is low. From the results of the model, we suggest adopting the proposed management strategies shown by the model to decrease the probability of a Suminoe oyster population becoming self-sustaining. Policy makers and fishery managers can use the model to predict potential outcomes of policy decisions, supporting the ability to make science-based policy decisions about the proposed introduction of triploid Suminoe oysters into the Chesapeake Bay.

Use of aerial survey and aerophotogrammetry methods in monitoring manatee populations

We evaluated the use of strip-transect survey methods for manatees through a series of replicate aerial surveys in the Banana River, Brevard County, Florida, during summer 1993 and summer 1994. Transect methods sample a representative portion of the total study area, thus allowing for statistical extrapolation to the total area. Other advantages of transect methods are less flight time and less cost than total coverage, ease of navigation, and reduced likelihood of double-counting. Our objectives were: (1) to identify visibility biases associated with the transect survey method and to adjust the counts accordingly; (2) to derive a population estimate with known variance for the Banana River during summer; and (3) to evaluate the potential value of this survey method for monitoring trends in manatee population size over time. (51 page document)

Automated, objective texture segmentation of multibeam echosounder data - Seafloor survey and substrate maps from James Island to Ozette Lake, Washington Outer Coast

Without knowledge of basic seafloor characteristics, the ability to address any number of critical marine and/or coastal management issues is diminished. For example, management and conservation of essential fish habitat (EFH), a requirement mandated by federally guided fishery management plans (FMPs), requires among other things a description of habitats for federally managed species. Although the list of attributes important to habitat are numerous, the ability to efficiently and effectively describe many, and especially at the scales required, does not exist with the tools currently available. However, several characteristics of seafloor morphology are readily obtainable at multiple scales and can serve as useful descriptors of habitat. Recent advancements in acoustic technology, such as multibeam echosounding (MBES), can provide remote indication of surficial sediment properties such as texture, hardness, or roughness, and further permit highly detailed renderings of seafloor morphology. With acoustic-based surveys providing a relatively efficient method for data acquisition, there exists a need for efficient and reproducible automated segmentation routines to process the data. Using MBES data collected by the Olympic Coast National Marine Sanctuary (OCNMS), and through a contracted seafloor survey, we expanded on the techniques of Cutter et al. (2003) to describe an objective repeatable process that uses parameterized local Fourier histogram (LFH) texture features to automate segmentation of surficial sediments from acoustic imagery using a maximum likelihood decision rule. Sonar signatures and classification performance were evaluated using video imagery obtained from a towed camera sled. Segmented raster images were converted to polygon features and attributed using a hierarchical deep-water marine benthic classification scheme (Greene et al. 1999) for use in a geographical information system (GIS). (PDF contains 41 pages.)

Potential Application of Adsorptive Media to Enhance Phosphorus Uptake in Stormwater Basins and Wetlands at Lake Tahoe : literature review. Report to the University of California Davis Tahoe Research Group

Phosphorus removal by wetlands and basins in Lake Tahoe may be improved through designing these systems to filter storm water through media having higher phosphorus removal capabilities than local parent material. Substrates rich in iron, aluminum and calcium oftentimes have enhanced phosphorus removal. These substrates can be naturally occurring, byproducts of industrial or water treatment processes, or engineered. Phosphorus removal fundamentally occurs through chemical adsorption and/or precipitation and much of the phosphorus can be irreversibly bound. In addition to these standard media, other engineered substrates are available to enhance P removal. One such substrate is locally available in Reno and uses lanthanum coated diatomaceous earth for arsenate removal. This material, which has a high positive surface charge, can also irreversibly remove phosphorus. Physical factors also affect P removal. Specifically, specific surface area and particle shape affect filtration capacity, contact area between water and the surface area, and likelihood of clogging and blinding. A number of substrates have been shown to effectively remove P in case studies. Based upon these studies, promising substrates include WTRs, blast furnace slag, steel furnace slag, OPC, calcite, marble Utelite and other LWAs, zeolite and shale. However, other nonperformance factors such as environmental considerations, application logistics, costs, and potential for cementification narrow the list of possible media for application at Tahoe. Industrial byproducts such as slags risk possible leaching of heavy metals and this potential cannot be easily predicted. Fly ash and other fine particle substrates would be more difficult to apply because they would need to be blended, making them less desirable and more costly to apply than larger diameter media. High transportation costs rule out non-local products. Finally, amorphous calcium products will eventually cementify reducing their effectiveness in filtration systems. Based upon these considerations, bauxite, LWAs and expanded shales/clays, iron-rich sands, activated alumina, marble and dolomite, and natural and lanthanum activated diatomaceous earth are the products most likely to be tested for application at Tahoe. These materials are typically iron, calcium or aluminum based; many have a high specific surface area; and all have low transportation costs. (PDF contains 21 pages)

Comparison of resident and tourist preferences for public beach access

As coastal destinations continue to grow, due to tourism and residential expansion, the demand for public beach access and related amenities will also increase. As a resultagencies that provide beach access and related amenities face challenges when considering both residents and visitors use beaches and likely possess different needs, as well as different preferences for management decisions. Being a resident of a coastal county provides more opportunity to use local beaches, but coastal tourism is an important and growing economic engine in coastal communities (Kriesel, Landry, & Keeler, 2005; Pogue & Lee, 1999). Therefore, providing agencies with a comprehensive assessment of the differences between these two groups will increase the likelihood of effective management programs and policies for the provision of public beach access and related amenities. The purpose of this paper was to use a stated preference choice method (SPCM) to identify the extent of both residents’ and visitors’ preferences for public beach management options. (PDF contains 4 pages)

The survival and growth of hybrid Clarias gariepinus and Heterotis longifilis fed with live and frozen Dapthnids

Six glass tanks, each containing dechlorinated tap water and stocked with 3-day-old Heterobranchus, larvae hybrid between, Heterobranchus longifilis (male and Clarias gariepinus (female) were administered two types of feeding trials viz: live and frozen Dapthnids. Each treatment was replicated three times. The larvae were each fed approximately 50 Dapthnids per feeding time for fifteen days. Morphometric measurements of weigh and total length were taken before and after the experiment, and water quality parameters were monitored throughout the experimental period. At the end of the experiment, fish larvae fed frozen, Dapthnids, showed higher survival rate than Heterobranchus fed live Dapthnids. Even though the statistical analysis revealed that there was no significant difference (P>0.05) in survival. However, the Heteroclarias fed live dapthnids performed better in terms of growth rate than Heteroclarias fed frozen dapthnids. The statistical analysis did not show significant difference (P>0.05) between the final length and weight of two treatments. Heteroclarias fed live dapthnids had higher length and weight than Heteroclarias fed frozen dapthnids. It was therefore concluded that based on this experiment there is the likelihood that frozen zooplankton (Daphnids) do not encourage growth of Heteroclarias, but improves its survival. However, it's suggested that frozen zooplankton can be used to supplement live zooplankton in situation of sacrcity

Subspecies, morphs and clines in the amphipod Gammarus duebeni from fresh and saline waters

The nature and extent of morphological variation within populations of Gammarus duebeni are examined. The exceptional tolerance of G. duebeni to salinities that encompass three orders of magnitude was known in the 19th Century, and has attracted considerable attention from physiologists and ecologists in the 20th Century, including the likelihood that populations in freshwater are distinct from those living in more saline environments. It is concluded that gradual evolution into discrete and readily distinguished subspecies is currently underway, eventually producing several new species in freshwater and saline habitats.

From the individual to the community and beyond: water quality, stress indicators and key species in coastal waters

This review examines water quality and stress indicators at levels of organisation from the individual to the community and beyond by means of three case studies concentrating on rocky shores within the north-east Atlantic. Responses of dogwhelks (Nucella) to tributyltin pollution from antifouling paints is examined as the main case study. There are effects at the individual level (development of male sexual characteristics in the female leading to effective sterility) and population level (reduction in juveniles, few females and eventual population disappearance of dogwhelks in badly contaminated areas) but information on community level effects of dogwhelk demise is sparse. Such effects were simulated by dogwhelk removal experiments on well studied, moderately exposed ledges on shores on the Isle of Man. The removal of dogwhelks reduced the size and longevity of newly established Fucus clumps that had escaped grazing. Removal of dogwhelks also increased the likelihood of algal escapes. In a factorial experiment dogwhelks were shown to be less important than limpets \{Patella) in structuring communities but still had a significant modifying effect by increasing the probability of algal escapes. Community level responses to stress on rocky shores are then explored by reference to catastrophic impacts such as oil spills, using the Torrey Canyon as a case study. Recovery of the system in response to this major perturbation took between 10-15 years through a series of damped oscillations. The final case study is that of indicators of ecosystem level change in response to climate fluctuations, using ratios of northern \{Semibalanus balanoides) and southern (Chthamalus spp.) barnacles. Indices derived from counts on the shore show good correlations with inshore sea-water temperatures after a 2-year lag phase. The use of barnacles to measure offshore changes is reviewed. The discussion considers the use of bioindicators at various levels of organisation.

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.

An investigation into the potential losses of fish on the filter plates at the Hodder Water Treatment Works Stocks Reservoir

Stocks Reservoir is situated amidst the Forest of Bowland in the upper reaches of the old river valley of the Hodder. The reservoir was built in 1927 for the Fylde Water Board who primarily supplied water to Blackpool. The objective of this study is to assess the degree and likelihood of fish ingress onto the fish plates at the present and proposed stocking densities. An additional aim is to evaluate the operational implications, and if necessary suggest methods of alleviating the problem. Three spheres of study have been undertaken to achieve these objectives, these being: 1. To selectively stock the reservoir and monitor the angling club catches in order to assess the total population, relating it to fish plate losses and proposed stocking densities. 2. To monitor the fish taken from the fish plates and assess the reasons for their ingress. 3. To study the draw off tower and fish plates, and suggest ways of ameliorating or halting the loss of fish and consequent operational problems.

Demographics by depth: spatially explicit life-history dynamics of a protogynous reef fish

Distribution and demographics of the hogfish (Lachnolaimus maximus) were investigated by using a combined approach of in situ observations and life history analyses. Presence, density, size, age, and size and age at sex change all varied with depth in the eastern Gulf of Mexico. Hogfish (64–774 mm fork length and 0–19 years old) were observed year-round and were most common over complex, natural hard bottom habitat. As depth increased, the presence and density of hogfish decreased, but mean size and age increased. Size at age was smaller nearshore (<30 m). Length and age at sex change of nearshore hogfish were half those of offshore hogfish and were coincident with the minimum legal size limit. Fishing pressure is presumably greater nearshore and presents a confounding source of increased mortality; however, a strong red tide occurred the year before this study began and likely also affected nearshore demographics. Nevertheless, these data indicate ontogenetic migration and escapement of fast-growing fish to offshore habitat, both of which should reduce the likelihood of fishing-induced evolution. Data regarding the hogfish fishery are limited and regionally dependent, which has confounded previous stock assessments; however, the spatially explicit vital rates reported herein can be applied to future monitoring efforts.

Integrating methods for determining length-at-age to improve growth estimates for two large scombrids

Fish growth is commonly estimated from length-at-age data obtained from otoliths. There are several techniques for estimating length-at-age from otoliths including 1) direct observed counts of annual increments; 2) age adjustment based on a categorization of otolith margins; 3) age adjustment based on known periods of spawning and annuli formation; 4) back-calculation to all annuli, and 5) back-calculation to the last annulus only. In this study we compared growth estimates (von Bertalanffy growth functions) obtained from the above five methods for estimating length-at-age from otoliths for two large scombrids: narrow-barred Spanish mackerel (Scomberomorus commerson) and broad-barred king mackerel (Scomberomorus semifasciatus). Likelihood ratio tests revealed that the largest differences in growth occurred between the back-calculation methods and the observed and adjusted methods for both species of mackerel. The pattern, however, was more pronounced for S. commerson than for S. semifasciatus, because of the pronounced effect of gear selectivity demonstrated for S. commerson. We propose a method of substituting length-at-age data from observed or adjusted methods with back-calculated length-at-age data to provide more appropriate estimates of population growth than those obtained with the individual methods alone, particularly when faster growing young fish are disproportionately selected for. Substitution of observed or adjusted length-at-age data with back-calculated length-at-age data provided more realistic estimates of length for younger ages than observed or adjusted methods as well as more realistic estimates of mean maximum length than those derived from backcalculation methods alone.

Spatial and temporal variability in the relative contribution of king mackerel (Scomberomorus cavalla) stocks to winter mixed fisheries off South Florida

King mackerel (Scomberomorus cavalla) are ecologically and economically important scombrids that inhabit U.S. waters of the Gulf of Mexico (GOM) and Atlantic Ocean (Atlantic). Separate migratory groups, or stocks, migrate from eastern GOM and southeastern U.S. Atlantic to south Florida waters where the stocks mix during winter. Currently, all winter landings from a management-defined south Florida mixing zone are attributed to the GOM stock. In this study, the stock composition of winter landings across three south Florida sampling zones was estimated by using stock-specific otolith morphological variables and Fourier harmonics. The mean accuracies of the jackknifed classifications from stepwise linear discriminant function analysis of otolith shape variables ranged from 66−76% for sex-specific models. Estimates of the contribution of the Atlantic stock to winter landings, derived from maximum likelihood stock mixing models, indicated the contribution was highest off southeastern Florida (as high as 82.8% for females in winter 2001−02) and lowest off southwestern Florida (as low as 14.5% for females in winter 2002−03). Overall, results provided evidence that the Atlantic stock contributes a certain, and perhaps a significant (i.e., ≥50%), percentage of landings taken in the management-defined winter mixing zone off south Florida, and the practice of assigning all winter mixing zone landings to the GOM stock should