Removing observational noise from fisheries-independent time series data using ARIMA models

Abundance indices derived from fishery-independent surveys typically exhibit much higher interannual variability than is consistent with the within-survey variance or the life history of a species. This extra variability is essentially observation noise (i.e. measurement error); it probably reflects environmentally driven factors that affect catchability over time. Unfortunately, high observation noise reduces the ability to detect important changes in the underlying population abundance. In our study, a noise-reduction technique for uncorrelated observation noise that is based on autoregressive integrated moving average (ARIMA) time series modeling is investigated. The approach is applied to 18 time series of finfish abundance, which were derived from trawl survey data from the U.S. northeast continental shelf. Although the a priori assumption of a random-walk-plus-uncorrelated-noise model generally yielded a smoothed result that is pleasing to the eye, we recommend that the most appropriate ARIMA model be identified for the observed time series if the smoothed time series will be used for further analysis of the population dynamics of a species.

Age validation, growth modeling, and mortality estimates for striped trumpeter (Latris lineata) from southeastern Australia: making the most of patchy data

Age estimates for striped trumpeter (Latris lineata) from Tasmanian waters were produced by counting annuli on the transverse section of sagittal otoliths and were validated by comparison of growth with known-age individuals and modal progression of a strong recruitment pulse. Estimated ages ranged from one to 43 years; fast growth rates were observed for the first five years. Minimal sexual dimorphism was shown to exist between length, weight, and growth characteristics of striped trumpeter. Seasonal growth variability was strong in individuals up to at least age four, and growth rates peaked approximately one month after the observed peak in sea surface temperature. A modified two-phase von Bertalanffy growth function was fitted to the length-at-age data, and the transition between growth phases was linked to apparent changes in physiological and life history traits, including offshore movement as fish approach maturity. The two-phase curve was found to represent the mean length at age in the data better than the standard von Bertalanffy growth function. Total mortality was estimated by using catch curve analysis based on the standard and two-phase von Bertalanffy growth functions, and estimates of natural mortality were calculated by using two empirical models, one based on longevity and the other based on the parameters L∞ and k from both growth functions. The interactions between an inshore gillnet fishery targeting predominately juveniles and an offshore hook fishery targeting predominately adults highlight the need to use a precautionary approach when developing harvest strategies.

The A Posteriori Aspects of Estuarine Modeling

This exercise is the application of an analytical method for systematically modeling ecosystems data to observations made on a naturally eutrophic, mesohaline planktonic microcosm. The theory and experimental design are briefly outlined and the particular steps in the acutal modeling process follow. Then there is a discussion as to how the whole endeavor can be refined to culminate in models with predictive capabilities. (PDF has 16 pages.)

Establishing National Ocean Service Priorities for Estuarine, Coastal, and Ocean Modeling: Capabilities, Gaps, and Preliminary Prioritization Factors

This report was developed to help establish National Ocean Service priorities and chart new directions for research and development of models for estuarine, coastal and ocean ecosystems based on user-driven requirements and supportive of sound coastal management, stewardship, and an ecosystem approach to management. (PDF contains 63 pages)

Panel discussion: U.S. EPA using modeling and ecosystem services to enhance coastal decision making

This panel will discuss the research being conducted, and the models being used in three current coastal EPA studies being conducted on ecosystem services in Tampa Bay, the Chesapeake Bay and the Coastal Carolinas. These studies are intended to provide a broader and more comprehensive approach to policy and decision-making affecting coastal ecosystems as well as provide an account of valued services that have heretofore been largely unrecognized. Interim research products, including updated and integrated spatial data, models and model frameworks, and interactive decision support systems will be demonstrated to engage potential users and to elicit feedback. It is anticipated that the near-term impact of the projects will be to increase the awareness by coastal communities and coastal managers of the implications of their actions and to foster partnerships for ecosystem services research and applications. (PDF contains 4 pages)

A modeling approach to identify optimal habitat and suitable fishing grounds for neon flying squid (Ommastrephes bartramii) in the Northwest Pacific Ocean

