Report of the special session of "Chesapeake Bay Ecological Forecasting: Moving Ecosystem Modeling from Research to Operation" of the 2010 Chesapeake Modeling Symposium

Moving ecosystem modeling from research to applications and operations has direct management relevance and will be integral to achieving the water quality and living resource goals of the 2010 Chesapeake Bay Executive Order. Yet despite decades of ecosystem modeling efforts of linking climate to water quality, plankton and fish, ecological models are rarely taken to the operational phase. In an effort to promote operational ecosystem modeling and ecological forecasting in Chesapeake Bay, a meeting was convened on this topic at the 2010 Chesapeake Modeling Symposium (May, 10-11). These presentations show that tremendous progress has been made over the last five years toward the development of operational ecological forecasting models, and that efforts in Chesapeake Bay are leading the way nationally. Ecological forecasts predict the impacts of chemical, biological, and physical changes on ecosystems, ecosystem components, and people. They have great potential to educate and inform not only ecosystem management, but also the outlook and opinion of the general public, for whom we manage coastal ecosystems. In the context of the Chesapeake Bay Executive Order, ecological forecasting can be used to identify favorable restoration sites, predict which sites and species will be viable under various climate scenarios, and predict the impact of a restoration project on water quality.

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.)

Modeling the Chesapeake Bay and Tributaries: a synopsis

The last decade has seen the development and application of a spectrum of physical and numerical hydrographic models of the Chesapeake Bay and its tributaries. The success of the James River Hydraulic Model has initiated the construction of an estuarine hydraulic model of the entire Chesapeake System. Numerical analogues for hydrographic behavior and contaminant dispersion in one-, two-, and three dimensional model estuaries exist for various regions of the Bay. From an engineering viewpoint, one dimensional models are sufficiently advanced to be routinely employed in aiding management decisions. Bay investigators are playing leading roles in the development of two- and three-dimensional models of estuarine flows.

The Empirical Modeling of an Ecosystem

The authors have endeavored to create a verified a-posteriori model of a planktonic ecosystem. Verification of an empirically derived set of first-order, quadratic differential equations proved elusive due to the sensitivity of the model system to changes in initial conditions. Efforts to verify a similarly derived set of linear differential equations were more encouraging, yielding reasonable behavior for half of the ten ecosystem compartments modeled. The well-behaved species models gave indications as to the rate-controlling processes in the ecosystem.

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.

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)

The ecological basis of fishery yield of the Puerto Rico-Virgin Islands Insular Shelf: 1987 Assessment

A literature review was conducted to locate information on the flow of energy from primary producers to the fishery stocks of the Puerto Rican-Virgin Islands insular shelf. This report uses site-specific information to describe the major ecological subsystems, or habitats, of the region, to identify the more common species and the subsystems in which they occur, to quantify productivity and biomass, and to outline trophic relationships. Discussions on each topic and subsystem vary in substance and detail, being limited by the availability and accessibility of information. (PDF contains 189 pages) Seven distinct subsystems are described: mangrove estuary, seagrass bed, coral reef, algal plain, sand/mud bottom, shelf break, and overlying pelagic. Over 50 tables provide lists of species found in each habitat on various surveys dating back to 1956. Estimates of density, relative abundance, and productivity are provided when possible. We evaluated whether sufficient information exists to support an analysis of the energy basis of fishery production in the area, beginning with the design and development of an ecosystem model. Data needs in three categories - species lists, biomass, and trophic relations - were examined for each subsystem and for each of three species groups - primary producers, invertebrates, and fish. We concluded that adequate data, sufficient for modeling purposes, are available in 16 (25%) of 64 categories; limited data, those requiring greater extrapolation, are available in 35 (55%) categories; and no data are available in 13 (20%) categories. The best-studied subsystems are seagrass beds and coral reefs, with at least limited data in all categories. Invertebrates, the intermediate link in the food web between primary producers and fishes, are the least quantified group in the region. Primary production and fishes, however, are relatively well-studied, providing sufficient data to support an ecosystem-level analysis and to initiate a modeling effort.

Stormwater runoff - modeling impacts of urbanization and climate change

