A Bayesian network approach to refining ecological risk assessments: Mercury and the Florida panther (Puma concolor coryi)

Development of a probabilistic risk model for pharmaceuticals in the environment under population and wastewater treatment scenarios

Preparing for future environmental pressures requires projections of how relevant risks will change over time. Current regulatory models of environmental risk assessment (ERA) of pollutants such as pharmaceuticals could be improved by considering the influence of global change factors (e.g., population growth) and by presenting uncertainty more transparently. In this article, we present the development of a prototype object-oriented Bayesian network (BN) for the prediction of environmental risk for six high-priority pharmaceuticals across 36 scenarios: current and three future population scenarios, combined with infrastructure scenarios, in three Norwegian counties. We compare the risk, characterized by probability distributions of risk quotients (RQs), across scenarios and pharmaceuticals. Our results suggest that RQs would be greatest in rural counties, due to the lower development of current wastewater treatment facilities, but that these areas consequently have the most potential for risk mitigation. This pattern intensifies under higher population growth scenarios. With this prototype, we developed a hierarchical probabilistic model and demonstrated its potential in forecasting the environmental risk of chemical stressors under plausible demographic and management scenarios, contributing to the further development of BNs for ERA. Integr Environ Assess Manag 2024;20:1715-1735. © 2024 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Key challenges and developments in wildlife ecological risk assessment: Problem formulation

Problem formulation (PF) is a critical initial step in planning risk assessments for chemical exposures to wildlife, used either explicitly or implicitly in various jurisdictions to include registration of new pesticides, evaluation of new and existing chemicals released to the environment, and characterization of impact when chemical releases have occurred. Despite improvements in our understanding of the environment, ecology, and biological sciences, few risk assessments have used this information to enhance their value and predictive capabilities. In addition to advances in organism-level mechanisms and methods, there have been substantive developments that focus on population- and systems-level processes. Although most of the advances have been recognized as being state-of-the-science for two decades or more, there is scant evidence that they have been incorporated into wildlife risk assessment or risk assessment in general. In this article, we identify opportunities to consider elevating the relevance of wildlife risk assessments by focusing on elements of the PF stage of risk assessment, especially in the construction of conceptual models and selection of assessment endpoints that target population- and system-level endpoints. Doing so will remain consistent with four established steps of existing guidance: (1) establish clear protection goals early in the process; (2) consider how data collection using new methods will affect decisions, given all possibilities, and develop a decision plan a priori; (3) engage all relevant stakeholders in creating a robust, holistic conceptual model that incorporates plausible stressors that could affect the targets defined in the protection goals; and (4) embrace the need for iteration throughout the PF steps (recognizing that multiple passes may be required before agreeing on a feasible plan for the rest of the risk assessment). Integr Environ Assess Manag 2024;20:658-673. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Incorporating climate projections in the environmental risk assessment of pesticides in aquatic ecosystems

Global climate change will significantly impact the biodiversity of freshwater ecosystems, both directly and indirectly via the exacerbation of impacts from other stressors. Pesticides form a prime example of chemical stressors that are expected to synergize with climate change. Aquatic exposures to pesticides might change in magnitude due to increased runoff from agricultural fields, and in composition, as application patterns will change due to changes in pest pressures and crop types. Any prospective chemical risk assessment that aims to capture the influence of climate change should properly and comprehensively account for the variabilities and uncertainties that are inherent to projections of future climate. This is only feasible if they probabilistically propagate extensive ensembles of climate model projections. However, current prospective risk assessments typically make use of process-based models of chemical fate that do not typically allow for such high-throughput applications. Here, we describe a Bayesian network model that does. It incorporates a two-step univariate regression model based on a 30-day antecedent precipitation index, circumventing the need for computationally laborious mechanistic models. We show its feasibility and application potential in a case study with two pesticides in a Norwegian stream: the fungicide trifloxystrobin and herbicide clopyralid. Our analysis showed that variations in pesticide application rates as well as precipitation intensity lead to variations in in-stream exposures. When relating to aquatic risks, the influence of these processes is reduced and distributions of risk are dominated by effect-related parameters. Predicted risks for clopyralid were negligible, but the probability of unacceptable future environmental risks due to exposure to trifloxystrobin (i.e., a risk quotient >1) was 8%-12%. This percentage further increased to 30%-35% when a more conservative precautionary factor of 100 instead of 30 was used. Integr Environ Assess Manag 2024;20:384-400. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

A novel approach to assessing natural resource injury with Bayesian networks

