Development of a species distribution model for fin whales (Balaenoptera physalus) within a Bayesian hierarchical framework: Implications for uncertainty

Integrated Statistics, Woods Hole, MA, USA
National Marine Mammal Laboratory, NOAA, National Marine Fisheries Service, Alaska Fisheries Science Center, Seattle, US
NOAA Southeast Fisheries Science Center, Miami, FL, US
NOAA Northeast Fisheries Science Center, Woods Hole, MA, US
DOI
10.7287/peerj.preprints.27424v1
Subject Areas
Conservation Biology, Ecology, Marine Biology, Natural Resource Management, Population Biology
Keywords
: Species Distribution Model (SDM), Bayesian Model, jagam, Generalized Additive Model (GAM), Tweedie distribution, Fin whales
Copyright
© 2018 Sigourney et al.
Licence
This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ Preprints) and either DOI or URL of the article must be cited.
Cite this article
Sigourney DB, Chavez-Rosales S, Conn P, Garrison L, Josephson E, Palka D. 2018. Development of a species distribution model for fin whales (Balaenoptera physalus) within a Bayesian hierarchical framework: Implications for uncertainty. PeerJ Preprints 6:e27424v1

Abstract

Species distribution models (SDMs) have proven to be an integral tool in the conservation and management of cetaceans. Many applications have adopted a two-step approach where a detection function is estimated using conventional distance sampling in the first step and subsequently used as an offset to a density-habitat model in the second step. A drawback to this approach, hereafter referred to as the conventional species distribution model (CSDM), is the difficulty in propagating the uncertainty from the first step to the final density estimates. We describe a Bayesian hierarchical species distribution model (BHSDM) which has the advantage of simultaneously propagating multiple sources of uncertainty. Our framework includes 1) a mark-recapture distance sampling observation model that can accommodate two team line transect data, 2) an informed prior for surface availability 3) spatial smoothers using spline-like bases and 4) a compound Poisson-gamma likelihood which is a special case of the Tweedie distribution. We compare our approach to the CSDM method using a simulation study and a case study of fin whales (Balaenoptera physalus) off the East Coast of the USA. Simulations showed that the BHSDM method produced estimates with lower precision but with confidence interval coverage closer to the nominal 95% rate (94% for the BSHDM vs 85% for the CSDM). Results from the fin whale analysis showed that density estimates and predicted distribution patterns were largely similar among methods. Abundance estimates were also similar though modestly higher for the CSDM (4700, CV=0.13) than the BHSDM (4526, CV=0.26). Estimated sampling error differed substantially among the two methods where the average CV for density estimates from BHSDM method was approximately 3.5 times greater than estimates from the CSDM method. Successful wildlife management hinges on the ability to properly quantify uncertainty. Underestimates of uncertainty can result in ill-informed management decisions. Our results highlight the additional sampling uncertainty that is propagated in a hierarchical framework. Future applications of SDMs should consider techniques that allow all sources of error to be fully represented in final density predictions.

Author Comment

This is a submission to PeerJ for review.

Supplemental Information

Additional information on the detection functions

DOI: 10.7287/peerj.preprints.27424v1/supp-1

Additional information on the details of the simulation study

DOI: 10.7287/peerj.preprints.27424v1/supp-2

Environmental Covariates

Description of environmental covariates included in the habitat models and transformations applied to normalize variables.

DOI: 10.7287/peerj.preprints.27424v1/supp-3

Description of the final MRDS models

Description of the final MRDS model used to model perception bias for each platform used in the AMAPPS surveys, species pooled in the analysis and sample size. The species pooled include humpback whales (huwh), fin whales (fiwh), minke whales (miwh), sei whales (sewh), unidentified fin-sei whale (fise), unidentified baleen whales (unbw), right whales (riwh) and sperm whales (spwh). The platform specific distance sampling (DS) model, truncation distance (W), key function (half-normal (HN) or hazard rate (HZ)), and total sample size (NTotal) for each analysis are included/shown/provided.

DOI: 10.7287/peerj.preprints.27424v1/supp-4

Code to run BHSDM model in JAGS

DOI: 10.7287/peerj.preprints.27424v1/supp-5

Raw sightings data of fin whales

DOI: 10.7287/peerj.preprints.27424v1/supp-6

Sightings and environmental covariate information for observed grid cells

DOI: 10.7287/peerj.preprints.27424v1/supp-7