Abstract. On 11 March 2011, an earthquake occurred about 130 km off the Pacific coast of Japan's main island Honshu, followed by a large tsunami. The resulting loss of electric power at the Fukushima Dai-ichi nuclear power plant developed into a disaster causing massive release of radioactivity into the atmosphere. In this study, we determine the emissions into the atmosphere of two isotopes, the noble gas xenon-133 (133Xe) and the aerosol-bound caesium-137 (137Cs), which have very different release characteristics as well as behavior in the atmosphere. To determine radionuclide emissions as a function of height and time until 20 April, we made a first guess of release rates based on fuel inventories and documented accident events at the site. This first guess was subsequently improved by inverse modeling, which combined it with the results of an atmospheric transport model, FLEXPART, and measurement data from several dozen stations in Japan, North America and other regions. We used both atmospheric activity concentration measurements as well as, for 137Cs, measurements of bulk deposition. Regarding 133Xe, we find a total release of 15.3 (uncertainty range 12.2–18.3) EBq, which is more than twice as high as the total release from Chernobyl and likely the largest radioactive noble gas release in history. The entire noble gas inventory of reactor units 1–3 was set free into the atmosphere between 11 and 15 March 2011. In fact, our release estimate is higher than the entire estimated 133Xe inventory of the Fukushima Dai-ichi nuclear power plant, which we explain with the decay of iodine-133 (half-life of 20.8 h) into 133Xe. There is strong evidence that the 133Xe release started before the first active venting was made, possibly indicating structural damage to reactor components and/or leaks due to overpressure which would have allowed early release of noble gases. For 137Cs, the inversion results give a total emission of 36.6 (20.1–53.1) PBq, or about 43% of the estimated Chernobyl emission. Our results indicate that 137Cs emissions peaked on 14–15 March but were generally high from 12 until 19 March, when they suddenly dropped by orders of magnitude at the time when spraying of water on the spent-fuel pool of unit 4 started. This indicates that emissions may not have originated only from the damaged reactor cores, but also from the spent-fuel pool of unit 4. This would also confirm that the spraying was an effective countermeasure. We explore the main dispersion and deposition patterns of the radioactive cloud, both regionally for Japan as well as for the entire Northern Hemisphere. While at first sight it seemed fortunate that westerly winds prevailed most of the time during the accident, a different picture emerges from our detailed analysis. Exactly during and following the period of the strongest 137Cs emissions on 14 and 15 March as well as after another period with strong emissions on 19 March, the radioactive plume was advected over Eastern Honshu Island, where precipitation deposited a large fraction of 137Cs on land surfaces. Radioactive clouds reached North America on 15 March and Europe on 22 March. By middle of April, 133Xe was fairly uniformly distributed in the middle latitudes of the entire Northern Hemisphere and was for the first time also measured in the Southern Hemisphere (Darwin station, Australia). In general, simulated and observed concentrations of 133Xe and 137Cs both at Japanese as well as at remote sites were in good quantitative agreement. Altogether, we estimate that 6.4 PBq of 137Cs, or 18% of the total fallout until 20 April, were deposited over Japanese land areas, while most of the rest fell over the North Pacific Ocean. Only 0.7 PBq, or 1.9% of the total fallout were deposited on land areas other than Japan.

The Lagrangian particle dispersion model FLEXPART was originally designed for calculating long-range and mesoscale dispersion of air pollutants from point sources, such that occurring after an accident in a nuclear power plant. In the meantime, FLEXPART has evolved into a comprehensive tool for atmospheric transport modeling and analysis at different scales. A need for further multiscale modeling and analysis has encouraged new developments in FLEXPART. In this paper, we present a FLEXPART version that works with the Weather Research and Forecasting (WRF) mesoscale meteorological model. We explain how to run this new model and present special options and features that differ from those of the preceding versions. For instance, a novel turbulence scheme for the convective boundary layer has been included that considers both the skewness of turbulence in the vertical velocity as well as the vertical gradient in the air density. To our knowledge, FLEXPART is the first model for which such a scheme has been developed. On a more technical level, FLEXPART-WRF now offers effective parallelization, and details on computational performance are presented here. FLEXPART-WRF output can either be in binary or Network Common Data Form (NetCDF) format, both of which have efficient data compression. In addition, test case data and the source code are provided to the reader as a Supplement. This material and future developments will be accessible at http://www.flexpart.eu.Peer ReviewedPostprint (published version

