Spatial Stochastic Simulation Offers Potential as a Quantitative Method for Pest Risk Analysis

Pest risk analysis represents an emerging field of risk analysis that evaluates the potential risks of the introduction and establishment of plant pests into a new geographic location and then assesses the management options to reduce those potential risks. Development of new and adapted methodology is required to answer questions concerning pest risk analysis of exotic plant pests. This research describes a new method for predicting the potential establishment and spread of a plant pest into new areas using a case study, Ralstonia solanacearum, a bacterial disease of potato. This method combines current quantitative methodologies, stochastic simulation, and geographic information systems with knowledge of pest biology and environmental data to derive new information about pest establishment potential in a geographical region where a pest had not been introduced. This proposed method extends an existing methodology for matching pest characteristics with environmental conditions by modeling and simulating dissemination behavior of a pest organism. Issues related to integrating spatial variables into risk analysis models are further discussed in this article.     

Allee Effects and the Risk of Biological Invasion

The Allee effect is a nonlinear phenomenon exhibited in the population dynamics of sparse populations in which the per capita population growth rate increases with increasing population density. In sufficiently sparse populations, the Allee effect may lead to extinction and is known to generate a threshold in the probability of establishment when presented as a function of introduced population size or density. As introduced populations are generally small, Allee effects are probably common in biological invasions and their consideration is necessary for accurately assessing the risk of invasion by many species, including all sexually reproducing species. Bythotrephes longimanus, an invasive, freshwater, cladoceran zooplankter from Europe, is one such species. Here, I review a previously published model of the Allee effect for continuously sexually reproducing species. Then, I develop a new model for seasonally parthenogenetic species such as Bythotrephes, and thereby demonstrate the potential consequences of Allee effects. This result underscores the importance of considering nonlinear phenomena, including thresholds, when conducting risk analysis for biological introductions.     

Improving Default Risk Prediction Using Bayesian Model Uncertainty Techniques

Credit risk is the potential exposure of a creditor to an obligor's failure or refusal to repay the debt in principal or interest. The potential of exposure is measured in terms of probability of default. Many models have been developed to estimate credit risk, with rating agencies dating back to the 19th century. They provide their assessment of probability of default and transition probabilities of various firms in their annual reports. Regulatory capital requirements for credit risk outlined by the Basel Committee on Banking Supervision have made it essential for banks and financial institutions to develop sophisticated models in an attempt to measure credit risk with higher accuracy. The Bayesian framework proposed in this article uses the techniques developed in physical sciences and engineering for dealing with model uncertainty and expert accuracy to obtain improved estimates of credit risk and associated uncertainties. The approach uses estimates from one or more rating agencies and incorporates their historical accuracy (past performance data) in estimating future default risk and transition probabilities. Several examples demonstrate that the proposed methodology can assess default probability with accuracy exceeding the estimations of all the individual models. Moreover, the methodology accounts for potentially significant departures from “nominal predictions” due to “upsetting events” such as the 2008 global banking crisis.     

Optimal Strategies for Interception,Detection, and Eradication in Plant Biosecurity

The introduction of invasive species causes damages from the economic and ecological point of view. Interception of plant pests and eradication of the established populations are two management options to prevent or limit the risk posed by an invasive species. Management options generate costs related to the interception at the point of entry, and the detection and eradication of established field populations. Risk managers have to decide how to allocate resources between interception, field detection, containment, and eradication minimizing the expected total costs. In this work is considered an optimization problem aiming at determining the optimal allocation of resources to minimize the expected total costs of the introduction of Bemisia tabaci‐transmitted viruses in Europe. The optimization problem takes into account a probabilistic model for the estimation of the percentage of viruliferous insect populations arriving through the trade of commodities, and a population dynamics model describing the process of the vector populations' establishment and spread. The time of field detection of viruliferous insect populations is considered as a random variable. The solution of the optimization problem allows to determine the optimal allocation of the search effort between interception and detection/eradication. The behavior of the search effort as a function of efficacy or search in interception and in detection is then analyzed. The importance of the vector population growth rate and the probability of virus establishment are also considered in the analysis of the optimization problem.     

Predictive Modeling and Categorizing Likelihoods of Quarantine Pest Introduction of Imported Propagative Commodities from Different Countries

