Relationship between water quality, watermilfoil frequency, and weevil distribution in the State of Washington

During the summer of 1997, we surveyed 50 waterbodies in Washington State to determine the distribution of the aquatic weevil Euhrychiopsis lecontei Dietz. We collected data on water quality and the frequency of occurrence of watermilfoil species within selected watermilfoil beds to compare the waterbodies and determine if they were related to the distribution E. lecontei . We found E. lecontei in 14 waterbodies, most of which were in eastern Washington. Only one lake with weevils was located in western Washington. Weevils were associated with both Eurasian ( Myriophyllum spicatum L.) and northern watermilfoil ( M. sibiricum K.). Waterbodies with E. lecontei had significantly higher ( P < 0.05) pH (8.7 ± 0.2) (mean ± 2SE), specific conductance (0.3 ± 0.08 mS cm -1 ) and total alkalinity (132.4 ± 30.8 mg CaCO 3 L -1 ). We also found that weevil presence was related to surface water temperature and waterbody location ( = 24.3, P ≤ 0.001) and of all the models tested, this model provided the best fit (Hosmer- Lemeshow goodness-of-fit = 4.0, P = 0.9). Our results suggest that in Washington State E. lecontei occurs primarily in eastern Washington in waterbodies with pH ≥ 8.2 and specific conductance ≥ 0.2 mS cm -1 . Furthermore, weevil distribution appears to be correlated with waterbody location (eastern versus western Washington) and surface water temperature.

Fishery resource surveys: beyond length-weight relationship models

Fisheries resource surveys are regular management tools for rational exploitation of commercial fisheries. In a growing number of cases, the use of these resource surveys has been largely restricted to assessment of the relative well being of fish stocks and the potential yields of such fisheries. This paper seeks to demonstrate that the data from such surveys can also be easily used to evaluate species diversity of such fisheries, both in terms of species richness and equitability of distribution. Using published data on two freshwater and two marine fisheries as case studies, Shannon-Wiener Diversity Function and Simpson's Index were computed for each of these fisheries. These biodiversity indices gave a deeper insight into the environmental status of each of these fisheries, beyond what the length-weight relationship models can reveal. Generally, while the marine fisheries showed more species richness, the freshwater fisheries apparently had more stable and equilibrated fish communities

Parameterizing probabilities for estimating age-composition distributions for mixture models

When estimating parameters that constitute a discrete probability distribution {pj}, it is difficult to determine how constraints should be made to guarantee that the estimated parameters { pˆj} constitute a probability distribution (i.e., pˆj>0, Σ pˆj =1). For age distributions estimated from mixtures of length-at-age distributions, the EM (expectationmaximization) algorithm (Hasselblad, 1966; Hoenig and Heisey, 1987; Kimura and Chikuni, 1987), restricted least squares (Clark, 1981), and weak quasisolutions (Troynikov, 2004) have all been used. Each of these methods appears to guarantee that the estimated distribution will be a true probability distribution with all categories greater than or equal to zero and with individual probabilities that sum to one. In addition, all these methods appear to provide a theoretical basis for solutions that will be either maximum-likelihood estimates or at least convergent to a probability distribut

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.

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.

Magnitude and distribution of sea turtle bycatch in the sea scallop (Placopecten magellanicus) dredge fishery in two areas of the northwestern Atlantic Ocean, 2001–2002

In May 2001, the National Marine Fisheries Service (NMFS) opened two areas in the northwestern Atlantic Ocean that had been previously closed to the U.S. sea scallop (Placopecten magellanicus) dredge fishery. Upon reopening these areas, termed the “Hudson Canyon Controlled Access Area” and the “Virginia Beach Controlled Access Area,” NMFS observers found that marine turtles were being caught incidentally in scallop dredges. This study uses the generalized linear model and the generalized additive model fitting techniques to identify environmental factors and gear characteristics that influence bycatch rates, and to predict total bycatch in these two areas during May-December 2001 and 2002 by incorporating environmental factors into the models. Significant factors affecting sea turtle bycatch were season, time-of-day, sea surface temperature, and depth zone. In estimating total bycatch, rates were stratified according to a combination of all these factors except time-of-day which was not available in fishing logbooks. Highest bycatch rates occurred during the summer season, in temperatures greater than 19°C, and in water depths from 49 to 57 m. Total estimated bycatch of sea turtles during May–December in 2001 and 2002 in both areas combined was 169 animals (CV=55.3), of which 164 (97%) animals were caught in the Hudson Canyon area. From these findings, it may be possible to predict hot spots for sea turtle bycatch in future years in the controlled access areas.

A general framework for integrating environmental time series into stock assessment models: model description, simulation testing, and example

We present a method to integrate environmental time series into stock assessment models and to test the significance of correlations between population processes and the environmental time series. Parameters that relate the environmental time series to population processes are included in the stock assessment model, and likelihood ratio tests are used to determine if the parameters improve the fit to the data significantly. Two approaches are considered to integrate the environmental relationship. In the environmental model, the population dynamics process (e.g. recruitment) is proportional to the environmental variable, whereas in the environmental model with process error it is proportional to the environmental variable, but the model allows an additional temporal variation (process error) constrained by a log-normal distribution. The methods are tested by using simulation analysis and compared to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. In the traditional method, the estimates of recruitment were provided by a model that allowed the recruitment only to have a temporal variation constrained by a log-normal distribution. We illustrate the methods by applying them to test the statistical significance of the correlation between sea-surface temperature (SST) and recruitment to the snapper (Pagrus auratus) stock in the Hauraki Gulf–Bay of Plenty, New Zealand. Simulation analyses indicated that the integrated approach with additional process error is superior to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. The results suggest that, for the snapper stock, recruitment is positively correlated with SST at the time of spawning.

Age, growth, mortality, and distribution of pinfish (Lagodon rhomboides) in Tampa Bay and adjacent Gulf of Mexico waters

Life-history dynamics of pinfish (Lagodon rhomboides) were examined from data derived from random station surveys conducted in Tampa Bay and adjacent Gulf of Mexico waters during 1993–97. In addition, patterns in spatial distribution and abundance in Gulf of Mexico waters were investigated. Ages determined from whole otoliths ranged from 0 to 7 years, and von Bertalanffy growth models for males and females were not significantly different. Von Bertalanffy growth model parameters were L∞=219.9 mm SL, k =0.33/yr, and t0 =–1.10 years for all fish combined. High gonadosomatic indices during October–December indicated that some spawning may occur in Tampa Bay. Estimated lengths at 50% maturity were 132 mm SL for males and 131 mm SL for females. Total instantaneous mortality rates derived from the Chapman-Robson estimator ranged from 0.88 to 1.08/yr, and natural mortality was estimated to be 0.78/yr. In Gulf of Mexico waters, pinfish catch rates declined with increasing depth, and most pinfish were caught in <17 m of water. Length distributions showed that pinfish segregate by size with increasing depth.