Quantifying the effects of environmental stressors on natural resources is problematic because of complex interactions among environmental factors that influence endpoints of interest. This complexity, coupled with data limitations, propagates uncertainty that can make it difficult to causally associate specific environmental stressors with injury endpoints. The Natural Resource Damage Assessment and Restoration (NRDAR) regulations under the Comprehensive Environmental Response, Compensation, and Liability Act and Oil Pollution Act aim to restore natural resources injured by oil spills and hazardous substances released into the environment; exploration of alternative statistical methods to evaluate effects could help address NRDAR legal claims. Bayesian networks (BNs) are statistical tools that can be used to estimate the influence and interrelatedness of abiotic and biotic environmental variables on environmental endpoints of interest. We investigated the application of a BN for injury assessment using a hypothetical case study by simulating data of acid mine drainage (AMD) affecting a fictional stream-dwelling bird species. We compared the BN-generated probability estimates for injury with a more traditional approach using toxicity thresholds for water and sediment chemistry. Bayesian networks offered several distinct advantages over traditional approaches, including formalizing the use of expert knowledge, probabilistic estimates of injury using intermediate direct and indirect effects, and the incorporation of a more nuanced and ecologically relevant representation of effects. Given the potential that BNs have for natural resource injury assessment, more research and field-based application are needed to determine their efficacy in NRDAR. We expect the resulting methods will be of interest to many US federal, state, and tribal programs devoted to the evaluation, mitigation, remediation, and/or restoration of natural resources injured by releases or spills of contaminants. Integr Environ Assess Manag 2024;20:562-573. Published 2023. This article is a U.S. Government work and is in the public domain in the USA. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Integrating climate model projections into environmental risk assessment: A probabilistic modeling approach

The Society of Environmental Toxicology and Chemistry (SETAC) convened a Pellston workshop in 2022 to examine how information on climate change could be better incorporated into the ecological risk assessment (ERA) process for chemicals as well as other environmental stressors. A major impetus for this workshop is that climate change can affect components of ecological risks in multiple direct and indirect ways, including the use patterns and environmental exposure pathways of chemical stressors such as pesticides, the toxicity of chemicals in receiving environments, and the vulnerability of species of concern related to habitat quality and use. This article explores a modeling approach for integrating climate model projections into the assessment of near- and long-term ecological risks, developed in collaboration with climate scientists. State-of-the-art global climate modeling and downscaling techniques may enable climate projections at scales appropriate for the study area. It is, however, also important to realize the limitations of individual global climate models and make use of climate model ensembles represented by statistical properties. Here, we present a probabilistic modeling approach aiming to combine projected climatic variables as well as the associated uncertainties from climate model ensembles in conjunction with ERA pathways. We draw upon three examples of ERA that utilized Bayesian networks for this purpose and that also represent methodological advancements for better prediction of future risks to ecosystems. We envision that the modeling approach developed from this international collaboration will contribute to better assessment and management of risks from chemical stressors in a changing climate. Integr Environ Assess Manag 2024;20:367-383. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

A Bayesian network-GIS probabilistic model for addressing human disturbance risk to ecological conservation redline areas

Population growth and associated ecological space occupation are posing great risks to regional ecological security and social stability. In China, "Ecological Conservation Redline" (ECR) that prohibited urbanization and industrial construction has been proposed as a national policy to resolve spatial mismatches and management contradictions. However, unfriendly human disturbance activities (e.g., cultivation, mining, and infrastructure construction) still exist within the ECR, posing a great threat to ecological stability and safety. In this article, a Bayesian network (BN)-GIS probabilistic model is proposed to spatially and quantitatively address the human disturbance risk to the ECR at the regional scale. The Bayesian models integrate multiple human activities, ecological receptors of the ECR, and their exposure relationships for calculating the human disturbance risk. The case learning method geographic information systems (GIS) is then introduced to train BN models based on the spatial attribute of variables to evaluate the spatial distribution and correlation of risks. This approach was applied to the human disturbance risk assessment for the ECR that was delineated in 2018 in Jiangsu Province, China. The results indicated that most of the ECRs were at a low or medium human disturbance risk level, while some drinking water sources and forest parks in Lianyungang City possessed the highest risk. The sensitivity analysis result showed the ECR vulnerability, especially for cropland, that contributed most to the human disturbance risk. This spatially probabilistic method can not only enhance model's prediction precision, but also help decision-makers to determine how to establish priorities for policy design and conservation interventions. Overall, it presents a foundation for later ECR adjustments as well as for human disturbance risk supervision and management at the regional scale.

Bayesian inference of heavy metals exposure in crayfish for assessing human non-carcinogenic health risk

