Interactive comment on “ Long-term assessment of airborne radiocesium after the Fukushima nuclear accident : Re-suspension from bare soil and forest ecosystems ”

This paper estimated significance of the resuspension/deposition processes in the Cs137 budget at the ground surface in a wide area in the northern part of Japan. Although the conclusion that the resuspension is insignificant in changing the contamination distribution is somewhat too obvious, the procedures and discussions that result in this conclusion are pertinent and informative. It is also interesting that the different sources were found to account for the air concentration variations in the different seasons. Since technical comments were already made in the previous reviewing process, the reviewer raises some points for discussion here.

In contrast, there have been few studies on the long-term (more than 1 year) quantitative assessment of radioactivity in the atmosphere associated with the Fukushima accident (Igarashi et al., 2015;Ishizuka et al., 2016;Kinase et al., 2016).More than 100 000 people were evacuated (METI, 2012), and most have still not been able to return to their homes and the public remains anxious about the safety of the affected areas.Radionuclides with long half-lives such as 134 Cs (2.07 years) and 137 Cs (30.1 years) are of particular concern (Evangeliou et al., 2014).Following the Chernobyl accident there were several studies on the re-suspension and long-term assessment of these radionuclides, such as Holländer and Garger (1996), Garger et al. (1998), Hatano and Hatano (2003), and Garger et al. (2012).For example, Garger et al. (2012) estimated the re-suspension "descending trend" as having a half-life of 300 days based on the surface activity concentration of 137 Cs.In the case of the Fukushima accident, Igarashi et al. (2015) estimated the half-reduction time by fitting multi-component exponential functions based on the 137 Cs concentration at the Meteorological Research Institute (Tsukuba) as 5.9 days, 16 days, and 1.1 years.These estimates were based on the trend in the observed surface air concentrations of 137 Cs, and thus the contributions from advection, diffusion, emission, and deposition terms were not quantified.
There are thousands of monitoring posts situated in the contaminated area in Fukushima prefecture to measure the ambient gamma dose rate.However, the data cannot be used for evaluating internal exposure: evaluation of internal exposure requires direct measurement of the surface air ac-tivity concentration.There are only a few observation sites that continuously measure the concentration of radiocesium (e.g., Igarashi et al., 2015;Ishizuka et al., 2016;Kinase et al., 2016).To assess the spatial distribution of the internal exposure hazard, 3-D numerical simulation is necessary to interpolate values for unmonitored locations.The numerical simulation requires emission flux as a boundary condition.However, the mechanism and thus the radioactivity flux associated with the re-suspension of 137 Cs were unknown, despite extensive efforts based on field observations (e.g., Igarashi et al., 2015;Ishizuka et al., 2016;Kinase et al., 2016).Garger et al. (2012) summarized the re-suspension sources following the Chernobyl accident as (1) dust emission, (2) human activity in fields as well as on roads and construction sites, (3) forest fires, and (4) emissions from the power plant (i.e., opening of the Chernobyl sarcophagus).Re-suspension sources (1) (Ishizuka et al., 2016) and (4), i.e., additional emissions from the reactor buildings of FDNPP (TEPCO, 2012(TEPCO, , 2013(TEPCO, , 2014a(TEPCO, , b, 2015)), were considered in the present study.With respect to source (2), since the town of Namie is located in the evacuation zone, human activity has been extremely limited except for decontaminationrelated work.Although the decontamination-related work could be a potential source of re-suspension, it is not considered in the current simulation as the re-suspension flux has been hardly quantified.As to source (3), there is a low risk of forest fires in Japan given the high humidity, although some open biomass burning remains possible.Kinase et al. (2016) found no increase in 137 Cs activities when the concentration of levoglucosan, a marker of biomass burning, was increased, and thus re-suspension due to biomass burning was not considered in the present study.In addition to the four sources of 137 Cs from the Chernobyl accident, re-suspension from terrestrial biota was considered as suggested by Kinase et al. (2016).They found substantial amounts of bioaerosols upon scanning electron microscopy samples collected in the summer, when 137 Cs concentration was high.
In the present study, the long-term effect of radiocesium re-suspension from contaminated soil and terrestrial biota was quantitatively assessed using 3-D numerical simulation, a field experiment on dust emission flux in a contaminated area (Namie, Fukushima prefecture), and air concentration measurements conducted inside (Namie) and outside (city of Tsukuba, Ibaraki prefecture) the contaminated area.Even though the re-suspension mechanisms still remain unknown, by utilizing the observational data both inside and outside, together with the transport model, we aimed to provide as robust an analysis as possible of the re-suspension, transport, and re-deposition of 137 Cs over the Tohoku and Kanto regions of Japan.

Numerical simulation
A brief description of the numerical method, such as the processes considered in the model and simulation settings, is presented in this section, and detailed model formulations are described in Appendix A. Because the schemes and assumptions regarding the emissions are key to the current study, they are described in detail in the following subsections.

Lagrangian model and simulation settings
Figure 1 shows the domain of the Lagrangian model (LM) with model terrestrial elevations, covering 138-143 • E and 34-39 • N. The model domain covers the southern part of the Tohoku region (the northern mountainous part of the domain, including Yamagata, Miyagi, and Fukushima prefectures), and includes the FDNPP and highly polluted areas such as the habitation-restricted zone (HRZ; 20-50 mSv yr −1 ) and difficult-to-return zone (DRZ; > 50 mSv yr −1 ; METI, 2012).It also covers the Kanto region (or Kanto Plain, the largest plain in Japan, approximately 120 km × 120 km), a highly populated region that includes low to moderately polluted areas such as Tokyo, Gunma, Tochigi, Ibaraki, Saitama, and Chiba prefectures.
The LM considers horizontal and vertical diffusion and advection, gravitational settling, dry and wet depositions, and radioactive decay.It uses simple parameterizations for dry and wet deposition schemes, and it can be driven by meteorological analysis datasets so that it does not require a meteorological model to predict detailed meteorological fields and variables.The model was designed to be easily handled and computationally efficient so that non-specialists of numerical simulations can conduct long-term assessments of atmospheric diffusion problems using their desktop or laptop computers.The LM was designed for rough budget estimates, as presented in the current study, or for sensitivity analyses using a large number of parameters (e.g., Groëll et al., 2014;Quérel et al., 2015;Girard et al., 2016), rather than for process-oriented analysis (e.g., Morino et al., 2013;Katata et al., 2015) or sensitivity analyses of the physical and chemical parameters of aerosols (Adachi et al., 2013).Details of each process and parameter are described in Appendix A. Statistical error of a Lagrangian simulation is inversely proportional to the square of the number of Lagrangian particles (LPs).The statistical accuracy of the current simulation setting is discussed in Appendix B.
The Grid Point Value-Mesoscale Model (GPV-MSM) of the Japan Meteorological Agency (JMA) was used for meteorological analysis to calculate the transport of LPs.It covers 120-150 • E and 23-47 • N and provides 3-hourly and 16 pressure levels of 3-D meteorological variables, from 1000 to 100 hPa, with a horizontal grid resolution of approximately 11 km ( longitude = 0.125 • and latitude = 0.1 • ) and surface variables at twice the resolution of that for the 3-D variables ( longitude = 0.0625 • and latitude = 0.05 • ).In the simulation, the whole model domain where LPs can travel is 138-143 • E, 34-39 • N and from ground surface to 500 hPa.
For output of the model results, LP fields are converted to Eulerian concentration (Bq m −3 ) and deposition (Bq m −2 ) fields in the same horizontal space as the 3-D variables but are vertically allocated from the ground surface to an altitude of 1 km at 100 m intervals.The observed surface air concentration was compared with the simulated mean concentration at 0-100 m above ground level (AGL).

Re-suspension from bare soil
Ishizuka et al. ( 2016) developed a re-suspension scheme for radiocesium from bare soil based on measurements on the schoolyard at Namie High School, Tsushima Campus (denoted as Namie (Tsushima) in Table 1 and Fig. 1), in the DRZ.
where F soil is the 137 Cs dust re-suspension flux from soil (Bq m −2 s −1 ), p 20µm is the surface area fraction of dust smaller than 20 µm in diameter against soil containing a maximum size of 2 mm particles and varies depending on soil texture (1.3 × 10 −8 for sand, 0.19 for loamy sand, 0.45 for sandy loam, and 0.80 for silt loam), F M is the total dust mass flux (kg m −2 s −1 ), f forest is the forest area fraction, and B 5 mm (t) is the specific radioactivity of surface soil (from the surface to a depth of 5 mm; Bq kg −1 ) as a function of time since March 2011.The formula is based on the assumption that dust particles smaller than 20 µm in diameter originating from the surface soil and to a depth of 5 mm were suspended and transported through the atmosphere.The surface area fraction was used for p 20µm based on the assumption that radiocesium is bound to the surface of soil particles (Evrard et al., 2015).F M is formulated as being proportional to the cube of the friction velocity u * (m s −1 ) as described by Loosmore and Hunt (2002) and was applied to the dust emission: Since u * is not available in GPV-MSM, u * was estimated using a wind speed at 10 m AGL by assuming neutral stratification conditions.B 5 mm (t) was derived from the combination of B obs , the observed horizontal distribution of 137 Cs deposition obtained from an airborne radiological survey (NRA, 2012; Bq m −2 ), and r 5 mm , the surface soil activity ratio of 0-5 mm to 0-5 cm obtained from a vertical profile measurement of 137 Cs in the ground soil at Namie High School (i.e., 0.57 Bq Bq −1 ), as where ρ b, soil is the bulk density of soil particles per unit volume in the ground space (kg m −3 ) obtained from the porosity (0.4 m 3 m −3 ) and the density of dust particles (2650 kg m −3 ).
For R decay , which is the decreasing ratio of activity in the ground to its initial state (R decay (t = 0) = 1), only radioactive decay was considered for the re-suspension calculation.The decreasing rate due to other processes such as land surface processes (e.g., runoff, erosion, percolation and all the processes resulting in migration of radiocesium in the soil and biota; Evrard et al., 2015;Matsuda et al., 2015) and decontamination were not considered here.Furthermore, suppression of dust emission due to soil moisture and snow cover was not considered.Therefore, it should be noted here that F soil in Eq. ( 1) is considered as the upper boundary of 137 Cs re-suspension flux from surface soil.Effects such as land surface processes, decontamination, and dust emission suppression due to snow cover are extensively discussed in Sect.5.3 using ambient gamma dose rate measurements obtained by the monitoring posts in Fukushima prefecture.Equation ( 1) is a function of soil texture.The areal fraction of soil texture of the model grid was obtained from the database of the advanced research Weather Research and Forecasting model version 3 (WRFV3; Skamarock et al., 2008).Sixteen categories of soil texture (Miller and White, 1998) with a 30 arcsec resolution dataset can be obtained from the web after subscription at http://www2.mmm.ucar.edu/wrf/users/download/get_sources_wps_geog.htmland were re-categorized into the above-mentioned four categories (i.e., sand, loamy sand, sandy loam, and silt loam) and interpolated to the LM resolution (∼ 11 km) as shown in Fig. 2a-c.Note that the loamy sand fraction is not presented because it is zero for the entire domain.The parameter f forest (Fig. 2d) was also obtained from the database of WRFV3 and was calculated based on the 24 United States Geological Survey (USGS) land use categories, which are constant over time.The land use category dataset can also be obtained from the above website.Ishizuka et al. (2016) validated their dust emission module by using a 1-D model and observed the surface air concen-tration of 137 Cs at Namie in the winter.After applying the module to our 3-D simulation, we found that the air concentration at Namie was underestimated by about factor of 5 for the same period.The module was formulated based on physical parameters (such as u * ) but contains parameters obtained at a single location (such as r 5 mm and ρ b, soil ) and under a fixed atmospheric condition (Ishizuka et al., 2016), whereas ideally parameters in Eqs. ( 1)-(3) should have considered variations among locations and atmospheric conditions for the 3-D simulation.We simply multiplied the dust emission flux by 5 after adjusting the simulation results against the observed concentration of 137 Cs at Namie in the winter.This is one of the simplest top-down approaches for adjusting the emission flux according to the air concentration.The module requires improvement in the future as more reliable parameters become available for various conditions and locations.

Re-suspension from the forest ecosystems
The re-suspension mechanism of radiocesium from land ecosystems remains unknown.Kinase et al. (2016) found substantial amounts of bioaerosols (rather than mineral dust particles) in samples collected for scanning electron microscopy in the summer, when the 137 Cs concentration was high.This does not prove that the bioaerosol was carrying radiocesium but that it could be a potential carrier.The behavior of Cs in the environment can be inferred by analogy with K, a congener of Cs.Potassium is a necessary and abundant element in plants and circulates between land ecosystems.The addition of potassium fertilizer to a rice field in Fukushima significantly reduced the Cs content of the rice (Ohmori et al., 2014).Substantial amounts of Ksalt-rich particles, possibly emitted by active biota such as plants and fungi, and coated with secondary organic aerosols, were observed in pristine Amazonian rainforest (Pöhlker et al., 2012).The major areal fraction of the contaminated area in Fukushima is covered by biota-rich mountain forests.Despite the differences in plant species and locations, it is plausible that water-soluble radiocesium circulating in the biota and soil in the forests was somehow re-emitted into the atmosphere and contributed to the surface air concentration.The re-suspension from the forest ecosystem was simply formulated as follows: where F forest is the 137 Cs re-suspension flux from forest (Bq m −2 s −1 ); f forest is the forest area fraction; f green is the monthly green area fraction; r const is the constant re-suspension coefficient (s −1 ); and, as in Eq. ( 3), B obs and R decay are the observed 137 Cs deposition (NRA, 2012; Bq m −2 ) and the decreasing ratio of activity in the ground, respectively.r const is a tunable parameter to adjust the simulated air concentration of 137 Cs to that observed.In the current study, r const is set to 10 −7 h −1 by adjusting the simulation data using the observed 137 Cs concentration at Namie in the summer, when the re-suspension from soil was negligible due to the higher soil moisture content (following considerable rain) and lower wind speed.As with re-suspension from bare soil, only radioactive decay was considered for R decay and the other processes were not considered.The parameter f green was obtained from the database of WRFV3 and was originally derived from satellite Advanced Very High Resolution Radiometer (AVHRR) normalized difference vegetation index (NDVI) data (Gutman and Ignatov, 1998).Whereas f forest remains constant, the monthly averaged f green was used in order to reflect seasonal changes in the activity of the biota.
2.4 Emission from FDNPP (primary emission, additional emissions from the reactor buildings, and unexpected re-suspension associated with debris removal operations) The Japan Atomic Energy Agency (JAEA)'s latest estimate of the primary emission from FDNPP, Katata et al. (2015), was applied for the emergency situation of March 2011 to evaluate the performance of the LM model against the hori- zontal distribution of 137 Cs deposition of the airborne radiological survey (NRA, 2012; as shown later in Fig. 4a) and surface air concentrations measured at Tsukuba (Fig. 1a).
We selected this inventory because it is JAEA's most up-todate version.Based on an integrated understanding of environmental radioactivity, atmospheric dispersion, and the nuclear reactors, the JAEA team has carefully established a series of inventories for about 5 years, starting with Chino et al. ( 2011), followed by Katata et al. (2012a, b) and Terada et al. (2012), and finally the current inventory (Katata et al., 2015), which is substantially improved compared to the previous versions.
Ongoing emissions during the study analysis period after the emergency situation, i.e., January to December 2013, were obtained from the Tokyo Electric Power Co. Inc. (TEPCO) monthly mean emission flux from the reactor buildings (TEPCO, 2012(TEPCO, , 2013(TEPCO, , 2014a(TEPCO, , b, 2015)).Because only the sum of 134 Cs and 137 Cs was provided, the fractions of these two isotopes were calculated based on their halflives and the assumption that their activities were equal in March 2011 (e.g., Katata et al., 2015), as shown in Fig. 3.The values range from 10 5 to 10 7 Bq h −1 ; however, for simplicity we set a constant value of 10 6 Bq h −1 in the current simulation.
In August 2013, unexpected re-suspension associated with debris removal operations was reported by TEPCO (2014c) and NRA ( 2014) and the gross amount was 10 10 -10 11 Bq of 137 Cs (TEPCO, 2014c;NRA, 2014;Steinhauser et al., 2015).The impact of this unexpected re-suspension is briefly discussed in Sect.5.2 along with an additional finding, but this emission was not considered in the present LM simulation.In this study we focused on the ongoing and continuous emission, mostly from the natural environment, that is difficult to control.

Field observations
Details of the surface air activity concentration measurement techniques can be found in Ishizuka et al. (2016) and Kinase et al. (2016) for Namie and Igarashi et al. (2015) for Tsukuba.At the both sites, 134 Cs and 137 Cs surface air concentrations were measured.The sampling intervals were 1-2 days at Namie and 1 week at Tsukuba for the analysis period of this study, the year 2013.The observations at Namie started on 14 December 2012, while those at Tsukuba started on 31 March 2003, before the FDNPP accident.In March 2011, the sampling interval was increased at Tsukuba to 6 h-1 day and the data for these periods were used for the validation of LM and its parameters, as presented in Sect.4.1.
The Namie site is located on the schoolyard of Namie High School, Tsushima Campus, in the Tsushima district of Namie in Fukushima prefecture, as shown in Table 1 and Fig. 1.The town of Namie extends from the Hamadori coastal area (denoted as C in Fig. 1) to the Abukuma highland area (B in Fig. 1).There are two sites in Namie used in the study.In order to distinguish the Tsushima Campus site in the highland area from the monitoring post site located in Omaru district in Namie in the coastal area, the Tsushima site is referred to as Namie (sometimes as Namie (Tsushima)) and the Omaru site is referred to as Omaru (sometimes as Namie (Omaru)) throughout the manuscript.The Tsukuba site is located on the premises of the Meteorological Research Institute (Table 1 and Fig. 1a).

Results
Section 4.1 presents a validation of the LM model and the optimization of the model deposition parameters by using airborne observations (NRA, 2012) and the emission inventory of Katata et al. (2015) for the emergency situation of March 2011.Using the optimized ranges of model parameters validated in Sect.4.1, the simulated re-suspension of 137 Cs from soil and forest, as well as emission from FDNPP, is presented in Sect.4.2, and the budgets for re-suspension, transport, and re-deposition are presented in Sect.4.3.
Figures 4 and 5 show the observed and simulated distribution of 137 Cs deposition in March 2011, and the scatterplot comparing the observational and simulation results.In the simulation shown in the figures, the "reference" sets used for dry and wet deposition parameters, namely the collection efficiency of aerosols using hydrometeors E c (Eq. A1) and the dry deposition velocity over land v d (Eq.A3), were 0.04 and 0.1 cm s −1 , respectively.
Since LM uses simple parameterizations for dry and wet deposition, as described in Appendix A, sensitivity tests were conducted for model validation, together with optimization of the deposition parameters.Table 2 summarizes the ranges of the deposition parameters for the sensitivity tests and the results of the 137 Cs budget and statistical measurements between the observations and the simulation.
The parameter E c was 0.05 for the JMA dispersion model (JMA, 1998) but the targeted species are different.For example, volcanic ash (particles larger than 1 µm in diameter) used for the JMA model was generally larger than the aerosols carrying 137 Cs (around 1 µm in diameter observed in the downwind area, Tsukuba; Kaneyasu et al., 2012;Adachi et al., 2013).Since the inertia of these smaller 137 Cs particles is likely to be smaller than that for volcanic ash, E c could be smaller.The range of E c was set as 0.02-0.06.The dry deposition velocity v d was selected as 0.1 cm s −1 for 137 Cs in Furuno et al. (1999).The range for v d was set as 0.05 to 0.15 in the present study.
The emission inventory of Katata et al. (2015) amounted to 14.1 PBq from 12 March to 1 April 2013.The simulated deposition over the model domain (138-143 • E, 34-39 • N) ranged from 3.4 to 4.7 PBq, which is approximately 24-33 % of the emission from the FDNPP.Sixty percent of the total deposition occurred over land, for a total of 2.0-2.8PBq, which is close to the observed value of 2.68 PBq, and the observed value is within the range of the sensitivity runs.Statistical measures such as the fractional bias FB, the correlation coefficient R (linear vs. linear), and FAx (fraction of the simulated values within a factor of x) are listed in Table 2. To find better combinations of (or to optimize) the dry and wet deposition parameters, sensitivity runs were screened based on the criteria FA10 > 0.9, FA5 > 0.7, R > 0.75, and an absolute value of FB < 10 %.Only one combination (E c , v d ) = (0.04, 0.1 cm s −1 ) satisfied the criteria, and thus this is referred to as the "reference" parameters.To evaluate the sensitivity (or uncertainty) of the re-suspension simulation for 2013 due to the deposition parameters, the range of the combination of parameters was set as (E c , v d ) = (0.03-0.05, 0.05-0.1 cm s −1 ) around the reference parameters (referred to as the "optimized range") by excluding the parameters with the worse performances.The ranges of the statistical measures of the optimized runs are listed in Table 2. FB, R, FA2, FA5, and FA10 after the optimization had the www.atmos-chem-phys.net/16/13149/2016/ranges −0.18 to −0.036, 0.74-0.77,0.26-0.30,0.68-0.74,and 0.91-0.92,respectively.These statistical measures were comparable to those reported in previous multi-model comparison studies (R: 0.27-0.85;FB: −0.84-0.56;and FA2: 0.14-0.57, in SCJ, 2014; Draxler et al., 2015).The current model is thus shown to be sufficiently credible for the budget analysis in this study, despite the simple parameterization and the low resolution in space (∼ 11 km) and time (3 h).
In agreement with many previous studies, the simulated contribution of wet deposition was larger than that of dry deposition: the ratio of the amount of dry to wet deposition ranged from 0.12 to 0.23 for the optimized parameter ranges, indicating that the results were less sensitive to the dry deposition parameter.Generally speaking, R became higher as E c became lower, whereas FAx became higher as E c became higher for the various ranges of the sensitivity tests.Therefore, lower E c did not meet the criteria of FAx and higher E c did not meet the criteria of R. Consequently, after the optimization, the maximum values of the statistical measures were lower but the minimum values became higher, indicating that the optimization was successful in excluding the parameters with the worse performances (rather than selecting  the best parameters).It should be noted here that the optimized deposition parameters are not necessarily physically valid but rather give simulation results consistent with the available observations.The results presented in this section indicate that the current LM simulation with these optimized parameters has the potential to reproduce consistent features of the radiocesium budget over the Tohoku and Kanto regions of Japan.
Figure 6 shows the temporal variation in simulated (red) and observed (black) 137 Cs concentrations at Tsukuba in March 2011.The model reproduced the three major plumes arriving at Tsukuba on 15-16, 20-21, and 28-30 March; a plume on 23 March was not observed but only appeared in the simulation.The red shaded areas indicate the range of 137 Cs concentrations obtained when the simulations were run using the optimized parameter ranges E c = 0.03-0.05and v d = 0.05-0.1 cm s −1 .Due to differences in the parameters used, the surface air concentration could vary by approximately 1 order of magnitude in transported plumes that experienced wet scavenging along their pathway.

Re-suspension in 2013
Figure 7 shows the observed and daily mean simulated (with the reference parameters) surface air concentrations at Namie and Tsukuba for the year 2013.The red, green, and blue lines indicate re-suspension from soil, re-suspension from forest, and emission from the FDNPP reactor buildings, respectively.Scatterplots between the observed and the simulated total concentrations are also shown in the figure and the statistical measures between them are listed in Table 3.The correlation coefficients (linear vs. linear) are low at the both sites, but the simulated and the observed medians are close to each other and FAx shows high values at the both sites.Consequently, the current simulation is credible enough for the quantitative assessment of the annual radiocesium budget.
As discussed in Sect.2.2, note that the re-suspension flux due to the dust emission module (Ishizuka et al., 2016) is multiplied by 5 in this study in order to adjust to the observation level at Namie in the cold season (January to March, October to December).Also note that the re-suspension coefficient r const in Eq. ( 4) was set as 10 −7 h −1 in order to adjust to the observation level at Namie in the warm season (May to September).The emission flux reported by TEPCO varied from 10 5 to 10 7 Bq h −1 during the study period but we set it to 10 6 Bq h −1 for simplicity.We did not try to precisely adjust r const to the observation by, for example, using inverse modeling, and instead we simply multiplied by power-of-10 values (except for dust emission, which is five), which are constant in time and space because (1) this facilitated straightforward interpretation of the simulation results by keeping the simulated variation solely originating from the variation in boundary conditions (i.e., emission and meteorology), and therefore (2) this provided simple but useful clues for understanding the re-suspension mechanisms, which remain unknown.
Using the dust emission module (which has a physical basis), 137 Cs flux re-suspended from soil could account for the level of the observed surface air concentration of 137 Cs at Namie in the cold season.Under the influence of the northwesterly winter monsoon, the surface wind speed is high over the contaminated area compared to the summer (as shown later in Fig. 10).Note that the flux might be a maximum estimate since it does not consider land surface processes (such as soil moisture, snow cover, or migration of 137 Cs in the soil and biota) and decontamination, which could reduce the 137 Cs re-suspension flux.In contrast, in the warm season, the estimated flux significantly underestimated the observation by 1 to 2 orders of magnitude due to the weak surface wind, indicating that the dust emission process may not be the sole process involved in sustaining the air concentration of 137 Cs during this period.Introducing the 137 Cs re-suspension component from forest with a resuspension coefficient of 10 −7 h −1 and a monthly variation in the green area fraction (derived from NDVI) could quantitatively account for the observed air concentration together with its seasonal variation at Namie.Even though both the simulated re-suspension from soil and forests reproduced the quantity and seasonal variation in the background concentration (in other words, concentrations originating from continuously presenting emissions, such as natural emissions, and not accidental ones) at Namie, sporadic peak events, such as the daily mean 137 Cs concentration exceeding 10 mBq m −3 as observed in June and August at Namie, were not simulated.Some specific re-suspension events might occur from the highly contaminated areas such as the premises of FD-NPP (e.g., debris removal operations) or very close to FD-NPP on these days, as indicated later in Sect.5.2.The simulated 137 Cs concentrations due to the monthly mean emission from the reactor buildings (i.e., 10 6 Bq h −1 ) significantly underestimated the observed concentration by more than 3 orders of magnitude at Namie and by 2 orders of magnitude at Tsukuba.Even the maximum estimate of 10 7 Bq h −1 does not reach the observed level.The emission from FDNPP may not have been the sole process sustaining the air concentration of 137 Cs in 2013, supporting the discussion in Igarashi et al. (2015), which concluded that direct emission from the FDNPP played a minor role in the observed atmospheric radiocesium concentrations over Tsukuba during 2013-2014.
The observed air concentration of 137 Cs at Tsukuba was about 1 to 2 orders of magnitude lower than that at Namie.The simulated difference between the two sites inside and outside of the contaminated areas was consistent with the observed difference.This finding indicates that the current LM simulation provided consistent features of re-suspension, transport, and re-deposition in the Tohoku and Kanto regions of Japan of 137 Cs originating from Fukushima.

Budget analysis
Figure 8 illustrates the simulated (with the reference parameters) annual amounts of total re-suspension and re-deposition of 137 Cs, together with their ratios to the observed deposition (Fig. 4a).The simulated areal total re-suspended amount was 1.28 TBq, which was equivalent to 0.048 % of the total deposited amount, 2.68 PBq.The areal total re-deposited amount (with the reference parameters) was 0.22 TBq (0.18-0.23 TBq for the optimized range of the deposition parameters), corresponding to approximately 17.2 (14.1-18.0)% of the re-suspended amount deposited mainly in the Tohoku region, with the remainder being transported out of the model domain.Therefore, the regional mean rate in the decrease of the land surface 137 Cs concentration due to re-suspension was estimated to be 0.040 (0.039-0.041) % yr −12 , equivalent to 1.1 (1.06-1.12)× 10 −6 day −1 .The spatial distribution of the re-suspension and re-deposition ratio to the primary deposition ranged from 0.01 to 0.3 and 0.001 to 0.03 %, respectively.The spatial distribution of the land surface 137 Cs deposition decay due to re-suspension ranged from 2.7 × 10 −7 to 8.2 × 10 −6 day −1 .Re-suspension therefore had a negligible effect on reducing land surface radioactive contamination.

Discussion
Seasonal variation in the surface activity concentration and its source contributions are extensively discussed in Sect.5.1.A possible source of the observed sporadic peak events, which could not be reproduced by the simulation, is dis- cussed in Sect.5.2.The effects of other processes that were not considered in the model, such as land surface processes and decontamination, are discussed based on the dose rate measurements from the monitoring posts in Fukushima in Sect.5.3.Future issues are summarized in Sect.5.4.

Seasonal variation and source contribution
The discussion in this section expands on that in Sect.4.2. Figure 9 shows the same temporal variation as Fig. 7 but for simulated (using the optimized ranges of parameters) results for 137 Cs from dust and FDNPP in winter (January to March) and from forest and FDNPP in summer (June to August).
In the winter, the simulated trend for dust agreed well with the observed trend (Fig. 9a), and the surface air concentration during this period was positively correlated with the surface wind speed in both the simulation (R = 0.88 on an hourly basis) and the observations (R = 0.32 on a daily basis).There was a sporadic peak in the observational data of 6.7 mBq m −3 from the 17 March at 13:00 local time (LT) to 18 March at 13:00 LT that could not be reproduced by the dust module, and this peak coincided with a plume arriving from FDNPP, as shown in Fig. 9c.Therefore, the observed peak could be accounted for by specific re-suspension events on the order of 10 9 Bq h −1 (the left axis divided by the right axis multiplied by 10 6 Bq h −1 in Fig. 9c) if they occurred on the premises of FDNPP or close to the area.There are also two events exceeding 2 mBq m −3 , one in January and another in February.It is unlikely that the two peaks originated from the direct emission from FDNPP and likely that they originated from the dust emission because the observed peaks coincided with the simulated dust peaks (Fig. 9a) and not with the simulated peaks due to the FDNPP emission (Fig. 9c).
In the summer, the simulated quantity as well as the variation in the forest data agreed well with the observed data (Fig. 9b).Because there is only monthly variation in the simulated emission, the simulated daily trend solely originated from variations in the meteorological parameters (wind field, turbulent mixing, and wet scavenging).A significant peak of 60.4 mBq m −3 is observed from 14 August at 13:00 LT to 15 August at 13:00 LT.This observed level was approximately 2 orders of magnitude larger than the simulated level and 1 to 2 orders of magnitude larger than the observed level for the other days in this period.Therefore, continuous emission such as re-suspension from forest is less likely to be the origin of the peak.Because the observed peak and the simulated peak of 137 Cs from FDNPP coincided (Fig. 9d), the observed level could be accounted for by specific re-suspension events on the order of 10 10 Bq h −1 either on the premises of FD-NPP or close to the area.There have been several arguments that the observed peaks in August 2013 were associated with debris removal operations at FDNPP, and this is discussed separately in Sect.5.2.
Figure 10 illustrates the seasonal mean surface wind vector and surface air 137 Cs concentration (simulated using the reference parameters) due to (a) dust re-suspension in the winter and (b) forest re-suspension in the summer.Due to the prevailing northwesterly winter monsoon, 137 Cs was carried southeastward in the winter.In the summer, under the influence of the Pacific high-pressure system, 137 Cs was carried inland.The monthly mean wind speed is high in winter and low in summer.The upper panels of Fig. 11 illustrate the observed and simulated (using the optimized ranges of parameters) total (from soil, forest, and FDNPP) 137 Cs concentration at Namie and Tsukuba.The time resolutions of the simulation are daily for Namie and weekly for Tsukuba to be consistent with the sampling intervals of the two respective sites.The simulation successfully reproduced the quantity and variation in the observed background concentration at Namie and Tsukuba but could not reproduce the sporadic peak events observed at Namie, as discussed above.The simulation also significantly underestimated the observations at Tsukuba from January to March, 2013.Due to the northwesterly monsoon (Fig. 10a), there was less air mass transported from FDNPP to Tsukuba in the winter and therefore this underestimation is probably due specifically to underestimation of the simulated re-suspension around Tsukuba.The lower panels of Fig. 11 show the relative contributions of 137 Cs from soil and forests at Namie and Tsukuba.The contribution from FDNPP was negligible throughout the year.At both sites, the contribution from dust was high (0.7-0.9) in the cold season and low (0.2-0.4) in the warm season due to higher surface wind speed in the cold season.
Figure 12a shows the observed and simulated (with the reference parameters) monthly Namie to Tsukuba 137 Cs concentration ratios.The mean concentration ratio exceeded 100 in June and 200 in August due to the sporadic peak events.The monthly median would be relevant for comparing the background observation with the simulation results by considering only continuous emission sources.The values of the simulated concentration ratio and its seasonal variation agreed fairly well with the observed monthly median ratio: the observed and simulated annual means were 38.9 and 31.1, respectively.Figure 12b shows the monthly mean simulated re-suspension source area contributions to the 137 Cs air concentration at Namie and Tsukuba.The re-suspension source area is defined as the model grid where the observed deposition amount exceeded 300 kBq m −2 (Fig. 4a) and includes the Namie grid (2300 kBq m −2 ).Eighty to 90 % of the 137 Cs air concentration at Namie originated from the source region, and there was no clear seasonal variation in the value.In contrast, Tsukuba is characterized as a downwind region and there was clear seasonal variation in the source contribution ratio: high in summer and low in winter, due to the summer and winter monsoons, as discussed above.Nonetheless, the highest value at Tsukuba was 0.4 in July, and so more than half of the 137 Cs concentration at Tsukuba originated locally or from areas other than the contaminated regions throughout the year.As shown in Fig. 11a and b, the variability in the simulated concentration at Tsukuba due to uncertainty in the deposition parameters was much larger than that at Namie.The differences in the variability indicated that the Namie and Tsukuba sites can be characterized as the source area and the downwind area, respectively: as the time required for the plume to move from the emission site to the observation site increases, the variability becomes larger due to the increased chance for the plume to experience dry and wet scavenging.

Possible source of sporadic peak events
There have been several scientific studies and governmental reports on the unexpected re-suspension from FDNPP in August 2013.The high-dose-rate alarm was activated on August 19 within the premises of FDNPP associated with the debris removal operation.Matsunami et al. (2016) related the radiocesium contamination of brown rice in Fukushima in 2013 to this operation, whereas MAFF (2015) denied any association.The NRA estimated the 137 Cs emission rate during the debris removal operation as 6.7 × 10 10 Bq h −1 and the cumulative amount as 1.1 × 10 11 Bq (NRA, 2014).TEPCO (2014c) estimated the emission rate during the operation as 5.8 × 10 10 -1.2 × 10 11 Bq h −1 and the cumulative amount as 1.3-2.6 × 10 11 Bq.Steinhauser et al. (2015) estimated the gross amount as 2.8 × 10 11 Bq using measurements of weekly air filter sampling and monthly deposition, and a numerical simulation.Their estimates are similar in magnitude to our estimate (10 10 Bq h −1 ; see Sect.5.1) but the dates are different: our observed peak was earlier than the reported removal operation.
Our daily sampling showed a peak concentration (60.4 mBq m −3 ) from 14 August at 13:00 LT to 15 August at 13:00 LT before the reported operation but did not detect high concentrations in the 19 August (0.33 mBq m −3 for August 18 at 13:00 LT to 19 August at 13:00 LT and 1.2 mBq m −3 for 19 August at 13:00 to 20 August at 13:00 LT). Figure 13 shows the forward trajectories predicted by the LM (statistical locations of LPs) starting from FDNPP on 14 August (left) and 19 August (right).The sky-blue lines and red dashed circles indicate areas containing approximately two-thirds of the LPs within 1 km AGL: the extent of the area reflects horizontal and vertical atmospheric diffusion.The highest dose rate peaks were observed from 13:50 to 14:10 LT on 19 August at 2.8-8.3 km north and north-northwest of FDNPP on the leeward side, as reported by Fukushima prefecture (https://www.pref.fukushima.lg.jp/ download/1/20130827moni.pdf.pdf).The forward trajectories on 12 August indicated that plumes during the debris removal operation traveled north-northwest to north of FDNPP (Fig. 13b and d), rather than towards the west-northwest, where the Namie site is located.On 14 August, on the other hand, plumes were transported towards the west (starting at 12:00 LT, Fig. 13a), and then to the north (starting at 15:00 LT, Fig. 13c) due to fast changes in wind direction, resulting in the simulated peak concentration shown in Fig. 9d during this period.Our simulation and observations together indicated that the same order of magnitude of 137 Cs emission occurred on August 14-15 and on 19 August.Alarm activation was not reported on August 14-15, but debris removal operation was also conducted on August 14 and 16 (MAFF, 2015).
www.atmos-chem-phys.net/16/13149/2016/5.3 Land surface processes, decontamination, and dust emission suppression due to snow cover Figure 14 presents the time series of ambient gamma dose rates measured at the monitoring posts in Fukushima prefecture indicated in Fig. 1b.The data were obtained from the Nuclear Regulation Authority (NRA), Japan, website (http://radioactivity.nsr.go.jp/map/ja/index.html).A total of six monitoring posts, two from each of the three geographical areas (Hamadori coastal area, Abukuma highland area, and Nakadori valley area), were selected.There are tens to hundreds of monitoring posts in each municipality (village, town, and city) in Fukushima prefecture.We selected a monitoring post in each municipality by applying the following conditions: a post showing the highest dose rate of all posts in the municipality at the time of downloading (around 11:00, 28 December 2015), data are available since 1 April 2012, and the instruments are situated 100 cm above the ground.
The dose rate significantly dropped when the ground was covered with snow, in January and December 2013 in Hamadori and Nakadori, and from January to early March and December 2013 in Abukuma (the elevation of Abukuma is 500-1000 m and higher than Hamadori and Nakadori).Snow cover suppresses re-suspension due to dust emission.Namie (Tsushima) is located in Abukuma and the ground was covered with snow until early March (Ishizuka et al., 2016).The observed air concentration of 137 Cs at Namie (Tsushima) was correlated with the wind speed in the winter from January to March, indicating that re-suspension during the period was mechanically induced.In the winter, dust resuspension from outside Abukuma, such as from Hamadori and Nakadori, or from land surface where the snow cover was partly melted due to solar radiation, might be the dominant source contributing to the 137 Cs surface air concentration at Namie (Tsushima) when the ground was covered with snow.
The first-order decreasing rates fitted by the least-squares approximation for the period without snow cover, May to October 2012 and 2013, are presented in Fig. 14.The rates ranged from 5.2 to 12.1 × 10 −4 day −1 .The monthly mean radioactive decay rates of total radiocesium ( 134 Cs + 137 Cs), determined by assuming that the activities of 134 Cs and 137 Cs were equivalent in March 2011, were 4.2 × 10 −4 and 3.0 × 10 −4 day −1 in April 2012 and March 2014, respectively, due to the difference in half-life of 134 Cs and 137 Cs (2.07 and 30.1 years, respectively).By assuming that the gamma dose rate primarily originated to radiation from the land surface radiocesium, the radioactive decay accounted for 35-50 % of the decreasing rate of total ground radioactivity; the exception was Shirakawa, where radioactive decay accounted for 55-80 % of the decrease.In other words, 50-65 % of the ground radioactivity decrease was likely due to land surface processes, decontamination, and re-suspension to air.As discussed in Sect.4.3, the estimated decreasing rate due to re-suspension was 2.7 × 10 −7 -8.2 × 10 −6 day −1 , which is 2 to 3 orders of magnitude smaller than the decreasing rates due to the other processes (10 −4 -10 −3 day −1 ).
It is difficult to distinguish the contributions of land surface processes and decontamination.By subtracting the radioactive decay rate (3.0-4.2 × 10 −4 day −1 ) and the decreasing rate due to re-suspension (2.7 × 10 −7 -8.2 × 10 −6 day −1 ) from the gross decreasing rate (5.2-12.1 × 10 −4 day −1 ), the estimated decreasing rates due to land surface processes and decontamination ranged from 1.0 to 7.9 × 10 −4 day −1 .Matsuda et al. (2015) summarized the depth profiles of radiocesium in soil at more than 80 locations in Fukushima, including Hamadori, Abukuma, and Nakadori.They found that the radiocesium levels have been slowly migrating downward, with rates ranging from 1.7 to 9.6 kg m −2 yr −1 (equivalent to 1.1-6.0mm yr −1 for a dust particle density of 2650 kg m −3 and a porosity of 0.4 m 3 m −3 , for example).Evrard et al. (2015) summarized that significant transfer of particulate-bound radiocesium occurs during major rainfall and runoff events (e.g., typhoons and spring snowmelt).Together with the relaxation depth-dose rate relationship provided by Saito and Petoussi-Hess (2014), the decreasing rate due to land surface processes such as downward migration, runoff, and erosion could be quantified and thus the decontamination effect could be separately extracted.

Future issues
Issues that remain to be resolved in future research are summarized in this section.First of all, re-suspension from biota could be predominant in the warm season, but the resuspension sources as well as mechanisms remain essentially unknown.Further study is needed to understand the mechanism based on field experiments and numerical simulations.The current estimation could account for the measured background concentration (0.1-1 mBq m −3 ) but could not reproduce the observed sporadic peak concentration (1-10 mBq m −3 ) at the Namie site.Further study is needed to identify the cause.The dust flux module has been validated at a single location.The module needs to be improved to be applicable to various land use and soil texture conditions.The decontamination may reduce resuspension afterward, whereas the resuspension may occur during decontamination-related work.This effect should be evalu-ated in the future.The current estimation was based on a single model simulation.Variability in multi-model simulations is rather large (SCJ, 2014;Draxler et al., 2015) and therefore multi-model assessment will be indispensable for long-term re-suspension analysis.

Conclusions
The long-term effect of 137 Cs re-suspension from contaminated soil and biota due to the Fukushima nuclear accident has been quantitatively assessed using a numerical simulation, a field experiment on dust emission in the contaminated area (Namie, Fukushima), and air concentration measurements inside (Namie) and outside (Tsukuba, Ibaraki) of the area.The re-suspension mechanism remains unknown.We therefore utilized the observational data obtained both inside and outside the contaminated area, together with a transport model, to provide a robust budget analysis of the re-suspension, transport, and re-deposition of 137 Cs in the eastern part (the Tohoku and Kanto regions) of Japan.Our findings are summarized as follows: 1. Optimization of the deposition parameters of the LM for simulating the emergency situation of March 2011, using aircraft observation data (NRA, 2012) and the prescribed emission inventory (Katata et al., 2015), provided 0.1 (0.05-0.1) cm s −1 for a dry deposition velocity over land and 0.04 (0.03-0.05) for a hydrometeor collection efficiency for aerosols.The optimized (or validated) ranges of the deposition parameters were applied to long-term re-suspension assessment for the year 2013.
2. Using the dust emission module (Ishizuka et al., 2016), which was developed based on physical parameters, simulated 137 Cs re-suspension from soil multiplied by 5 accounted for the observed 137 Cs surface air concentration measured at Namie during only the cold season; the module underestimated the 137 Cs concentration by 1 to 2 orders of magnitude in the warm season.
3. Introducing re-suspension from forest using a constant re-suspension coefficient of 10 −7 h −1 and monthly green area fraction could quantitatively account for the observed concentration together with its seasonal variation.
4. The contribution from additional emission from the reactor buildings of FDNPP (10 6 Bq h −1 ) was negligible throughout the year and underestimated the observed air concentration by 2 to 3 orders of magnitude at both observation sites.
5. At Namie and Tsukuba, the simulated contribution of re-suspension from soil was high (0.7-0.9) in the cold season and low (0.2-0.4) in the warm season; the remaining contribution was from forest and was low in winter and high in summer.The contribution of the resuspension from the source area (where the aircraftobserved deposition exceeded 300 kBq m −2 ) to the air concentration at Namie was 0.8-0.9throughout the year, while that at Tsukuba varied from 0.1 to 0.4, and was high in the summer and low in the winter.
6.The simulated annual total re-suspended amount for the whole region was 1.28 TBq, equivalent to 0.048 % of the aircraft-observed total deposited amount of 2.68 PBq.The total re-deposition was 0.18-0.23TBq, equivalent to 14.1-18.0% of the total re-suspended amount: the rest of the 137 Cs was transported out of the model domain.The spatial distribution of the decreasing rate of land surface 137 Cs due to re-suspension ranged from 2.7 × 10 −7 -8.2 × 10 −6 day −1 .
7. The first-order decrease rate of the ambient gamma dose rate in Fukushima prefecture ranged from 5.2 to 12.1 × 10 −4 day −1 .By subtracting the radioactive decay rate of 3.0-4.2× 10 −4 day −1 , the ground radioactivity decay due to land surface processes, decontamination, and re-suspension was found to range from 1.0 to 7.9 × 10 −4 day −1 .The estimated re-suspension rate was 2 to 3 orders of magnitude lower than the decrease in rate due to the other processes, showing that re-suspension contributed negligibly towards reducing ground radioactivity.

Data availability
All data are available upon request to the authors.
www.atmos-chem-phys.net/16/13149/2016/Atmos.Chem.Phys., 16, 13149-13172, 2016 with the land surface (soil, water, and vegetation), and gravitational settling for aerosols (e.g., Wesely and Hicks, 2000).Therefore, v d is a function of height as well as of turbulent flux and surface conditions.Nevertheless, v d is set as constant in the simulation, but the height dependency of v d is considered in the dry scavenging rate dry (s −1 ), following Furuno et al. (1999) as where z is the height of the LP (m AGL) and z srf is the surface layer height set as 100 m AGL in the study.Instead of explicitly predicting v d , its range was set at around 0.1 cm s −1 , a typical speed for a range of aerosols around 1 µm in diameter, for the sensitivity tests, as listed in Table 2.The value of v d is applied over land, whereas v d over the ocean is multiplied by 0.1, because v d over a flat surface is approximately 1 and 2 orders of magnitude smaller than v d over short vegetation such as grass and tall vegetation such as forest, respectively (e.g., Petroff and Zhang, 2010).

Appendix B: Statistical accuracy of the current simulation setting
Because the statistical error of Lagrangian simulation is inversely proportional to the square of the number of LPs, the statistical accuracy of the current simulation setting was evaluated using the following measures (relative errors of quantities of the sensitivity runs to those of the reference run):

B1 Point source case
The number emission rate of LPs, N LP , was set as 32 000 h −1 (i.e., N LP_ref ) for a point source emission case such as the primary emission in March 2011 and additional emission from the reactor buildings in 2013.The median values together with the 25th and 75th percentile values of E con and E dep of the sensitivity runs (sensitivity to deposition parameters and sensitivity to N LP ) against the reference run are listed on the rows of the top half of Table B1.Both E con and E dep of N LP_ref × 4 were significantly lower than those for the deposition parameters sensitivity run.This result indicates that 32 000 h −1 for N LP was sufficient to allow a statistically significant simulation for the purpose of this study, as shown in Figs. 4 and 5: the difference in concentration and deposition due to the deposition parameters was much larger than the difference due to model uncertainty in N LP .E con and E dep of N LP_ref × 0.25 (which are also smaller than those of the deposition parameters sensitivity run) are larger than those of N LP_ref × 4, indicating fairly good accuracy convergence of the LM model.

B2 Areal emission case
N LP_ref was 16 h −1 grid −1 for the areal emission case simulating re-suspension from soil and forests in 2013.The lower half of Table B1 is the same as the upper half except for the areal emission case (re-suspension from forest).Both E con and E dep of N LP_ref × 4 were much lower than those for the deposition parameters sensitivity run, indicating that 16 h −1 grid −1 supports a statistically significant simulation for the purpose of this study, as shown in Figs.7-12.Usually, Lagrangian-type models are not appropriate for solving areal emission problems, because they require many more LPs for areal emission cases compared to point source cases and thus become computationally too expensive to obtain statistical accuracy (i.e., to set large number of LPs).In the case of this simulation, especially for the concentration, sensitivity to deposition parameters was much more significant than sensitivity to model uncertainty in N LP using the sufficiently small number of N LP_ref = 16 h −1 grid −1 .E con and E dep of N LP_ref × 4 are smaller than those of N LP_ref × 0.25, indicating fair accuracy convergence of the LM model.

Figure 1 .
Figure 1.(a) The model domain showing the model terrestrial elevation, observation sites, and other locations described in the study.The linear distances from FDNPP to Namie (Tsushima district, town of Namie), Tsukuba and Tokyo are approximately 30, 170, and 220 km, respectively.The numbers denote prefectures: 1, Miyagi; 2, Fukushima; 3, Ibaraki; 4, Tochigi; 5, Gunma; 6, Saitama; 7, Chiba; and 8, Tokyo.(b) Fukushima prefecture and (colored circles) the locations (village, town, or city name) of monitoring posts used in this study (see Fig. 14).The red letters in both (a) and (b) denote the name of the area based on geographical features: A, Nakadori valley; B, Abukuma highland; and C, Hamadori coastal area.

Figure 2 .
Figure 2. The areal fractions of (a)-(c) soil texture and (d)-(f) land use category used for the boundary conditions of the simulation.

Figure 3 .
Figure 3. Monthly mean emission flux of radiocesium released from the reactor buildings of FDNPP from October 2011 to August 2015 as estimated by TEPCO (TEPCO, 2012-2015).
(A1) b Dry deposition velocity over land; see Eq. (A3).c Simulated deposited amount over the whole model domain.d Simulated deposited amount only over land.e Fractional bias between D land and D obs .f Correlation coefficient between each grid cell of the observed and simulated deposition (linear vs. linear).g Fraction of simulated values within a factor of 2 of the observed values.h Fraction of simulated values within a factor of 5 of the observed values.f, g, h Compared only at grids where the observed values are greater than 10 kBq m −2 .

Figure 4 .
Figure 4. (a) Aircraft observation and (b, c) simulation of 137 Cs deposition depicted (b) only over land and (c) for the whole domain (kBq m −2 ).The observation was interpolated to the model grid ( longitude = 0.125 • and latitude = 0.1 • ).A decay correction for the observation was made for March-May 2012, varying depending on the region.The simulation provided a cumulative value from 12 March to 1 April 2011.The total activities are presented as numbers.The color scales are the same for (a)-(c) but observed values are not depicted below the detection limit, 10 kBq m −2 .The deposition parameters of the simulation in the figure are E c = 0.04 and v d = 0.1 cm s −1 .

Figure 5 .
Figure 5. Scatterplot between the observational data and the simulation results for 137 Cs deposition (kBq m −2 ).The deposition parameters of the simulation in the figure are E c = 0.04 and v d = 0.1 cm s −1 .N indicates the number of samples, and the statistical measures R, FA2, and FA5 are described in Table2.

Figure 6 .
Figure 6.Time series of the surface activity concentration of (black) observed and (red) simulated 137 Cs concentrations at Tsukuba from 12 March to 1 April 2011 (Bq m −3 ).The red shaded areas indicate the range of 137 Cs concentrations obtained when the simulations were run using the optimized parameter ranges E c = 0.03-0.05and v d = 0.05-0.1 cm s −1 .

Table 3 .Figure 7 .
Figure 7. Time series of the surface air concentration of (black) observed 137 Cs and (colors) simulated daily 137 Cs levels on the left and scatterplots between observed and simulated total (i.e., sum of colors) 137 Cs at (a) Namie and (b) Tsukuba.The colored lines indicate simulated 137 Cs concentrations due to (red) re-suspension from soil using the scheme given in Ishizuka et al. (2016) (five times), (green) re-suspension from forest with a re-suspension rate of 10 −7 h −1 , and (blue) emission from the FDNPP reactor buildings with a constant emission rate of 10 6 Bq h −1 .The deposition parameters of the simulation in the figure are E c = 0.04 and v d = 0.1 cm s −1 .

Figure 8 .
Figure 8. Annual total amounts of (a) the re-suspension and (b) the re-deposition of 137 Cs.The total activities are presented as numbers.Fractions of the re-suspension and the re-deposition to the observed 137 Cs deposited amounts are also shown in (c) and (d), respectively.The deposition parameters of the simulation in the figure are E c = 0.04 and v d = 0.1 cm s −1 .

Figure 9 .
Figure 9.Time series of observed (black) and simulated (colors) 137 Cs surface air concentration at Namie (a, c) in winter from January to March and (b, d) in summer from June to August 2013.The colors of the lines are the same as in Fig. 7 but the time variation is hourly.The shaded areas indicate the range of 137 Cs concentrations obtained when the simulations were run using the optimized parameter ranges E c = 0.03-0.05and v d = 0.05-0.1 cm s −1 .

Figure 10 .
Figure 10.Seasonal mean surface (10 m above ground level) wind vector and 137 Cs surface concentration (a) due to dust re-suspension in winter from January to March and (b) due to forest re-suspension in summer from June to August.The deposition parameters of the simulation in the figure are E c = 0.04 and v d = 0.1 cm s −1 .

Figure 11 .
Figure 11.Time series of (a, b; black) the observed and (purple) the simulated surface 137 Cs concentration due to total re-suspension and (c, d) the relative contribution of (red) dust and (green) forest re-suspension to the 137 Cs concentration at (a, c) Namie and (b, d) Tsukuba.The shaded areas of (a) and (b) indicate the range of 137 Cs concentrations obtained when the simulations were run using the optimized parameter ranges E c = 0.03-0.05and v d = 0.05-0.1 cm s −1 .

Figure 12 .
Figure 12.Monthly mean (a) observed mean, observed median and simulated Namie to Tsukuba 137 Cs concentration ratio and (b) simulated re-suspension source area (where the observed deposition amount > 300 kBq m −2 ) contributions to 137 Cs air concentration at Namie and Tsukuba.The deposition parameters of the simulation in the figure are E c = 0.04 and v d = 0.1 cm s −1 .

Figure 13 .
Figure 13.48 h forward trajectory (statistical locations of LPs within 1 km AGL) predicted by the LM starting at (a) 12:00 LT and (c) 15:00 LT on 14 August and (b) 09:00 LT and (d) 15:00 LT on 19 August.Blue lines indicate median locations of LPs at 1 h and 1-day intervals as blue and red dots, respectively.Sky-blue lines indicate 17th and 83rd percentile locations of LPs and red dashed circles indicate areas containing 66 % of the LPs.

Figure 14 .
Figure 14.Daily mean gamma dose rate at the six monitoring sites: (a; red) Okuma and (orange) Namie (Omaru district) in the coastal area, (b; blue) village of Iitate and (sky blue) town of Kawamata in the Abukuma highland area, and (c; green) cities of Fukushima and (greenish yellow) Shirakawa in the Nakadori valley area of Fukushima prefecture as presented in Fig. 1b for two fiscal years (from April 2012 to March 2014).The first-order decreasing rates of the least-squares approximation over the period of no snow cover (May to October, 2012 and 2013) is also presented, along with the radioactive decay ( 134 Cs plus 137 Cs, by assuming the same amount of activity on 12 March 2011) in April 2012 and March 2013.

E
con (x, y) = |C sens (x, y) − C ref (x, y)| C ref (x, y) , (B1) E dep (x, y) = |D sens (x, y) − D ref (x, y)| D ref (x, y) , (B2) where x and y indicate grid points on the longitudinal and latitudinal axes, respectively.C sens and D sens indicate temporal mean surface concentrations (Bq m −3 ) and temporal cumulative depositions (Bq m −2 ) of the sensitivity runs, respectively.C ref and D ref are the same as C sens and D sens but for the reference run.E con and E dep were sampled only at grids where C ref (x, y) and D ref (x, y) are greater than their areal mean values, respectively.

Table 1 .
The observation sites and monitoring posts used to provide data for this study.Original location (now moved to city of Nihonmatsu). 2 Ottozawa San-ku Chiku Shukaijo (in Japanese).3OmaruTamokutekiShukaijo(in Japanese).4Originallocation (now moved to city of Fukushima). 5 Yamakiya Otsu Hachi-ku Community Shoubou Center (in Japanese). 1

Table 2 .
137Cs budget and statistical analysis for the comparison of observed and simulated deposition data for March 2011.

Table B1 .
Statistical measures of temporal mean 137 Cs surface concentration (E con ) and cumulative deposition (E dep ) of the sensitivity runs against the reference run for (top) the point source case and (bottom) the areal emission case.Dry deposition velocity over land; see Eq. (A3).c Relative errors of temporal mean surface concentration at each grid cell of the sensitivity run to that of the reference run; see Eq. (B1).d Same as E con but for cumulative deposition; see Eq. (B2).