RESUMO
Density surface models (DSMs) are an important tool in the conservation and management of cetaceans. Most previous applications of DSMs have adopted a two-step approach to model fitting (hereafter referred to as the Two-Stage Method), whereby detection probabilities are first estimated using distance sampling detection functions and subsequently used as an offset when fitting a density-habitat model. Although variance propagation techniques have recently become available for the Two-Stage Method, most previous applications have not propagated detection probability uncertainty into final density estimates. In this paper, we describe an alternative approach for fitting DSMs based on Bayesian hierarchical inference (hereafter referred to as the Bayesian Method), which is a natural framework for simultaneously propagating multiple sources of uncertainty into final estimates. Our framework includes (1) a mark-recapture distance sampling observation model that can accommodate two team line transect data, (2) an informed prior for the probability a group of animals is at the surface and available for detection (i.e. surface availability) (3) a density-habitat model incorporating spatial smoothers and (4) a flexible compound Poisson-gamma model for count data that incorporates overdispersion and zero-inflation. We evaluate our method and compare its performance to the Two-Stage Method with simulations and an application to line transect data of fin whales (Balaenoptera physalus) off the east coast of the USA. Simulations showed that both methods had low bias (<1.5%) and confidence interval coverage close to the nominal 95% rate when variance was propagated from the first step. Results from the fin whale analysis showed that density estimates and predicted distribution patterns were largely similar among methods; however, the coefficient of variation of the final abundance estimate more than doubled (0.14 vs 0.31) when detection variance was correctly propagated into final estimates. An analysis of the variance components demonstrated that overall detectability as well as surface availability contributed substantial amounts of variance in the final abundance estimates whereas uncertainty in mean group size contributed a negligible amount. Our method provides a Bayesian alternative to DSMs that incorporates much of the flexibility available in the Two-Stage Method. In addition, these results demonstrate the degree to which uncertainty can be underestimated if certain components of a DSM are assumed fixed.
RESUMO
The objective of this study was to identify the main environmental covariates related to the abundance of 17 cetacean species/groups in the western North Atlantic Ocean based on generalized additive models, to establish a current habitat suitability baseline, and to estimate abundance that incorporates habitat characteristics. Habitat models were developed from dedicated sighting survey data collected by NOAA- Northeast and Southeast Fisheries Science Centers during July 2010 to August 2013. A group of 7 static physiographic characteristics and 9 dynamic environmental covariates were included in the models. For the small cetacean models, the explained deviance ranged from 16% to 69%. For the large whale models, the explained deviance ranged from 32% to 52.5%. Latitude, sea surface temperature, bottom temperature, primary productivity and distance to the coast were the most common covariates included and their individual contribution to the deviance explained ranged from 5.9% to 18.5%. The habitat-density models were used to produce seasonal average abundance estimates and habitat suitability maps that provided a good correspondence with observed sighting locations and historical sightings for each species in the study area. Thus, these models, maps and abundance estimates established a current habitat characterization of cetacean species in these waters and have the potential to be used to support management decisions and conservation measures in a marine spatial planning context.