Contaminant concentrations often presented left censorship that below the limit of detection (LOD), which may contain true zero values because of no residue. In this study, we analyzed the concentrations of lead (Pb), cadmium (Cd), mercury (Hg), arsenic (As), and chromium (Cr) in 391 crayfish samples collected from 24 provinces of China, modelled the concentrations with two different models in a Bayesian framework, and evaluated health risk by estimating the target hazard quotient (THQ). The highest concentration was associated with Cr in Northwest China (0.1753 ± 0.1520 mg/kg w. w.), and the minimum concentration was Cd in Southwest China (0.0052 ± 0.0144 mg/kg w. w.), all heavy metal concentrations were below their safety limits. The posterior means of not detect rates P₀ of Pb, Cd, Hg, As, and Cr obtained with two models were both nearly equal to the observed not detect rates (51.15%, 36.83%, 27.37%, 64.71%, and 43.73%, respectively). The posterior probabilistic density lines for Pb, Hg, and Cd concentrations obtained with two models were similar, and fitted the empirical distributions well. The posterior density lines of THQs showed that the non-carcinogenic risk of As and Hg were significant high. Moreover, Bayesian approach presented a better understanding of the percentage of population exposed to potential risk.

Bayesian Network Applications for Sustainable Holistic Water Resources Management: Modeling Opportunities for South Africa

Anthropogenic transformation of land globally is threatening water resources in terms of quality and availability. Managing water resources to ensure sustainable utilization is important for a semiarid country such as South Africa. Bayesian networks (BNs) are probabilistic graphical models that have been applied globally to a range of water resources management studies; however, there has been very limited application of BNs to similar studies in South Africa. This article explores the benefits and challenges of BN application in the context of water resources management, specifically in relation to South Africa. A brief overview describes BNs, followed by details of some of the possible opportunities for BNs to benefit water resources management. These include the ability to use quantitative and qualitative information, data, and expert knowledge. BN models can be integrated into geographic information systems and predict impact of ecosystem services and sustainability indicators. With additional data and information, BNs can be updated, allowing for integration into an adaptive management process. Challenges in the application of BNs include oversimplification of complex systems, constraints of BNs with categorical nodes for continuous variables, unclear use of expert knowledge, and treatment of uncertainty. BNs have tremendous potential to guide decision making by providing a holistic approach to water resources management.

Development of a Bayesian network for probabilistic risk assessment of pesticides

Conventional environmental risk assessment of chemicals is based on a calculated risk quotient, representing the ratio of exposure to effects of the chemical, in combination with assessment factors to account for uncertainty. Probabilistic risk assessment approaches can offer more transparency by using probability distributions for exposure and/or effects to account for variability and uncertainty. In this study, a probabilistic approach using Bayesian network modeling is explored as an alternative to traditional risk calculation. Bayesian networks can serve as meta-models that link information from several sources and offer a transparent way of incorporating the required characterization of uncertainty for environmental risk assessment. To this end, a Bayesian network has been developed and parameterized for the pesticides azoxystrobin, metribuzin, and imidacloprid. We illustrate the development from deterministic (traditional) risk calculation, via intermediate versions, to fully probabilistic risk characterization using azoxystrobin as an example. We also demonstrate the seasonal risk calculation for the three pesticides. Integr Environ Assess Manag 2022;18:1072-1087. © 2021 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

A Bayesian Approach to Incorporating Spatiotemporal Variation and Uncertainty Limits into Modeling of Predicted Environmental Concentrations from Chemical Monitoring Campaigns

Environmental monitoring studies provide key information to assess ecosystem health. Results of chemical monitoring campaigns can be used to identify the exposure scenarios of regulatory concern. In environmental risk assessment (ERA), measured concentrations of chemicals can be used to model predicted environmental concentrations (PECs). As the PEC is, by definition, a predicted variable, it is highly dependent on the underlying modeling approach from which it is derived. We demonstrate the use of Bayesian distributional regression models to derive PECs by incorporating spatiotemporal conditional variances, and limits of quantification (LOQ) and detection (LOD) as de facto data censoring. Model accuracies increase when incorporating spatiotemporal conditional variances, and the inclusion of LOQ and LOD results in potentially more robust PEC distributions. The methodology is flexible, credibly quantifies uncertainty, and can be adjusted to different scientific and regulatory needs. Posterior sampling allows to express PECs as distributions, which makes this modeling procedure directly compatible with other Bayesian ERA approaches. We recommend the use of Bayesian modeling approaches with chemical monitoring data to make realistic and robust PEC estimations and encourage the scientific debate about the benefits and challenges of Bayesian methodologies in the context of ERA.

