skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Dynamic Data-Driven Event Reconstruction for Atmospheric Releases

Abstract

Accidental or terrorist releases of hazardous materials into the atmosphere can impact large populations and cause significant loss of life or property damage. Plume predictions have been shown to be extremely valuable in guiding an effective and timely response. The two greatest sources of uncertainty in the prediction of the consequences of hazardous atmospheric releases result from poorly characterized source terms and lack of knowledge about the state of the atmosphere as reflected in the available meteorological data. In this report, we discuss the development of a new event reconstruction methodology that provides probabilistic source term estimates from field measurement data for both accidental and clandestine releases. Accurate plume dispersion prediction requires the following questions to be answered: What was released? When was it released? How much material was released? Where was it released? We have developed a dynamic data-driven event reconstruction capability which couples data and predictive models through Bayesian inference to obtain a solution to this inverse problem. The solution consists of a probability distribution of unknown source term parameters. For consequence assessment, we then use this probability distribution to construct a ''''composite'' forward plume prediction which accounts for the uncertainties in the source term. Since in mostmore » cases of practical significance it is impossible to find a closed form solution, Bayesian inference is accomplished by utilizing stochastic sampling methods. This approach takes into consideration both measurement and forward model errors and thus incorporates all the sources of uncertainty in the solution to the inverse problem. Stochastic sampling methods have the additional advantage of being suitable for problems characterized by a non-Gaussian distribution of source term parameters and for cases in which the underlying dynamical system is non-linear. We initially developed a Markov Chain Monte Carlo (MCMC) stochastic methodology and demonstrated its effectiveness by reconstructing a wide range of release scenarios, using synthetic as well as real-world data. Data for evaluation of our event reconstruction capability were drawn from the short-range Prairie Grass, Copenhagen, and Joint Urban 2003 field experiments and a continental-scale real-world accidental release in Algeciras, Spain. The method was tested using a variety of forward models, including a Gaussian puff dispersion model INPUFF, the regional-to-continental scale Lagrangian dispersion model LODI (the work-horse real-time operational dispersion model used by the National Atmospheric Release Advisory Center), the empirical urban model UDM, and the building-scale computational computational fluid dynamics code FEM3MP. The robustness of the Bayesian methodology was demonstrated via the use of subsets of the available concentration data and by introducing error into some of the measurements. These tests showed that the Bayesian approach is capable of providing reliable estimates of source characteristics even in cases of limited or significantly corrupted data. For more effective treatment of strongly time-dependent problems, we developed a Sequential Monte Carlo (SMC) approach. To achieve the best performance under a wide range of conditions we combined SMC and MCMC sampling into a hybrid methodology. We compared the effectiveness and advantages of this approach relative to MCMC using a set of synthetic data examples. Our dynamic data-driven event reconstruction capability seamlessly integrates observational data streams with predictive models, in order to provide the best possible estimates of unknown source term parameters, as well as optimal and timely situation analyses consistent with both models and data.« less

Authors:
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
908126
Report Number(s):
UCRL-TR-229417
TRN: US200722%%486
DOE Contract Number:
W-7405-ENG-48
Resource Type:
Technical Report
Country of Publication:
United States
Language:
English
Subject:
54 ENVIRONMENTAL SCIENCES; COMPUTERIZED SIMULATION; DISTRIBUTION; FLUID MECHANICS; HAZARDOUS MATERIALS; LAGRANGIAN FUNCTION; SAMPLING; SOURCE TERMS; ATMOSPHERICS

Citation Formats

Kosovic, B, Belles, R, Chow, F K, Monache, L D, Dyer, K, Glascoe, L, Hanley, W, Johannesson, G, Larsen, S, Loosmore, G, Lundquist, J K, Mirin, A, Neuman, S, Nitao, J, Serban, R, Sugiyama, G, and Aines, R. Dynamic Data-Driven Event Reconstruction for Atmospheric Releases. United States: N. p., 2007. Web. doi:10.2172/908126.
Kosovic, B, Belles, R, Chow, F K, Monache, L D, Dyer, K, Glascoe, L, Hanley, W, Johannesson, G, Larsen, S, Loosmore, G, Lundquist, J K, Mirin, A, Neuman, S, Nitao, J, Serban, R, Sugiyama, G, & Aines, R. Dynamic Data-Driven Event Reconstruction for Atmospheric Releases. United States. doi:10.2172/908126.
Kosovic, B, Belles, R, Chow, F K, Monache, L D, Dyer, K, Glascoe, L, Hanley, W, Johannesson, G, Larsen, S, Loosmore, G, Lundquist, J K, Mirin, A, Neuman, S, Nitao, J, Serban, R, Sugiyama, G, and Aines, R. Thu . "Dynamic Data-Driven Event Reconstruction for Atmospheric Releases". United States. doi:10.2172/908126. https://www.osti.gov/servlets/purl/908126.
@article{osti_908126,
title = {Dynamic Data-Driven Event Reconstruction for Atmospheric Releases},
author = {Kosovic, B and Belles, R and Chow, F K and Monache, L D and Dyer, K and Glascoe, L and Hanley, W and Johannesson, G and Larsen, S and Loosmore, G and Lundquist, J K and Mirin, A and Neuman, S and Nitao, J and Serban, R and Sugiyama, G and Aines, R},
abstractNote = {Accidental or terrorist releases of hazardous materials into the atmosphere can impact large populations and cause significant loss of life or property damage. Plume predictions have been shown to be extremely valuable in guiding an effective and timely response. The two greatest sources of uncertainty in the prediction of the consequences of hazardous atmospheric releases result from poorly characterized source terms and lack of knowledge about the state of the atmosphere as reflected in the available meteorological data. In this report, we discuss the development of a new event reconstruction methodology that provides probabilistic source term estimates from field measurement data for both accidental and clandestine releases. Accurate plume dispersion prediction requires the following questions to be answered: What was released? When was it released? How much material was released? Where was it released? We have developed a dynamic data-driven event reconstruction capability which couples data and predictive models through Bayesian inference to obtain a solution to this inverse problem. The solution consists of a probability distribution of unknown source term parameters. For consequence assessment, we then use this probability distribution to construct a ''''composite'' forward plume prediction which accounts for the uncertainties in the source term. Since in most cases of practical significance it is impossible to find a closed form solution, Bayesian inference is accomplished by utilizing stochastic sampling methods. This approach takes into consideration both measurement and forward model errors and thus incorporates all the sources of uncertainty in the solution to the inverse problem. Stochastic sampling methods have the additional advantage of being suitable for problems characterized by a non-Gaussian distribution of source term parameters and for cases in which the underlying dynamical system is non-linear. We initially developed a Markov Chain Monte Carlo (MCMC) stochastic methodology and demonstrated its effectiveness by reconstructing a wide range of release scenarios, using synthetic as well as real-world data. Data for evaluation of our event reconstruction capability were drawn from the short-range Prairie Grass, Copenhagen, and Joint Urban 2003 field experiments and a continental-scale real-world accidental release in Algeciras, Spain. The method was tested using a variety of forward models, including a Gaussian puff dispersion model INPUFF, the regional-to-continental scale Lagrangian dispersion model LODI (the work-horse real-time operational dispersion model used by the National Atmospheric Release Advisory Center), the empirical urban model UDM, and the building-scale computational computational fluid dynamics code FEM3MP. The robustness of the Bayesian methodology was demonstrated via the use of subsets of the available concentration data and by introducing error into some of the measurements. These tests showed that the Bayesian approach is capable of providing reliable estimates of source characteristics even in cases of limited or significantly corrupted data. For more effective treatment of strongly time-dependent problems, we developed a Sequential Monte Carlo (SMC) approach. To achieve the best performance under a wide range of conditions we combined SMC and MCMC sampling into a hybrid methodology. We compared the effectiveness and advantages of this approach relative to MCMC using a set of synthetic data examples. Our dynamic data-driven event reconstruction capability seamlessly integrates observational data streams with predictive models, in order to provide the best possible estimates of unknown source term parameters, as well as optimal and timely situation analyses consistent with both models and data.},
doi = {10.2172/908126},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Thu Feb 22 00:00:00 EST 2007},
month = {Thu Feb 22 00:00:00 EST 2007}
}

