A spatial approach to combatting wildlife crime

Poaching can have devastating impacts on animal and plant numbers, and in many countries has reached crisis levels, with illegal hunters employing increasingly sophisticated techniques. We used data from an 8‐year study in Savé Valley Conservancy, Zimbabwe, to show how geographic profiling—a mathematical technique originally developed in criminology and recently applied to animal foraging and epidemiology—can be adapted for use in investigations of wildlife crime. The data set contained information on over 10,000 incidents of illegal hunting and the deaths of 6,454 wild animals. We used a subset of data for which the illegal hunters’ identities were known. Our model identified the illegal hunters’ home villages based on the spatial locations of the hunting incidences (e.g., snares). Identification of the villages was improved by manipulating the probability surface inside the conservancy to reflect the fact that although the illegal hunters mostly live outside the conservancy, the majority of hunting occurs inside the conservancy (in criminology terms, commuter crime). These results combined with rigorous simulations showed for the first time how geographic profiling can be combined with GIS data and applied to situations with more complex spatial patterns, for example, where landscape heterogeneity means some parts of the study area are less likely to be used (e.g., aquatic areas for terrestrial animals) or where landscape permeability differs (e.g., forest bats tend not to fly over open areas). More broadly, these results show how geographic profiling can be used to target antipoaching interventions more effectively and more efficiently and to develop management strategies and conservation plans in a range of conservation scenarios.     

Antipoaching standards in onshore hydrocarbon concessions drawn from a Central African case study

Unsustainable hunting outside protected areas is threatening tropical biodiversity worldwide and requires conservationists to engage increasingly in antipoaching activities. Following the example of ecocertified logging companies, we argue that other extractive industries managing large concessions should engage in antipoaching activities as part of their environmental management plans. Onshore hydrocarbon concessions should also adopt antipoaching protocols as a standard because they represent a biodiversity threat comparable to logging. We examined the spatiotemporal patterns of small‐ and large‐mammal poaching in an onshore oil concession in Gabon, Central Africa, with a Bayesian occupancy model based on signs of poaching collected from 2010 to 2015 on antipoaching patrols. Patrol locations were initially determined based on local intelligence and past patrol successes (adaptive management) and subsequently with a systematic sampling of the concession. We generated maps of poaching probability in the concession and determined the temporal trends of this threat over 5 years. The spatiotemporal patterns of large‐ and small‐mammal poaching differed throughout the concession, and likely these groups will need different management strategies. By elucidating the relationship between site‐specific sampling effort and detection probability, the Bayesian method allowed us to set goals for future antipoaching patrols. Our results indicate that a combination of systematic sampling and adaptive management data is necessary to infer spatiotemporal patterns with the statistical method we used. On the basis of our case study, we recommend hydrocarbon companies interested in implementing efficient antipoaching activities in their onshore concessions to lay the foundation of long‐needed industry standards by: adequately measuring antipoaching effort; mixing adaptive management and balanced sampling; setting goals for antipoaching effort; pairing patrols with large‐mammal monitoring; supporting antipoaching patrols across the landscape; restricting access to their concessions; performing random searches for bushmeat and mammal products at points of entry; controlling urban and agricultural expansion; supporting bushmeat alternatives; and supporting land‐use planning.     

Testing Decision Rules for Categorizing Species’ Extinction Risk to Help Develop Quantitative Listing Criteria for the U.S. Endangered Species Act

