In search of climate effects on Atlantic Croaker (Micropogonias undulatus) stock off the U.S. Atlantic coast with Bayesian state-space biomass dynamic models

Atlantic Croaker (Micropogonias undulatus) production dynamics along the U.S. Atlantic coast are regulated by fishing and winter water temperature. Stakeholders for this resource have recommended investigating the effects of climate covariates in assessment models. This study used state-space biomass dynamic models without (model 1) and with (model 2) the minimum winter estuarine temperature (MWET) to examine MWET effects on Atlantic Croaker population dynamics during 1972–2008. In model 2, MWET was introduced into the intrinsic rate of population increase (r). For both models, a prior probability distribution (prior) was constructed for r or a scaling parameter (r0); imputs were the fishery removals, and fall biomass indices developed by using data from the Multispecies Bottom Trawl Survey of the Northeast Fisheries Science Center, National Marine Fisheries Service, and the Coastal Trawl Survey of the Southeast Area Monitoring and Assessment Program. Model sensitivity runs incorporated a uniform (0.01,1.5) prior for r or r0 and bycatch data from the shrimp-trawl fishery. All model variants produced similar results and therefore supported the conclusion of low risk of overfishing for the Atlantic Croaker stock in the 2000s. However, the data statistically supported only model 1 and its configuration that included the shrimp-trawl fishery bycatch. The process errors of these models showed slightly positive and significant correlations with MWET, indicating that warmer winters would enhance Atlantic Croaker biomass production. Inconclusive, somewhat conflicting results indicate that biomass dynamic models should not integrate MWET, pending, perhaps, accumulation of longer time series of the variables controlling the production dynamics of Atlantic Croaker, preferably including winter-induced estimates of Atlantic Croaker kills.

Structural dynamic eutrophication models

This article discusses problems of modelling the seasonal succession of algal species in lakes and reservoirs, and the adaptive selection of certain groups of algae in response to changes in the inputs and relative concentrations of nutrients and other environmental variables. A new generation of quantitative models is being developed which attempts to translate some important biological properties of species (survival, variation, inheritance, reproductive rates and population growth) into predictions about the survival of the fittest, where ”fitness” is measured or estimated in thermodynamic terms. The concept of ”exergy” and its calculation is explored to examine maximal exergy as a measure of fitness in ecosystems, and its use for calculating changes in species composition by means of structural dynamic models. These models accomodate short-term changes in parameters that affect the adaptive responses (species selection) of algae.

Modeling Environmental Stress

The word stress when applied to ecosystems is ambiguous. Stress may be low-level, with accompanying near-linear strain, or it may be of finite magnitude, with nonlinear response and possible disintegration of the system. Since there are practically no widely accepted definitions of ecosystem strain, classification of models of stressed systems is tenuous. Despite appearances, most ecosystem models seem to fall into the low-level linear response category. Although they sometimes simulate systems behavior well, they do not provide necessary and sufficient information about sudden structural changes nor structure after transition. Dynamic models of finiteamplitude response to stress are rare because of analytical difficulties. Some idea as to future transition states can be obtained by regarding the behavior of unperturbed functions under limiting strain conditions. Preliminary work shows that, since community variables do respond in a coherent manner to stress, macroscopic analyses of stressed ecosystems offer possible alternatives to compartmental models.

An examination of the usefulness of single-species biomass models for the management of Lake Victoria fisheries

Nile perch (Lates niloticus), tilapia (Oreochromis spp), dagaa (Rastrineobola argentea, silver cyprinid), and haplochromines (Tribe Haplochromini) form the backbone of the commercial fishery on Lake Victoria. These fish stocks account for about 70% of the total catch in the three riparian states Uganda, Kenya, and Tanzania. The lake fisheries have been poorly managed, in part due to inadequate scientific analysis and management advice. The overall objective of this project was to model the stocks of the commercial fisheries of Lake Victoria with the view of determining reference points and current stock status. The Schaefer biomass model was fitted to available data for each stock (starting in the 1960s or later) in the form of landings, catch per unit effort, acoustic survey indices, and trawl survey indices. In most cases, the Schaefer model did not fit all data components very well, but attempts were made to find the best model for each stock. When the model was fitted to the Nile perch data starting from 1996, the estimated current biomass is 654 kt (95% CI 466–763); below the optimum of 692 kt and current harvest rate is 38% (33–73%), close to the optimum of 35%. At best, these can be used as tentative guidelines for the management of these fisheries. The results indicate that there have been strong multispecies interactions in the lake ecosystem. The findings from our study can be used as a baseline reference for future studies using more complex models, which could take these multispecies interactions into account.

