Meta-analysis of radiocesium contamination data in Japanese cedar and cypress forests over the period 2011-2017

Long-term behaviour of Cs-137, Cs-133 and K in beech trees of French forests

In the long-term after atmospheric deposit onto a forest ecosystem, Cs-137 becomes incorporated into the biogeochemical cycle of stable elements and progressively reaches a quasi-equilibrium state. This study aimed at determining to what extent Cs-137 activity distribution in tree vegetation could be predicted from that of stable caesium (Cs-133) and potassium (K), which are known to be stable chemical analogues and competitors for Cs-137 intake in tree organs. Field campaigns that focused on beech trees (Fagus sylvatica L.) were conducted in 2021 in three French forest stands with contrasted characteristics regarding either the contribution of global vs. Chornobyl fallouts, soil or climatic conditions. Decades after Cs-137 fallouts, it was found that more than 80% of the total radioactive inventory in the system remained confined in the top 20 cm mineral layers, while organic layers and beech vegetation (including roots) contributed each to less than 1.5%. The enhanced downward migration of Cs-137 in cambisol than podzol forest sites was presumably due to migration of clay particles and bioturbation. The distribution of Cs-137 and Cs-133 inventories in beech trees was very similar among sites but differed from that of K due a higher accumulation of Cs isotopes in roots (40-50% vs. < 25% for K). The aggregated transfer factor (Tag) of Cs-137 calculated for aerial beech organs were all lower than those reported in literature more than 20 years ago, this suggesting a decrease of bioavailability in soil due to ageing processes. Regarding their variability, Tags were generally lower by a factor 5 at the cambisol site, which was fairly well explained by a much higher value of RIP (radiocesium immobilisation potential). Cs-137 concentrations in trees organs normalized by the soil exchangeable fractions were linearly correlated to those of Cs-133 and the best fit was found for the linear regression model without intercept indicating that no more contribution of the foliar uptake could be observed on long term. Provided that the vertical distribution of caesium concentrations and fine root density are properly measured or estimated, Cs-133 was shown to be a much better proxy than K to estimate the root transfer of Cs-137.

Estimation of spatio-temporal distribution of 137Cs concentrations in litter layer of forest ecosystems in Fukushima using FoRothCs model

The nuclear power plant accident in Fukushima had led to pollution of forest ecosystems with ¹³⁷Cs in 2011. In this study, we simulated the spatiotemporal distribution of ¹³⁷Cs concentrations of litter layer in the contaminated forest ecosystems in two decades from 2011, which is one of the key environmental components of ¹³⁷Cs migration in the environment due to the high bioavailability of ¹³⁷Cs in the litter. Our simulations showed that ¹³⁷Cs deposition is the most important factor in the degree of contamination of the litter layer but vegetation type (evergreen coniferous/deciduous broadleaf) and mean annual temperature are also important for changes over time. Deciduous broadleaf trees had higher initial concentrations in the litter layer due to the direct initial deposition on the forest floor. However, the concentrations remained higher than those in evergreen conifers after 10 years due to redistribution of ¹³⁷Cs by vegetation. Moreover, areas with lower average annual temperatures and lower litter decomposition activity retained higher ¹³⁷Cs concentrations in the litter layer. The results of the spatiotemporal distribution estimation of the radioecological model suggest that, in addition to ¹³⁷Cs deposition, elevation and vegetation distribution should also be considered in the long-term management of contaminated watersheds, which can be informative in identifying hotspots of ¹³⁷Cs contamination on a long-term scale.

Mode of Atmospheric Deposition in Forests Demonstrates Notable Differences in Initial Radiocesium Behavior