Development pressure throughout the coastal areas of the United States continues to build, particularly in the southeast (Allen and Lu 2003, Crossett et al. 2004). It is well known that development alters watershed hydrology: as land becomes covered with surfaces impervious to rain, water is redirected from groundwater recharge and evapotranspiration to stormwater runoff, and as the area of impervious cover increases, so does the volume and rate of runoff (Schueler 1994, Corbett et al. 1997). Pollutants accumulate on impervious surfaces, and the increased runoff with urbanization is a leading cause of nonpoint source pollution (USEPA 2002). Sediment, chemicals, bacteria, viruses, and other pollutants are carried into receiving water bodies, resulting in degraded water quality (Holland et al. 2004, Sanger et al. 2008). (PDF contains 5 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)

Watershed monitoring and the TMDL modeling to assess bacterial loading in estuarine environments to improve management

Shellfish bed closures along the North Carolina coast have increased over the years seemingly concurrent with increases in population (Mallin 2000). More and faster flowing storm water has come to mean more bacteria, and fecal indicator bacterial (FIB) standards for shellfish harvesting are often exceeded when no source of contamination is readily apparent (Kator and Rhodes, 1994). Could management reduce bacterial loads if the source of the bacteria where known? Several potentially useful methods for differentiating human versus animal pollution sources have emerged including Ribotyping and Multiple Antibiotic Resistance (MAR) (US EPA, 2005). Total Maximum Daily Load (TMDL) studies on bacterial sources have been conducted for streams in NC mountain and Piedmont areas (U.S. EPA, 1991 and 2005) and are likely to be mandated for coastal waters. TMDL analysis estimates allowable pollutant loads and allocates them to known sources so management actions may be taken to restore water to its intended uses (U.S. EPA, 1991 and 2005). This project sought first to quantify and compare fecal contamination levels for three different types of land use on the coast, and second, to apply MAR and ribotyping techniques and assess their effectiveness for indentifying bacterial sources. Third, results from these studies would be applied to one watershed to develop a case study coastal TMDL. All three watershed study areas are within Carteret County, North Carolina. Jumping Run Creek and Pettiford Creek are within the White Oak River Basin management unit whereas the South River falls within the Neuse River Basin. Jumping Run Creek watershed encompasses approximately 320 ha. Its watershed was a dense, coastal pocosin on sandy, relic dune ridges, but current land uses are primarily medium density residential. Pettiford Creek is in the Croatan National Forest, is 1133 ha. and is basically undeveloped. The third study area is on Open Grounds Farm in the South River watershed. Half of the 630 ha. watershed is under cultivation with most under active water control (flashboard risers). The remaining portion is forested silviculture.(PDF contains 4 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.

Age validation, growth, mortality, and demographic modeling of spotted gully shark (Triakis megalopterus) from the southeast coast of South Africa

This study documents validation of vertebral band-pair formation in spotted gully shark (Triakis megalopterus) with the use of fluorochrome injection and tagging of captive and wild sharks over a 21-year period. Growth and mortality rates of T. megalopterus were also estimated and a demographic analysis of the species was conducted. Of the 23 OTC (oxytetracycline) -marked vertebrae examined (12 from captive and 11 from wild sharks), seven vertebrae (three from captive and four from wild sharks) exhibited chelation of the OTC and fluoresced under ultraviolet light. It was concluded that a single opaque and translucent band pair was deposited annually up to at least 25 years of age, the maximum age recorded. Reader precision was assessed by using an index of average percent error calculated at 5%. No significant differences were found between male and female growth patterns (P>0.05), and von Bertalanffy growth model parameters for combined sexes were estimated to be L∞=1711.07 mm TL, k=0.11/yr and t0=–2.43 yr (n=86). Natural mortality was estimated at 0.17/yr. Age at maturity was estimated at 11 years for males and 15 years for females. Results of the demographic analysis showed that the population, in the absence of fishing mortality, was stable and not significantly different from zero and particularly sensitive to overfishing. At the current age at first capture and natural mortality rate, the fishing mortality rate required to result in negative population growth was low at F>0.004/ yr. Elasticity analysis revealed that juvenile survival was the principal factor in explaining variability in population growth rate.

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.

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

Ecological modeling for all

A brief review of the status of the ECOPATH modeling approach and software is presented, with emphasis on the recent release of a Windows version (ECOPATH 3.0), which enables consideration of uncertainties, and sets the stage for simulation modeling using ECOSIM. Modeling of coral reefs is emphasized.