The present study investigates U.S. Department of Agriculture inspection records in the Agricultural Quarantine Activity System database to estimate the probability of quarantine pests on propagative plant materials imported from various countries of origin and to develop a methodology ranking the risk of country–commodity combinations based on quarantine pest interceptions. Data collected from October 2014 to January 2016 were used for developing predictive models and validation study. A generalized linear model with Bayesian inference and a generalized linear mixed effects model were used to compare the interception rates of quarantine pests on different country–commodity combinations. Prediction ability of generalized linear mixed effects models was greater than that of generalized linear models. The estimated pest interception probability and confidence interval for each country–commodity combination was categorized into one of four compliance levels: “High,” “Medium,” “Low,” and “Poor/Unacceptable,” Using K‐means clustering analysis. This study presents risk‐based categorization for each country–commodity combination based on the probability of quarantine pest interceptions and the uncertainty in that assessment.     

Risk Analysis for Biological Hazards: What We Need to Know about Invasive Species

Risk analysis for biological invasions is similar to other types of natural and human hazards. For example, risk analysis for chemical spills requires the evaluation of basic information on where a spill occurs; exposure level and toxicity of the chemical agent; knowledge of the physical processes involved in its rate and direction of spread; and potential impacts to the environment, economy, and human health relative to containment costs. Unlike typical chemical spills, biological invasions can have long lag times from introduction and establishment to successful invasion, they reproduce, and they can spread rapidly by physical and biological processes. We use a risk analysis framework to suggest a general strategy for risk analysis for invasive species and invaded habitats. It requires: (1) problem formation (scoping the problem, defining assessment endpoints); (2) analysis (information on species traits, matching species traits to suitable habitats, estimating exposure, surveys of current distribution and abundance); (3) risk characterization (understanding of data completeness, estimates of the "potential" distribution and abundance; estimates of the potential rate of spread; and probable risks, impacts, and costs); and (4) risk management (containment potential, costs, and opportunity costs; legal mandates and social considerations and information science and technology needs).     

Quantile Uncertainty and Value‐at‐Risk Model Risk

This article develops a methodology for quantifying model risk in quantile risk estimates. The application of quantile estimates to risk assessment has become common practice in many disciplines, including hydrology, climate change, statistical process control, insurance and actuarial science, and the uncertainty surrounding these estimates has long been recognized. Our work is particularly important in finance, where quantile estimates (called Value‐at‐Risk) have been the cornerstone of banking risk management since the mid 1980s. A recent amendment to the Basel II Accord recommends additional market risk capital to cover all sources of “model risk” in the estimation of these quantiles. We provide a novel and elegant framework whereby quantile estimates are adjusted for model risk, relative to a benchmark which represents the state of knowledge of the authority that is responsible for model risk. A simulation experiment in which the degree of model risk is controlled illustrates how to quantify Value‐at‐Risk model risk and compute the required regulatory capital add‐on for banks. An empirical example based on real data shows how the methodology can be put into practice, using only two time series (daily Value‐at‐Risk and daily profit and loss) from a large bank. We conclude with a discussion of potential applications to nonfinancial risks.     

Probabilistic Framework for the Estimation of the Adult and Child Toxicokinetic Intraspecies Uncertainty Factors

Human health risk assessments use point values to develop risk estimates and thus impart a deterministic character to risk, which, by definition, is a probability phenomenon. The risk estimates are calculated based on individuals and then, using uncertainty factors (UFs), are extrapolated to the population that is characterized by variability. Regulatory agencies have recommended the quantification of the impact of variability in risk assessments through the application of probabilistic methods. In the present study, a framework that deals with the quantitative analysis of uncertainty (U) and variability (V) in target tissue dose in the population was developed by applying probabilistic analysis to physiologically-based toxicokinetic models. The mechanistic parameters that determine kinetics were described with probability density functions (PDFs). Since each PDF depicts the frequency of occurrence of all expected values of each parameter in the population, the combined effects of multiple sources of U/V were accounted for in the estimated distribution of tissue dose in the population, and a unified (adult and child) intraspecies toxicokinetic uncertainty factor UFH-TK was determined. The results show that the proposed framework accounts effectively for U/V in population toxicokinetics. The ratio of the 95th percentile to the 50th percentile of the annual average concentration of the chemical at the target tissue organ (i.e., the UFH-TK) varies with age. The ratio is equivalent to a unified intraspecies toxicokinetic UF, and it is one of the UFs by which the NOAEL can be divided to obtain the RfC/RfD. The 10-fold intraspecies UF is intended to account for uncertainty and variability in toxicokinetics (3.2x) and toxicodynamics (3.2x). This article deals exclusively with toxicokinetic component of UF. The framework provides an alternative to the default methodology and is advantageous in that the evaluation of toxicokinetic variability is based on the distribution of the effective target tissue dose, rather than applied dose. It allows for the replacement of the default adult and children intraspecies UF with toxicokinetic data-derived values and provides accurate chemical-specific estimates for their magnitude. It shows that proper application of probability and toxicokinetic theories can reduce uncertainties when establishing exposure limits for specific compounds and provide better assurance that established limits are adequately protective. It contributes to the development of a probabilistic noncancer risk assessment framework and will ultimately lead to the unification of cancer and noncancer risk assessment methodologies.     