Stock assessment of protogynous fish: evaluating measures of spawning biomass used to estimate biological reference points

In stock assessments, recruitment is typically modeled as a function of females only. For protogynous stocks, however, disproportionate fishing on males increases the possibility of reduced fertilization rates. To incorporate the importance of males in protogynous stocks, assessment models have been used to predict recruitment not just from female spawning biomass (Sf), but also from that of males (Sm) or both sexes (Sb). We conducted a simulation study to evaluate the ability of these three measures to estimate biological reference points used in fishery management. Of the three, Sf provides best estimates if the potential for decreased fertilization is weak, whereas Sm is best only if the potential is very strong. In general, Sb estimates the true reference points most closely, which indicates that if the potential for decreased fertilization is moderate or unknown, Sb should be used in assessments of protogynous stocks. Moreover, for a broad range of scenarios, relative errors from Sf and Sb occur in opposite directions, indicating that estimates from these measures could be used to bound uncertainty.

Active reservoir management: a model solution

Steady-state procedures, of their very nature, cannot deal with dynamic situations. Statistical models require extensive calibration, and predictions often have to be made for environmental conditions which are often outside the original calibration conditions. In addition, the calibration requirement makes them difficult to transfer to other lakes. To date, no computer programs have been developed which will successfully predict changes in species of algae. The obvious solution to these limitations is to apply our limnological knowledge to the problem and develop functional models, so reducing the requirement for such rigorous calibration. Reynolds has proposed a model, based on fundamental principles of algal response to environmental events, which has successfully recreated the maximum observed biomass, the timing of events and a fair simulation of the species succession in several lakes. A forerunner of this model was developed jointly with Welsh Water under contract to Messrs. Wallace Evans and Partners, for use in the Cardiff Bay Barrage study. In this paper the authors test a much developed form of this original model against a more complex data-set and, using a simple example, show how it can be applied as an aid in the choice of management strategy for the reduction of problems caused by eutrophication. Some further developments of the model are indicated.

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.

Comparison of habitat-based indices of abundance with fishery-independent biomass estimates from bottom trawl surveys

Rockfish species are notoriously difficult to sample with multispecies bottom trawl survey methods. Typically, biomass estimates have high coefficients of variation and can fluctuate outside the bounds of biological reality from year to year. This variation may be due in part to their patchy distribution related to very specific habitat preferences. We successfully modeled the distribution of five commercially important and abundant rockf ish species. A two-stage modeling method (modeling both presence-absence and abundance) and a collection of important habitat variables were used to predict bottom trawl survey catch per unit of effort. The resulting models explained between 22% and 66% of the variation in rockfish distribution. The models were largely driven by depth, local slope, bottom temperature, abundance of coral and sponge, and measures of water column productivity (i.e., phytoplankton and zooplankton). A year-effect in the models was back-transformed and used as an index of the time series of abundance. The abundance index trajectories of three of five species were similar to the existing estimates of their biomass. In the majority of cases the habitat-based indices exhibited less interannual variability and similar precision when compared with stratified survey-based biomass estimates. These indices may provide for stock assessment models a more stable alternative to current biomass estimates produced by the multispecies bottom trawl survey in the Gulf of Alaska.

Life-stage recruitment models for Atlantic cod (Gadus morhua) and haddock (Melanogrammus aeglefinus) on Georges Bank