We developed a habitat suitability index (HSI) model to understand and identify the optimal habitat and potential fishing grounds for neon f lying squid (Ommastrephes bartramii) in the Northwest Pacific Ocean. Remote sensing data, including sea surface temperature, sea surface salinity, sea surface height, and chlorophyll-a concentrations, as well as fishery data from Chinese mainland squid f leets in the main fishing ground (150–165°E longitude) from August to October, from 1999 to 2004, were used. The HSI model was validated by using fishery data from 2005. The arithmetic mean modeling with three of the environmental variables—sea surface temperature, sea surface height anomaly, and chlorophyll- a concentrations—was defined as the most parsimonious HSI model. In 2005, monthly HSI values >0.6 coincided with productive fishing grounds and high fishing effort from August to October. This result implies that the model can reliably predict potential f ishing grounds for O. bartramii. Because spatially explicit fisheries and environmental data are becoming readily available, it is feasible to develop a dynamic, near real-time habitat model for improving the process of identifying potential fishing areas for and optimal habitats of neon flying squid.

Use of genetic data to assess the uncertainty in stock assessments due to the assumed stock structure: the case of albacore (Thunnus alalunga) from the Atlantic Ocean

Stock assessments can be problematic because of uncertainties associated with the data or because of simplified assumptions made when modeling biological processes (Rosenberg and Restrepo, 1995). For example, the common assumption in stock assessments that stocks are homogeneous and discrete (i.e., there is no migration between the stocks) is not necessarily true (Kell et al., 2004a, 2004b).

A production modeling approach to the assessment of the horseshoe crab (Limulus polyphemus) population in Delaware Bay

