Effects of species and geo-information on the 137Cs concentrations in edible wild mushrooms and plants collected by residents after the Fukushima nuclear accident

After the accident at the Fukushima Daiichi Nuclear Power Plant (FDNPP), much of the wild and edible mushrooms and plants in the surrounding areas were contaminated with radiocesium (137Cs). To elucidate their concentration characteristics, we analyzed 137Cs radioactivity data in edible forest products brought in for food inspection by the residents of Kawauchi Village, 12–30 km away from the FDNPP, from 2012 to 2019. A Bayesian model to estimate 137Cs concentration was constructed. Parameters of the normalized concentration of species (NCsp) for mushrooms were similar to those of the same species obtained in a previous study. Although NCsp values were highly varied among species, mycorrhizal mushrooms tended to have high NCsp values, followed by saprotrophic mushrooms, and wild edible plants values were low. Also, half of mycorrhizal mushroom species (8 of 16) showed an increasing trend in concentration with time; however, saprotrophic mushrooms and wild plants generally demonstrated a decreasing trend (22 of 24). The model considering the sub-village location information decreased the error of individual samples by 40% compared to the model not considering any location information, indicating that the detailed geo-information improved estimation accuracy. Our results indicate that the radioactivity data from samples collected by local residents can be used to accurately assess internal exposure to radiation due to self-consumption of contaminated wild mushrooms and plants.


Results
Estimation results of the two models. From the food inspection results in Kawauchi Village,137 Cs concentration data of wild mushrooms and plants for a total of 4534 specimens of 40 species taken from May 2012 to March 2020 were obtained (Table 1). Among these data, 1834 specimens had measured values below the detection limit. Tables S1 and S2 show the data of specimens and sub-village, respectively. In the analysis, two models ("sub-village model" and "whole village model") were used. Both models considered the species, date of collection, and 137 Cs deposition amount as explanatory variables. NC sp and λ sp are the concentration parameters of species standardized by deposition and slope parameter with time, respectively. Furthermore, the sub-village model hypothesized that each sub-village shown in Fig. 1b would have effects on sample 137 Cs concentration (defined as r loc ) dependent to the deposition (D loc ). The whole village model did not consider sample location in the village and used geometric mean of 137 Cs deposition for whole village (D village ) in the analysis.  26 was shown for each sub-village. The sub-village boundaries were created from a hand-drawn map obtained from the administrative office. The 137 Cs deposition data in the plotted areas were based on the results of the third airborne monitoring survey performed by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) from May 31 to 2 July 2011 28 . These maps were created using R version 4.1.0 29  We obtained posterior results of each parameter by Bayesian estimation. Using these parameters, we estimated the 137 Cs concentration of each specimen based on models and compared it with the measured value (Fig. 2). In Fig. 2a, the estimated median from the sub-village model was plotted, and results of whole village model were shown in Fig. 2b. For additional comparison, the plots were separately shown for mycorrhizal mushrooms, saprotrophic mushrooms, and wild plants. Among mushroom species, mycorrhizal fungi form symbiotic relationships with plants and uptake carbon via plant roots, while saprotrophic fungi decompose dead or dying organic matter for their nutrition. In both models, the estimated and measured values were distributed around the 1:1 line. The deviation from the 1:1 line was slightly smaller in the sub-village model. Among detected samples (blue symbols), 71% (1926/2700) and 67% (1801/2700) of the estimated values from the sub-village model and the whole village model were within ten times the range of the measured values (yellow belts in Fig. 2), respectively. The difference was significantly different (p < 0.001, chi-squared test, Table S3). Table 2 shows the posterior distributions of individual error (σ) and the hierarchical parameters. The sigma (σ), which indicates the magnitude of the deviation from the mean estimated by the model (μ sample ), was 0.51 for the sub-village model and 0.57 for the whole village model. The 95% confidence intervals for σ did not overlap between the two models (0.50-0.53 in the sub-village model and 0.55-0.59 in the whole village model); this indicates that the variability of the concentration becomes smaller when the sampling points can be specified at the sub-village level. The variability parameters of the concentration characteristics among species (σ sp ) were 0.81 and 0.83, which were larger than the individual error.