Ichthyoplankton surveys have been used to provide an independent estimate of adult spawning biomass of commercially exploited species and to further our understanding of the recruitment processes in the early life stages. However, predicting recruitment has been difficult because of the complex interaction of physical and biological processes operating at different spatial and temporal scales that can occur at the different life stages. A model of first-year life-stage recruitment was applied to Georges Bank Atlantic cod (Gadus morhua) and haddock (Melanogrammus aeglefinus) stocks over the years 1977–2004 by using environmental and densitydependent relationships. The best lifestage mortality relationships for eggs, larvae, pelagic juveniles, and demersal juveniles were first determined by hindcasting recruitment estimates based on egg and larval abundance and mortality rates derived from two intensive sampling periods, 1977–87 and 1995–99. A wind-driven egg mortality relationship was used to estimate losses due to transport off the bank, and a wind-stress larval mortality relationship was derived from feeding and survival studies. A simple metric for the density-dependent effects of Atlantic cod was used for both Atlantic cod and haddock. These life stage proxies were then applied to the virtual population analysis (VPA) derived annual egg abundances to predict age-1 recruitment. Best models were determined from the correlation of predicted and VPA-derived age-1 abundance. The larval stage was the most quantifiable of any stage from surveys, whereas abundance estimates of the demersal juvenile stage were not available because of undersampling. Attempts to forecast recruitment from spawning stock biomass or egg abundance, however, will always be poor because of variable egg survival.

Interactions of age-dependent mortality and selectivity functions in age-based stock assessment models

The natural mortality rate (M) of fish varies with size and age, although it is often assumed to be constant in stock assessments. Misspecification of M may bias important assessment quantities. We simulated fishery data, using an age-based population model, and then conducted stock assessments on the simulated data. Results were compared to known values. Misspecification of M had a negligible effect on the estimation of relative stock depletion; however, misspecification of M had a large effect on the estimation of parameters describing the stock recruitment relationship, age-specific selectivity, and catchability. If high M occurs in juvenile and old fish, but is misspecified in the assessment model, virgin biomass and catchability are often poorly estimated. In addition, stock recruitment relationships are often very difficult to estimate, and steepness values are commonly estimated at the upper bound (1.0) and overfishing limits tend to be biased low. Natural mortality can be estimated in assessment models if M is constant across ages or if selectivity is asymptotic. However if M is higher in old fish and selectivity is dome-shaped, M and the selectivity cannot both be adequately estimated because of strong interactions between M and selectivity.

Measurement errors in body size of sea scallops (Placopecten magellanicus) and their effect on stock assessment models

Body-size measurement errors are usually ignored in stock assessments, but may be important when body-size data (e.g., from visual sur veys) are imprecise. We used experiments and models to quantify measurement errors and their effects on assessment models for sea scallops (Placopecten magellanicus). Errors in size data obscured modes from strong year classes and increased frequency and size of the largest and smallest sizes, potentially biasing growth, mortality, and biomass estimates. Modeling techniques for errors in age data proved useful for errors in size data. In terms of a goodness of model fit to the assessment data, it was more important to accommodate variance than bias. Models that accommodated size errors fitted size data substantially better. We recommend experimental quantification of errors along with a modeling approach that accommodates measurement errors because a direct algebraic approach was not robust and because error parameters were diff icult to estimate in our assessment model. The importance of measurement errors depends on many factors and should be evaluated on a case by case basis.

Deriving acceptable biological catch from the overfishing limit: implications for assessment models

The recently revised Magnuson–Stevens Fishery Conservation and Management Act requires that U.S. fishery management councils avoid overfishing by setting annual catch limits (ACLs) not exceeding recommendations of the councils’ scientific advisers. To meet that requirement, the scientific advisers will need to know the overfishing limit (OFL) estimated in each stock assessment, with OFL being the catch available from applying the limit fishing mortality rate to current or projected stock biomass. The advisers then will derive ‘‘acceptable biological catch’’ (ABC) from OFL by reducing OFL to allow for scientific uncertainty, and ABC becomes their recommendation to the council. We suggest methodology based on simple probability theory by which scientific advisers can compute ABC from OFL and the statistical distribution of OFL as estimated by a stock assessment. Our method includes approximations to the distribution of OFL if it is not known from the assessment; however, we find it preferable to have the assessment model estimate the distribution of OFL directly. Probability-based methods such as this one provide well-defined approaches to setting ABC and may be helpful to scientific advisers as they translate the new legal requirement into concrete advice.

