Behavior of yellowfin (<Emphasis Type="Italic">Thunnus albacares</Emphasis>) and bigeye (<Emphasis Type="Italic">T. obesus</Emphasis>) tuna in a network of fish aggregating devices (FADs)

The influence of multiple anchored fish aggregating devices (FADs) on the spatial behavior of yellowfin (Thunnus albacares) and bigeye tuna (T. obesus) was investigated by equipping all thirteen FADs surrounding the island of Oahu (HI, USA) with automated sonic receivers (“listening stations”) and intra-peritoneally implanting individually coded acoustic transmitters in 45 yellowfin and 12 bigeye tuna. Thus, the FAD network became a multi-element passive observatory of the residence and movement characteristics of tuna within the array. Yellowfin tuna were detected within the FAD array for up to 150 days, while bigeye tuna were only observed up to a maximum of 10 days after tagging. Only eight yellowfin tuna (out of 45) and one bigeye tuna (out of 12) visited FADs other than their FAD of release. Those nine fish tended to visit nearest neighboring FADs and, in general, spent more time at their FAD of release than at the others. Fish visiting the same FAD several times or visiting other FADs tended to stay longer in the FAD network. A majority of tagged fish exhibited some synchronicity when departing the FADs but not all tagged fish departed a FAD at the same time: small groups of tagged fish left together while others remained. We hypothesize that tuna (at an individual or collective level) consider local conditions around any given FAD to be representative of the environment on a larger scale (e.g., the entire island) and when those conditions become unfavorable the tuna move to a completely different area. Thus, while the anchored FADs surrounding the island of Oahu might concentrate fish and make them more vulnerable to fishing, at a meso-scale they might not entrain fish longer than if there were no (or very few) FADs in the area. At the existing FAD density, the ‘island effect’ is more likely to be responsible for the general presence of fish around the island than the FADs. We recommend further investigation of this hypothesis.

---

Statistics in disease ecology: introduction to a special issue

The three papers included in this special issue represent a set of presentations in an invited session on disease ecology at the 2005 Spring Meeting of the Eastern North American Region of the International Biometric Society. The papers each address statistical estimation and inference for particular components of different disease processes and, taken together, illustrate the breadth of statistical issues arising in the study of the ecology and public health impact of disease. As an introduction, we provide a very brief overview of the area of “disease ecology”, a variety of synonyms addressing different aspects of disease ecology, and present a schematic structure illustrating general components of the underlying disease process, data collection issues, and different disciplinary perspectives ranging from microbiology to public health surveillance.

---

Testing homogeneity of capture probabilities of animals in a removal experiment

Consider the removal experiment used to estimate population sizes. Statistical methods towards testing the homogeneity of capture probabilities of animals, including a graphical diagnostic and a formal test, are presented and illustrated by real biological examples. Simulation is used to assess the test and compare it with the χ² test.

Trap mortality in mark-recapture studies

When animals die in traps in a mark-recapture study, straightforward likelihood inferences are possible in a class of models. The class includes M₀, M_t, and M_b as reported by White et al. (Los Alamos National Laboratory, LA-8787-NERP, pp 235, 1982), those that do not involve heterogeneity. We include three Markov chain “persistence” models and show that they provide good fits in a trapping study of deer mice in the Cascade-Siskiyou National Monument of Southern Oregon where trapping mortality was high.

Bootstrap methods for simultaneous benchmark analysis with quantal response data

A primary objective in quantitative risk assessment is the characterization of risk which is defined to be the likelihood of an adverse effect caused by an environmental toxin or chemcial agent. In modern risk-benchmark analysis, attention centers on the “benchmark dose” at which a fixed benchmark level of risk is achieved, with a lower confidence limits on this dose being of primary interest. In practice, a range of benchmark risks may be under study, so that the individual lower confidence limits on benchmark dose must be corrected for simultaneity in order to maintain a specified overall level of confidence. For the case of quantal data, simultaneous methods have been constructed that appeal to the large sample normality of parameter estimates. The suitability of these methods for use with small sample sizes will be considered. A new bootstrap technique is proposed as an alternative to the large sample methodology. This technique is evaluated via a simulation study and examples from environmental toxicology.

Replace the NOAEL and LOAEL with the BMDL01 and BMDL10

Although benchmark-dose methodology has existed for more than 20 years, benchmark doses (BMDs) still have not fully supplanted the no-observed-adverse-effect level (NOAEL) and lowest-observed-adverse-effect level (LOAEL) as points of departure from the experimental dose–response range for setting acceptable exposure levels of toxic substances. Among the issues involved in replacing the NOAEL (LOAEL) with a BMD are (1) which added risk level(s) above background risk should be targeted as benchmark responses (BMRs), (2) whether to apply the BMD methodology to both carcinogenic and noncarcinogenic toxic effects, and (3) how to model continuous health effects that aren’t observed in a natural risk-based context like dichotomous health effects. This paper addresses these issues and recommends specific BMDs to replace the NOAEL and LOAEL.

Confidence limits on one-stage model parameters in benchmark risk assessment