The Lagrangian particle dispersion model FLEX-PART in its original version in the mid-1990s was designed for calculating the long-range and mesoscale dispersion of hazardous substances from point sources, such as those released after an accident in a nuclear power plant. Over the past decades, the model has evolved into a comprehensive tool for multi-scale atmospheric transport modeling and analysis and has attracted a global user community. Its application fields have been extended to a large range of atmospheric gases and aerosols, e.g., greenhouse gases, short-lived climate forcers like black carbon and volcanic ash, and it has also been used to study the atmospheric branch of the water cycle. Given suitable meteorological input data, it can be used for scales from dozens of meters to global. In particular, inverse modeling based on source-receptor relationships from FLEXPART has become widely used. In this paper, we present FLEXPART version 10.4, which works with meteorological input data from the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS) and data from the United States National Cen-ters of Environmental Prediction (NCEP) Global Forecast System (GFS). Since the last publication of a detailed FLEX-PART description (version 6.2), the model has been improved in different aspects such as performance, physicochemical parameterizations, input/output formats, and available preprocessing and post-processing software. The model code has also been parallelized using the Message Passing Interface (MPI). We demonstrate that the model scales well up to using 256 processors, with a parallel efficiency greater than 75 % for up to 64 processes on multiple nodes in runs with very large numbers of particles. The deviation from 100 % efficiency is almost entirely due to the remaining nonparallelized parts of the code, suggesting large potential for further speedup. A new turbulence scheme for the convective boundary layer has been developed that considers the skewness in the vertical velocity distribution (updrafts and downdrafts) and vertical gradients in air density. FLEXPART is the only model available considering both effects, making it highly accurate for small-scale applications, e.g., to quantify dispersion in the vicinity of a point source. The wet deposi-Published by Copernicus Publications on behalf of the European Geosciences Union. 4956 I. Pisso et al.: FLEXPART version 10.4tion scheme for aerosols has been completely rewritten and a new, more detailed gravitational settling parameterization for aerosols has also been implemented. FLEXPART has had the option of running backward in time from atmospheric concentrations at receptor locations for many years, but this has now been extended to also work for deposition values and may become useful, for instance, for the interpretation of ice core measurements. To our knowledge, to date FLEX-PART is the only model with that capability. Furthermore, the temporal variation and temperature dependence of chemical reactions with ...

Five different atmospheric transport and dispersion model's (ATDM) deposition and air concentration results for atmospheric releases from the Fukushima Daiichi nuclear power plant accident were evaluated over Japan using regional (137)Cs deposition measurements and (137)Cs and (131)I air concentration time series at one location about 110 km from the plant. Some of the ATDMs used the same and others different meteorological data consistent with their normal operating practices. There were four global meteorological analyses data sets available and two regional high-resolution analyses. Not all of the ATDMs were able to use all of the meteorological data combinations. The ATDMs were configured identically as much as possible with respect to the release duration, release height, concentration grid size, and averaging time. However, each ATDM retained its unique treatment of the vertical velocity field and the wet and dry deposition, one of the largest uncertainties in these calculations. There were 18 ATDM-meteorology combinations available for evaluation. The deposition results showed that even when using the same meteorological analysis, each ATDM can produce quite different deposition patterns. The better calculations in terms of both deposition and air concentration were associated with the smoother ATDM deposition patterns. The best model with respect to the deposition was not always the best model with respect to air concentrations. The use of high-resolution mesoscale analyses improved ATDM performance; however, high-resolution precipitation analyses did not improve ATDM predictions. Although some ATDMs could be identified as better performers for either deposition or air concentration calculations, overall, the ensemble mean of a subset of better performing members provided more consistent results for both types of calculations.

On 11 March 2011, an earthquake occurred about 130 km off the Pacific coast of Japan's main island Honshu, followed by a large tsunami. The resulting loss of electric power at the Fukushima Dai-ichi nuclear power plant (FD-NPP) developed into a disaster causing massive release of radioactivity into the atmosphere. In this study, we determine the emissions of two isotopes, the noble gas xenon-133 (¹³³Xe) and the aerosol-bound caesium-137 (¹³⁷Cs), which have very different release characteristics as well as behavior in the atmosphere. To determine radionuclide emissions as a function of height and time until 20 April, we made a first guess of release rates based on fuel inventories and documented accident events at the site. This first guess was subsequently improved by inverse modeling, which combined the first guess with the results of an atmospheric transport model, FLEXPART, and measurement data from several dozen stations in Japan, North America and other regions. We used both atmospheric activity concentration measurements as well as, for ¹³⁷Cs, measurements of bulk deposition. Regarding ¹³³Xe, we find a total release of 16.7 (uncertainty range 13.4–20.0) EBq, which is the largest radioactive noble gas release in history not associated with nuclear bomb testing. There is strong evidence that the first strong ¹³³Xe release started very early, possibly immediately after the earthquake and the emergency shutdown on 11 March at 06:00 UTC. The entire noble gas inventory of reactor units 1–3 was set free into the atmosphere between 11 and 15 March 2011. For ¹³⁷Cs, the inversion results give a total emission of 35.8 (23.3–50.1) PBq, or about 42% of the estimated Chernobyl emission. Our results indicate that ¹³⁷Cs emissions peaked on 14–15 March but were generally high from 12 until 19 March, when they suddenly dropped by orders of magnitude exactly when spraying of water on the spent-fuel pool of unit 4 started. This indicates that emissions were not only coming from the damaged reactor cores, but also from the spent-fuel pool of unit 4 and confirms that the spraying was an effective countermeasure. We also explore the main dispersion and deposition patterns of the radioactive cloud, both regionally for Japan as well as for the entire Northern Hemisphere. While at first sight it seemed fortunate that westerly winds prevailed most of the time during the accident, a different picture emerges from our detailed analysis. Exactly during and following the period of the strongest ¹³⁷Cs emissions on 14 and 15 March as well as after another period with strong emissions on 19 March, the radioactive plume was advected over Eastern Honshu Island, where precipitation deposited a large fraction of ¹³⁷Cs on land surfaces. The plume was also dispersed quickly over the entire Northern Hemisphere, first reaching North America on 1...