Statistical inference about the relative efficiency of a new survey protocol, based on paired-tow survey calibration data

Paired-tow calibration studies provide information on changes in survey catchability that may occur because of some necessary change in protocols (e.g., change in vessel or vessel gear) in a fish stock survey. This information is important to ensure the continuity of annual time-series of survey indices of stock size that provide the basis for fish stock assessments. There are several statistical models used to analyze the paired-catch data from calibration studies. Our main contributions are results from simulation experiments designed to measure the accuracy of statistical inferences derived from some of these models. Our results show that a model commonly used to analyze calibration data can provide unreliable statistical results when there is between-tow spatial variation in the stock densities at each paired-tow site. However, a generalized linear mixed-effects model gave very reliable results over a wide range of spatial variations in densities and we recommend it for the analysis of paired-tow survey calibration data. This conclusion also applies if there is between-tow variation in catchability.

Modified sorting technique to mitigate the collateral mortality of trawled school prawns (Metapenaeus macleayi)

The potential for changes to onboard handling practices in order to improve the fate of juvenile school prawns (Metapenaeus macleayi) discarded during trawling were investigated in two Australian rivers (Clarence and Hunter) by comparing a purpose-built, water-filled sorting tray against a conventional dry tray across various conditions, including the range of typical delays before the start of sorting the catch (2 min vs. 15 min). Juvenile school prawns (n= 5760), caught during 32 and 16 deployments in each river, were caged and sacrificed at four times: immediately (T0), and at 24 (T24), 72 (T72), and 120 (T12 0) hours after having been discarded. In both rivers, most mortalities occurred between T0 and T24 and, after adjusting for control deaths (<12%), were greatest for the 15-min conventional treatment (up to 41% at T120). Mixed-effects logistic models revealed that in addition to the sampling time, method of sorting, and delay in sorting, the weight of the catch, salinity, and percentage cloud cover were significant predictors of mortality. Although trawling caused some mortalities and comparable stress (measured as L -lactate) in all school prawns, use of the water tray lessened the negative impacts of some of the above factors across both the 2-min and 15-min delays in sorting so that the overall discard mortality was reduced by more than a third. When used in conjunction with selective trawls, widespread application of the water tray should help to improve the sustainability of trawling for school prawns.

Effects of a large fishing closure on benthic communities in the western Gulf of Maine: recovery from the effects of gillnets and otter trawls

The recovery of benthic communities inside the western Gulf of Maine fishing closure area was evaluated by comparing invertebrate assemblages at sites inside and outside of the closure four to six years after the closure was established. The major restriction imposed by the closure was a year-round prohibition of bottom gillnets and otter trawls. A total of 163 seafloor sites (~half inside and half outside the closure) within a 515-km2 study area were sampled with some combination of Shipek grab, Wildco box corer, or underwater video. Bottom types ranged from mud (silt and clay) to boulders, and the effects of the closure on univariate measures (total density, biomass, taxonomic richness) of benthos varied widely among sediment types. For sites with predominantly mud sediments, there were mixed effects on inside and outside infauna and no effect on epifauna. For sites with mainly sand sediments, there were higher density, biomass, and taxonomic richness for infauna inside the closure, but no significant effects on epifauna. For sites dominated by gravel (which included boulders in some areas), there were no effects on infauna but strong effects on epifaunal density and taxonomic richness. For fishing gear, the data indicated that infauna recovered in sand from the impacts of otter trawls operated inside the closure but that they did not recover in mud, and that epifauna recovered on gravel bottoms from the impact of gillnets used inside the closure. The magnitudes of impact and recovery, however, cannot be inferred directly from our data because of a confounding factor of different fishing intensities outside the closure for a direct comparison of preclosure and postclosure data. The overall negative impact of trawls is likely underestimated by our data, whereas the negative impact of gillnets is likely overestimated.

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.

Assessment of Atlantic Red Drum for 1999: Northern and southern regions

