Emulation and Sobol' sensitivity analysis of an atmospheric dispersion model applied to the Fukushima nuclear accident. Issue 7 (6th April 2016)

• Record Type: Journal Article
• Authors: Girard, Sylvain; Mallet, Vivien; Korsakissok, Irène; Mathieu, Anne
• Is Part Of: Journal of geophysical research. Volume 121:Issue 7(2016)
• Page Start: 3484
• Page End: 3496
• Publication Date: 2016-04-06
• Subjects: Fukushima -- sensitivity analysis -- Gaussian process -- atmospheric dispersion -- Sobol' -- emulation; Atmospheric physics -- Periodicals; Geophysics -- Periodicals; 551.5
• DOI: 10.1002/2015JD023993
• Languages: English
• ISSNs: 2169-897X

Abstract

Abstract: Simulations of the atmospheric dispersion of radionuclides involve large uncertainties originating from the limited knowledge of meteorological input data, composition, amount and timing of emissions, and some model parameters. The estimation of these uncertainties is an essential complement to modeling for decision making in case of an accidental release. We have studied the relative influence of a set of uncertain inputs on several outputs from the Eulerian model Polyphemus/Polair3D on the Fukushima case. We chose to use the variance‐based sensitivity analysis method of Sobol'. This method requires a large number of model evaluations which was not achievable directly due to the high computational cost of Polyphemus/Polair3D. To circumvent this issue, we built a mathematical approximation of the model using Gaussian process emulation. We observed that aggregated outputs are mainly driven by the amount of emitted radionuclides, while local outputs are mostly sensitive to wind perturbations. The release height is notably influential, but only in the vicinity of the source. Finally, averaging either spatially or temporally tends to cancel out interactions between uncertain inputs. Key Points: We performed a Sobol' sensitivity analysis of an atmospheric dispersion model on the Fukushima case The computational cost was drastically reduced using a Gaussian process emulator of the model Aggregated outputs are controlled by emitted amounts and local outputs by wind perturbations