In modern environmental risk analysis, inferences are often desired on those low dose levels at which a fixed benchmark risk is achieved. In this paper, we study the use of confidence limits on parameters from a simple one-stage model of risk historically popular in benchmark analysis with quantal data. Based on these confidence bounds, we present methods for deriving upper confidence limits on extra risk and lower bounds on the benchmark dose. The methods are seen to extend automatically to the case where simultaneous inferences are desired at multiple doses. Monte Carlo evaluations explore characteristics of the parameter estimates and the confidence limits under this setting.

Multi-parameter renewal processes and simulation of forest fires

The concept of the renewal property is extended to processes indexed by a multidimensional time parameter. The definition given includes not only partial sum processes, but also Poisson processes and many other point processes whose jump points are not totally ordered. Various properties of renewal processes are discussed. Renewal processes are proposed as a basis for modelling the spread of a forest fire under a prevailing wind.

A note on data augmentation and iterative simulation for survival rate estimation

In this paper we examine the use of data augmentation techniques for simplifying iterative simulation in the context of both Bayesian and classical statistical inference for survival rate estimation. We examine two distinct model families common in population ecology to illustrate our ideas, ring-recovery models and capture–recapture models, and we present the computational advantage of this approach. We discuss also the fact that problems associated with identifiability in the classical framework can be overcome using data augmentation, but highlight the dangers in doing so under both inferential paradigms.

Hierarchical spatial point process analysis for a plant community with high biodiversity

A complex multivariate spatial point pattern of a plant community with high biodiversity is modelled using a hierarchical multivariate point process model. In the model, interactions between plants with different post-fire regeneration strategies are of key interest. We consider initially a maximum likelihood approach to inference where problems arise due to unknown interaction radii for the plants. We next demonstrate that a Bayesian approach provides a flexible framework for incorporating prior information concerning the interaction radii. From an ecological perspective, we are able both to confirm existing knowledge on species’ interactions and to generate new biological questions and hypotheses on species’ interactions.

Distribution and abundance of penaeid shrimps in a hypersaline lagoon in northwestern Mexico,emphasizing the brown shrimp <Emphasis Type="Italic">Farfantepenaeus californiensis</Emphasis> life cycle

Life histories of penaeid shrimp have been classified according to the preferred habitats of postlarval, juvenile, and adult stages, ranging from exclusively estuarine to exclusively offshore waters. Brown shrimp Farfantepenaeus californiensis migrate to an offshore habitat at the juvenile stage or even a smaller body size. This paper presents results of monthly samplings from 24 stations over 1 year in the Agiabampo Lagoon complex, a hypersaline lagoon in northwestern Mexico. Five species of penaeid shrimp were identified, with brown shrimp the most abundant during the year of sampling. Results suggest that residency of brown shrimp inside this lagoon is longer than reported in previous studies. An interaction between length and environmental variables (near-surface temperature, salinity, and rainfall) appear to be cues concerning migration.

---

Fiat boundaries: some implications for interpretation,decision-support,and multi-temporal analysis

Polygon-based thematic maps can be composed of boundaries that exist by definition—i.e., bona fide boundaries—or those that exist relative to a specific interpretation of a spatial phenomenon—i.e., fiat boundaries. The construction of maps composed of fiat boundaries is usually based on a subjective interpretive methodology that is affected by the data used to construct the map and the minimum mapping unit employed. That fiat boundaries are not the same as bona fide boundaries affects their use in computer-based spatial decision support tools. This is discussed both in terms of an analysis conducted at one specific moment, and in respect to increasingly common multi-temporal analysis.

A stochastic forest fire growth model

We consider a stochastic fire growth model, with the aim of predicting the behaviour of large forest fires. Such a model can describe not only average growth, but also the variability of the growth. Implementing such a model in a computing environment allows one to obtain probability contour plots, burn size distributions, and distributions of time to specified events. Such a model also allows the incorporation of a stochastic spotting mechanism.

Social network analysis and valid Markov chain Monte Carlo tests of null models

Analyses of animal social networks derived from group-based associations often rely on randomisation methods developed in ecology (Manly, Ecology 76:1109–1115, 1995) and made available to the animal behaviour community through implementation of a pair-wise swapping algorithm by Bejder et al. (Anim Behav 56:719–725, 1998). We report a correctable flaw in this method and point the reader to a wider literature on the subject of null models in the ecology literature. We illustrate the importance of correcting the method using a toy network and use it to make a preliminary analysis of a network of associations among eagle rays.

Macrozoobenthos diversity in an oxygen minimum zone off northern Namibia