A Framework for Understanding Uncertainty in Seismic Risk Assessment

A better understanding of the uncertainty that exists in models used for seismic risk assessment is critical to improving risk-based decisions pertaining to earthquake safety. Current models estimating the probability of collapse of a building do not consider comprehensively the nature and impact of uncertainty. This article presents a model framework to enhance seismic risk assessment and thus gives decisionmakers a fuller understanding of the nature and limitations of the estimates. This can help ensure that risks are not over- or underestimated and the value of acquiring accurate data is appreciated fully. The methodology presented provides a novel treatment of uncertainties in input variables, their propagation through the model, and their effect on the results. The study presents ranges of possible annual collapse probabilities for different case studies on buildings in different parts of the world, exposed to different levels of seismicity, and with different vulnerabilities. A global sensitivity analysis was conducted to determine the significance of uncertain variables. Two key outcomes are (1) that the uncertainty in ground-motion conversion equations has the largest effect on the uncertainty in the calculation of annual collapse probability; and (2) the vulnerability of a building appears to have an effect on the range of annual collapse probabilities produced, i.e., the level of uncertainty in the estimate of annual collapse probability, with less vulnerable buildings having a smaller uncertainty.     

Domino Effect Analysis Using Bayesian Networks

A new methodology is introduced based on Bayesian network both to model domino effect propagation patterns and to estimate the domino effect probability at different levels. The flexible structure and the unique modeling techniques offered by Bayesian network make it possible to analyze domino effects through a probabilistic framework, considering synergistic effects, noisy probabilities, and common cause failures. Further, the uncertainties and the complex interactions among the domino effect components are captured using Bayesian network. The probabilities of events are updated in the light of new information, and the most probable path of the domino effect is determined on the basis of the new data gathered. This study shows how probability updating helps to update the domino effect model either qualitatively or quantitatively. The methodology is applied to a hypothetical example and also to an earlier‐studied case study. These examples accentuate the effectiveness of Bayesian network in modeling domino effects in processing facility.     

Bayesian Quantile Impairment Threshold Benchmark Dose Estimation for Continuous Endpoints

Quantitative risk assessment often begins with an estimate of the exposure or dose associated with a particular risk level from which exposure levels posing low risk to populations can be extrapolated. For continuous exposures, this value, the benchmark dose, is often defined by a specified increase (or decrease) from the median or mean response at no exposure. This method of calculating the benchmark dose does not take into account the response distribution and, consequently, cannot be interpreted based upon probability statements of the target population. We investigate quantile regression as an alternative to the use of the median or mean regression. By defining the dose–response quantile relationship and an impairment threshold, we specify a benchmark dose as the dose associated with a specified probability that the population will have a response equal to or more extreme than the specified impairment threshold. In addition, in an effort to minimize model uncertainty, we use Bayesian monotonic semiparametric regression to define the exposure–response quantile relationship, which gives the model flexibility to estimate the quantal dose–response function. We describe this methodology and apply it to both epidemiology and toxicology data.     

An Analytical Framework for Relating Dose, Risk, and Incidence: An Application to Occupational Tuberculosis Infection

