Assessing conservation status with extensive but low‐resolution data: Application of frequentist and Bayesian models to endangered Athabasca River rainbow trout

• Published in: Conservation Biology
• Main Authors: Post, John R., Ward, Hillary G. M., Wilson, Kyle L., Sterling, George L., Cantin, Ariane, Taylor, Eric B.
• Format: Article in Journal/Newspaper
• Language: English
• Published: Wiley 2021
• Subjects: Athabasca River; Canada
• Online Access: http://dx.doi.org/10.1111/cobi.13783; https://onlinelibrary.wiley.com/doi/pdf/10.1111/cobi.13783; https://onlinelibrary.wiley.com/doi/full-xml/10.1111/cobi.13783

Abstract Use of extensive but low‐resolution abundance data is common in the assessment of species at‐risk status based on quantitative decline criteria under International Union for Conservation of Nature (IUCN) and national endangered species legislation. Such data can be problematic for 3 reasons. First, statistical power to reject the null hypothesis of no change is often low because of small sample size and high sampling uncertainty leading to a high frequency of type II errors. Second, range‐wide assessments composed of multiple site‐specific observations do not effectively weight site‐specific trends into global trends. Third, uncertainty in site‐specific temporal trends and relative abundance are not propagated at the appropriate spatial scale. A common result is the propensity to underestimate the magnitude of declines and therefore fail to identify the appropriate at‐risk status for a species. We used 3 statistical approaches, from simple to more complex, to estimate temporal decline rates for a designatable unit (DU) of rainbow trout in the Athabasca River watershed in western Canada. This DU is considered a native species for purposes of listing because of its genetic composition characterized as >0.95 indigenous origin in the face of continuing introgressive hybridization with introduced populations in the watershed. Analysis of abundance trends from 57 time series with a fixed‐effects model identified 33 sites with negative trends, but only 2 were statistically significant. By contrast, a hierarchical linear mixed model weighted by site‐specific abundance provided a DU‐wide decline estimate of 16.4% per year and a 3‐generation decline of 93.2%. A hierarchical Bayesian mixed model yielded a similar 3‐generation decline trend of 91.3% and the posterior distribution showed that the estimate had a >99% probability of exceeding thresholds for an endangered listing. We conclude that the Bayesian approach was the most useful because it provided a probabilistic statement of threshold exceedance in ...