Quantifying and predicting ecological and human health risks for binary heavy metal pollution accidents at the watershed scale using Bayesian Networks

The accidental leakage of industrial wastewater containing heavy metals from enterprises poses great risks to resident health, social instability, and ecological safety. During 2005-2018, heavy metal mixed pollution accidents comprised approximately 33% of the major environmental ones in China. A Bayesian Networks-based probabilistic approach is developed to quantitatively predict ecological and human health risks for heavy metal mixed pollution accidents at the watershed scale. To estimate the probability distributions of joint ecological exposure once a heavy metal mixed pollution accident occurs, a Copula-based joint exposure calculation method, comprised of a hydro-dynamic model, emergent heavy metal pollution transport model, and the Copula functions, is embedded. This approach was applied to the risk assessment of acute Cr⁶⁺-Hg²⁺ mixed pollution accidents at 76 electroplating enterprises in 24 risk sub-watersheds of the Dongjiang River downstream watershed. The results indicated that nine sub-watersheds created high ecological risks, while only five created high human health risks. In addition, the ecological and human health risk levels were highest in the tributary (the Xizhijiang River), while the ecological risk was more critical in the river network, and the human health risk was more serious in the mainstream of the Dongjiang River. The quantitative risk assessment provides a substantial support to incident prevention and control, risk management, as well as regulatory decision making for electroplating enterprises.

Bayesian Networks in Environmental Risk Assessment: A Review

Human activities both depend upon and have consequences on the environment. Environmental risk assessment (ERA) is a process of estimating the probability and consequences of the adverse effects of human activities and other stressors on the environment. Bayesian networks (BNs) can synthesize different types of knowledge and explicitly account for the probabilities of different scenarios, therefore offering a useful tool for ERA. Their use in formal ERA practice has not been evaluated, however, despite their increasing popularity in environmental modeling. This paper reviews the use of BNs in ERA based on peer-reviewed publications. Following a systematic mapping protocol, we identified studies in which BNs have been used in an environmental risk context and evaluated the scope, technical aspects, and use of the models and their results. The review shows that BNs have been applied in ERA, particularly in recent years, and that there is room to develop both the model implementation and participatory modeling practices. Based on this review and the authors' experience, we outline general guidelines and development ideas for using BNs in ERA. Integr Environ Assess Manag 2021;17:62-78. © 2020 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

"This Is What We Don't Know": Treating Epistemic Uncertainty in Bayesian Networks for Risk Assessment

Failing to communicate current knowledge limitations, that is, epistemic uncertainty, in environmental risk assessment (ERA) may have severe consequences for decision making. Bayesian networks (BNs) have gained popularity in ERA, primarily because they can combine variables from different models and integrate data and expert judgment. This paper highlights potential gaps in the treatment of uncertainty when using BNs for ERA and proposes a consistent framework (and a set of methods) for treating epistemic uncertainty to help close these gaps. The proposed framework describes the treatment of epistemic uncertainty about the model structure, parameters, expert judgment, data, management scenarios, and the assessment's output. We identify issues related to the differentiation between aleatory and epistemic uncertainty and the importance of communicating both uncertainties associated with the assessment predictions (direct uncertainty) and the strength of knowledge supporting the assessment (indirect uncertainty). Probabilities, intervals, or scenarios are expressions of direct epistemic uncertainty. The type of BN determines the treatment of parameter uncertainty: epistemic, aleatory, or predictive. Epistemic BNs are useful for probabilistic reasoning about states of the world in light of evidence. Aleatory BNs are the most relevant for ERA, but they are not sufficient to treat epistemic uncertainty alone because they do not explicitly express parameter uncertainty. For uncertainty analysis, we recommend embedding an aleatory BN into a model for parameter uncertainty. Bayesian networks do not contain information about uncertainty in the model structure, which requires several models. Statistical models (e.g., hierarchical modeling outside the BNs) are required to consider uncertainties and variability associated with data. We highlight the importance of being open about things one does not know and carefully choosing a method to precisely communicate both direct and indirect uncertainty in ERA. Integr Environ Assess Manag 2021;17:221-232. © 2020 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Increased Use of Bayesian Network Models Has Improved Environmental Risk Assessments