Technical Report:

Save / Share:
  • The role of an event reconstruction capability in a case of an atmospheric release is to characterize the source by answering the critical questions--How much material was released? When? Where? and What are the potential consequences? Accurate estimation of the source term is essential to accurately predict plume dispersion, effectively manage the emergency response, and mitigate consequences in a case of an atmospheric release of hazardous material. We are developing a capability that seamlessly integrates observational data streams with predictive models in order to provide probabilistic estimates of unknown source term parameters consistent with both data and model predictions. Ourmore » approach utilizes Bayesian inference with stochastic sampling using Markov Chain and Sequential Monte Carlo methodology. The inverse dispersion problem is reformulated into a solution based on efficient sampling of an ensemble of predictive simulations, guided by statistical comparisons with data. We are developing a flexible and adaptable data-driven event-reconstruction capability for atmospheric releases that provides (1) quantitative probabilistic estimates of the principal source-term parameters (e.g., the time-varying release rate and location); (2) predictions of increasing fidelity as an event progresses and additional data become available; and (3) analysis tools for sensor network design and uncertainty studies. Our computational framework incorporates multiple stochastic algorithms, operates with a range and variety of atmospheric models, and runs on multiple computer platforms, from workstations to large-scale computing resources. Our final goal is a multi-resolution capability for both real-time operational response and high fidelity multi-scale applications.« less
  • Accidental or terrorist releases of hazardous materials into the atmosphere can impact large populations and cause significant loss of life or property damage. Plume predictions have been shown to be extremely valuable in guiding an effective and timely response. The two greatest sources of uncertainty in the prediction of the consequences of hazardous atmospheric releases result from poorly characterized source terms and lack of knowledge about the state of the atmosphere as reflected in the available meteorological data. We have developed a new event reconstruction methodology that provides probabilistic source term estimates from field measurement data for both accidental andmore » clandestine releases. Accurate plume dispersion prediction requires the following questions to be answered: What was released? When was it released? How much material was released? Where was it released? We have developed a dynamic-data-driven event reconstruction capability that couples data and predictive methods through Bayesian inference to obtain a solution to this inverse problem. The solution consists of a probability distribution of unknown source term parameters. For consequence assessment, we then use this probability distribution to construct a 'composite' forward plume prediction that accounts for the uncertainties in the source term. Since in most cases of practical significance it is impossible to find a closed form solution, Bayesian inference is accomplished by utilizing stochastic sampling methods. This approach takes into consideration both measurement and forward model errors and thus incorporates all the sources of uncertainty in the solution to the inverse problem. Stochastic sampling methods have the additional advantage of being suitable for problems characterized by a non-Gaussian distribution of source term parameters and for cases in which the underlying dynamical system is nonlinear. We initially developed a Markov Chain Monte Carlo (MCMC) stochastic methodology and demonstrated its effectiveness by reconstructing a wide range of release scenarios, using synthetic as well as real-world data. Data for evaluation of our event reconstruction capability were drawn from the short-range Prairie Grass, Copenhagen, and Joint Urban 2003 field experiments and a continental-scale real-world accidental release in Algeciras, Spain. The method was tested using a variety of forward models, including a Gaussian puff dispersion model INPUFF, the regional-to-continental scale Lagrangian dispersion model LODI (the work-horse real-time operational dispersion model used by the National Atmospheric Release Advisory Center), the empirical urban model UDM, and the building-scale computational fluid dynamics code FEM3MP. The robustness of the Bayesian methodology was demonstrated via the use of subsets of the available concentration data and by introducing error into some of the measurements (Fig. 1). These tests showed that the Bayesian approach is capable of providing reliable estimates of source characteristics even in cases of limited or significantly corrupted data. An example of an urban release scenario is shown in Fig. 2. For more effective treatment of strongly time-dependent problems, we developed a Sequential Monte Carlo (SMC) approach. To achieve the best performance under a wide range of conditions we combined SMC and MCMC sampling into a hybrid methodology. We compared the effectiveness and advantages of this approach relative to MCMC using a set of synthetic data examples. We created a modular, scalable computational framework to accommodate the full set of stochastic methodologies (e.g., MCMC, SMC, hybrid stochastic algorithms, 'Green's function', 'reciprocal' methods), as well as a selection of key classes of dispersion models. This design provides a clear separation of stochastic algorithms from predictive models and supports parallelization at both the stochastic algorithm and individual model level. In other words, it supports a parallel stochastic algorithm (e.g., SMC) that invokes parallel forward models. The framework is written in Python and utilizes pyMPI. It invokes forward models either through system calls or as shared objects. Our dynamic-data-driven event reconstruction capability seamlessly integrates observational data streams with predictive models, in order to provide the best possible estimates of unknown source-term parameters, as well as optimal and timely situation analyses consistent with both models and data. This new methodology is shown to be both flexible and robust, adaptable for use with any atmospheric dispersion model, and suitable for use in operational emergency response applications.« less
  • This is a collaborative LDRD Exploratory Research project involving four directorates--Energy & Environment, Engineering, NAI and Computation. The project seeks to answer the following critical questions regarding atmospheric releases--''How much material was released? When? Where? and What are the potential consequences?'' Inaccurate estimation of the source term can lead to gross errors, time delays during a crisis, and even fatalities. We are developing a capability that seamlessly integrates observational data streams with predictive models in order to provide the best possible estimates of unknown source term parameters, as well as optimal and timely situation analyses consistent with both models andmore » data. Our approach utilizes Bayesian inference and stochastic sampling methods (Markov Chain and Sequential Monte Carlo) to reformulate the inverse problem into a solution based on efficient sampling of an ensemble of predictive simulations, guided by statistical comparisons with data. We are developing a flexible and adaptable data-driven event-reconstruction capability for atmospheric releases that provides (1) quantitative probabilistic estimates of the principal source-term parameters (e.g., the time-varying release rate and location); (2) predictions of increasing fidelity as an event progresses and additional data become available; and (3) analysis tools for sensor network design and uncertainty studies. Our computational framework incorporates multiple stochastic algorithms, operates with a range and variety of atmospheric models, and runs on multiple computer platforms, from workstations to large-scale computing resources. Our final goal is a multi-resolution capability for both real-time operational response and high fidelity multi-scale applications.« less
  • For atmospheric releases, event reconstruction answers the critical questions - How much material was released? When? Where? and What are the potential consequences? Inaccurate estimation of the source term can lead to gross errors, time delays during a crisis, and even fatalities. We are developing a capability that seamlessly integrates observational data streams with predictive models in order to provide the best possible estimates of unknown source term parameters, as well as optimal and timely situation analyses consistent with both models and data. Our approach utilizes Bayesian inference and stochastic sampling methods (Markov Chain and Sequential Monte Carlo) to reformulatemore » the inverse problem into a solution based on efficient sampling of an ensemble of predictive simulations, guided by statistical comparisons with data.« less
  • Atmospheric releases of hazardous materials are highly effective means to impact large populations. We propose an atmospheric event reconstruction framework that couples observed data and predictive computer-intensive dispersion models via Bayesian methodology. Due to the complexity of the model framework, a sampling-based approach is taken for posterior inference that combines Markov chain Monte Carlo (MCMC) and sequential Monte Carlo (SMC) strategies.