An assessment of the status of the Atlantic stock of red drum is conducted using recreational and commercial data from 1986 through 1998. This assessment updates data and analyses from the 1989, 1991, 1992 and 1995 stock assessments on Atlantic coast red drum (Vaughan and Helser, 1990; Vaughan 1992; 1993; 1996). Since 1981, coastwide recreational catches ranged between 762,300 pounds in 1980 and 2,623,900 pounds in 1984, while commercial landings ranged between 60,900 pounds in 1997 and 422,500 pounds in 1984. In weight of fish caught, Atlantic red drum constitute predominantly a recreational fishery (ranging between 85 and 95% during the 1990s). Commercially, red drum continue to be harvested as part of mixed species fisheries. Using available length-frequency distributions and age-length keys, recreational and commercial catches are converted to catch in numbers at age. Separable and tuned virtual population analyses are conducted on the catch in numbers at age to obtain estimates of fishing mortality rates and population size (including recruitment to age 1). In tum, these estimates of fishing mortality rates combined with estimates of growth (length and weight), sex ratios, sexual maturity and fecundity are used to estimate yield per recruit, escapement to age 4, and static (or equilibrium) spawning potential ratio (static SPR, based on both female biomass and egg production). Three virtual analysis approaches (separable, spreadsheet, and FADAPT) were applied to catch matrices for two time periods (early: 1986-1991, and late: 1992-1998) and two regions (Northern: North Carolina and north, and Southern: South Carolina through east coast of Florida). Additional catch matrices were developed based on different treatments for the catch-and-release recreationally-caught red drum (B2-type). These approaches included assuming 0% mortality (BASEO) versus 10% mortality for B2 fish. For the 10% mortality on B2 fish, sizes were assumed the same as caught fish (BASEl), or positive difference in size distribution between the early period and the later period (DELTA), or intermediate (PROP). Hence, a total of 8 catch matrices were developed (2 regions, and 4 B2 assumptions for 1986-1998) to which the three VPA approaches were applied. The question of when offshore emigration or reduced availability begins (during or after age 3) continues to be a source of bias that tends to result in overestimates of fishing mortality. Additionally, the continued assumption (Vaughan and Helser, 1990; Vaughan 1992; 1993; 1996) of no fishing mortality on adults (ages 6 and older), causes a bias that results in underestimates of fishing mortality for adult ages (0 versus some positive value). Because of emigration and the effect of the slot limit for the later period, a range in relative exploitations of age 3 to age 2 red drum was considered. Tuning indices were developed from the MRFSS, and state indices for use in the spreadsheet and FADAPT VPAs. The SAFMC Red Drum Assessment Group (Appendix A) favored the FADAPT approach with catch matrix based on DELTA and a selectivity for age 3 relative to age 2 of 0.70 for the northern region and 0.87 for the southern region. In the northern region, estimates of static SPR increased from about 1.3% for the period 1987-1991 to approximately 18% (15% and 20%) for the period 1992-1998. For the southern region, estimates of static SPR increased from about 0.5% for the period 1988-1991 to approximately 15% for the period 1992-1998. Population models used in this assessment (specifically yield per recruit and static spawning potential ratio) are based on equilibrium assumptions: because no direct estimates are available as to the current status of the adult stock, model results imply potential longer term, equilibrium effects. Because current status of the adult stock is unknown, a specific rebuilding schedule cannot be determined. However, the duration of a rebuilding schedule should reflect, in part, a measure of the generation time of the fish species under consideration. For a long-lived, but relatively early spawning, species as red drum, mean generation time would be on the order of 15 to 20 years based on age-specific egg production. Maximum age is 50 to 60 years for the northern region, and about 40 years for the southern region. The ASMFC Red Drum Board's first phase recovery goal of increasing %SPR to at least 10% appears to have been met. (PDF contains 79 pages)

Effects of macrobenthos on sediment-water oxygen and ammonium fluxes. Final Report.

Sediments are an important location in determining the fate of nutrients entering the estuary. Role of sediments needs to be incorporated into water quality models. Purpose of this study was to estimate the portion of sediment oxygen consumption (SOC) and sediment ammonium (NH4+) release directly attributable to benthic invertebrates via the respiratory use of oxygen and catabolic release of ammonium. Samples were collected at 8 locations from August 1985 through November 1988. (PDF contains 45 pages)

Predicting the effects of climate change on sea turtle nesting habitat in Florida

Rising global temperatures threaten the survival of many plant and animal species. Having already risen at an unprecedented rate in the past century, temperatures are predicted to rise between 0.3 and 7.5C in North America over the next 100 years (Hawkes et al. 2007). Studies have documented the effects of climate warming on phenology (timing of seasonal activities), with observations of early arrival at breeding grounds, earlier ends to the reproductive season, and delayed autumnal migrations (Pike et al. 2006). In addition, for species not suited to the physiological demands of cold winter temperatures, increasing temperatures could shift tolerable habitats to higher latitudes (Hawkes et al. 2007). More directly, climate warming will impact thermally sensitive species like sea turtles, who exhibit temperature-dependent sexual determination. Temperatures in the middle third of the incubation period determine the sex of sea turtle offspring, with higher temperatures resulting in a greater abundance of female offspring. Consequently, increasing temperatures from climate warming would drastically change the offspring sex ratio (Hawkes et al. 2007). Of the seven extant species of sea turtles, three (leatherback, Kemp’s ridley, and hawksbill) are critically endangered, two (olive ridley and green) are endangered, and one (loggerhead) is threatened. Considering the predicted scenarios of climate warming and the already tenuous status of sea turtle populations, it is essential that efforts are made to understand how increasing temperatures may affect sea turtle populations and how these species might adapt in the face of such changes. In this analysis, I seek to identify the impact of changing climate conditions over the next 50 years on the availability of sea turtle nesting habitat in Florida given predicted changes in temperature and precipitation. I predict that future conditions in Florida will be less suitable for sea turtle nesting during the historic nesting season. This may imply that sea turtles will nest at a different time of year, in more northern latitudes, to a lesser extent, or possibly not at all. It seems likely that changes in temperature and precipitation patterns will alter the distribution of sea turtle nesting locations worldwide, provided that beaches where the conditions are suitable for nesting still exist. Hijmans and Graham (2006) evaluate a range of climate envelope models in terms of their ability to predict species distributions under climate change scenarios. Their results suggested that the choice of species distribution model is dependent on the specifics of each individual study. Fuller et al. (2008) used a maximum entropy approach to model the potential distribution of 11 species in the Arctic Coastal Plain of Alaska under a series of projected climate scenarios. Recently, Pike (in press) developed Maxent models to investigate the impacts of climate change on green sea turtle nest distribution and timing. In each of these studies, a set of environmental predictor variables (including climate variables), for which ‘current’ conditions are available and ‘future’ conditions have been projected, is used in conjunction with species occurrence data to map potential species distribution under the projected conditions. In this study, I will take a similar approach in mapping the potential sea turtle nesting habitat in Florida by developing a Maxent model based on environmental and climate data and projecting the model for future climate data. (PDF contains 5 pages)