The March 2011 Fukushima Dai-ichi Nuclear Power Plant accident in Japan released 520 PBq of radionuclides compared to a total release of 5300 PBq from the Chornobyl Nuclear Power Plant accident. Both nuclear accidents resulted in deposition of radiocesium throughout the northern hemisphere, and a plethora of studies have been performed regarding radiocesium (¹³⁷Cs) behavior. However, few studies have assessed the impact of precipitation on ¹³⁷Cs deposition in forests. Wide-scale environmental measurements from 2011 and 2016 were used to determine the differences in ¹³⁷Cs deposition because of precipitation following the Fukushima accident. In areas where wet deposition processes were dominant, dense forests generally had lower ambient dose rates and levels of contamination on forest floors than other stands with fewer stems per hectare in 2011. Similar tendencies were not observed in areas that were primarily subject to dry deposition nor were any trends observed in 2016. ¹³⁷Cs was retained in dense forest canopies for an extended period regardless of the deposition mode. Additionally, it was found that the initial retention of radionuclides by forest canopies is in general higher for areas with predominantly dry deposition. Incorporation of radiocesium into wood tissues was the same for both wet and dry deposition.

Late phase radiocesium dynamics in Fukushima forests post deposition

The long term dynamics of radiocesium in typical forest ecosystems was studied in the radioactive contaminated areas in Fukushima Prefecture. Six observations sites located in Yamakiya Village (Kawamata Town; since 2014), Tsushima Village (Namie Town, since 2015), and Tomioka Town (since 2017) were monitored. The forests consisted of artificial plantations of Japanese cedar (Cryptomeria japonica) at Yamakiya Village, Tsushima Village, and Tomioka Town. Tsushima Village also had a natural mixed forest dominated by Japanese red pine (Pinus densiflora), and Tomioka Town had a young and a mature artificial plantation of Japanese cypress (Chamaecyparis obtuse). Concentrations of ¹³⁷Cs were monitored in the samples collected from the main aboveground biomass compartments, fresh litterfall, forest litter, and soil. Concentrations of exchangeable forms of ¹³⁷Cs and stable K were measured in soil samples. During the observation period, the litter radiocesium inventories at all sites decreased significantly to approximately 1% or less of the total ground deposition. Approximately 80% of the total radiocesium inventory is localized in the upper 5-cm layer of soil and there is little downward migration of radiocesium. At the sites with the longest monitoring series (Yamakiya and Tsushima), the radiocesium expectation depths and expectation mass depths were relatively constant at 2-3 cm and 5-6 kg m^-2, respectively. Aboveground biomass compartments showed similar decreasing trends in radiocesium aggregated transfer factors, T_ag, in the compartments that were exposed to atmospheric fallout in March 2011 (old foliage, small branches, and outer bark). The mean T_ag in cedar stand compartments currently are in the range of 10^-3-10^-2 m² kg^-1 dw. However, the mean T_ag and their dynamic trend significantly differed in the wood compartments of the cedar stands, which may indicate root uptake differences of orders of magnitude between observation sites. The difference in radiocesium concentration in wood between the sites becomes less pronounced when normalized by the ratio of exchangeable ¹³⁷Cs/K in the soils.

Contamination processes of tree components in Japanese forest ecosystems affected by the Fukushima Daiichi Nuclear Power Plant accident 137Cs fallout

