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.

Growth, mortality, and hatchdate distributions of larval and juvenile spotted seatrout (Cynoscion nebulosus) in Florida Bay, Everglades National Park

Life history aspects of larval and, mainly, juvenile spotted seatrout (Cynoscion nebulosus) were studied in Florida Bay, Everglades National Park, Florida. Collections were made in 1994−97, although the majority of juveniles were collected in 1995. The main objective was to obtain life history data to eventually develop a spatially explicit model and provide baseline data to understand how Everglades restoration plans (i.e. increased freshwater flows) could influence spotted seatrout vital rates. Growth of larvae and juveniles (<80 mm SL) was best described by the equation loge standard length = –1.31 + 1.2162 (loge age). Growth in length of juveniles (12–80 mm SL) was best described by the equation standard length = –7.50 + 0.8417 (age). Growth in wet weight of juveniles (15–69 mm SL) was best described by the equation loge wet-weight = –4.44 + 0.0748 (age). There were no significant differences in juvenile growth in length of spotted seatrout in 1995 between three geographical subdivisions of Florida Bay: central, western, and waters adjacent to the Gulf of Mexico. We found a significant difference in wet-weight for one of six cohorts categorized by month of hatchdate in 1995, and a significant difference in length for another cohort. Juveniles (i.e. survivors) used to calculate weekly hatchdate distributions during 1995 had estimated spawning times that were cyclical and protracted, and there was no correlation between spawning and moon phase. Temperature influenced otolith increment widths during certain growth periods in 1995. There was no evidence of a relationship between otolith growth rate and temperature for the first 21 increments. For increments 22–60, otolith growth rates decreased with increasing age and the extent of the decrease depended strongly in a quadratic fashion on the temperature to which the fish was exposed. For temperatures at the lower and higher range, increment growth rates were highest. We suggest that this quadratic relationship might be influenced by an environmental factor other than temperature. There was insufficient information to obtain reliable inferences on the relationship of increment growth rate to salinity.

Recruitment as an evolving random process of aggregation and mortality

The dynamics of the survival of recruiting fish are analyzed as evolving random processes of aggregation and mortality. The analyses draw on recent advances in the physics of complex networks and, in particular, the scale-free degree distribution arising from growing random networks with preferential attachment of links to nodes. In this study simulations were conducted in which recruiting fish 1) were subjected to mortality by using alternative mortality encounter models and 2) aggregated according to random encounters (two schools randomly encountering one another join into a single school) or preferential attachment (the probability of a successful aggregation of two schools is proportional to the school sizes). The simulations started from either a “disaggregated” (all schools comprised a single fish) or an aggregated initial condition. Results showed the transition of the school-size distribution with preferential attachment evolving toward a scale-free school size distribution, whereas random attachment evolved toward an exponential distribution. Preferential attachment strategies performed better than random attachment strategies in terms of recruitment survival at time when mortality encounters were weighted toward schools rather than to individual fish. Mathematical models were developed whose solutions (either analytic or numerical) mimicked the simulation results. The resulting models included both Beverton-Holt and Ricker-like recruitment, which predict recruitment as a function of initial mean school size as well as initial stock size. Results suggest that school-size distributions during recruitment may provide information on recruitment processes. The models also provide a template for expanding both theoretical and empirical recruitment research.

Capture probability of a survey trawl for red king crab (Paralithodes camtschaticus)

The relative abundance of Bristol Bay red king crab (Paralithodes camtschaticus) is estimated each year for stock assessment by using catch-per-swept-area data collected on the Alaska Fisheries Science Center’s annual eastern Bering Sea bottom trawl survey. To estimate survey trawl capture efficiency for red king crab, an experiment was conducted with an auxiliary net (fitted with its own heavy chain-link footrope) that was attached beneath the trawl to capture crabs escaping under the survey trawl footrope. Capture probability was then estimated by fitting a model to the proportion of crabs captured and crab size data. For males, mean capture probability was 72% at 95 mm (carapace length), the size at which full vulnerability to the survey trawl is assigned in the current management model; 84.1% at 135 mm, the legal size for the fishery; and 93% at 184 mm, the maximum size observed in this study. For females, mean capture probability was 70% at 90 mm, the size at which full vulnerability to the survey trawl is assigned in the current management model, and 77% at 162 mm, the maximum size observed in this study. The precision of our estimates for each sex decreased for juveniles under 60 mm and for the largest crab because of small sample sizes. In situ data collected from trawl-mounted video cameras were used to determine the importance of various factors associated with the capture of individual crabs. Capture probability was significantly higher when a crab was standing when struck by the footrope, rather than crouching, and higher when a crab was hit along its body axis, rather than from the side. Capture probability also increased as a function of increasing crab size but decreased with increasing footrope distance from the bottom and when artificial light was provided for the video camera.