Group
Subgroup Species  Table 1. Posterior distribution of normalized concentration (NC sp ), temporal change parameter (λ sp ), and effective annual change ratio (R eff_sp ) of each species. 1) P. squarrosa is reported as poisonous. Also, it may contain similar edible species. 2) Among bamboo shoots, P. heterocycla (moso-bamboo) is the most popular in Japan but may contain multiple species. 3   www.nature.com/scientificreports/ Characteristics of radiocesium concentration in mushroom and edible plant species. Table 1 shows the calculated posterior results for concentration (NC sp ) and temporal trends (λ sp ) for each species. Here, the effective annual change ratio (R eff_sp ) and half-lives (doubling time) were calculated using λ sp . The median and mean values of NC sp differed by more than 100-fold (up to 845 times between Cortinarius salor (C. iodes) and Petasites japonicus (leaf stalk)) among species. NC sp of mushroom species in this study was compared with NC sp of previous research 22 in Fig. 3, and with Cs ess from Tagami and Uchida 31 in Supplementary Fig. S1. Both comparisons showed a positive and high correlation coefficient (R 2 = 0.713, 0.764, respectively). Figure 4 shows the median posterior NC sp and R eff_sp values by type presented as boxplots. The NC sp values of mycorrhizal mushrooms (Fig. 4a) are the highest, followed by saprotrophic mushrooms, and the lowest values are for wild plants. However, even among the same types, there was a large variation in the NC sp values. For example, in mycorrhizal fungi, the median NC sp values ranged 189-fold between C. salor and hygroscopic earthstar (Astraeus hygrometricus). Most wild plants had NC sp values of < 10 −3 , but only koshiabura (Eleutherococcus sciadophylloides) had significantly higher values (3.0 × 10 −3 m 2 /kg, FM). In addition, stalk and scape of Japanese butterbur (Petasites japonicus) were separately analyzed, and the NC sp of butterbur scape was 3.2 times higher than the stalk.
The R eff_sp values of mycorrhizal mushrooms (Fig. 4b) are higher than the other types, with half of medians (8 of 16) greater than one. In particular, the R eff_sp values of Sarcodon aspratus ("Kotake" in Japanese) are clearly greater than 1 (the median and lower 95% confidence interval value are 1.21 and 1.16, respectively, Table 1). Meanwhile, the R eff_sp values of most saprotrophic mushrooms and all wild plants are less than 1 (22 of 24), indicating a decreasing trend of 137 Cs concentrations with time.
The characteristics of the regional parameter. The spatial distribution of the median r loc of the regional parameters affecting the radioactivity of samples is shown in Fig. 5a. The r loc tended to be higher in the village's central and western part. However, the area bordering the northwestern village boundary and the south-  www.nature.com/scientificreports/ eastern area had lower r loc . Because the model was expressed as a sum of normalized logarithmic parameters, the value of r loc to the 10th power (10 rloc ), which ranged 3.61 times as maximum, had the effect of adjusting the concentration characteristics of each sub-village from the amount of 137 Cs deposition (D loc ). When estimated r loc and the geometric mean of 137 Cs deposition of each sub-village was plotted, there was a significant negative correlation (R 2 = 0.253, p < 0.001, Fig. 5b). In Fig S2, we show a map of the product of 137 Cs deposition (D loc ) and r loc to the 10th power (D loc × 10 rloc ) for each sub-village. The multiplied values indicate the adjusted effect on the sample 137 Cs concentration at the sub-villages. Although this figure is similar to Fig. 1b, the values increased at the central area, and the difference of the colored value is slightly reduced between sub-villages; the maximum ranges of the adjusted concentration effect (D loc × 10 rloc , Fig. S2) and 137 Cs deposition (D loc , Fig. 1b) are 10.5 and 11.7, respectively.