Horseshoe crab (Limulus polyphemus) is harvested commercially, used by the biomedical industry, and provides food for migrating shorebirds, particularly in Delaware Bay. Recently, decreasing crab population trends in this region have raised concerns that the stock may be insufficient to fulfill the needs of these diverse user groups. To assess the Delaware Bay horseshoe crab population, we used surplus production models (programmed in ASPIC), which incorporated data from fishery-independent surveys, fishery-dependent catch-per-unit-of-effort data, and regional harvest. Results showed a depleted population (B2003/=0.03−0.71) BMSY and high relative fishing mortality /FMSY=0.9−9.5). Future harvest (F2002strategies for a 15-year period were evaluated by using population projections with ASPICP software. Under 2003 harvest levels (1356 t), population recovery to BMSY would take at least four years, and four of the seven models predicted that the population would not reach BMSY within the 15year period. Production models for horseshoe crab assessment provided management benchmarks for a species with limited data and no prior stock assessment

Quantifying intrapopulation variability in stable isotope data for Spotted Seatrout (Cynoscion nebulosus)

Stable isotope (SI) values of carbon (δ13C) and nitrogen (δ15N) are useful for determining the trophic connectivity between species within an ecosystem, but interpretation of these data involves important assumptions about sources of intrapopulation variability. We compared intrapopulation variability in δ13C and δ15N for an estuarine omnivore, Spotted Seatrout (Cynoscion nebulosus), to test assumptions and assess the utility of SI analysis for delineation of the connectivity of this species with other species in estuarine food webs. Both δ13C and δ15N values showed patterns of enrichment in fish caught from coastal to offshore sites and as a function of fish size. Results for δ13C were consistent in liver and muscle tissue, but liver δ15N showed a negative bias when compared with muscle that increased with absolute δ15N value. Natural variability in both isotopes was 5–10 times higher than that observed in laboratory populations, indicating that environmentally driven intrapopulation variability is detectable particularly after individual bias is removed through sample pooling. These results corroborate the utility of SI analysis for examination of the position of Spotted Seatrout in an estuarine food web. On the basis of these results, we conclude that interpretation of SI data in fishes should account for measurable and ecologically relevant intrapopulation variability for each species and system on a case by case basis.

Modeling seabird group size: implications for ecological impact assessments

The purpose of this project is to model seabird flock size data to provide recommendations to the Bureau of Ocean and Energy Management for offshore wind turbine placement. Our hypothesis is that ecological characteristics influence which statistical distribution will provide the best fit to seabird flock size data. To test this, seabird species can be grouped based on shared ecological traits, such as foraging mechanism or diet.

An investigation of survey technologies and modeling techniques for improving deepwater surveys of nonindigenous species in the Northwestern Hawaiian Islands

Nonindigenous species (NIS) are a major threat to marine ecosystems, with possible dramatic effects on biodiversity, biological productivity, habitat structure and fisheries. The Papahānaumokuākea Marine National Monument (PMNM) has taken active steps to mitigate the threats of NIS in Northwestern Hawaiian Islands (NWHI). Of particular concern are the 13 NIS already detected in NWHI and two invasive species found among the main Hawaiian Islands, snowflake coral (Carijoa riseii) and a red alga (Hypnea musciformis). Much of the information regarding NIS in NWHI has been collected or informed by surveys using conventional SCUBA or fishing gear. These technologies have significant drawbacks. SCUBA is generally constrained to depths shallower than 40 m and several NIS of concern have been detected well below this limit (e.g., L. kasmira – 256 m) and fishing gear is highly selective. Consequently, not all habitats or species can be properly represented. Effective management of NIS requires knowledge of their spatial distribution and abundance over their entire range. Surveys which provide this requisite information can be expensive, especially in the marine environment and even more so in deepwater. Technologies which minimize costs, increase the probability of detection and are capable of satisfying multiple objectives simultaneously are desired. This report examines survey technologies, with a focus on towed camera systems (TCSs), and modeling techniques which can increase NIS detection and sampling efficiency in deepwater habitats of NWHI; thus filling a critical data gap in present datasets. A pilot study conducted in 2008 at French Frigate Shoals and Brooks Banks was used to investigate the application of TCSs for surveying NIS in habitats deeper than 40 m. Cost and data quality were assessed. Over 100 hours of video was collected, in which 124 sightings of NIS were made among benthic habitats from 20 to 250 m. Most sightings were of a single cosmopolitan species, Lutjanus kasmira, but Cephalopholis argus, and Lutjanus fulvus, were also detected. The data expand the spatial distributions of observed NIS into deepwater habitats, identify algal plain as an important habitat and complement existing data collected using SCUBA and fishing gear. The technology’s principal drawback was its inability to identify organisms of particular concern, such as Carijoa riseii and Hypnea musciformis due to inadequate camera resolution and inability to thoroughly inspect sites. To solve this issue we recommend incorporating high-resolution cameras into TCSs, or using alternative technologies, such as technical SCUBA diving or remotely operated vehicles, in place of TCSs. We compared several different survey technologies by cost and their ability to detect NIS and these results are summarized in Table 3.

Characterizing Jobos Bay, Puerto Rico: a watershed modeling analysis and monitoring plan

Land-based pollution is commonly identified as a major contributor to the observed deterioration of shallow-water coral reef ecosystem health. Human activity on the coastal landscape often induces nutrient enrichment, hypoxia, harmful algal blooms, toxic contamination and other stressors that have degraded the quality of coastal waters. Coral reef ecosystems throughout Puerto Rico, including Jobos Bay, are under threat from coastal land uses such as urban development, industry and agriculture. The objectives of this report were two-fold: 1. To identify potentially harmful land use activities to the benthic habitats of Jobos Bay, and 2. To describe a monitoring plan for Jobos Bay designed to assess the impacts of conservation practices implemented on the watershed. This characterization is a component of the partnership between the U.S. Department of Agriculture (USDA) and the National Oceanic and Atmospheric Administration (NOAA) established by the Conservation Effects Assessment Project (CEAP) in Jobos Bay. CEAP is a multi-agency effort to quantify the environmental benefits of conservation practices used by private landowners participating in USDA programs. The Jobos Bay watershed, located in southeastern Puerto Rico, was selected as the first tropical CEAP Special Emphasis Watershed (SEW). Both USDA and NOAA use their respective expertise in terrestrial and marine environments to model and monitor Jobos Bay resources. This report documents NOAA activities conducted in the first year of the three-year CEAP effort in Jobos Bay. Chapter 1 provides a brief overview of the project and background information on Jobos Bay and its watershed. Chapter 2 implements NOAA’s Summit to Sea approach to summarize the existing resource conditions on the watershed and in the estuary. Summit to Sea uses a GIS-based procedure that links patterns of land use in coastal watersheds to sediment and pollutant loading predictions at the interface between terrestrial and marine environments. The outcome of Summit to Sea analysis is an inventory of coastal land use and predicted pollution threats, consisting of spatial data and descriptive statistics, which allows for better management of coral reef ecosystems. Chapters 3 and 4 describe the monitoring plan to assess the ecological response to conservation practices established by USDA on the watershed. Jobos Bay is the second largest estuary in Puerto Rico, but has more than three times the shoreline of any other estuarine area on the island. It is a natural harbor protected from offshore wind and waves by a series of mangrove islands and the Punta Pozuelo peninsula. The Jobos Bay marine ecosystem includes 48 km² of mangrove, seagrass, coral reef and other habitat types that span both intertidal and subtidal areas. Mapping of Jobos Bay revealed 10 different benthic habitats of varying prevalence, and a large area of unknown bottom type covering 38% of the entire bay. Of the known benthic habitats, submerged aquatic vegetation, primarily seagrass, is the most common bottom type, covering slightly less than 30% of the bay. Mangroves are the dominant shoreline feature, while coral reefs comprise only 4% of the total benthic habitat. However, coral reefs are some of the most productive habitats found in Jobos Bay, and provide important habitat and nursery grounds for fish and invertebrates of commercial and recreational value.

Coupling ecology and economy: modeling optimal release scenarios for summer flounder (Paralichthys dentatus) stock enhancement

Increasing interest in the use of stock enhancement as a management tool necessitates a better understanding of the relative costs and benefits of alternative release strategies. We present a relatively simple model coupling ecology and economic costs to make inferences about optimal release scenarios for summer flounder (Paralichthys dentatus), a subject of stock enhancement interest in North Carolina. The model, parameterized from mark-recapture experiments, predicts optimal release scenarios from both survival and economic standpoints for varyious dates-of-release, sizes-at-release, and numbers of fish released. Although most stock enhancement efforts involve the release of relatively small fish, the model suggests that optimal results (maximum survival and minimum costs) will be obtained when relatively large fish (75–80 mm total length) are released early in the nursery season (April). We investigated the sensitivity of model predictions to violations of the assumption of density-independent mortality by including density-mortality relationships based on weak and strong type-2 and type-3 predator functional responses (resulting in depensatory mortality at elevated densities). Depending on postrelease density, density-mortality relationships included in the model considerably affect predicted postrelease survival and economic costs associated with enhancement efforts, but do not alter the release scenario (i.e. combination of release variables) that produces optimal results. Predicted (from model output) declines in flounder over time most closely match declines observed in replicate field sites when mortality in the model is density-independent or governed by a weak type-3 functional response. The model provides an example of a relatively easy-to-develop predictive tool with which to make inferences about the ecological and economic potential of stock enhancement of summer flounder and provides a template for model creation for additional species that are subjects of stock enhancement interest, but for which limited empirical data exist.

Modeling the effect of striped bass (Morone saxatilis) on the population viability of Sacramento River winter-run chinook salmon (Onchorhynchus tshawytscha)

We estimated the impact of striped bass (Morone saxatilis) predation on winter-run chinook salmon (Oncorhynchus tshawytscha) with a Bayesian population dynamics model using striped bass and winter-run chinook salmon population abundance data. Winter-run chinook salmon extinction and recovery probabilities under different future striped bass abundance levels were estimated by simulating from the posterior distribution of model parameters. The model predicts that if the striped bass population declines to 512,000 adults as expected in the absence of stocking, winter-run chinook salmon will have about a 28% chance of quasi-extinction (defined as three consecutive spawning runs of fewer than 200 adults) within 50 years. If stocking stabilizes the striped bass population at 700,000 adults, the predicted quasi-extinction probability is 30%. A more ambitious stocking program that maintains a population of 3 million adult striped bass would increase the predicted quasi-extinction probability to 55%. Extinction probability, but not recovery probability, was fairly insensitive to assumptions about density dependence. We conclude that winter-run chinook salmon face a serious extinction risk without augmentation of the striped bass population and that substantial increases in striped bass abundance could significantly increase the threat to winter-run chi-nook salmon if not mitigated by increasing winter chinook salmon survival in some other way.

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.