Estimates of survival probabilities for oceanic-stage loggerhead sea turtles (Caretta caretta) in the North Atlantic

Estimates of instantaneous mortality rates (Z) and annual apparent survival probabilities (Φ) were generated from catch-curve analyses for oceanic-stage juvenile loggerheads (Caretta caretta) in the waters of the Azores. Two age distributions were analyzed: the “total sample” of 1600 loggerheads primarily captured by sighting and dipnetting from a variety of vessels in the Azores between 1984 and 1995 and the “tuna sample” of 733 loggerheads (a subset of the total sample) captured by sighting and dipnetting from vessels in the commercial tuna fleet in the Azores between 1990 and 1992. Because loggerhead sea turtles begin to emigrate from oceanic to neritic habitats at age 7, the best estimates of instantaneous mortality rate (0.094) and annual survival probability (0.911) not confounded with permanent emigration were generated for age classes 2 through 6. These estimates must be interpreted with caution because of the assumptions upon which catch-curve analyses are based. However, these are the first directly derived estimates of mortality and survival probabilities for oceanic-stage sea turtles. Estimation of survival probabilities was identified as “an immediate and critical requirement” in 2000 by the Turtle Expert Working Group of the U.S. National Marine Fisheries Service.

Assessing the precision of frequency distributions estimated from trawl-survey samples

In trawl surveys a cluster of fish are caught at each station, and fish caught together tend to have more similar characteristics, such as length, age, stomach contents etc., than those in the entire population. When this is the case, the effective sample size for estimates of the frequency distribution of a population characteristic can, therefore, be much smaller than the number of fish sampled during a survey. As examples, it is shown that the effective sample size for estimates of length-frequency distributions generated by trawl surveys conducted in the Barents Sea, off Namibia, and off South Africa is on average approximately one fish per tow. Thus many more fish than necessary are measured at each station (location). One way to increase the effective sample size for these surveys and, hence, increase the precision of the length-frequency estimates, is to reduce tow duration and use the time saved to collect samples at more stations.

Stock-rebuilding time isopleths and constant-F stock-rebuilding plans for overfished stocks

Stock-rebuilding time isopleths relate constant levels of fishing mortality (F), stock biomass, and management goals to rebuilding times for overfished stocks. We used simulation models with uncertainty about FMSY and variability in annual intrinsic growth rates (ry) to calculate rebuilding time isopleths for Georges Bank yellowtail flounder, Limanda ferruginea, and cowcod rockfish, Sebastes levis, in the Southern California Bight. Stock-rebuilding time distributions from stochastic models were variable and right-skewed, indicating that rebuilding may take less or substantially more time than expected. The probability of long rebuilding times increased with lower biomass, higher F, uncertainty about FMSY, and autocorrelation in ry values. Uncertainty about FMSY had the greatest effect on rebuilding times. Median recovery times from simulations were insensitive to model assumptions about uncertainty and variability, suggesting that median recovery times should be considered in rebuilding plans. Isopleths calculated in previous studies by deterministic models approximate median, rather than mean, rebuilding times. Stochastic models allow managers to specify and evaluate the risk (measured as a probability) of not achieving a rebuilding goal according to schedule. Rebuilding time isopleths can be used for stocks with a range of life histories and can be based on any type of population dynamics model. They are directly applicable with constant F rebuilding plans but are also useful in other cases. We used new algorithms for simulating autocorrelated process errors from a gamma distribution and evaluated sensitivity to statistical distributions assumed for ry. Uncertainty about current biomass and fishing mortality rates can be considered with rebuilding time isopleths in evaluating and designing constant-F rebuilding plans.