Discussion
The radiocesium levels of mushroom species in our study agree with those of previous studies (Fig. 3, Fig. S1). Komatsu et al. 22 analyzed the food monitoring data of wild mushrooms obtained at various municipalities in eastern Japan after the FDNPP accident and evaluated NC sp as in this research. Tagami and Uchida 31 analyzed radiocesium concentration of wild mushrooms from previous literature and monitoring data and evaluated 137 Cs in mushrooms originating before the FDNPP accident, such as nuclear bomb tests. Each of the two results is a generalization of radiocesium concentration characteristics in wild mushrooms, but it was necessary to confirm  www.nature.com/scientificreports/ the results' universality using another data set. Because the results obtained in one small area (Kawauchi Village) demonstrated a high positive correlation with the results of data collected over a wide area, it was confirmed that each wild mushroom species' concentration characteristics exist systematically in Japan. Although the concentration characteristics of wild plants were not validated sufficiently in this research, koshiabura, which has the highest NC sp among wild mushrooms, has the largest number of municipalities with shipping restrictions (113 in March 2021), while taranome (Aralia elata) and Asian royal fern (Osmunda japonica), which have a moderately higher NC sp of 10 -3 , have the second largest number of restricted municipalities after koshiabra (44 and 22, respectively) 10 . Therefore, there is a correspondence between the concentration characteristics and the number of cities, towns, and villages with shipping restrictions. As with wild mushrooms, each species of wild plants' concentration characteristics is universal. A comparison of the trends of NC sp and R eff_sp values of the ecological groups shows clear group differences in these parameters (Fig. 4, Table 1). The NC sp values reflect the 137 Cs concentration trend for each species, normalized for the effects of deposition, sampling location (sub-village in this research), and sampling date (year). The fact that the radiocesium contents in mycorrhizal mushrooms are generally higher than saprotrophic mushrooms has been widely accepted by research conducted after both Fukushima and Chernobyl nuclear power plant accidents 22,32,33 . Yoshida and Muramatsu 34 hypothesizes that differences in their mycelium distribution zones would explain differences in radiocesium concentrations among ecological types of mushroom species. Mycorrhizal fungi generally develop their mycelium in soil, while the saprotrophic family grows mycelium on dead wood and fallen leaves for nutrient absorption. It is known that most of the radiocesium is transferred to the surface layer of mineral soil with time after the accident. In a survey conducted by Imamura et al. 5 in a cedar forest in Kawauchi Village, 30.7%, 31.0% and 24.7% of 137 Cs was distributed in trees, organic layer, and the mineral soil layer in the summer of 2011, respectively. However, the 137 Cs proportions in the trees and organic layer decreased to 4.0% and 7.8%, respectively, whereas those of mineral soil increased significantly to 88.2%. In mineral soils of forest ecosystems, 137 Cs has been demonstrated to remain in the surface (0-5 cm depth) mineral surface layer (65%-91% of the total 137 Cs in the mineral soils). The hypothesis that radiocesium concentration in mycelium's habitat affects the mushrooms' concentration could explain the differences in concentration between types. Future research is required to quantitatively evaluate the relationship between radiocesium concentration and the mycelium distribution among species and types.
The higher radiocesium contents in mushrooms than in wild plants are thought to be due to differences in their ability to absorb nutrients. Yoshida and Muramatsu 35 compared the concentration of various metal elements in fungi and plants with the radiocesium concentration, and they found that mushrooms had higher potassium (K) concentrations and radioactive and stable cesium ( 137 Cs and 133 Cs) concentrations than plant leaves. However, radiocesium concentrations are variable even among wild plant species, and for example, Eleutherococcus sciadophylloides ("Koshiabura" in Japanese) is close to the values of mycorrhizal mushrooms (Table 1). Kiyono and Akama have shown that radiocesium concentration in koshiabura is related to radiocesium deposition in the litter 19,36 . These results indicate that koshiabura grows its roots at the surface and can absorb radiocesium more efficiently. Yamaji et al. 37 suggested that siderophores produced by endophytic root bacteria in koshiabra promote absorption by desorption of 137 Cs in the rhizosphere. It would be desirable to clarify the mechanisms causing intraspecific variation in concentration referencing the present results.
Among species, the time-dependent parameter, R eff_sp , also showed different trends between ecological groups (Fig. 4b, Table 1). The 137 Cs concentrations in samples were decay-corrected at the time of the accident (March 11th, 2011) and the physical decay of 137 Cs (30.17 years) was not accounted for in the λ sp calculation. On the other hand, R eff_sp involves physical decay of 137 Cs, and R eff_sp values greater or less than one indicate increasing or decreasing trends, respectively. Therefore, the trends indicate that the 137 Cs concentrations in the half of mycorrhizal species are increasing, while those of saprotrophic mushrooms and wild plants are decreasing. Some of the differences in the annual trends among these groups can be explained by the temporal changes in radiocesium concentrations in their respective substrates. The increasing 137 Cs concentrations of mycorrhizal fungi corresponds to the increasing concentration of 137 Cs in the soil, and the decreasing 137 Cs concentrations of saprotrophic fungi corresponds to the decreasing concentrations of 137 Cs in litter and dead wood 5 . In particular, kotake mushrooms demonstrated a clear increasing trend. Yamada et al. reported the possibility of accumulation of 137 Cs in mycelium because many of the mycorrhizal mushrooms with high concentrations retained a large amount of 137 Cs originating from events prior to the Fukushima accident 38 . It is necessary to confirm whether such retention of 137 Cs in mycelium is actually occurring.
On the other hand, the decreasing 137 Cs concentration in wild plants that uptake nutrients from mineral soils cannot be explained by the increasing levels of 137 Cs in their substrates. Manaka et al. 39 shows that the proportions of exchangeable radiocesium in organic layer and mineral soil layer were high in 2011 (6-10%), immediately after the accident 2011 and then decreased (2-4%). The distribution of radiocesium in the forest changes over time, and the accessibility of radiocesium changes too. Experiments demonstrated that mycelium has the ability to solubilize Cs adsorbed in soil 40 . However, plants are not able to solubilize such Cs adsorbed in soil; therefore, the decrease in accessibility of radiocesium over time in forest soils was considered to be responsible for decrease in the concentration of wild plants.
The concentration characteristics of the species varied greatly ( Table 1), confirming that it is extremely important to consider the effect of species in estimating the concentration of wild mushrooms and plants. On the other hand, from comparison results, more detailed information of sampling points would improve the estimation accuracy of radioactivity in wild mushrooms and plants when collection points could be obtained ( Fig. 2 and Table 2). In Fig. S4, the normalized concentration, which was calculated by dividing the 137 Cs concentration of wild mushrooms and plants by the amount of 137 Cs deposition, is colored according to the r loc values of the sampled sub-villages. This visualization shows that samples collected at sub-villages with higher r loc values tended to exhibit higher normalized concentrations within a given species. The 95% confidence intervals for www.nature.com/scientificreports/ 137 Cs concentration in the samples from the normal distribution and the median of posterior σ ( Table 2) was estimated to be a factor of 103 (= 10 0.51×2×1.96 ) for the sub-village model and a factor of 171 (= 10 0.57×2×1.96 ) for the whole-village model. This result simply indicates that information of sample collection site in each village would reduce sample variability by 40%, assuming that the concentration varies among sub-villages according to the deposition (D loc ) adjusted by r loc . Such an improvement in estimation accuracy is expected to contribute to the evaluation of the internal dose according to ingestion and collection behavior. However, the sub-village model, which takes into account species and site effects, still showed large errors. This result suggests variation in concentration among specimens of the same species collected within a few square kilometers of a sub-village. Previous studies have shown that radiocesium concentrations in wild mushrooms of a given species can fluctuate within an even smaller area. Dahlberg et al. 41 investigated intraspecific variability of 137 Cs concentrations of one mushroom species (Suillus variegatus) from seven sites within 1 km 2 area. Almost 40% of the variation in 137 Cs levels was accounted for by the sites, and the other 60% was accounted for within sites. Komatsu et al. 42 investigated 137 Cs concentration of Boletus hiratsukae at two plots about 300 m apart and showed that the geometric mean of the mushroom concentrations at the two sites differed by a factor of 3.5, and that even within a site of about 200 m 2 the mushroom concentrations differed by a factor of nine. Location information of greater accuracy is expected to reduce the error in estimating the concentration of specimens. However, there may be a limit to the use of location information alone because it is affected by variation in the amount of radiocesium deposited in the forest and the environment. The factors governing variation in concentration within a given forest must be explained.
The location parameters affecting concentration (r loc ) are spatially biased (Fig. 5a). Radiocesium concentration in wild mushrooms and trees is known to be affected by deposition and by soil properties such as exchangeable K concentration, pH, and clay mineral content [42][43][44][45] . Therefore, one possibility is that the characteristics of the soil may have caused the deviation between the r loc concentration and the amount of deposition. However, r loc was higher in the central area where deposition was lower, and there was a significant negative correlation between log 10 D loc and r loc (Fig. 5b). One possible explanation for this relationship is that the deposition across the forest was not sufficiently evaluated by aircraft monitoring. The central area is relatively flat along valley floors, and it has farmland and residential areas. In general, radioactive cesium deposited on farmland and residential areas is more likely to move than forests 46 . Therefore, the deposition of 137 Cs on farmlands and residential areas is usually lower than that on nearby forests. The 137 Cs deposition data used in this research have a 250-m resolution, having been measured from helicopters piloted at an altitude of approximately 300 m 47 . Therefore, the amount of cesium deposited in forests in the central area may have been underestimated due to the mixing effect of other land covers, and this may have been recognized as a difference in r loc .
This study has data quality limitation because of the collection of the data by residents. The data may include some ambiguous reporting; for example, some locations of collection sites may actually be the collector's residence. This is because it is sometimes difficult for the collectors to report the exact location of collection within the vast forested areas. In addition, for expensive mushrooms such as matsutake (Tricholoma matsutake), which can fetch some hundreds of dollars per kilogram, there is a risk that collection locations will not be accurately reported. If samples with high radiocesium activity were collected at the outer side of the village (with high amounts of radiocesium deposition) but were reported as having been collected in the central area (with rates of lower deposition), then the r loc value of the central area will be overestimated as shown in Fig. 5b. However, Bayesian analysis of the data set that excluded matsutake shows little change in the negative correlation between r loc and D loc (Supplementary Fig. S3). Thus, uncertainty in the sampling location did not have a significant effect on r loc . There was a large number and variety of samples distributed widely throughout the village, indicating that our results are robust. Therefore, the data originating from samples obtained by the residents are worth analyzing and are likely to provide valuable research conclusions.

