Knowing what amount of radioactive material was released from Fukushima in March 2011 and at what time instants is crucial to assess the risk, the pollution, and to understand the scope of the consequences. Moreover, it could be used in forward simulations to obtain accurate maps of deposition. But these data are often not publicly available. We propose to estimate the emission waveforms by solving an inverse problem. Previous approaches have relied on a detailed expert guess of how the releases appeared, and they produce a solution strongly biased by this guess. If we plant a nonexistent peak in the guess, the solution also exhibits a nonexistent peak. We pro- pose a method that solves the Fukushima inverse problem blindly. Using atmospheric dispersion models and worldwide radioactivity measurements together with sparse regularization, the method correctly reconstructs the times of major events during the accident, and gives plausible estimates of the released quantities of Xenon.
This code and all associated files are the supplementary material to the paper The Fukushima inverse problem, ICASSP 2013, Vancouver, 2013.
The original paper is available online and attacks the problem of estimating the magnitude and timing of radioactive release at the Fukushima Dai-Ichi power plant during the accident in March 2011.
Andreas Stohl is with the Norwegian Air Research Institute (NILU)
If necessary, the location of the CVX installation directory can be changed directly in makeFigures.m
All the parameters of the algorithms and simulations are described
in a detailed fashion in the comments in
If you need more detail, do not hesitate to contact us.
To recreate all four figures from the paper, follow these steps.
Note: We could unfortunately not release the exact position of the CTBTO stations. As a consequence, figure 2C is omitted.
matrixCleaning.m: Reduce the transport matrix until a given condition number is obtained.
reconstructFromSyntheticData_Xe.m: Simulate reconstruction error for synthetic data for different reconstruction algorithm.
reconstructSourceL1Pos.m: Reconstruct the emission source using L1 minimization with positivity constraints.
opt_routines: Directory containing all the optimization routines we use.
Data/matrixGFSXe.mat : the transport matrix
Data/measXe.mat : the measurement vector
Data/aPrioriSource.mat : the a priori source used in 
Copyright (c) 2013, Marta Martinez-Camara, Ivan Dokmanic, Juri Ranieri, Robin Scheibler, Martin Vetterli, Andreas Stohl
This code is free to reuse for non-commercial purpose such as academic or educational. For any other use, please contact the authors.
The Fukushima Inverse Problem by Marta Martinez-Camara, Ivan Dokmanic, Juri Ranieri, Robin Scheibler, Martin Vetterli, Andreas Stohl is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.
Based on a work at https://github.com/LCAV/FukushimaInverse-ICASSP2013.