An adverse health impact is often treated as a binary variable (response vs. no response), in which case the risk of response is defined as a monotonically increasing function R of the dose received D. For a population of size N , specifying the forms of R(D) and of the probability density function (pdf) for D allows determination of the pdf for risk, and computation of the mean and variance of the distribution of incidence, where the latter parameters are denoted E[S _N] and Var[ S _N], respectively. The distribution of S _N describes uncertainty in the future incidence value. Given variability in dose (and risk) among population members, the distribution of incidence is Poisson-binomial. However, depending on the value of E[S _N], the distribution of incidence is adequately approximated by a Poisson distribution with parameter μ= E[S _N], or by a normal distribution with mean and variance equal to E[S _N] and Var[ S _N]. The general analytical framework is applied to occupational infection by Mycobacterium tuberculosis (M. tb). Tuberculosis is transmitted by inhalation of 1–5 μm particles carrying viable M. tb bacilli. Infection risk has traditionally been modeled by the expression: R(D) = 1 – exp(– D ), where D is the expected number of bacilli that deposit in the pulmonary region. This model assumes that the infectious dose is one bacillus. The beta pdf and the gamma pdf are shown to be reasonable and especially convenient forms for modeling the distribution of the expected cumulative dose across a large healthcare worker cohort. Use of the the analytical framework is illustrated by estimating the efficacy of different respiratory protective devices in reducing healthcare worker infection risk.     

Evaluating Critical Uncertainty Thresholds in a Spatial Model of Forest Pest Invasion Risk

Pest risk maps can provide useful decision support in invasive species management, but most do not adequately consider the uncertainty associated with predicted risk values. This study explores how increased uncertainty in a risk model's numeric assumptions might affect the resultant risk map. We used a spatial stochastic model, integrating components for entry, establishment, and spread, to estimate the risks of invasion and their variation across a two-dimensional landscape for Sirex noctilio , a nonnative woodwasp recently detected in the United States and Canada. Here, we present a sensitivity analysis of the mapped risk estimates to variation in key model parameters. The tested parameter values were sampled from symmetric uniform distributions defined by a series of nested bounds (±5%, … , ±40%) around the parameters' initial values. The results suggest that the maximum annual spread distance, which governs long-distance dispersal, was by far the most sensitive parameter. At ±15% or larger variability bound increments for this parameter, there were noteworthy shifts in map risk values, but no other parameter had a major effect, even at wider bounds of variation. The methodology presented here is generic and can be used to assess the impact of uncertainties on the stability of pest risk maps as well as to identify geographic areas for which management decisions can be made confidently, regardless of uncertainty.     

General Probability‐Time Tradeoff and Intertemporal Risk‐Value Model

This article proposes an intertemporal risk‐value (IRV) model that integrates probability‐time tradeoff, time‐value tradeoff, and risk‐value tradeoff into one unified framework. We obtain a general probability‐time tradeoff, which yields a formal representation form to reflect the psychological distance of a decisionmaker in evaluating a temporal lottery. This intuition of probability‐time tradeoff is supported by robust empirical findings as well as by psychological theory. Through an explicit formalization of probability‐time tradeoff, an IRV model taking into account three fundamental dimensions, namely, value, probability, and time, is established. The object of evaluation in our framework is a complex lottery. We also give some insights into the structure of the IRV model using a wildcatter problem.     

A Scale of Risk

This article proposes a conceptual framework for ranking the relative gravity of diverse risks. This framework identifies the moral considerations that should inform the evaluation and comparison of diverse risks. A common definition of risk includes two dimensions: the probability of occurrence and the associated consequences of a set of hazardous scenarios. This article first expands this definition to include a third dimension: the source of a risk. The source of a risk refers to the agents involved in the creation or maintenance of a risk and captures a central moral concern about risks. Then, a scale of risk is proposed to categorize risks along a multidimensional ranking, based on a comparative evaluation of the consequences, probability, and source of a given risk. A risk is ranked higher on the scale the larger the consequences, the greater the probability, and the more morally culpable the source. The information from the proposed comparative evaluation of risks can inform the selection of priorities for risk mitigation.     

Ecological Risk Assessment Framework for Low-Altitude Aircraft Overflights: I. Planning the Analysis and Estimating Exposure

An ecological risk assessment framework for low-altitude aircraft overflights was developed, with special emphasis on military applications. The problem formulation and exposure analysis phases are presented in this article; an analysis of effects and risk characterization is presented in a companion article. The intent of this article is threefold: (1) to illustrate the development of a generic framework for the ecological risk assessment of an activity, (2) to show how the U.S. Environmental Protection Agency's ecological risk assessment paradigm can be applied to an activity other than the release of a chemical, and (3) to provide guidance for the assessment of ecological risks from low-altitude aircraft overflights. The key stressor for low-altitude aircraft overflights is usually sound, although visual and physical (collision) stressors may also be significant. Susceptible and regulated wildlife populations are the major assessment endpoint entities, although plant communities may be impacted by takeoffs and landings. The exposure analysis utilizes measurements of wildlife locations, measurements of sound levels at the wildlife locations, measurements of slant distances from aircraft to wildlife, models that extrapolate sound from the source aircraft to the ground, and bird-strike probability models. Some of the challenges to conducting a risk assessment for aircraft overflights include prioritizing potential stressors and endpoints, choosing exposure metrics that relate to wildlife responses, obtaining good estimates of sound or distance, and estimating wildlife locations.     