Conclusion
We showed the radiocesium concentration's trend in wild mushrooms and wild edible plants using inspection data of food. For mushrooms, since the values were close to those of previous studies, the estimated values of species were considered to be widely applicable, at least in Japan. In addition, this study shows that the data from samples collected by residents can be used to analyze the concentration characteristics of species and regions. In the future, it will be possible to use such data for internal exposure assessment by collecting wild food consumption data that takes into account regional and individual differences. Furthermore, such results may also encourage sustained sampling efforts by citizens. The inspection of food brought in by residents is already widely conducted at municipalities in Fukushima Prefecture. Should this continue, then we expect further development of participatory research by citizens.
It is said that radiocesium in forests approaches an equilibrium state with time 13 . However, the radiocesium concentrations in some species, such as Sarcodon aspratus (one of the popular and economically valuable edible mushrooms in the area) are still trending upward and are not in equilibrium. It is necessary to continue monitoring to show the prospect of concentration changes.

Materials and methods
Site information. We collected radioactivity data of wild mushrooms and wild edible plants from inspection results of specimens brought in by residents in Kawauchi Village, which is located 12-30 km away from the FDNPP (Fig. 1). Kawauchi Village is considered small, with an area of 197.4 km 2 , and a population of about 2500 (2820 in 2010 and 2518 in 2021) 48 . It is located in the middle of the Abukuma Highlands, where the elevation ranges from 270 to 1,192 m above the sea level. It has a forest coverage of 89.0%, which is higher than the average for Fukushima Prefecture (71%) and Japan as a whole (69%) 49  www.nature.com/scientificreports/ 960 kBq/m 2 in 2011, estimated from an aircraft monitoring 28 . Before the accident, its residents were accustomed to gathering wild foods, such as wild edible mushrooms, plants, mammals, and wild honey 50 ; many have been brought in for inspection. Information on collection areas of sub-village levels, called "Ko-aza" in Japanese, is also recorded. For these reasons, we thought that the data of the brought in inspection in Kawauchi Village would possess high value as data for inter-species and inter-region analysis on the wild mushrooms and edible plants' radioactivity concentrations.
Radioactivity data of mushrooms and wild plants. Fukushima Table S4 shows the specifications of the machines 51-53 . All instruments have been confirmed to meet the radiocesium screening method requirements for food 53 . Among these machines, FD-08Cs1000-1 can measure radioactivity non-destructively, and the others conduct destructive measurements. The sample weight is approximately 500 g, and the counting time is 30 min. FD-08Cs1000-1 outputs the summed concentration of the two radiocesium nuclides ( 134 Cs and 137 Cs), and its detection limit is 10 Bq/kg (for total 134 Cs + 137 Cs). Each of the other three machines separately outputs the concentrations of 134 Cs and 137 Cs, and the detection limit is 10 Bq/kg for each radionuclide. Energy calibrations and background checks were performed daily, and the accuracy was periodically verified with brown rice whose radiocesium concentration was verified by calibrated high-purity Germanium (HPGe) detectors installed in the Fukushima Nuclear Center 49 . Table S4 shows the results of quality control using brown rice.
Data preparation of radioactivity of samples. From the radioactivity data of wild mushrooms and plants, we picked up data that met the following criteria; • Data have information of sampling location at sub-village levels • Items that are not confirmed to be cooked products such as "boiled" or "dried. " • Species with more than ten samples in which radiocesium was detected.
In cases where mushrooms and wild plants were given in dialects, we confirmed the species' names with residents. The names of the species were determined from the Japanese names of the items, but in some cases, it was not possible to distinguish between Cortinarius salor ("Murasakiaburashimejimodoki" in Japanese) and C. iodes ("Murasakiaburashimeji"), considered to be closely related species, so the two were mixed for analysis. The leaf stalk and scape of Petasites japonicus (Japanese butterbur) are called "Fuki" and "Fukinotou" in Japanese, respectively, and are registered separately. Therefore, despite being the same species, they were distinguished in the analysis. In this data, there were not sampling date but measurement date. Therefore, the date of measurement and sample collection were assumed to be the same.
The 137 Cs concentration results were used in the model analysis. The reason for not using the 134 Cs concentration among the measured values is explained in the subsection of "Bayesian estimation". 137 Cs concentrations were decay-corrected to March 11th, 2011 for comparison with Komatsu et al. (2019). Based on the assumption that the 134 Cs/ 137 Cs ratio at the time of the accident was one 54 , the summed concentration of 134 Cs and 137 Cs concentration taken by FD08-Cs1000-1 was converted to a 137 Cs concentration, which was decay-corrected to March 11th, 2011, using the following equation; In this equation, dy indicates the period from March 11th, 2011, to the date of measuring, and it is expressed by decimal years.