Lack of guidance for interpreting the definitions of endangered and threatened in the U.S. Endangered Species Act (ESA) has resulted in case‐by‐case decision making leaving the process vulnerable to being considered arbitrary or capricious. Adopting quantitative decision rules would remedy this but requires the agency to specify the relative urgency concerning extinction events over time, cutoff risk values corresponding to different levels of protection, and the importance given to different types of listing errors. We tested the performance of 3 sets of decision rules that use alternative functions for weighting the relative urgency of future extinction events: a threshold rule set, which uses a decision rule of x% probability of extinction over y years; a concave rule set, where the relative importance of future extinction events declines exponentially over time; and a shoulder rule set that uses a sigmoid shape function, where relative importance declines slowly at first and then more rapidly. We obtained decision cutoffs by interviewing several biologists and then emulated the listing process with simulations that covered a range of extinction risks typical of ESA listing decisions. We evaluated performance of the decision rules under different data quantities and qualities on the basis of the relative importance of misclassification errors. Although there was little difference between the performance of alternative decision rules for correct listings, the distribution of misclassifications differed depending on the function used. Misclassifications for the threshold and concave listing criteria resulted in more overprotection errors, particularly as uncertainty increased, whereas errors for the shoulder listing criteria were more symmetrical. We developed and tested the framework for quantitative decision rules for listing species under the U.S. ESA. If policy values can be agreed on, use of this framework would improve the implementation of the ESA by increasing transparency and consistency. Evaluando Reglas de Decisión para Categorizar el Riesgo de Extinción de Especies con el Fin de Desarrollar de Criterios Cuantitativos de Alistamiento en el Acta de Especies en Peligro de los EE. UU.     

The Importance of Being Spatial (and Reserved): Assessing Northern Spotted Owl Habitat Relationships with Hierarchical Bayesian Models

Abstract:  Regional conservation planning increasingly draws on habitat suitability models to support decisions regarding land allocation and management. Nevertheless, statistical techniques commonly used for developing such models may give misleading results because they fail to account for 3 factors common in data sets of species distribution: spatial autocorrelation, the large number of sites where the species is absent (zero inflation), and uneven survey effort. We used spatial autoregressive models fit with Bayesian Markov Chain Monte Carlo techniques to assess the relationship between older coniferous forest and the abundance of Northern Spotted Owl nest and activity sites throughout the species' range. The spatial random-effect term incorporated in the autoregressive models successfully accounted for zero inflation and reduced the effect of survey bias on estimates of species–habitat associations. Our results support the hypothesis that the relationship between owl distribution and older forest varies with latitude. A quadratic relationship between owl abundance and older forest was evident in the southern portion of the range, and a pseudothreshold relationship was evident in the northern portion of the range. Our results suggest that proposed changes to the network of owl habitat reserves would reduce the proportion of the population protected by up to one-third, and that proposed guidelines for forest management within reserves underestimate the proportion of older forest associated with maximum owl abundance and inappropriately generalize threshold relationships among subregions. Bayesian spatial models can greatly enhance the utility of habitat analysis for conservation planning because they add the statistical flexibility necessary for analyzing regional survey data while retaining the interpretability of simpler models.     

Estimating freshwater turtle mortality rates and population declines following hook ingestion

Freshwater turtle populations are susceptible to declines following small increases in the mortality of adults, making it essential to identify and understand potential threats. Freshwater turtles ingest fish hooks associated with recreational angling, and this is likely a problem because hook ingestion is a source of additive mortality for sea turtles. We used a Bayesian‐modeling framework, observed rates of hook ingestion by freshwater turtles, and mortality of sea turtles from hook ingestion to examine the probability that a freshwater turtle in a given population ingests a hook and subsequently dies from it. We used the results of these analyses and previously published life‐history data to simulate the effects of hook ingestion on population growth for 3 species of freshwater turtle. In our simulation, the probability that an individual turtle ingests a hook and dies as a result was 1.2–11%. Our simulation results suggest that this rate of mortality from hook ingestion is sufficient to cause population declines. We believe we have identified fish‐hook ingestion as a serious yet generally overlooked threat to the viability of freshwater turtle populations.     

A Comparison of Bayesian Methods for Uncertainty Analysis in Hydraulic and Hydrodynamic Modeling

We evaluate and compare the performance of Bayesian Monte Carlo (BMC), Markov chain Monte Carlo (MCMC), and the Generalized Likelihood Uncertainty Estimation (GLUE) for uncertainty analysis in hydraulic and hydrodynamic modeling (HHM) studies. The methods are evaluated in a synthetic 1D wave routing exercise based on the diffusion wave model, and in a multidimensional hydrodynamic study based on the Environmental Fluid Dynamics Code to simulate estuarine circulation processes in Weeks Bay, Alabama. Results show that BMC and MCMC provide similar estimates of uncertainty. The posterior parameter densities computed by both methods are highly consistent, as well as the calibrated parameter estimates and uncertainty bounds. Although some studies suggest that MCMC is more efficient than BMC, our results did not show a clear difference between the performance of the two methods. This seems to be due to the low number of model parameters typically involved in HHM studies, and the use of the same likelihood function. In fact, for these studies, the implementation of BMC results simpler and provides similar results to MCMC. The results of GLUE are, on the other hand, less consistent to the results of BMC and MCMC in both applications. The posterior probability densities tend to be flat and similar to the uniform priors, which can result in calibrated parameter estimates centered in the parametric space.     

