Bayesian transdimensional inverse reconstruction of the Fukushima Daiichi caesium 137 release
<p>The accident at the Fukushima Daiichi nuclear power plant (NPP) yielded massive and rapidly varying atmospheric radionuclide releases. The assessment of these releases and of the corresponding uncertainties can be performed using inverse modelling methods that combine an atmospheric transpo...
| Main Authors: | , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Copernicus Publications
2023-02-01
|
| Series: | Geoscientific Model Development |
| Online Access: | https://gmd.copernicus.org/articles/16/1039/2023/gmd-16-1039-2023.pdf |
| _version_ | 1811168370575278080 |
|---|---|
| author | J. Dumont Le Brazidec J. Dumont Le Brazidec M. Bocquet O. Saunier Y. Roustan |
| author_facet | J. Dumont Le Brazidec J. Dumont Le Brazidec M. Bocquet O. Saunier Y. Roustan |
| author_sort | J. Dumont Le Brazidec |
| collection | DOAJ |
| description | <p>The accident at the Fukushima Daiichi nuclear power plant (NPP) yielded massive and rapidly varying atmospheric radionuclide releases.
The assessment of these releases and of the corresponding uncertainties can be performed using inverse modelling methods that combine an atmospheric transport model with a set of observations and have proven to be very effective for this type of problem.
In the case of the Fukushima Daiichi NPP, a Bayesian inversion is particularly suitable because it allows errors to be modelled rigorously and a large number of observations of different natures to be assimilated at the same time.
More specifically, one of the major sources of uncertainty in the source assessment of the Fukushima Daiichi NPP releases stems from the temporal representation of the source.
To obtain a well-time-resolved estimate, we implement a sampling algorithm within a Bayesian framework – the reversible-jump Markov chain Monte Carlo – in order to retrieve the distributions of the magnitude of the Fukushima Daiichi NPP caesium 137 (<span class="inline-formula"><sup>137</sup></span>Cs) source as well as its temporal discretization.
In addition, we develop Bayesian methods that allow us to combine air concentration and deposition measurements as well as to assess the spatio-temporal information of the air concentration observations in the definition of the observation error matrix.</p>
<p>These methods are applied to the reconstruction of the posterior distributions of the magnitude and temporal evolution of the <span class="inline-formula"><sup>137</sup></span>Cs release. They yield a source estimate between 11 and 24 March as well as an assessment of the uncertainties associated with the observations, the model, and the source estimate.
The total reconstructed release activity is estimated to be between 10 and 20 PBq, although it increases when the deposition measurements are taken into account.
Finally, the variable discretization of the source term yields an almost hourly profile over certain intervals of high temporal variability, signalling identifiable portions of the source term.</p> |
| first_indexed | 2024-04-10T16:25:10Z |
| format | Article |
| id | doaj.art-dc5d3dcd4b05465c9b6b8ed74c312bfc |
| institution | Directory Open Access Journal |
| issn | 1991-959X 1991-9603 |
| language | English |
| last_indexed | 2024-04-10T16:25:10Z |
| publishDate | 2023-02-01 |
| publisher | Copernicus Publications |
| record_format | Article |
| series | Geoscientific Model Development |
| spelling | doaj.art-dc5d3dcd4b05465c9b6b8ed74c312bfc2023-02-09T07:17:09ZengCopernicus PublicationsGeoscientific Model Development1991-959X1991-96032023-02-01161039105210.5194/gmd-16-1039-2023Bayesian transdimensional inverse reconstruction of the Fukushima Daiichi caesium 137 releaseJ. Dumont Le Brazidec0J. Dumont Le Brazidec1M. Bocquet2O. Saunier3Y. Roustan4IRSN, PSE-SANTE, SESUC, BMCA, Fontenay-aux-Roses, FranceCEREA, École des Ponts and EDF R&D, Île-de-France, FranceCEREA, École des Ponts and EDF R&D, Île-de-France, FranceIRSN, PSE-SANTE, SESUC, BMCA, Fontenay-aux-Roses, FranceCEREA, École des Ponts and EDF R&D, Île-de-France, France<p>The accident at the Fukushima Daiichi nuclear power plant (NPP) yielded massive and rapidly varying atmospheric radionuclide releases. The assessment of these releases and of the corresponding uncertainties can be performed using inverse modelling methods that combine an atmospheric transport model with a set of observations and have proven to be very effective for this type of problem. In the case of the Fukushima Daiichi NPP, a Bayesian inversion is particularly suitable because it allows errors to be modelled rigorously and a large number of observations of different natures to be assimilated at the same time. More specifically, one of the major sources of uncertainty in the source assessment of the Fukushima Daiichi NPP releases stems from the temporal representation of the source. To obtain a well-time-resolved estimate, we implement a sampling algorithm within a Bayesian framework – the reversible-jump Markov chain Monte Carlo – in order to retrieve the distributions of the magnitude of the Fukushima Daiichi NPP caesium 137 (<span class="inline-formula"><sup>137</sup></span>Cs) source as well as its temporal discretization. In addition, we develop Bayesian methods that allow us to combine air concentration and deposition measurements as well as to assess the spatio-temporal information of the air concentration observations in the definition of the observation error matrix.</p> <p>These methods are applied to the reconstruction of the posterior distributions of the magnitude and temporal evolution of the <span class="inline-formula"><sup>137</sup></span>Cs release. They yield a source estimate between 11 and 24 March as well as an assessment of the uncertainties associated with the observations, the model, and the source estimate. The total reconstructed release activity is estimated to be between 10 and 20 PBq, although it increases when the deposition measurements are taken into account. Finally, the variable discretization of the source term yields an almost hourly profile over certain intervals of high temporal variability, signalling identifiable portions of the source term.</p>https://gmd.copernicus.org/articles/16/1039/2023/gmd-16-1039-2023.pdf |
| spellingShingle | J. Dumont Le Brazidec J. Dumont Le Brazidec M. Bocquet O. Saunier Y. Roustan Bayesian transdimensional inverse reconstruction of the Fukushima Daiichi caesium 137 release Geoscientific Model Development |
| title | Bayesian transdimensional inverse reconstruction of the Fukushima Daiichi caesium 137 release |
| title_full | Bayesian transdimensional inverse reconstruction of the Fukushima Daiichi caesium 137 release |
| title_fullStr | Bayesian transdimensional inverse reconstruction of the Fukushima Daiichi caesium 137 release |
| title_full_unstemmed | Bayesian transdimensional inverse reconstruction of the Fukushima Daiichi caesium 137 release |
| title_short | Bayesian transdimensional inverse reconstruction of the Fukushima Daiichi caesium 137 release |
| title_sort | bayesian transdimensional inverse reconstruction of the fukushima daiichi caesium 137 release |
| url | https://gmd.copernicus.org/articles/16/1039/2023/gmd-16-1039-2023.pdf |
| work_keys_str_mv | AT jdumontlebrazidec bayesiantransdimensionalinversereconstructionofthefukushimadaiichicaesium137release AT jdumontlebrazidec bayesiantransdimensionalinversereconstructionofthefukushimadaiichicaesium137release AT mbocquet bayesiantransdimensionalinversereconstructionofthefukushimadaiichicaesium137release AT osaunier bayesiantransdimensionalinversereconstructionofthefukushimadaiichicaesium137release AT yroustan bayesiantransdimensionalinversereconstructionofthefukushimadaiichicaesium137release |