Sub-village ("Ko-aza") boundary map of Kawauchi Village. Kawauchi Village comprises eight
administrative communities (called "Oh-aza" in Japanese), which are further subdivided into small administrative units known as "Ko-aza". Here, we refer to these small administrative units as sub-villages. We obtained a sub-village map from the administrative office. The printed map was originally drawn by hand and had been used for village administration. To create a polygon shapefile of the map, we digitized it by scanning, geo-rectifying, and digitizing using GIS software in TNTmips v2014 (MicroImages, Inc, NE) and ArcGIS  where D indicates the radiocesium deposition amount based on the aircraft monitoring. Then the above equation is transformed and logarithmized to yield; In this expression of the model equation, we further assumed that the logartihm of NC encompassed the summed effects of species identity, collection date, and collection site, and that the logarithm of NC was normally distributed around the estimated mean as per the following equations; where NC sp , λ sp , D loc and r loc indicate characteristics of concentration of species, temporal trends of species, 137 Cs deposition of each sub-village area and effects of sub-village on concentration, respectively. r loc is a parameter with zero mean that represents the deviation of the concentration effect from the expected value based on the deposition (D loc ) value at the point of collection. These parameters except D loc were obtained from hierarchically sampled from normal distribution with hierarchical parameters (μ sp , σ sp , μ λsp , σ λsp , σ loc ). Additionally, r loc was sampled using the Intrinsic Conditional Auto-Regressive (Intrinsic CAR) model 56 , which is one of the models considering spatial auto-correlation. For samples whose measured radiocesium concentrations were below the detection limit, radiocesium concentration values were estimated by a censoring distribution in which the detection limit was treated as the upper bound 57 . This model was defined as the "sub-village model" for this research. This model is similar to model 6 in Komatsu et al. 22 but differs in that their previous model takes into account 134 Cs values and differences between 134 and 137 Cs values. Komatsu et al. evaluated the regional trend in the difference between 134 Cs and 137 Cs concentrations across eastern Japan because 137 Cs originating from nuclear bomb tests before the FDNPP accident was detected in wild mushrooms sampled in the northern and southern parts of eastern Japan, which are far from the FDNPP and received less deposition from the accident (< 10 Bq/ m 2 according to aircraft monitoring). However, in Kawauchi Village, the amount of 137 Cs deposition because of the FDNPP accident (42-960 Bq/m 2 in Kawauchi Village 26 ) was larger than the amount of 137 Cs deposition from nuclear tests (< 10 kBq/m 2 in Japan 58 ), and it is difficult to evaluate differences between 137 Cs and 134 Cs concentrations of each specimen. Therefore, our analysis used only 137 Cs values given the longer half-life of that isotope.
Alternatively, for the comparison, we also used the "whole village model, " in which r loc was not taken into account and the geometric mean of the 137 Cs deposition in the village (D village ) was used for the calculation.
For the Bayesian estimation, we used WinBUGS ver 1.4.3 57 . In the calculation, we set burn-in, burn-out, chain, and thin as 10,000, 20,000, 3, and 10, respectively. As a result, we gained 3000 posterior results for each parameter.
The effective annual change ratio according to species (R eff_sp ) was calculated based on the radiological decay of 137 Cs (half-life: 30.17 years) and using λ sp as follows: Limitation. The date of collection of each sample was not registered, and the substitution of the measurement date for the collection date may cause uncertainty in the analysis of temporal trends. However, the measurement dates were reported in a concentrated manner for each species (Fig. S5), and we assume that the collection and measurement dates were generally close. Considering the small annual rate of change (< 20%) and the length of the sampling period, it is unlikely that a substantial bias affected the conclusions.
It has been pointed out that non-destructive testing machines cannot fill samples homogeneously and can sometimes introduce a bias of > 20% vs. HPGe measurements 52 . However, in the case of wild mushrooms and plants, the variation among species and samples is substantial, on the one to two orders of magnitude. Therefore, we consider the influence of errors because of the measurement equipment to be dwarfed by other factors.