In forests affected by the Fukushima Daiichi Nuclear Power Plant accident, trees became contaminated with ¹³⁷Cs. However, ¹³⁷Cs transfer processes determining tree contamination (particularly for stem wood, a prominent commercial resource) remain insufficiently understood. We propose a model for simulating dynamic behavior of ¹³⁷Cs in a forest tree-litter-soil system and applied it to contaminated forests of cedar plantation and natural oak stand in Fukushima to elucidate relative impact of distinct ¹³⁷Cs transfer processes determining the tree contamination. The transfer of ¹³⁷Cs to the trees occurred mostly (>99%) through surface uptake of ¹³⁷Cs trapped by needles and bark during the fallout. Root uptake of soil ¹³⁷Cs was several orders of magnitude lower than the surface uptake over a 50-year period following the accident. As a result, internal contamination of the trees proceeded through an enduring recycling (translocation) of ¹³⁷Cs absorbed on the tree surface. A significant surface uptake of ¹³⁷Cs through bark was suggested, contributing to 100% (leafless oak tree) and 30% (foliated cedar tree; the remaining uptake occurred at needles) of the total uptake by the trees, although that pathway still needs to be evaluated by experimental evidence. It was suggested that the activity concentration of ¹³⁷Cs in stem wood of the trees at these sites are currently (as of 2021) decreasing by ~3% per year, mainly through radioactive decay of ¹³⁷Cs and partly through dilution effect from tree growth. Although further refinement of the model is recommended, for example by including tree species specific ¹³⁷Cs transportation in stem, these findings provide vital information for planning of forestry reactivation in Fukushima; e.g., removal of forest floor organic layer will not reduce the tree contamination for a long term because of the ¹³⁷Cs absorption via the tree surface substantially greater than root uptake of ¹³⁷Cs deposited to the floor.

Evaluation of contribution rate of the infiltrated water collected using zero-tension lysimeter to the downward migration of 137Cs derived from the FDNPP accident in a cedar forest soil

The vertical distribution of ¹³⁷Cs in forest soil is important for predicting air dose rates and future cycling in forest ecosystems. However, there are many unexplained questions about the mechanisms of its downward migration. In this study, the ¹³⁷Cs flux by rainfall infiltration was observed for three years from August 2017 using zero-tension lysimeters in a mature cedar forest where monitoring of the vertical distribution of ¹³⁷Cs has been conducted since 2011. As a result, the ¹³⁷Cs concentration in infiltrated water through the litter layer, 5 cm and 10 cm showed a tendency to be high in summer, but no such seasonal variation was found at 20 cm. Although the ¹³⁷Cs inventory in the litter layer has been exponentially decreasing, the annual ¹³⁷Cs fluxes in infiltrated water through the litter layer were almost the same in three years, and about 0.14-0.17% of the deposition density of ¹³⁷Cs. Comparing these ¹³⁷Cs fluxes with the apparent amounts of downward migration of ¹³⁷Cs estimated from the change in the vertical distribution of ¹³⁷Cs, the contribution rate of the infiltrated water to downward migration of ¹³⁷Cs from litter to soil was calculated to be 8.5-17.7%. Similarly, the contribution rate in mineral soil layers was calculated to be 0.6-0.8% on a measured basis and estimated to be 3.0 ± 0.2% after correcting the amount of collected water, which is a problem with zero-tension lysimeter. It indicates that rainfall infiltration can explain a small part of the downward migration of ¹³⁷Cs, thus further studies are required to clarify the contribution rate of remaining mechanisms such as advection-diffusion, colloidal transport, physical mixing, bioturbation, and growth and death of plant roots.

Machine learning analysis of 137Cs contamination of terrestrial plants after the Fukushima accident using the random forest algorithm

Radioactive contamination of terrestrial plants was extensively investigated and quantitatively modeled after the Fukushima nuclear power plant accident. This phenomenon, which is important for ecosystem functioning and protection of human health, is influenced by multiple factors, including plant species, time after the accident, and climate. Machine learning algorithms such as random forests (RF) have a record of strong performance on large multi-dimensional data sets, but, to our knowledge, combined data on post-Fukushima plant contamination with radionuclides were not yet subjected to a machine learning analysis. Here we performed such analysis on two large published data sets: (1) ¹³⁷Cs activity concentrations in four common Japanese forest tree species. (2) Plant/soil ¹³⁷Cs concentration ratios in multiple perennial plant species. The goal was to show the usefulness of machine learning for identifying and quantifying the main trends of ¹³⁷Cs contamination in terrestrial plants. Each data set was split randomly into training and testing parts, RF was fitted and tuned on the training parts, and its performance was assessed on the testing parts by three metrics: coefficient of determination (R²), root mean squared error, and mean absolute error. Synthetic noise variables and the Boruta algorithm were used in a customized procedure to identify the most important predictor variables, which consistently outperformed random noise. Good agreement between observations and RF predictions (e.g. R²∼0.9 on testing data) was obtained on both data sets. The effects of the most important predictors (e.g. time after the accident, ¹³⁷Cs land contamination level, and plant species) and interactions between them were quantified by partial dependence plots. These results of machine learning analyses of large data collections can help to complement previous modeling efforts, and to clarify the patterns of ¹³⁷Cs contamination of plants after the Fukushima accident.