A benthological survey in the Benguela upwelling area off northern Namibia (located at 17.3°S and water depth ranging between 26 and 117 m) showed the concentration of dissolved oxygen and the accumulation of organic-rich sediments to control macrozoobenthic community patterns. In contrast to highly biodiverse nearshore areas with well-structured shell deposits of the brachiopod Discinisca tenuis (Sowerby 1847), the benthic community in the oxygen minimum zone (OMZ) decreased strongly in species numbers. Nevertheless, a well-established community ranging from 13 to 31 species persisted. Species densities (300–3,350 ind m⁻²) and biomass (4–109 g afdw/m²) were surprisingly high for areas with near bottom oxygen concentrations from 0.06 to 0.88 ml l⁻¹. In contrast to OMZ’s of other upwelling areas, where the benthic macrofauna is generally dominated by small-bodied polychaetes, off Namibia larger key organisms like the bivalve Nuculana bicuspidata (Gould 1845) and the snail Nassarius vinctus (Marrett 1877) accounted for a large proportion of the macrozoobenthos >1 mm. This is supposed to have a distinct effect on the functional properties of the sediments.

---

Flexible hierarchical mark-recapture modeling for open populations using WinBUGS

Hierarchical mark-recapture models offer three advantages over classical mark-recapture models: (i) they allow expression of complicated models in terms of simple components; (ii) they provide a convenient way of modeling missing data and latent variables in a way that allows expression of relationships involving latent variables in the model; (iii) they provide a convenient way of introducing parsimony into models involving many nuisance parameters. Expressing models using the complete data likelihood we show how many of the standard mark-recapture models for open populations can be readily fitted using the software WinBUGS. We include examples that illustrate fitting the Cormack–Jolly–Seber model, multi-state and multi-event models, models including auxiliary data, and models including density dependence.

A note on non-parametric and semi-parametric modeling of wildfire hazard in Los Angeles County, California

This paper explores the use of, and problems that arise in, kernel smoothing and parametric estimation of the relationships between wildfire incidence and various meteorological variables. Such relationships may be treated as components in separable point process models for wildfire activity. The resulting models can be used for comparative purposes in order to assess the predictive performance of the Burning Index.

The swimming crab Charybdis smithii: distribution, biology and trophic role in the pelagic ecosystem of the western Indian Ocean

Surface “swarms” of the swimming crabs Charybdis smithii are still considered as an unusual phenomenon in the open Indian Ocean, although their dense pelagic aggregations were already reported in waters off the Indian coast and in the northern Arabian Sea. Based on an extensive large-scale data series taken over 45 years, we demonstrate that C. smithii is common in the pelagic provinces of the western Indian Ocean driven by the wind monsoon regime. Swimming crabs are dispersed by the monsoon currents throughout the equatorial Indian Ocean. They aggregate at night in the upper 150-m layer, where their estimated biomass derived from pelagic trawling data can exceed 130 kg km⁻². Abundance of C. smithii can reach >15,000 ind. km⁻² in July (i.e. the peak of the south-west monsoon), declines by 50-fold in March and is negligible in May. C. smithii is an important prey for more than 30 species of abundant epipelagic top predators. In turn, it feeds on mesopelagic species. This swimming crab is a major species of the intermediate trophic levels and represents a crucial seasonal trophic link in the open ocean ecosystem of the western Indian Ocean. Outbursts in pelagic waters of huge biomasses of ordinarily benthic crustaceans (C. smithii and Natosquilla investigatoris) are a remarkable feature of the Indian Ocean, although similar, but smaller, events are reported in the Pacific and Atlantic Oceans.

---

Vertical movement and habitat of opah (Lampris guttatus) in the central North Pacific recorded with pop-up archival tags

Data from 11 pop-up archival transmitting tags attached to opah (Lampris guttatus, F. Lampridae) in the central North Pacific between November 2003 and March 2005 were used to describe their vertical movement and habitat. In the subtropical gyre northwest of the Hawaiian Islands, opah generally inhabited a 50–400 m depth range and 8–22°C temperatures. They were frequently found in depths of 50–150 m at night and in greater depths (100–400 m) during the day, but were constantly moving vertically within this broad range. At night, excursions below 200 m were not uncommon and during the day they were very likely to spend some time at depths <175 m. Their vertical speeds were generally <25 cm s⁻¹ but on one occasion an opah descended at a burst speed of 4 m s⁻¹. Vertical habitat use by individual opah apparently varied with local oceanographic conditions, but over a 24-h period the average temperature experienced was always in the narrow range of 14.7 to 16.5°C.

---

Comparing model averaging with other model selection strategies for benchmark dose estimation

Model averaging (MA) has been proposed as a method of accommodating model uncertainty when estimating risk. Although the use of MA is inherently appealing, little is known about its performance using general modeling conditions. We investigate the use of MA for estimating excess risk using a Monte Carlo simulation. Dichotomous response data are simulated under various assumed underlying dose–response curves, and nine dose–response models (from the USEPA Benchmark dose model suite) are fit to obtain both model specific and MA risk estimates. The benchmark dose estimates (BMDs) from the MA method, as well as estimates from other commonly selected models, e.g., best fitting model or the model resulting in the smallest BMD, are compared to the true benchmark dose value to better understand both bias and coverage behavior in the estimation procedure. The MA method has a small bias when estimating the BMD that is similar to the bias of BMD estimates derived from the assumed model. Further, when a broader range of models are included in the family of models considered in the MA process, the lower bound estimate provided coverage close to the nominal level, which is superior to the other strategies considered. This approach provides an alternative method for risk managers to estimate risk while incorporating model uncertainty.