Testable models of aquatic systems: A NERC special topic for integrating experiments and modelling

A major part of the support for fundamental research on aquatic ecosystems continues to be provided by the Natural Environment Research Council (NERC). Funds are released for ”thematic” studies in a selected special topic or programme. ”Testable Models of Aquatic Ecosystems” was a Special Topic of the NERC, initiated in 1995, the aim of which was to promote ecological modelling by making new links between experimental aquatic biologists and state-of-the-art modellers. The Topic covered both marine and freshwater systems. This paper summarises projects on aspects of the responses of individual organisms to the effects of environmental variability, on the assembly, permanence and resilience of communities, and on aspects of spatial models. The authors conclude that the NERC Special Topic has been highly successful in promoting the development and application of models, most particularly through the interplay between experimental ecologists and formal modellers.

The predicting power of models for eutrophication

Restoration of water-bodies from eutrophication has proved to be extremely difficult. Mathematical models have been used extensively to provide guidance for management decisions. The aim of this paper is to elucidate important problems of using models for predicting environmental changes. First, the necessity for a proper uncertainty assessment of the model, upon calibration, has not been widely recognized. Predictions must not be a single time trajectory; they should be a band, expressing system uncertainty and natural variability. Availability of this information may alter the decision to be taken. Second, even with well-calibrated models, there is no guarantee they will give correct projections in situations where the model is used to predict the effects of measures designed to bring the system into an entirely different ”operating point”, as is typically the case in eutrophication abatement. The concept of educated speculation is introduced to partially overcome this difficulty. Lake Veluwe is used as a case to illustrate the point. Third, as questions become more detailed, such as ”what about expected algal composition”, there is a greater probability of running into fundamental problems that are associated with predicting the behaviour of complex non-linear systems. Some of these systems show extreme initial condition sensitivity and even, perhaps, chaotic behaviour, and are therefore fundamentally unpredictable.

Mathematical models for the prediction of temperature distributions resulting from the discharge of heated water into large bodies of water

Mathematical models for heated water outfalls were developed for three flow regions. Near the source, the subsurface discharge into a stratified ambient water issuing from a row of buoyant jets was solved with the jet interference included in the analysis. The analysis of the flow zone close to and at intermediate distances from a surface buoyant jet was developed for the two-dimensional and axisymmetric cases. Far away from the source, a passive dispersion model was solved for a two dimensional situation taking into consideration the effects of shear current and vertical changes in diffusivity. A significant result from the surface buoyant jet analysis is the ability to predict the onset and location of an internal hydraulic jump. Prediction can be made simply from the knowledge of the source Froude number and a dimensionless surface exchange coefficient. Parametric computer programs of the above models are also developed as a part of this study. This report was submitted in fulfillment of Contract No. 14-12-570 under the sponsorship of the Federal Water Quality Administration.

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.

The effects of temperature on hatching and survival of northern rock sole larvae (Lepidopsetta polyxystra)

Northern rock sole (Lepidopsetta polyxystra) is a commercially important flatfish in Alaska and was recently classified as a distinct species from southern rock sole (L. bilineata). Taxonomic and vital rate data for northern rock sole are still not fully described, notably at early egg and larval stages. In this study, we provide new taxonomic descriptions of late-stage eggs and newly hatched larvae, as well as temperature-response models of hatching (timing, duration, success), and larval size-at-hatch and posthatch survival at four temperatures (2°, 5°, 9°, and 12°C). Time-to-first-hatch, hatch cycle duration, and overall hatching success showed a negative relationship with temperature. Early hatching larvae within each temperature treatment were smaller and had larger yolk sacs, but larvae incubated at higher temperatures (9° and 12°C) had the largest yolk reserves overall. Despite having smaller yolks, size-at-hatch and the maximum size achieved during the hatching cycle was highest for larvae reared at cold temperatures (2° and 5°C), indicating that endogenous reserves are more efficiently used for growth at these temperatures. In addition, larvae reared at high temperatures died more rapidly in the absence of food despite having more yolk reserves than cold-incubated larvae. Overall, northern rock sole eggs and larvae display early life history traits consistent with coldwater adaptation for winter spawning in the North Pacific.

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.

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