Pinniped diet composition: a comparison of estimation models

Along the west coast of the United States, the potential impact of increasing pinniped populations on declining salmonid (Oncorhynchus spp.) stocks has become an issue of concern. Fisheries managers need species-specific estimates of consumption by pinnipeds to evaluate their impact on salmonid stocks. To estimate consumption, we developed a model that estimates diet composition by reconstructing prey biomass from fecal samples. We applied the model to data collected from harbor seals (Phoca vitulina) that are present year-round in the lower Columbia River where endangered stocks of salmonids pass as returning adults and as seaward-migrating smolts. Using the same data, we applied the split-sample frequency of occurrence model, which avoids reconstructing biomass by assuming that each fecal sample represents an equal volume of consumption and that within each sample each prey item represents an equal proportion of the volume. The two models for estimating diet composition yielded size-specific differences in consumption estimates that were as large as tenfold for the smallest and largest prey. Conclusions about the impact of harbor seal predation on adult salmonids, some of their largest prey species, remain uncertain without some appropriate rationale or further information (e.g. empirical captive studies) to discriminate between these models.

The study of biological characteristic and population dynamic of Caspian vimba (Vimba vimba persa), in coastal waters of Caspian Sea (Mazandaran Province)

This study was conducted to determine reproduction characteristics, diet regime, age structure and population dynamics parameters of the vimba vimba persa (Pallas, 1811) in Mazandaran waters of the Caspian Sea, from October 2008 to September 2009. A total of 994 specimens were monthly collected by beach seine and cast net from six fish landings of Ramsar, Tonekabon, Chaloos, Mahmood Abad, Sari and Behshahr. Biometric characters were measured for each specimen at the laboratory. Scales were used for age determination. Sex determination and fecundity were determined. Population dynamic parameters as well as stock assessment including cohort analysis were estimated using FISAT software. The finding showed that the mean of fork length and body weight of the Caspian Vimba were 168.4±2.6 mm and 71.94±32.24 g respectively. Strong correlation was found between these two variables (a= 0.012; b = 3.047; r2 = 0.955). 92 specimens were studied from the fecundity point of view. This species was found to have more abundance in spring (esp. Apr-May). The samples composed of 397(42.6%) male, 537(57.4%) female; Overall sex ratio (M: F =1: 1.35) was significantly different from the expected 1:1 ratio (p ≤0.05). The advanced stages of maturity (4th & 5th) were found in April and May. The highest Gonadosomatic Index in female was in May and the lowest one was in July. This fish is therefore a spring spawner. The maximum absolute and relative fecundities were 34640 and 260.9, respectively; the minimum absolute and relative fecundities were 5400 and 94.5 respectively. The averages of absolute and relative fecundities were 17198±7710 and 171.85±48.8, respectively. Coefficient vacuity index was 59.2% which indicates that this fish is mesophagous. Among of living creature consumes by Caspian Vimba mollusks, 76 arthropods, worms, plants, detritus and fishes were found 32.9% , 26.7% , 13.4% , 17% , 4.4% and 1.6% respectively. The infinite fork lengths were 261 mm for females, 25mm for males and 261 mm for both sexes respectively. For population growth and mortality parameters; K ( 0.28 per year for both sexes, 0.3 per year for males, 0.33 per year for females); t0 ( -0.65 year for both sexes, -0.23 year in females, -0.51 year in males ); Φ' ( 2.28 ); Z ( 0.98 per year ); M ( 0.59 per year); F ( 0.39 per year) and Exploitation coefficient was 0.4. The analysis showed that total biomass and MSY were 1336 and 528.8 tonnes respectively.