Environmental and ecological risk assessments are defined as the process for evaluating the likelihood that the environment may be impacted as a result of exposure to stressors. Although this definition implies the calculation of probabilities, risk assessments traditionally rely on nonprobabilistic methods such as calculation of a risk quotient. Bayesian network (BN) models are a tool for probabilistic and causal modeling, increasingly used in many fields of environmental science. Bayesian networks are defined as directed acyclic graphs where the causal relationships and the associated uncertainty are quantified in conditional probability tables. Bayesian networks inherently incorporate uncertainty and can integrate a variety of information types, including expert elicitation. During the last 2 decades, there has been a steady increase in reports on BN applications in environmental risk assessment and management. At recent annual meetings of the Society of Environmental Toxicology and Chemistry (SETAC) North America and SETAC Europe, a number of applications of BN models were presented along with new theoretical developments. Likewise, recent meetings of the European Geosciences Union (EGU) have dedicated sessions to Bayesian modeling in relation to water quality. This special series contains 10 articles based on presentations in these sessions, reflecting a range of BN applications to systems, ranging from cells and populations to watersheds and national scale. The articles report on recent progress in many topics, including climate and management scenarios, ecological and socioeconomic endpoints, machine learning, diagnostic inference, and model evaluation. They demonstrate that BNs can be adapted to established conceptual frameworks used to support environmental risk assessments, such as adverse outcome pathways and the relative risk model. The contributions from EGU demonstrate recent advancements in areas such as spatial (geographic information system [GIS]-based) and temporal (dynamic) BN modeling. In conclusion, this special series supports the prediction that increased use of Bayesian network models will improve environmental risk assessments. Integr Environ Assess Manag 2021;17:53-61. © 2020 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Quantification of an Adverse Outcome Pathway Network by Bayesian Regression and Bayesian Network Modeling

The adverse outcome pathway (AOP) framework has gained international recognition as a systematic approach linking mechanistic processes to toxicity endpoints. Nevertheless, successful implementation into risk assessments is still limited by the lack of quantitative AOP models (qAOPs) and assessment of uncertainties. The few published qAOP models so far are typically based on data-demanding systems biology models. Here, we propose a less data-demanding approach for quantification of AOPs and AOP networks, based on regression modeling and Bayesian networks (BNs). We demonstrate this approach with the proposed AOP #245, "Uncoupling of photophosphorylation leading to reduced ATP production associated growth inhibition," using a small experimental data set from exposure of Lemna minor to the pesticide 3,5-dichlorophenol. The AOP-BN reflects the network structure of AOP #245 containing 2 molecular initiating events (MIEs), 3 key events (KEs), and 1 adverse outcome (AO). First, for each dose-response and response-response (KE) relationship, we quantify the causal relationship by Bayesian regression modeling. The regression models correspond to dose-response functions commonly applied in ecotoxicology. Secondly, we apply the fitted regression models with associated uncertainty to simulate 10 000 response values along the predictor gradient. Thirdly, we use the simulated values to parameterize the conditional probability tables of the BN model. The quantified AOP-BN model can be run in several directions: 1) prognostic inference, run forward from the stressor node to predict the AO level; 2) diagnostic inference, run backward from the AO node; and 3) omnidirectionally, run from the intermediate MIEs and/or KEs. Internal validation shows that the AOP-BN can obtain a high accuracy rate, when run is from intermediate nodes and when a low resolution is acceptable for the AO. Although the performance of this AOP-BN is limited by the small data set, our study demonstrates a proof-of-concept: the combined use of Bayesian regression modeling and Bayesian network modeling for quantifying AOPs. Integr Environ Assess Manag 2021;17:147-164. © 2020 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Multi-pathway human exposure risk assessment using Bayesian modeling at the historically largest mercury mining district

The largest mercury (Hg) mining district in the world is located in Almadén (Spain), with well-known environmental impacts in the surrounding ecosystem. However, the impact of mercury on the health of the inhabitants of this area has not been documented accordingly. This study aims to carry out a probabilistic human health risk assessment using Bayesian modeling to estimate the non-carcinogenic risk related to Hg through multiple exposure pathways. Samples of vegetables, wild mushrooms, fish, soil, water, and air were analyzed, and adult residents were randomly surveyed to adjust the risk models to the specific population data. On the one hand, the results for the non-carcinogenic risk based on Hazard Quotient (HQ) showed unacceptable risk levels through ingestion of Hg-contaminated vegetables and fish, with HQ values 20 and 3 times higher, respectively, than the safe exposure threshold of 1 for the 97.5th percentile. On the other hand, ingestion of mushrooms, dermal contact with soil, ingestion of water, dermal contact with water and inhalation of air, were below the safety limit for the 97.5th percentile, and did not represent a risk to the health of residents. In addition, the probabilistic approach was compared with the conservative deterministic approach, and similar results were obtained. This is the first study conducted in Almadén, which clearly reveals the high levels of human health risk to which the population is exposed due to the legacy of two millennia of Hg mining.