Dynamic failure assessment of an ammonia storage unit: A case study

Chemical Process Industries usually contain a diverse inventory of hazardous chemicals and complex systems required to perform process operations such as storage, separation, reaction, compression etc. The complex interactions between the equipment make them vulnerable to catastrophic accidents. Risk and failure assessment provide engineers with an intuitive tool for decision making in the operation of such plants. Abnormal events and near-miss situations occur regularly during the operation of a system. Accident Sequence Precursors (ASP) can be used to demonstrate the real-time operating condition of a plant. Dynamic Failure Assessment (DFA) methodology is based on Bayesian statistical methods incorporates ASP data to revise the generic failure probabilities of the systems during its operational lifetime.In this paper, DFA methodology is applied on an ammonia storage unit in a specialized chemical industry. Ammonia is stored in cold storage tanks as liquefied gas at atmospheric pressure. These tanks are susceptible to failures due to various abnormal conditions arising due process failures.Tank failures due to three such abnormal conditions are considered. Variation of the failure probability of the safety systems is demonstrated. The authors use ASP data collected from plant specific sources and safety expert judgement. The failure probabilities of some safety systems concerned show considerable deviation from the generic values. The method helps to locate the components which have undergone more degradation over the period and hence must be paid attention to. In addition, a Bayesian predictive model has been used to predict the number of abnormal events in the next time interval. The user-friendly and intuitive nature of the tool makes it appropriate for application in safety assessment reports in process industries.     

Operational risk assessment: A case of the Bhopal disaster

Accidental releases of hazardous chemicals from process facilities can cause catastrophic consequences. The Bhopal disaster resulting from a combination of inherently unsafe designs and poorly managed operations is a well-known case. Effective risk modeling approaches that provide early warnings are helpful to prevent and control such rare but catastrophic events. Probability estimation of these events is a constant challenge due to the scarcity of directly relevant data. Therefore, precursor-based methods that adopt the Bayesian theorem to update prior judgments on event probabilities using empirical data have been proposed. The updated probabilities are then integrated with consequences of varying severity to produce the risk profile.This paper proposes an operational risk assessment framework, in which a precursor-based Bayesian network approach is used for probability estimation, and loss functions are applied for consequence assessment. The estimated risk profile can be updated continuously given real-time operational data. As process facilities operate, this method integrates a failure-updating mechanism with potential consequences to generate a real-time operational risk profile. The real time risk profile is valuable in activating accident prevention and control strategies. The approach is applied to the Bhopal accident to demonstrate its applicability and effectiveness.     

The theoretical value of encounters with parasitized hosts for parasitoids

A female parasitoid searching for hosts in a patch experiences a diminishing encounter rate with unparasitized and thus suitable hosts. To use the available time most efficiently, it constantly has to decide whether to stay in the patch and continue to search for hosts or to search for and travel to another patch in the habitat. Several informational cues can be used to optimize the searching success. Theoretically, encounters with unparasitized hosts should lead to a prolonged search in a given patch if hosts are distributed contagiously. The results of empirical studies strongly support this hypothesis. However, it has, to date, not been investigated theoretically whether encounters with already parasitized hosts (which usually entail time costs) provide a parasitoid with valuable information for the optimization of its search in depletable patches, although the empirical studies concerning this question so far have produced ambiguous results. Building on recent advances in Bayesian foraging strategies, we approached this problem by modeling a priori searching strategies (which differ in the amount of information considered) and then testing them in computer simulations. By comparing the strategies, we were able to determine whether and how encounters with already parasitized hosts can yield information that can be used to enhance a parasitoid’s searching success.