Dynamics of radiocaesium within forests in Fukushima-results and analysis of a model inter-comparison

Forests cover approximately 70% of the area contaminated by the Fukushima Daiichi Nuclear Power Plant accident in 2011. Following this severe contamination event, radiocaesium (¹³⁷Cs) is anticipated to circulate within these forest ecosystems for several decades. Since the accident, a number of models have been constructed to evaluate the past and future dynamics of ¹³⁷Cs in these forests. To explore the performance and uncertainties of these models we conducted a model inter-comparison exercise using Fukushima data. The main scenario addressed an evergreen needleleaf forest (cedar/cypress), which is the most common and commercially important forest type in Japan. We also tested the models with two forest management scenarios (decontamination by removal of soil surface litter and forest regeneration) and, furthermore, a deciduous broadleaf forest (konara oak) scenario as a preliminary modelling study of this type of forest. After appropriate calibration, the models reproduced the observed data reliably and the ranges of calculated trajectories were narrow in the early phase after the fallout. Successful model performances in the early phase were probably attributable to the availability of comprehensive data characterizing radiocaesium partitioning in the early phase. However, the envelope of the calculated model end points enlarged in long-term simulations over 50 years after the fallout. It is essential to continue repetitive verification/validation processes using decadal data for various forest types to improve the models and to update the forecasting capacity of the models.

Influence of non-equilibrium sorption on the vertical migration of 137Cs in forest mineral soils of Fukushima Prefecture

Sorption hypotheses and models are required for the prediction of ¹³⁷Cs migration in soils contaminated after nuclear reactor accidents and nuclear weapons tests. In assessment models, the K_d (distribution coefficient) hypothesis for sorption, which assumes that sorption is instantaneous, linear and reversible, has often been coupled with the convection-diffusion equation (CDE) to model ¹³⁷Cs migration. However, it fails to describe ¹³⁷Cs migration velocities which often decrease with time. Alternative equilibrium-kinetic (EK) hypotheses of ¹³⁷Cs sorption/desorption have been suggested by laboratory experiments but have not been fully validated in field conditions. This work addressed the influence and magnitude of non-equilibrium ¹³⁷Cs sorption in field conditions by reinterpreting, with an inverse approach, series of ¹³⁷Cs profiles measured in mineral soils of forest plots located in Fukushima Prefecture (2013-2018). Our results show that the inclusion of non-equilibrium sorption significantly improves, compared to the equilibrium hypothesis, the realism of simulated ¹³⁷Cs profiles. Fitted sorption parameters suggest a fast sorption kinetic (half-time of 1-7 h) and a pseudo-irreversible desorption rate (half-time of 3.2 × 10⁰-3.4 × 10⁶ years), whereas equilibrium sorption (4.0 × 10^-3 L kg^-1 on average) only affects a negligible portion of ¹³⁷Cs inventory. By June 2011, such EK parameters fitted on our plots realistically reproduced profiles measured in the same forest study site (Takahashi et al., 2015). Predictive modeling of ¹³⁷Cs profiles in soil suggests a strong persistence of the surface ¹³⁷Cs contamination by 2030, with exponential profiles consistent with those reported after the Chernobyl accident. This study demonstrates that hypotheses and parameters of ¹³⁷Cs sorption can be partially inferred from in situ measurements. However, further experiments in controlled conditions are required to better estimate the sorption parameters and to identify the processes behind non-equilibrium sorption.