Expert judgement and uncertainty quantification for climate change
|
|
- Magdalene Poole
- 5 years ago
- Views:
Transcription
1 Expert judgement and uncertainty quantification for climate change 1. Probabilistic inversion as applied to atmospheric dispersion Atmospheric dispersion is mathematically simple, yet involves many features of more complex environmental models. It has been extensively studied in connection with radionuclide dispersion after an accident with nuclear power plants (for references see 1,2 ) and represents a poster child for uncertainty quantification (UQ) with probabilistic inversion (PI) for parametrized models. Consider a simple Gaussian plume dispersion model for concentration C(x,y,z) in point (x,y,z) of a neutrally buoyant contaminant released at height H in a constant wind field with wind speed u in the x direction: Q exp(-y 2 /(2σ 2 y(x))) exp[-(z-h) 2 /(2σ 2 z(x)) + exp[-(z+h) 2 /(2σ 2 z(x)) + mixing layer terms C(x,y,z) = [g/m 3 ] (1) uσ y (x) (2π) σ z (x) (2π) Q = release rate [g/s] σ y (x) = A y x By ; σ z (x) = A z x Bz [m] (2) The term in (z+h) captures reflection from the ground. Similar terms in mixing layer height L capture reflections from the mixing layer, to and fro. Without the parameters H and L, (1) emerges as a solution of the Navier Stokes laws of motion under ideal (anisotropic Fickian) conditions. The standard deviations σ y (x), σ z (x) are given by σ i (x) 2 = 2K i x/u where the eddy diffusivities K i [m 2 /s] are constant, i = y,z. Since each standard deviation grows with the square root of downwind distance x, B y = B z = ½ and the dimensions on both sides of (1) agree. These ideal conditions are known not to apply. The windfield is not constant, there is an inversion layer which reflects the spreading plume back to Earth and ground reflects the plume skyward. There is surface roughness, vertical wind profile, plume meander, Coriolis forces, plume depletion and so on. Fixes are added to (1) in an attempt to account for these features while keeping the data appetite manageable. One such fix involves categorizing the stability of the atmosphere and expressing the diffusion coefficients as functions of downwind distance σ y (x) = A y x By ; σ z (x) = A z x Bz (2) where the parameters A y, B y,a z, B z depend on atmospheric stability. B y, and, B z thereby differ from their theoretical value of ½. Stability is classified by any of a variety of schemes, commonly depending on time of day, insolation, cloud cover and release height. Adding such fixes is sometimes called a "parameterization". (2) is not derived from physical laws. In the best case NATURE CLIMATE CHANGE 1
2 these coefficients are based on regression of experimental results. The meaning of these coefficients is dependent on all the other fixes that have been applied, and no one actually believes these models. They are used because they are judged fit to purpose given data and computing constraints. Independent experts are reluctant to quantify their uncertainty on the diffusion coefficients A y, B y, A z, B z directly, for a variety of reasons: (i) these coefficients' meaning depends on the whole set of parametric fixes (ii) the experts may not subscribe to this model cum fixes, and (iii) the uncertainty on these coefficients is strongly dependent both on each other (if A is large, then B must be small) and on other parameters. To get independent expert input, experts were asked to quantify their uncertainty on measurable quantities such as cross wind dispersion at various downwind distances after a release under stipulated conditions. Focusing the elicitation on measurable quantities has two advantages; first it renders the uncertainty quantification model independent and second, in enables comparison of expert uncertainties with actual measured values. The expert distributions are combined using weights derived from their performance based on comparisons with measured values. Supplementary Table 1 compares the statistical accuracy and informativeness i scores of the participating experts with the scores from the performance weighted and equal weighted combinations of experts, and also with the modelers' uncertainty distributions ("MUSEMET"). Dispersion Statistical accuracy (P-value) Relative Information # Calibration variables Expert Expert Expert Expert Expert Expert Expert Expert Expert Expert Expert Equal weighted Performance weighted MUSEMET Supplementary Table 1 Expert and combination performance scores for dispersion panel (from Cooke 1994) 4. 2 NATURE CLIMATE CHANGE
3 After combining the experts' distributions, something like Figure 1 (main text) for σ y (x) emerges. The dotted curves represent experts' combined uncertainty for σ y at four downwind distances. The five solid lines represent the function σ y (x) for five choices of (A y, B y ). Probabilistic inversion (PI) draws a large sample of values (A y, B y ) from a diffuse distribution which covers the support of the dotted curves. Weights are then assigned to each solid curve such that when the initial distribution of (A y, B y ) is re-sampled using these weights, the dotted distributions are optimally recovered. If the PI problem is feasible, an optimal set of weights, in the sense of minimal departure from the starting distribution, is quickly found. This distribution typically introduces complex dependencies between A y and B y, and is given numerically. Summary of Results Modelers have placed distributions on the coefficients in (2) in an attempt to capture uncertainty in the use of model (1) with (2). Figure S1 illustrates the differences between the modelers' and the independent experts' uncertainties. Eleven experts quantified their uncertainties on the ratio of the centerline concentration to the release rate (a proxy for σ y σ z at four downwind distances 4 ). The ratio of the 95th to the 5th percentiles is depicted, indicating the spread in the uncertainty distribution. Whereas the uncertainty of the experts, and of the equal weighted (eqdm) and performance weighted (dm) combinations of experts' distributions grow considerably at downwind distances of 10km and 30km, this is scarcely the case for the modeler's distribution ("musemet"). Indeed, the modelers' uncertainty is generally smaller than that of any of the experts. Figure S1: Ratio of 95/5 percentiles for 11 experts, optimized performance based combinations (dm), equal weight combinations (eqdm) and the modelers' in house uncertainty assessment, for the centerline concentration ratio in stability class B at downwind distances 1937m, 4250m, 10km, 30km (from Cooke 1994) 4. NATURE CLIMATE CHANGE 3
4 If a probabilistic inverse exists, it is generally not unique and we seek a preferred inverse. If no inverse exists, the problem is infeasible, and we seek a distribution which minimizes some appropriate error function. The mathematics of PI are embedded in algorithms for iteratively finding the weights for weighted re-sampling the original distribution. The principal tool, the iterative proportional fitting algorithm (IPF), was introduced by 5 and re discovered by 6 and many others. In case of feasibility, reference 7 showed that IPF converges to the solution which is the minimally informative distribution with respect to the starting distribution satisfying the constraints. In case of infeasibility, little is known about the behavior of IPF. However, variations of this algorithm will always converge to a distribution which minimizes an appropriate error functional 8. For more information on probabilistic inversion see references 3, 9, and 10. It may happen that no distribution over model parameters adequately captures the distributions over the target variables and this is an important element in model criticism. 2. Spreadsheet Example for Ice Sheets PI (Supplementary Table 2) BA13 11 performed a structured expert elicitation to assess the uncertainty in sea level rise in 2100, and found 5%, 50% and 95% values which differ significantly from the values emerging from the simplified Process Based Model in L Here, we "invert the PBM model " at the BA13 11 distributions, to obtain distributions over D(EAIS), D(WAIS), SMB(EAIS) and SMB(WAIS) such that when pushed through the PBM, the resulting 5%, 50% and 95% values for SLR(EAIS) and SLR(WAIS) agree with those in BA The operation of inverting a model at a set of joint distributions is termed "probabilistic inversion" (PI), and the primary tool is the iterative proportional fitting algorithm (IPF). A formal definition of a probabilistic inversion is as follows: Let X R M, and let functions H i : X R, i = 1 n be given. Suppose distributions F i are assigned to H i, i = 1 n. We seek a distribution over X such that the functions H i take the distributions F i under this distribution. Such a distribution is a probabilistic inverse of H 1 H n at F 1 F n. We may assign a set C i of distributions to H i, for example the set of distributions satisfying some quantile constraints. A probabilistic inverse at C 1,..C n is defined analoguously. For more detail, see Kurowicka and Cooke (2006) 3. The Ice sheet PI Example (Supplementary Table 2) gives detailed calculations of the probabilistic inversion using IPF applied to Antarctic Sea Level rise described in the paper. Given an initial sample from a prior distribution over D(EAIS), D(WAIS), SMB(EAIS) and SMB(WAIS), IPF iteratively computes weights for re-weighting the original sample. When the original sample is re-sampled using the IPF weights, the resulting distribution complies with the BA13 11 percentile values. All calculations are spreadsheet functions so the reader can follow all details. In the real application, 100,000 samples were drawn from the prior distribution for discharge, surface mass balance and sea level rise for the East and West Antarctica Ice Sheets. In this toy calculation 4 NATURE CLIMATE CHANGE
5 only 30 samples are used for illustrative purposes. Abstract formulations of the IPF and similar algorithms can be found in the literature 3. i. Experts assessed their 5, 50 and 95th percentiles of their subjective uncertainty for all variables, including calibration variables. Statistical accuracy is measured as the p-value at which the hypothesis that the expert's probability statements were accurate would be falsely rejected, assuming the calibration variables are independently drawn from a distribution complying with the expert's quantile assessments. Information is measured as Shannon relative information with respect to a background measure determined by all experts' assessed quantiles. References 1. Radiation Protection Dosimetry Expert Judgement and Accident consequence Uncertainty Analysis, (2000) special issue L.H.J. Goossens and G.N. Kelly (eds) Vol 90, No,3, 2. Harper FT, Hora SC, Young ML, Miller LA, Lui CH, McKay MD, Helton JC, Goossens, LHJ, Cooke RM Pasler-Sauer J, Kraan B, Jones JA. (1995) Probabilistic accident consequence uncertainty study: Dispersion and deposition uncertainty assessment. Prepared for U.S. Nuclear Regulatory Commission and Commission of European Communities NUREG/CR-6244, EUR EN, SAND , Washington/USA, and Brussels-Luxembourg, November 1994, published January Volume I: Main report, Volume II: Appendices A and B, Volume III: Appendices C, D, E, F, G, H; 3. Kurowicka, D., and Cooke, Roger M. (2006) Uncertainty Analysis and High Dimensional Dependence Modeling Wiley, New York. ISBN-13: , ISBN-10: Cooke, R.M. (1994) Uncertainty in dispersion and deposition in accident consequence modeling assessed with performance-based expert judgment, Reliability Engineering & System Safety, Volume 45, Issues 1 2, 1994, Pages 35-46, ISSN , ( 5. Kruithof, J. (1937). Telefoonverkeersrekening. De Ingenieur, 52(8):15-25, Deming, E.W. and Stephan, F.D.(1940). On a least squares adjustment of a sampled frequency table when the expected marginal totals are known. Annals of Mathematical Statistics, 11(4): , Csiszar I.(1975) I-divergence geometry of probability distributions and minimization problems. The Annals of Probability, 3(1):146(158, Matus, F. (2007) On iterated averages of I-projections, Statistiek und Informatik, Universit at Bielefeld, Bielefeld, Germany, matus@utia.cas.cz. 9. Kraan, B.C.P and Bedford. T.J (2005) Probabilistic inversion of expert judgements in the quantification of model uncertainty. Management Science, 51(6): , Du C., Kurowicka D., Cooke R.M.(2006) Techniques for generic probabilistic inversion, Comp. Stat. & Data Analysis (50) 2006, Bamber, J. L. & Aspinall, W. P. An expert judgement assessment of future sea level rise from the ice sheets. Nature Climate Change 3, (2013). NATURE CLIMATE CHANGE 5
6 12. Little, C. M., M. Oppenheimer, and N. M. Urban (2013b), Upper bounds on twenty-firstcentury Antarctic ice loss assessed using a probabilistic framework, Nature Climate Change, 3, , doi: /nclimate TU Delft TU Delft TU Delft TU Delft (/ TU Delft TU Delft TU Delft TU Delft TU Delft TU Delft TNO TNO TNO NO TNO TNO Virtual experts TUD OP.DM TNO OP.DM" TUD + TNO OP.DM TUD EQ.DM TNO EQ.DM" TUD + TNO EQ.DM MUSEMET b NATURE CLIMATE CHANGE
Techniques for Generic Probabilistic Inversion
Techniques for Generic Probabilistic Inversion D.Kurowicka 1, R.M. Cooke 2 1 Department of Mathematics, Delft University of Technology, Delft, The Netherlands; d.kurowicka@its.tudelft.nl 2 Department of
More informationUncertainty in Compartmental Models for Hazardous Materials - A Case Study
Appearing in Journal of Hazardous Materials 71 (2000) 253-268 Uncertainty in Compartmental Models for Hazardous Materials - A Case Study Roger Cooke Bernd Kraan Delft University of Technology 1 Abstract
More informationProcessing Expert Judgements in Accident Consequence Modelling
Appearing in Radiation Protection Dosimetry (Expert Judgement and Accident Consequence Uncertainty Analysis; special issue) Vol.90 No.3 2000 311-315. Processing Expert Judgements in Accident Consequence
More informationTechniques for Generic Probabilistic Inversion
Techniques for Generic Probabilistic Inversion C. Du 1, D. Kurowicka 2, R.M.Cooke 3 ABSTRACT Probabilistic inversion problems are defined, existing algorithms for solving such problems are discussed and
More informationProbabilistic Inversion for Chicken Processing Lines
Probabilistic Inversion for Chicken Processing Lines Roger M. Cooke a Maarten Nauta b, Arie H. Havelaar b, Ine van der Fels. b a Department of Mathematics Delft University of Technology Delft, The Netherlands
More informationEXPERT AGGREGATION WITH DEPENDENCE
EXPERT AGGREGATION WITH DEPENDENCE M. J. Kallen, R.M. Cooke 2 Department of Mathematics, Delft University of Technology, Delft, The Netherlands 2 Department of Mathematics, Delft University of Technology,
More informationVines in Overview Invited Paper Third Brazilian Conference on Statistical Modelling in Insurance and Finance Maresias, March 25-30, 2007
Vines in Overview Invited Paper Third Brazilian Conference on Statistical Modelling in Insurance and Finance Maresias, March 25-30, 2007 R.M. Cooke, O. Morales, D. Kurowicka Resources for the Future and
More informationEliciting Conditional and Unconditional Rank Correlations from Conditional Probabilities.
Eliciting Conditional and Unconditional Rank Correlations from Conditional Probabilities. O.Morales, D.Kurowicka Delft Institute of Applied Mathematics Delft University of Technology The Netherlands A.Roelen
More informationDiscrete Choice with Probabilistic Inversion: Application to energy policy choice and wiring failure
Discrete Choice with Probabilistic Inversion: Application to energy policy choice and wiring failure Roger Cooke Resources for the Future Dept. Math. TU Delft Jolanta Misiewicz Dept. Math. T.U.Delft Dept
More informationModeling and Validating Stakeholder Preferences with Probabilistic Inversion
Noname manuscript No. (will be inserted by the editor) Modeling and Validating Stakeholder Preferences with Probabilistic Inversion R.E.J Neslo, R.M. Cooke 6-9-2010 Abstract Despite the vast number of
More informationGraphical Methods for Uncertainty and Sensitivity Analysis
Graphical Methods for Uncertainty and Sensitivity Analysis Roger M. Cooke Delft University of Technology, The Netherlands Jan M. van Noortwijk HKV Consultants, Lelystad, The Netherlands Contents Introduction
More informationUSING METEOROLOGICAL ENSEMBLES FOR ATMOSPHERIC DISPERSION MODELLING OF THE FUKUSHIMA NUCLEAR ACCIDENT
USING METEOROLOGICAL ENSEMBLES FOR ATMOSPHERIC DISPERSION MODELLING OF THE FUKUSHIMA NUCLEAR ACCIDENT 17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory
More informationEnsembles, Uncertainty and Climate Projections. Chris Brierley (Room 117)
Ensembles, Uncertainty and Climate Projections Chris Brierley (Room 117) Definitions Ensemble: group of model simulations Uncertainty: doubt and ambiguity about future conditions Climate Projection: modeled
More informationStudy of the Effect of Weather, Topography and Radionuclide on the TEDE in a Fire Scenario Involving a Dispersion of a Plume in the Atmosphere
Middle-East Journal of Scientific Research 18 (8): 1192-1198, 2013 ISSN 1990-9233 IDOSI Publications, 2013 DOI: 10.5829/idosi.mejsr.2013.18.8.12405 Study of the Effect of Weather, Topography and Radionuclide
More informationNUCLEAR EMERGENCY RESPONSE ACTIVITIES COORDINATION GROUP Original: ENGLISH EC-JRC ENSEMBLE APPROACH TO LONG RANGE ATMOSPHERIC DISPERSION
WORLD METEOROLOGICAL ORGANIZATION COMMISSION FOR BASIC SYSTEMS CBS/CG-NERA/Doc. (number) (26.04.2006) OPAG on DPFS NUCLEAR EMERGENCY RESPONSE ACTIVITIES COORDINATION GROUP Original: ENGLISH VIENNA, AUSTRIA,
More informationEliciting conditional and unconditional rank correlations from conditional probabilities
Reliability Engineering and System Safety 93 (2008) 699 710 www.elsevier.com/locate/ress Eliciting conditional and unconditional rank correlations from conditional probabilities O. Morales a,, D. Kurowicka
More informationExpert Judgment in the Uncertainty Analysis of Dike Ring Failure Frequency
Expert Judgment in the Uncertainty Analysis of Dike Ring Failure Frequency Roger M. Cooke Department of Mathematics Delft University of Technology Delft, The Netherlands Karen A. Slijkhuis Ministerie Verkeer
More informationCFD calculations of the test 2-4 experiments. Author: G. de With
CFD calculations of the test 2-4 experiments Author: G. de With 34. Model setup and boundary conditions Dimensions CFD model: x=1000m / y=100m / z=2000m. CFD Model: Transient simulation, with steady-state
More informationProbabilistic Inversion in Priority Setting of Emerging Zoonoses
Risk Analysis, Vol. 30, No. 5, 2010 DOI: 10.1111/j.1539-6924.2010.01378.x Probabilistic Inversion in Priority Setting of Emerging Zoonoses Dorota Kurowicka, 1, Catalin Bucura, 1 Roger Cooke, 1,2 and Arie
More informationTHE VINE COPULA METHOD FOR REPRESENTING HIGH DIMENSIONAL DEPENDENT DISTRIBUTIONS: APPLICATION TO CONTINUOUS BELIEF NETS
Proceedings of the 00 Winter Simulation Conference E. Yücesan, C.-H. Chen, J. L. Snowdon, and J. M. Charnes, eds. THE VINE COPULA METHOD FOR REPRESENTING HIGH DIMENSIONAL DEPENDENT DISTRIBUTIONS: APPLICATION
More information9/12/17. Types of learning. Modeling data. Supervised learning: Classification. Supervised learning: Regression. Unsupervised learning: Clustering
Types of learning Modeling data Supervised: we know input and targets Goal is to learn a model that, given input data, accurately predicts target data Unsupervised: we know the input only and want to make
More informationGamma process model for time-dependent structural reliability analysis
Gamma process model for time-dependent structural reliability analysis M.D. Pandey Department of Civil Engineering, University of Waterloo, Waterloo, Ontario, Canada J.M. van Noortwijk HKV Consultants,
More informationExperimental results about the assessments of conditional rank correlations by experts: Example with air pollution estimates
Safety, Reliability and Risk Analysis: Beyond the Horizon Steenbergen et al. (Eds) 24 Taylor & Francis Group, London, ISBN 978--38-23-7 Experimental results about the assessments of conditional rank correlations
More informationFuture sea level rise through 2100 and beyond
Future sea level rise through 2100 and beyond Emma Stone Uncertain World Summit (Tuesday 20 October 2015) 1 Why are we interested in future sea level rise? How certain are we about the worst case scenario
More informationPATTERN RECOGNITION AND MACHINE LEARNING
PATTERN RECOGNITION AND MACHINE LEARNING Chapter 1. Introduction Shuai Huang April 21, 2014 Outline 1 What is Machine Learning? 2 Curve Fitting 3 Probability Theory 4 Model Selection 5 The curse of dimensionality
More informationMODELING FOR ENVIRONMENTAL RADIATION DOSE RECONSTRUCTION. Bruce Napier 23 May 2011
MODELING FOR ENVIRONMENTAL RADIATION DOSE RECONSTRUCTION Bruce Napier 23 May 2011 1 Topics NCRP Report No. 163 Atmospheric dispersion modeling Reconstruction of dose from releases of iodines and noble
More informationTitle: Quantifying Response Time in Carbon Monoxide Monitors to Produce High-
Title: Quantifying Response Time in Carbon Monoxide Monitors to Produce High- Resolution Measurements Authors: D.B. Young, R.M. Jones Division of Environmental and Occupational Health Sciences School of
More informationModule 01 Lecture - 06 Pollution modeling I
Health, Safety and Environmental Management in Offshore and Petroleum Engineering Prof. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institution of Technology, Madras Module 01 Lecture
More informationAnalysis on the formation process of high dose rate zone in the northwest direction of the Fukushima Daiichi nuclear power plant
Analysis on the formation process of high dose rate zone in the northwest direction of the Fukushima Daiichi nuclear power plant Jun. 13, 2011 Research group of Japan Atomic Energy Agency (JAEA) has analyzed
More informationWeighted space-filling designs
Weighted space-filling designs Dave Woods D.Woods@southampton.ac.uk www.southampton.ac.uk/ davew Joint work with Veronica Bowman (Dstl) Supported by the Defence Threat Reduction Agency 1 / 34 Outline Motivation
More informationStochastic Hydrology. a) Data Mining for Evolution of Association Rules for Droughts and Floods in India using Climate Inputs
Stochastic Hydrology a) Data Mining for Evolution of Association Rules for Droughts and Floods in India using Climate Inputs An accurate prediction of extreme rainfall events can significantly aid in policy
More informationLecture : Probabilistic Machine Learning
Lecture : Probabilistic Machine Learning Riashat Islam Reasoning and Learning Lab McGill University September 11, 2018 ML : Many Methods with Many Links Modelling Views of Machine Learning Machine Learning
More informationNaka-Gun, Ibaraki, , Japan
Examination of Atmospheric Dispersion Model s Performance - Comparison with the Monitoring Data under the Normal Operation of the Tokai Reprocessing Plant - M. Takeyasu 1, M. Nakano 1, N. Miyagawa 1, M.
More informationExpert Judgment study for Placement Ladder Bowtie
Expert Judgment study for Placement Ladder Bowtie D. Kurowicka, R.Cooke, L. Goossens, B. Ale Delft University of Technology, The Netherlands ABSTRACT: Statistics of occupational accidents motivated a large
More informationA VALIDATION EXERCISE ON THE SAFE-AIR VIEW SOFTWARE. Joint Research Centre NDFM Ispra, Italy 2
A VALIDATION EXERCISE ON THE SAFE-AIR VIEW SOFTWARE F. D Alberti 1, F. d Amati 1, E. Canepa 2, G. Triacchini 3 1 Joint Research Centre NDFM Ispra, Italy 2 CNR INFM CNISM Department of Physics, University
More informationBasic Verification Concepts
Basic Verification Concepts Barbara Brown National Center for Atmospheric Research Boulder Colorado USA bgb@ucar.edu Basic concepts - outline What is verification? Why verify? Identifying verification
More informationDepartment of Meteorology University of Nairobi. Laboratory Manual. Micrometeorology and Air pollution SMR 407. Prof. Nzioka John Muthama
Department of Meteorology University of Nairobi Laboratory Manual Micrometeorology and Air pollution SMR 407 Prof. Nioka John Muthama Signature Date December 04 Version Lab : Introduction to the operations
More informationDispersion for point sources CE 524 February
Dispersion for point sources CE 524 February 2011 1 Concentration Air pollution law in most industrial countries based on concentration of contaminants NAAQS in US Need method dto predict concentrations
More informationMeteorological Data Collection, X/Q and D/Q, Critical Receptors
Meteorological Data Collection, X/Q and D/Q, Critical Receptors Ken Sejkora Entergy Nuclear Northeast Pilgrim Station Presented at the 23rd Annual RETS-REMP Workshop Training Session Westminster, CO /
More information17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary
17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary CFD MODELLING OF DISPERSION IN NEUTRAL AND STABLE ATMOSPHERIC
More informationarxiv: v1 [physics.flu-dyn] 18 Mar 2018
APS Persistent incomplete mixing in reactive flows Alexandre M. Tartakovsky 1 and David Barajas-Solano 1 1 Pacific Northwest National Laboratory, Richland, WA 99352, USA arxiv:1803.06693v1 [physics.flu-dyn]
More information(Directions for Excel Mac: 2011) Most of the global average warming over the past 50 years is very likely due to anthropogenic GHG increases
(Directions for Excel Mac: 2011) Most of the global average warming over the past 50 years is very likely due to anthropogenic GHG increases How does the IPCC know whether the statement about global warming
More informationEnsemble Data Assimilation and Uncertainty Quantification
Ensemble Data Assimilation and Uncertainty Quantification Jeff Anderson National Center for Atmospheric Research pg 1 What is Data Assimilation? Observations combined with a Model forecast + to produce
More informationINTER-COMPARISON AND VALIDATION OF RANS AND LES COMPUTATIONAL APPROACHES FOR ATMOSPHERIC DISPERSION AROUND A CUBIC OBSTACLE. Resources, Kozani, Greece
INTER-COMPARISON AND VALIDATION OF AND LES COMPUTATIONAL APPROACHES FOR ATMOSPHERIC DISPERSION AROUND A CUBIC OBSTACLE S. Andronopoulos 1, D.G.E. Grigoriadis 1, I. Mavroidis 2, R.F. Griffiths 3 and J.G.
More informationSRI Briefing Note Series No.8 Communicating uncertainty in seasonal and interannual climate forecasts in Europe: organisational needs and preferences
ISSN 2056-8843 Sustainability Research Institute SCHOOL OF EARTH AND ENVIRONMENT SRI Briefing Note Series No.8 Communicating uncertainty in seasonal and interannual climate forecasts in Europe: organisational
More informationFINNISH EXPERIMENTS ON LEVEL 3 PRA. VTT Technical Research Centre of Finland Ltd.: P.O. Box 1000, Espoo, Finland, 02044, and
FINNISH EXPERIMENTS ON LEVEL 3 PRA Tero Tyrväinen 1, Ilkka Karanta 1, Jukka Rossi 1 1 VTT Technical Research Centre of Finland Ltd.: P.O. Box 1000, Espoo, Finland, 02044, and tero.tyrvainen@vtt.fi Level
More informationStochastic methods for representing atmospheric model uncertainties in ECMWF's IFS model
Stochastic methods for representing atmospheric model uncertainties in ECMWF's IFS model Sarah-Jane Lock Model Uncertainty, Research Department, ECMWF With thanks to Martin Leutbecher, Simon Lang, Pirkka
More informationModelling and data assimilation of hazardous volcanic ash plumes in the chemical-transport model MOCAGE
Modelling and data assimilation of hazardous volcanic ash plumes in the chemical-transport model MOCAGE Bojan Sic, Laaziz El Amraoui, Matthieu Plu CNRM/Météo-France 2 Introduction Model MOCAGE of Météo-France
More informationSupplementary Material for: Scientific Uncertainty and Climate Change: Part I. Uncertainty and Unabated Emissions
Noname manuscript No. (will be inserted by the editor) Supplementary Material for: Scientific Uncertainty and Climate Change: Part I. Uncertainty and Unabated Emissions Stephan Lewandowsky James S. Risbey
More informationAir Pollution Meteorology
Air Pollution Meteorology Government Pilots Utilities Public Farmers Severe Weather Storm / Hurricane Frost / Freeze Significant Weather Fog / Haze / Cloud Precipitation High Resolution Weather & Dispersion
More informationTheory and Practice of Data Assimilation in Ocean Modeling
Theory and Practice of Data Assimilation in Ocean Modeling Robert N. Miller College of Oceanic and Atmospheric Sciences Oregon State University Oceanography Admin. Bldg. 104 Corvallis, OR 97331-5503 Phone:
More informationThis document is a preview generated by EVS
INTERNATIONAL STANDARD ISO 13168 First edition 2015-07-01 Water quality Simultaneous determination of tritium and carbon 14 activities Test method using liquid scintillation counting Qualité de l eau Détermination
More informationUncertainty Management and Quantification in Industrial Analysis and Design
Uncertainty Management and Quantification in Industrial Analysis and Design www.numeca.com Charles Hirsch Professor, em. Vrije Universiteit Brussel President, NUMECA International The Role of Uncertainties
More informationDropping Ice Shelves onto an Ocean Model and Moving Grounding Lines. Robert Hallberg NOAA / GFDL
Dropping Ice Shelves onto an Ocean Model and Moving Grounding Lines Robert Hallberg NOAA / GFDL Projected Global Mean Sea Level Rise Sources of uncertainty in 2100 global mean sea level projections: Forcing
More informationGaussian Processes for Computer Experiments
Gaussian Processes for Computer Experiments Jeremy Oakley School of Mathematics and Statistics, University of Sheffield www.jeremy-oakley.staff.shef.ac.uk 1 / 43 Computer models Computer model represented
More informationLecture 7. Testing for Poisson cdf - Poisson regression - Random points in space 1
Lecture 7. Testing for Poisson cdf - Poisson regression - Random points in space 1 Igor Rychlik Chalmers Department of Mathematical Sciences Probability, Statistics and Risk, MVE300 Chalmers April 2010
More informationDover- Sherborn High School Mathematics Curriculum Probability and Statistics
Mathematics Curriculum A. DESCRIPTION This is a full year courses designed to introduce students to the basic elements of statistics and probability. Emphasis is placed on understanding terminology and
More informationAnswer keys for Assignment 10: Measurement of study variables (The correct answer is underlined in bold text)
Answer keys for Assignment 10: Measurement of study variables (The correct answer is underlined in bold text) 1. A quick and easy indicator of dispersion is a. Arithmetic mean b. Variance c. Standard deviation
More informationMultiscale integral analysis of a HT leakage in a fusion nuclear power plant
Journal of Physics: Conference Series PAPER OPEN ACCESS Multiscale integral analysis of a HT leakage in a fusion nuclear power plant To cite this article: M Velarde et al 2016 J. Phys.: Conf. Ser. 717
More informationComparison between test field data and Gaussian plume model
Environmental Modelling for RAdiation Safety II Working group 9 Comparison between test field data and Gaussian plume model Laura Urso Helmholtz Zentrum München Institut für Strahlenschutz AG Radioecological
More informationCUSUM(t) D data. Supplementary Figure 1: Examples of changes in linear trend and CUSUM. (a) The existence of
Supplementary Figures a b a(t) a(t) c t d t d(t) d(t) t t e f CUSUM(t) D data t CUSUM(t) D data t Supplementary Figure 1: Examples of changes in linear trend and CUSUM. (a) The existence of an abrupt jump
More informationA New Reliability Allocation Method Based on FTA and AHP for Nuclear Power Plant!
A New Reliability Allocation Method Based on FTA and AHP for Nuclear Power Plant! Presented by Rongxiang Hu Contributed by FDS Team Institute of Nuclear Energy Safety Technology (INEST) Chinese Academy
More informationSea level over glacial cycles.
Sea level over glacial cycles. T1: 1 m/century With short periods of with rates up to 4m/century LIG Sea level higher than today (+9m) Siddall et al. 2003 Very stable sea level in human history. Image
More informationHow to predict the probability of a major nuclear accident after Fukushima Da
How to predict the probability of a major nuclear accident after Fukushima Dai-ichi? CERNA Mines ParisTech March 14, 2012 1 2 Frequentist approach Bayesian approach Issues 3 Allowing safety progress Dealing
More information14: Correlation. Introduction Scatter Plot The Correlational Coefficient Hypothesis Test Assumptions An Additional Example
14: Correlation Introduction Scatter Plot The Correlational Coefficient Hypothesis Test Assumptions An Additional Example Introduction Correlation quantifies the extent to which two quantitative variables,
More informationImproved Consequence Modelling for Hypothetical Accidents in Fusion Power Plants
Improved Consequence Modelling for Hypothetical Accidents in Fusion Power Plants W. E. Han Euratom/UKAEA Fusion Association, Culham Science Centre, Abingdon, Oxfordshire, OX14 3DB, Tel.: +44 135 588 e-mail:
More informationMeteorology and Dispersion Forecast in Nuclear Emergency in Argentina
Meteorology and Dispersion Forecast in Nuclear Emergency in Argentina Kunst, J.J.; Boutet, L.I.; Jordán, O.D.; Hernández, D.G.; Guichandut, M.E. and Chiappesoni, H. Presentado en: 12 th International Congress
More informationPart I: Dry Convection
Turbulent dispersion and chemical transformation in the atmospheric boundary layer: Part I: Dry Convection DISPERSION Thanks: Alessandro Dosio Jordi Vilà-Guerau de Arellano WA G E N I N G E N U N I VE
More informationResearch Collection. Basics of structural reliability and links with structural design codes FBH Herbsttagung November 22nd, 2013.
Research Collection Presentation Basics of structural reliability and links with structural design codes FBH Herbsttagung November 22nd, 2013 Author(s): Sudret, Bruno Publication Date: 2013 Permanent Link:
More information17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary
17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary USING METEOROLOGICAL ENSEMBLES FOR ATMOSPHERIC DISPERSION
More informationAnalytical Model for Dispersion of Rocket Exhaust Source Size Impact Assessment
American Journal of Environmental Engineering 28, 8(4): 35-39 DOI:.5923/j.ajee.2884.8 Analytical Model for Dispersion of Rocket Exhaust Source Size Impact Assessment Bruno K. Bainy,*, Bardo E. J. Bodmann
More informationBasic Verification Concepts
Basic Verification Concepts Barbara Brown National Center for Atmospheric Research Boulder Colorado USA bgb@ucar.edu May 2017 Berlin, Germany Basic concepts - outline What is verification? Why verify?
More informationCHOOSING THE RIGHT SAMPLING TECHNIQUE FOR YOUR RESEARCH. Awanis Ku Ishak, PhD SBM
CHOOSING THE RIGHT SAMPLING TECHNIQUE FOR YOUR RESEARCH Awanis Ku Ishak, PhD SBM Sampling The process of selecting a number of individuals for a study in such a way that the individuals represent the larger
More informationUniversity of Bristol - Explore Bristol Research
Pausata, F., LeGrande, A., & Roberts, W. H. G. (2016). How Will Sea Ice Loss Affect the Greenland Ice Sheet?: On the Puzzling Features of Greenland Ice-Core Isotopic Composition; Copenhagen, Denmark, 26
More informationOn the lack of information in hydraulic engineering models
On the lack of information in hydraulic engineering models Karen A.H. Slijkhuis Risk Analysis Department, Ministry of Public Works and Water Management, Utrecht, the Netherlands Pieter H.A.J.M. van Gelder
More informationProbabilities, Uncertainties & Units Used to Quantify Climate Change
Name Class Probabilities, Uncertainties & Units Used to Quantify Climate Change Most of the global average warming over the past 50 years is very likely due to anthropogenic GHG increases How does the
More informationRECORD OF REVISIONS. Page 2 of 17 GEO. DCPP.TR.14.06, Rev. 0
Page 2 of 17 RECORD OF REVISIONS Rev. No. Reason for Revision Revision Date 0 Initial Report - this work is being tracked under Notification SAPN 50638425-1 8/6/2014 Page 3 of 17 TABLE OF CONTENTS Page
More informationR. Hallberg, A. Adcroft, J. P. Dunne, J. P. Krasting and R. J. Stouffer NOAA/GFDL & Princeton University
Sensitivity of 21st Century Steric Sea Level Rise to Ocean Model Formulation R. Hallberg, A. Adcroft, J. P. Dunne, J. P. Krasting and R. J. Stouffer NOAA/GFDL & Princeton University Hallberg, R., A. Adcroft,
More informationGlobal Climate Change
Why? Global Climate Change What evidence have scientists collected to support global climate change? Global climate change is a topic that is frequently discussed but often misunderstood due to the complexities
More informationAnalysis of Covariance. The following example illustrates a case where the covariate is affected by the treatments.
Analysis of Covariance In some experiments, the experimental units (subjects) are nonhomogeneous or there is variation in the experimental conditions that are not due to the treatments. For example, a
More information2. REGIONAL DISPERSION
Real-time Transport and Dispersion from Illinois Nuclear Power Plants Thomas E. Bellinger, CCM Illinois Emergency Management Agency Springfield, Illinois 1. INTRODUCTION Meteorological data routinely used
More information5. Advection and Diffusion of an Instantaneous, Point Source
1 5. Advection and Diffusion of an Instantaneous, Point Source In this chapter consider the combined transport by advection and diffusion for an instantaneous point release. We neglect source and sink
More informationINDEX. (The index refers to the continuous pagination)
(The index refers to the continuous pagination) Accuracy in physical models methods for assessing overall assessment acquisition of information acrylonitrile hazards polymerisation toxic effects toxic
More informationProblem Set 2. MAS 622J/1.126J: Pattern Recognition and Analysis. Due: 5:00 p.m. on September 30
Problem Set 2 MAS 622J/1.126J: Pattern Recognition and Analysis Due: 5:00 p.m. on September 30 [Note: All instructions to plot data or write a program should be carried out using Matlab. In order to maintain
More informationThe Canadian approach to ensemble prediction
The Canadian approach to ensemble prediction ECMWF 2017 Annual seminar: Ensemble prediction : past, present and future. Pieter Houtekamer Montreal, Canada Overview. The Canadian approach. What are the
More informationA METHOD FOR EVALUATING HAZARDS OF LOW VOLATILITY TOXIC LIQUIDS
A METHOD FOR EVALUANG HAZARDS OF LOW VOLALITY TOXIC LIQUIDS Timothy A. Melton and David W. Johnson Presented At 1998 Plant Safety Symposium Houston, Texas October 26-27, 1998 Presented By Quest Consultants
More informationAppendix 1: UK climate projections
Appendix 1: UK climate projections The UK Climate Projections 2009 provide the most up-to-date estimates of how the climate may change over the next 100 years. They are an invaluable source of information
More informationLower Bound of Optimization for the Public Considering Dose Distribution of Radiation due to Natural Background Radiation
Lower Bound of Optimization for the Public Considering Dose Distribution of Radiation due to Natural Background Radiation Takatoshi Hattori a* a Central Research Institute of Electric Power Industry, Radiation
More informationModeling Radiological Consequences of Sever Accidents in BWRs: Review of Models Development, Verification and Validation
Modeling Radiological Consequences of Sever Accidents in BWRs: Review of Models Development, Verification and Validation Mohamed GAHEEN Department of Research on Nuclear Safety and Radiological Emergencies,
More informationIMPACT OF WEATHER CHANGES ON TVA NUCLEAR PLANT CHI/Q (χ/q) Kenneth G. Wastrack Doyle E. Pittman Jennifer M. Call Tennessee Valley Authority
IMPACT OF WEATHER CHANGES ON TVA NUCLEAR PLANT CHI/Q (χ/q) Kenneth G. Wastrack Doyle E. Pittman Jennifer M. Call Tennessee Valley Authority The TVA nuclear plants, like most others in the United States,
More informationSIMULATION OF KAMENNA-EXPERIMENTS TEST 1 4 WITH THE DECISION SUPPORT MODEL LASAIR
SIMULATION OF KAMENNA-EXPERIMENTS TEST 1 4 WITH THE DECISION SUPPORT MODEL LASAIR IAEA, EMRAS II Working Group 9 Urban Areas Sevilla, Spain, 08. - 10.06. 2010 Federal Office for Radiation Protection Ingolstaedter
More informationVERFICATION OF OCEAN WAVE ENSEMBLE FORECAST AT NCEP 1. Degui Cao, H.S. Chen and Hendrik Tolman
VERFICATION OF OCEAN WAVE ENSEMBLE FORECAST AT NCEP Degui Cao, H.S. Chen and Hendrik Tolman NOAA /National Centers for Environmental Prediction Environmental Modeling Center Marine Modeling and Analysis
More information1990 Intergovernmental Panel on Climate Change Impacts Assessment
1990 Intergovernmental Panel on Climate Change Impacts Assessment Although the variability of weather and associated shifts in the frequency and magnitude of climate events were not available from the
More informationWe Prediction of Geological Characteristic Using Gaussian Mixture Model
We-07-06 Prediction of Geological Characteristic Using Gaussian Mixture Model L. Li* (BGP,CNPC), Z.H. Wan (BGP,CNPC), S.F. Zhan (BGP,CNPC), C.F. Tao (BGP,CNPC) & X.H. Ran (BGP,CNPC) SUMMARY The multi-attribute
More informationRELIABILITY ASSESSMENT OF EXISTING STRUCTURES
1. Introduction RELIABILITY ASSESSMENT OF EXISTING STRUCTURES EN Eurocodes, presently being implemented into the system of national standards nearly in the whole Europe, are particularly intended for the
More informationAssessment of atmospheric dispersion and radiological consequences for the Fukushima Dai-ichi Nuclear Power Plant accident
Assessment of atmospheric dispersion and radiological consequences for the Fukushima Dai-ichi Nuclear Power Plant accident IRPA 13 - Fukushima session A.Mathieu, I.Korsakissok, D.Quélo, J.Groëll, M.Tombette,
More informationST495: Survival Analysis: Hypothesis testing and confidence intervals
ST495: Survival Analysis: Hypothesis testing and confidence intervals Eric B. Laber Department of Statistics, North Carolina State University April 3, 2014 I remember that one fateful day when Coach took
More informationUsing Fuzzy Logic as a Complement to Probabilistic Radioactive Waste Disposal Facilities Safety Assessment -8450
Using Fuzzy Logic as a Complement to Probabilistic Radioactive Waste Disposal Facilities Safety Assessment -8450 F. L. De Lemos CNEN- National Nuclear Energy Commission; Rua Prof. Mario Werneck, s/n, BH
More informationStructure Learning in Sequential Data
Structure Learning in Sequential Data Liam Stewart liam@cs.toronto.edu Richard Zemel zemel@cs.toronto.edu 2005.09.19 Motivation. Cau, R. Kuiper, and W.-P. de Roever. Formalising Dijkstra's development
More informationSteve Colwell. British Antarctic Survey
Global Climate Observing System (GCOS) Steve Colwell British Antarctic Survey Goal and Structure of GCOS The Goal of GCOS is to provide continuous, reliable, comprehensive data and information on the state
More information