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.

Wave forecasts associated with hurricane landfalls in the vicinity of Cape Lookout, North Carolina

Hurricanes can cause extensive damage to the coastline and coastal communities due to wind-generated waves and storm surge. While extensive modeling efforts have been conducted regarding storm surge, there is far less information about the effects of waves on these communities and ecosystems as storms make landfall. This report describes a preliminary use of NCCOS’ WEMo (Wave Exposure Model; Fonseca and Malhotra 2010) to compute the wind wave exposure within an area of approximately 25 miles radius from Beaufort, North Carolina for estuarine waters encompassing Bogue Sound, Back Sound and Core Sound during three hurricane landfall scenarios. The wind wave heights and energy of a site was a computation based on wind speed, direction, fetch and local bathymetry. We used our local area (Beaufort, North Carolina) as a test bed for this product because it is frequently impacted by hurricanes and we had confidence in the bathymetry data. Our test bed conditions were based on two recent Hurricanes that strongly affected this area. First, we used hurricane Isabel which made landfall near Beaufort in September 2003. Two hurricane simulations were run first by passing hurricane Isabel along its actual path (east of Beaufort) and second by passing the same storm to the west of Beaufort to show the potential effect of the reversed wind field. We then simulated impacts by a hurricane (Ophelia) with a different landfall track, which occurred in September of 2005. The simulations produced a geographic description of wave heights revealing the changing wind and wave exposure of the region as a consequence of landfall location and storm intensity. This highly conservative simulation (water levels were that of low tide) revealed that many inhabited and developed shorelines would receive wind waves for prolonged periods of time at heights far above that found during even the top few percent of non-hurricane events. The simulations also provided a sense for how rapidly conditions could transition from moderate to highly threatening; wave heights were shown to far exceed normal conditions often long before the main body of the storm arrived and importantly, at many locations that could impede and endanger late-fleeing vessels seeking safe harbor. When joined with other factors, such as storm surge and event duration, we anticipate that the WEMo forecasting tool will have significant use by local emergency agencies and the public to anticipate the relative exposure of their property arising as a function of storm location and may also be used by resource managers to examine the effects of storms in a quantitative fashion on local living marine resources.

Estimation of bycatch in shrimp trawl fisheries: a comparison of estimation methods using field data and simulated data

Bycatch, or the incidental catch of nontarget organisms during fi shing operations, is a major issue in U.S. shrimp trawl fisheries. Because bycatch is typically discarded at sea, total bycatch is usually estimated by extrapolating from an observed bycatch sample to the entire fleet with either mean-per-unit or ratio estimators. Using both field observations of commercial shrimp trawlers and computer simulations, I compared five methods for generating bycatch estimates that were used in past studies, a mean-per-unit estimator and four forms of the ratio estimator, respectively: 1) the mean fish catch per unit of effort, where unit effort was a proxy for sample size, 2) the mean of the individual fish to shrimp ratios, 3) the ratio of mean fish catch to mean shrimp catch, 4) the mean of the ratios of fish catch per time fished (a variable measure of effort), and 5) the ratio of mean fish catch per mean time fished. For field data, different methods used to estimate bycatch of Atlantic croaker, spot, and weakfish yielded extremely different results, with no discernible pattern in the estimates by method, geographic region, or species. Simulated fishing fleets were used to compare bycatch estimated by the fi ve methods with “actual” (simulated) bycatch. Simulations were conducted by using both normal and delta lognormal distributions of fish and shrimp and employed a range of values for several parameters, including mean catches of fish and shrimp, variability in the catches of fish and shrimp, variability in fishing effort, number of observations, and correlations between fish and shrimp catches. Results indicated that only the mean per unit estimators provided statistically unbiased estimates, while all other methods overestimated bycatch. The mean of the individual fish to shrimp ratios, the method used in the South Atlantic Bight before the 1990s, gave the most biased estimates. Because of the statistically significant two- and 3-way interactions among parameters, it is unlikely that estimates generated by one method can be converted or corrected to estimates made by another method: therefore bycatch estimates obtained with different methods should not be compared directly.

An evaluation of back-calculation methodology using simulated otolith data

I simulated somatic growth and accompanying otolith growth using an individual-based bioenergetics model in order to examine the performance of several back-calculation methods. Four shapes of otolith radius-total length relations (OR-TL) were simulated. Ten different back-calculation equations, two different regression models of radius length, and two schemes of annulus selection were examined for a total of 20 different methods to estimate size at age from simulated data sets of length and annulus measurements. The accuracy of each of the twenty methods was evaluated by comparing the back-calculated length-at-age and the true length-at-age. The best back-calculation technique was directly related to how well the OR-TL model fitted. When the OR-TL was sigmoid shaped and all annuli were used, employing a least squares linear regression coupled with a log-transformed Lee back-calculation equation (y-intercept corrected) resulted in the least error; when only the last annulus was used, employing a direct proportionality back-calculation equation resulted in the least error. When the OR-TL was linear, employing a functional regression coupled with the Lee back-calculation equation resulted in the least error when all annuli were used, and also when only the last annulus was used. If the OR-TL was exponentially shaped, direct substitution into the fitted quadratic equation resulted in the least error when all annuli were used, and when only the last annulus was used. Finally, an asymptotically shaped OR-TL was best modeled by the individually corrected Weibull cumulative distribution function when all annuli were used, and when only the last annulus was used.

Tracing sediment transport and bed regime in Nha Trang Bay

Three spatial structure groups of radionuclides in U and Th series, 210Pb-excess and 137Cs, and 40K were found based on analyzing temporal and spatial datum of their content by factor analysis with oblique rotation in Nhatrang bay. U and Th spatial structure with their contours decreased toward the offshore, ran longshore and divided seawater of bay into two parts with strong gradient on both sides. Inside part located from center of Nhatrang bay toward the seashore with three main deposit centers of their contents higher than 23 Bq/kg.dry for 238U and 40 Bq/kg.dry for 232Th, indicated unstability of shoreline. Almost sediments coming from river extended toward the offshore, were stopped and transported toward southeastern. The outside part was less than above mentioned content. The boundary line between two parts superposed with the constantly limit line of turbid plume in the rainy season. Direct influence of the continental runoff was limited by the 9 Bq/kg.dry contour of 238U, 19 Bq/kg.dry contour of 232Th. Longshore current was a predominant process whereas lateral transport as sifting and winnowing process of finer grains in sediments of Nhatrang bay. Areas that had very low content of 137Cs and 210 Pb-excess adjoining shoreline showed areas being eroded. Accumulation of 137Cs and 210 Pbexcess nearby river mouth characterized for fine compositions of sediments controlled by seasonal plumes and sites further toward the south indicated finer materials transported from river and accumulated in lack of hydrodynamic process. Near shore accumulation of 40K revealed the sediments there originated from bed erosion.

Viability of Penaeus monodon eggs after simulated transport conditions

P. monodon spawners, transported from maturation pens suffer from stress which in turn may lead to lowered spawning rate or fertility. Spawning the females in the maturation site and transporting the eggs to the hatchery site is being considered as an alternative. Egg transport costs may be reduced to a minimum by using eggs from ablated spawners, transported at high density with no aeration. Experiments on higher egg densities as well as on transport of nauplii should, however, be undertaken.