A Probabilistic Framework for Risk Analysis of Widespread Flood Events: A Proof‐of‐Concept Study

This article presents a flood risk analysis model that considers the spatially heterogeneous nature of flood events. The basic concept of this approach is to generate a large sample of flood events that can be regarded as temporal extrapolation of flood events. These are combined with cumulative flood impact indicators, such as building damages, to finally derive time series of damages for risk estimation. Therefore, a multivariate modeling procedure that is able to take into account the spatial characteristics of flooding, the regionalization method top‐kriging, and three different impact indicators are combined in a model chain. Eventually, the expected annual flood impact (e.g., expected annual damages) and the flood impact associated with a low probability of occurrence are determined for a study area. The risk model has the potential to augment the understanding of flood risk in a region and thereby contribute to enhanced risk management of, for example, risk analysts and policymakers or insurance companies. The modeling framework was successfully applied in a proof‐of‐concept exercise in Vorarlberg (Austria). The results of the case study show that risk analysis has to be based on spatially heterogeneous flood events in order to estimate flood risk adequately.     

Risk Assessment for Aflatoxin: An Evaluation Based on the Multistage Model

Lifetime cancer potency of alfatoxin was assessed based on the Yeh et al. study from China in which both aflatoxin exposure and hepatitis B prevalence were measured. This study provides the best available information for estimating the carcinogenic risk posed by aflatoxin to the U.S. population. Cancer potency of aflatoxin was estimated using a biologically motivated risk assessment model. The best estimate of aflatoxin potency was 9 (mg/kg/day)⁻¹ for individuals negative for hepatitis B and 230 (mg/kg/day)⁻¹ for individuals positive for hepatitis B.     

Risk filtering, ranking, and management framework using hierarchical holographic modeling.

This paper contributes a methodological framework to identify, prioritize, assess, and manage risk scenarios of a large-scale system. Qualitative screening of scenarios and classes of scenarios is appropriate initially, while quantitative assessments may be applied once the set of all scenarios (hundreds) has been prioritized in several phases. The eight-phase methodology is described in detail and is applied to operations other than war. The eight phases are as follows: Phase I, Scenario Identification-A hierarchical holographic model (HHM) is developed to describe the system's "as planned" or "success" scenario. Phase II, Scenario Filtering-The risk scenarios identified in Phase I are filtered according to the responsibilities and interests of the current system user. Phase III, Bi-Criteria Filtering and Ranking. Phase IV, Multi-Criteria Evaluation. Phase V, Quantitative Ranking-We continue to filter and rank scenarios based on quantitative and qualitative matrix scales of likelihood and consequence; and ordinal response to system resiliency, robustness, redundancy. Phase VI, Risk Management is performed, involving identification of management options for dealing with the filtered scenarios, and estimating the cost, performance benefits, and risk reduction of each. Phase VII, Safeguarding Against Missing Critical Items--We examine the performance of the options selected in Phase VI against the scenarios previously filtered out during Phases II to V. Phase VIII, Operational Feedback-We use the experience and information gained during application to refine the scenario filtering and decision processes in earlier phases. These eight phases reflect a philosophical approach rather than a mechanical methodology. In this philosophy, the filtering and ranking of discrete scenarios is viewed as a precursor to, rather than a substitute for, consideration of the totality of all risk scenarios.     

概率预期下在线报童问题的最小风险策略

报童问题是库存管理中一个基本模型。已有的报童模型主要利用均值-方差方法和期望效用目标方法进行风险的描述和度量。这些方法假设需求分布信息已知,而实际中需求分布信息往往难以完全刻画。本文使用概率预期作为刻画不完全需求分布信息的格式,基于在线风险补偿的思想,为需求分布信息不完全的报童问题建立了最小风险模型。使用该模型设计了最小风险策略,使报童可以根据自己设定的不同收益和未来概率预期选择最优订购量。