Modeling dynamics of (137)Cs in forest surface environments: Application to a contaminated forest site near Fukushima and assessment of potential impacts of soil organic matter interactions

Adaptation of FDMT to include process-based models and consider regional aspects following radionuclide deposition events

Experience from earlier nuclear accidents has clearly shown the need for maintaining and developing appropriate modelling capabilities. Dealing with complex issues such as human exposure following a nuclear accident necessitates the implementation of a set of interconnected models such as FDMT. FDMT is an integrated module within the two main European decision support systems for radiological emergency preparedness, ARGOS and JRODOS, to simulate the transfer of radionuclides along terrestrial food chains and to predict their activity concentrations in foodstuffs. In order to make the module more fit-for-purpose, FDMT has been implemented in a new modelling platform (ECOLEGO) which provides a high degree of flexibility with regard to conducting developmental work. This paper presents improvements in FDMT further through either the incorporation of new models or further elaboration of existing ones, as well as updates in default parameters. Models have also been made more fit-for-purpose through consideration of regional-specific parameters. Specific improvements include modelling developments related to dry deposition, radioactive particle weathering, radiocaesium transfer influenced by soil characteristics and, for a region-specific case, animal uptake. In addition, the paper presents new pathways and parameters (and updated values) to be considered for making FDMT more adapted for Norwegian conditions. Overall, the improvements made in the present work should significantly reduce the uncertainties associated with the outputs of the FDMT models.

Investigation of dense plumes in porous media using CsNO3 solute and phenolphthalein tracer

We present results from the laboratory experiments performed to study the behaviour of dense plumes in porous media using Caesium Nitrate (CsNO3) as the solute. Specifically, we consider the case of fixed volume release of a dense fluid into a saturated porous medium of lighter density. Our experiments consider the injection volume and concentration of CsNO3 and the porous medium permeability as experimental variables. Our measurements show that the plume length and volume increase with time approximately as t0.5. The mean concentration of the plume decreases with time as plume monotonically dilutes due to continuous mixing with the surrounding. From our measurements and with the help of dimensionless scaling, we present empirical models for the length, volume and mean concentration of the plume. Our results may be used to gain preliminary insights into the spatio-temporal evolution of plumes formed in the subsurface environment.

Exploring simple ways to avoid collecting highly 137Cs-contaminated Aralia elata buds for the revival of local wild vegetable cultures

Collection and cooking of wild vegetables have provided seasonal enjoyments for Japanese local people as provisioning and cultural ecosystem services. However, the Fukushima Daiichi Nuclear Power Plant accident in March 2011 caused extensive radiocesium contamination of wild vegetables. Restrictions on commercial shipments of wild vegetables have been in place for the last 10 years. Some species, including buds of Aralia elata, are currently showing radiocesium concentrations both above and below the Japanese reference level for food (100 Bq/kg), implying that there are factors decreasing and increasing the 137Cs concentration. Here, we evaluated easy-to-measure environmental variables (dose rate at the soil surface, organic soil layer thickness, slope steepness, and presence/absence of decontamination practices) and the 137Cs concentrations of 40 A. elata buds at 38 locations in Fukushima Prefecture to provide helpful information on avoiding collecting highly contaminated buds. The 137Cs concentrations in A. elata buds ranged from 1 to 6,280 Bq/kg fresh weight and increased significantly with increases in the dose rate at the soil surface (0.10-6.50 μSv/h). Meanwhile, the 137Cs concentration in A. elata buds were not reduced by decontamination practices. These findings suggest that measuring the latest dose rate at the soil surface at the base of A. elata plants is a helpful way to avoid collecting buds with higher 137Cs concentrations and aid in the management of species in polluted regions.

Effect of soil organic matter on the fate of 137Cs vertical distribution in forest soils

Understanding the fate of the vertical distribution of radiocesium (¹³⁷Cs) in Japanese forest soils is key to assessing the radioecological consequences of the Fukushima Daiichi Nuclear Power Plant (FDNPP) accident. The ¹³⁷Cs behavior in mineral soil is known to be primarily governed by interaction with clay minerals; however, some observations suggest the role of soil organic matter (SOM) in enhancing the mobility of ¹³⁷Cs. Here we hypothesized that soil organic carbon (SOC) concentration profile determines the ultimate vertical pattern of ¹³⁷Cs distribution in Japanese forest soils. In testing this hypothesis, we obtained soil samples that were collected before the FDNPP accident at four forest sites with varying SOC concentration profiles and quantified the detailed vertical profile of ¹³⁷Cs inventory in the soils roughly half a century after global fallout in the early 1960 s. Results showed that the higher the SOC concentration in the soil profile, the deeper the ¹³⁷Cs downward penetration. On the basis of the data for surface soils (0-10 cm), the ¹³⁷Cs retention ratio for each of the 2-cm thick layers was evaluated as the ratio of ¹³⁷Cs inventory in the target soil layer to the total ¹³⁷Cs inventory in and below the soil layer. A negative correlation was found between the ratio and SOC concentration of the layer across all soils and depths. This indicates that the ultimate fate of ¹³⁷Cs vertical distribution can be predicted as a function of SOC concentration for Japanese forest soils, and provides further evidence for SOM effects on the mobility and bioavailability of ¹³⁷Cs in soils.

Soil dust and bioaerosols as potential sources for resuspended 137Cs occurring near the Fukushima Dai-ichi nuclear power plant

One of the current pathways to radiation exposure, caused by the radionuclides discharged during the Fukushima Dai-ichi Nuclear Power Plant (FDNPP) accident, is the inhalation of resuspended ¹³⁷Cs present in the air. Although wind-induced soil particle resuspension is recognized as a primary resuspension mechanism, studies regarding the aftermath of the FDNPP accident have suggested that bioaerosols can also be a potential source of atmospheric ¹³⁷Cs in rural areas, although the quantitative impact on the atmospheric ¹³⁷Cs concentration is still largely unknown. We propose a model for simulating the ¹³⁷Cs resuspension as soil particles and bioaerosols in the form of fungal spores, which are regarded as a potential candidate for the source of ¹³⁷Cs-bearing bioaerosol emission into the air. We apply the model to the difficult-to-return zone (DRZ) near the FDNPP to characterize the relative importance of the two resuspension mechanisms. Our model calculations show that soil particle resuspension is responsible for the surface-air ¹³⁷Cs observed during winter-spring but could not account for the higher ¹³⁷Cs concentrations observed in summer-autumn. Higher ¹³⁷Cs concentrations are reproduced by the emission of ¹³⁷Cs-bearing bioaerosols (fungal spores) that replenishes the low-level soil particle resuspension in summer-autumn. Our model results show that the accumulation of ¹³⁷Cs in fungal spores and large emissions of spores characteristic of the rural environment are likely responsible for the presence of biogenic ¹³⁷Cs in the air, although the former must be experimentally validated. These findings provide vital information for the assessment of the atmospheric ¹³⁷Cs concentration in the DRZ, as applying the resuspension factor (m^-1) from urban areas, where soil particle resuspension would dominate, can lead to a biased estimate of the surface-air ¹³⁷Cs concentration. Moreover, the influence of bioaerosol ¹³⁷Cs on the atmospheric ¹³⁷Cs concentration would last longer, because undecontaminated forests commonly exist within the DRZ.

How dynamic transfer models can complement an equilibrium-based approach: Case studies of radiocesium transfer to forest trees following accidental atmospheric release

Accidental release of radionuclides caused by nuclear accidents like those in Fukushima and Chernobyl can result in pulses of radioactivity entering the forest environment. Due to intense recycling in the forest, equilibrium between radioactivity concentrations in trees and in soil may not be reached during the period of severe short-term radionuclide transport following the accident. Another question arises as to whether the equilibrium hypothesis using empirical concentration ratios (CRs) can be applied to the long-term period. Using two atmospheric ¹³⁷Cs fallout scenarios in the Fukushima and Chernobyl sites, this study investigated whether the CR approach could provide conservative predictions of ¹³⁷Cs levels in trees following ¹³⁷Cs fallout events by comparing predictions from the CR approach using data gathered for trees by the IAEA to those from dynamic transfer models and actual measured data. The inter-comparisons also aimed to investigate whether the CR approach could account for the variability of ¹³⁷Cs levels across different tree organs. The results showed that caution may be necessary when using the CR approach, which relies on the IAEA dataset, to estimate ¹³⁷Cs accumulation in forest trees in the short - and long term following atmospheric ¹³⁷Cs fallout events. A calculation by TRIPS 2.0 demonstrated the importance of considering the distribution within tree organs for in-depth analysis of radiological impact of forest trees. Our findings suggest that it may be preferable to use CR values based on site-specific data rather than generic data collected from various sites. This is particularly relevant when studying the sites where the bioavailability of ¹³⁷Cs for trees and thus possible exposures are higher. This study also showed that dynamic modeling approaches could offer an alternative means of estimating CR values of the entire tree or specific tree organs in situations where empirically derived values are not available.

Influence of non-equilibrium and nonlinear sorption of 137Cs in soils. Study with stirred flow-through reactor experiments and quantification with a nonlinear equilibrium-kinetic model

This paper addresses the modelling of cesium sorption in non-equilibrium and nonlinear conditions with a two-site model. Compared to the classical K_d approach, the proposed model better reproduced the breakthrough curves observed during continuous-flow stirred tank reactor experiments conducted on two contrasted soils. Fitted parameters suggested contrasted conditions of cesium sorption between 1) equilibrium sites, with low affinity and high sorption capacity comparable to CEC and 2) non-equilibrium sites, with a fast sorption rate (half-time of 0.2-0.3 h), a slow desorption rate (half-time of 3-9 days) and a very low sorption capacity (0.02-0.04% of CEC). Comparison of EK sites densities with sorption capacities derived from the literature suggests that the EK equilibrium and kinetic sites might correspond to ion exchange and surface complexation of soil clay minerals respectively. This work stresses the limits of the K_d model to predict ¹³⁷Cs sorption in reactive transport conditions and supports an alternative non-equilibrium nonlinear approach.

Analysis of the vertical distribution of 137Cs and 239,240Pu in waterlogged and non-boggy soils by the sequential extraction method

Mobile capabilities of ¹³⁷Cs, ^239,240Pu and some stable element physicochemical forms were studied in soil layers at a depth of 6-8 cm (maximum concentrations of radionuclides) of non-boggy and waterlogged soil cores sampled on the shores of Lake Bedugnis (Lithuania). Soil samples were acidic with small amount of clay (muscovite). The radionuclide activity concentrations were about 2.5 times higher in non-boggy soil. It was explained by different conditions of deposition of radioactive fallout at these sampling sites and density of the soil samples. The value of the exchangeable fraction of radionuclides is shown to be an indicator of their migration capabilities in the soil. Under anaerobic conditions in waterlogged soil (stagnant pore water conditions) and due to the presence of clay admixture, ∼90% of ¹³⁷Cs was concentrated in the residual fraction and its mobility was low. ^239,240Pu was concentrated mainly in organic matter with rather large exchangeable fractions (∼9.6-∼13.9%). Under oxidizing conditions in non-boggy soil (dynamic pore water conditions in the case of rain), ^239,240Pu was mostly concentrated in the oxide fraction. Its exchangeable fractions were less than those in waterlogged soil and, respectively, its mobility was lower. In non-boggy soils, exchangeable fractions of ¹³⁷Cs were large and varied in the range of 10.1-12.2%, which indicated its high mobility. In the case, ¹³⁷Cs adsorption by clay materials was reduced and its residual fraction did not exceed 71.3%. The obtained data show that in the area of Lake Bedugnis, the migration capabilities of ¹³⁷Cs and ^239,240Pu change in antiphase.

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.

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.

Untangling radiocesium dynamics of forest-stream ecosystems: A review of Fukushima studies in the decade after the accident

Forest-stream ecosystems are widespread and biodiverse terrestrial landscapes with physical and social connections to downstream human activities. After radiocesium is introduced into these ecosystems, various material flows cause its accumulation or dispersal. We review studies conducted in the decade after the Fukushima nuclear accident to clarify the mechanisms of radiocesium transfer within ecosystems and to downstream areas through biological, hydrological, and geomorphological processes. After its introduction, radiocesium is heavily deposited in the organic soil layer, leading to persistent circulation due to biological activities in soils. Some radiocesium in soils, litter, and organisms is transported to stream ecosystems, forming contamination spots in depositional habitats. While reservoir dams function as effective traps, radiocesium leaching from sediments is a continual phenomenon causing re-contamination downstream. Integration of data regarding radiocesium dynamics and contamination sites, as proposed here, is essential for contamination management in societies depending on nuclear power to address the climate crisis.

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.

Leaching characteristics of 137Cs for forest floor affected by the Fukushima nuclear accident: A litterbag experiment

In forest ecosystems, forest litter is considered an active medium for radiocesium (¹³⁷Cs). To understand discharge mechanisms of highly bioavailable dissolved ¹³⁷Cs from forests to river systems, we investigated the characteristics of ¹³⁷Cs leaching from forest litter as observed from litterbag experiments. Leaching experiments with conifer needle and deciduous broadleaf litters were then conducted. After soaking conifer needles and broadleaf litters for 20 min, 140 min, and 1 day, the mean values of the ¹³⁷Cs leaching ratios were 0.13-2.0% and 0.81-6.6%, respectively, indicating that ¹³⁷Cs leaching ratios are different between forest litter types. To elucidate the factors affecting ¹³⁷Cs leaching from forest litter, a multi-regression analysis of ¹³⁷Cs leaching ratios was conducted against antecedent mean precipitation and temperature before sampling the litterbag and accumulated temperature during the litterbag experiments. The ¹³⁷Cs leaching ratios showed a negative correlation to the antecedent mean precipitation for both litters and the accumulated temperature for broadleaf litters, whereas it exhibited a positive correlation with the antecedent mean temperature for both litters and the accumulated temperature for conifer needle litters. It was proposed that the fraction of ¹³⁷Cs in labile sites in forest litter increased/decreased due to litter decomposition by antecedent/accumulated temperature, and that this fraction can be washed off by the antecedent precipitation. The different effects of accumulated temperature on ¹³⁷Cs leaching from conifer needles and broadleaf litters could be due to their different decomposition rates. Our results contribute further the understanding of the mechanisms associated with dissolved ¹³⁷Cs discharge from forested catchments.

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

Since Fukushima accident, dozens of field studies have been conducted in order to quantify and understand the behaviour of atmospheric radiocesium (¹³⁷Cs) fallouts in contaminated forests of Fukushima and neighbouring prefectures. In this paper, we carry out a detailed review of data acquired over 2011-2017 in Japanese cedar and cypress plantations, focusing on aerial tree organs, soil layers and tree-to-soil depuration fluxes. To enable comparison and reinforce the consistency between sites, radiological measurements were normalized by the deposit and interpolated onto the same spatio-temporal frame. Despite some (poorly explained) residual variability, we derived a "mean" pattern by log-averaging data among sites. These "mean" results were analysed with the help of a simple mass-balance approach and discussed in the light of post-Fukushima literature. We demonstrated that the activity levels and dynamics in all compartments were consistent and generally well reproduced by the mass balance approach, for values of the interception fraction between 0.7 and 0.85. The analysis indicated that about 5% of the initial deposit remained in the aerial vegetation after 6 years, more than two thirds of intercepted ¹³⁷Cs being transferred to the soil due to throughfall. The simulations indicated that foliar uptake might have contributed between 40% and 100% to the activity transferred to stem wood. The activity concentration in canopy organs rapidly decreased in the first few months then more slowly, according to an effective half-life of about 1.6 years. The activity level in the organic layer peaked in summer 2011 then decreased according to an effective half-life of 2.2 years. After a rapid increase in 2011, the contamination of mineral horizons continued to increase more slowly, 85% of ¹³⁷Cs incoming through the organic layer being retained in the 0-5 cm layer according to a mean residence time longer than in the upper layer (7 against 1.5 years).

Chemical sequential extraction of O horizon samples from Fukushima forests: Assessment for degradability and radiocesium retention capacity of organic matters

To investigate how radiocesium (¹³⁷Cs) is retained in the O horizon via interactions with organic matter, we collected O horizon samples in Japanese cedar (Cryptomeria japonica) and konara oak (Quercus serrata) forest sites in Fukushima during the 8 years following the Fukushima Dai-ichi Nuclear Power Plant accident. To assess degradability and ¹³⁷Cs retention capacity of organic matter, we conducted chemical sequential extraction with organic solvent and sulfuric acid, collecting the following fractions: organic solvent extractives (Fraction 1), acid-soluble carbohydrates (Fraction 3), and acid-insoluble residue (Fraction 4). In all samples, across sampling years and sites, ¹³⁷Cs content in Fractions 1, 3, and 4, as a proportion of the total ¹³⁷Cs content, was 0.0-23.6%, 18.4-42.9%, and 44.8-76.0%, respectively. Generally, ¹³⁷Cs is considered to be electrostatically bound to organic matter and relatively mobile, making it easily extractable by sulfuric acid treatment. However, we observed a relatively high proportion of ¹³⁷Cs in Fraction 4, suggesting strong retention of ¹³⁷Cs and their immobility in the O horizon. Complex organic matter such as lignin or tannin may contribute this retention. We also noted that some part of ¹³⁷Cs may be also retained by clay minerals in the O horizon. Although organic matter in Fractions 1 and 3 is considered to decompose faster than that in Fraction 4, over the observation period the ¹³⁷Cs proportion and net rate of decrease in ¹³⁷Cs content (in total and in each fraction) remained nearly constant. This result implies that decomposition of organic matter and the consequent release of bound ¹³⁷Cs may be partly compensated by additional input of ¹³⁷Cs from the canopy and ¹³⁷Cs recycling by soil microorganisms. Our study highlights the potential role of organic matter in the O horizon as a temporary reservoir of ¹³⁷Cs and a driver of the ¹³⁷Cs cycle in forest ecosystems.

Processes affecting land-surface dynamics of 129I impacted by atmospheric 129I releases from a spent nuclear fuel reprocessing plant

Terrestrial environments impacted by atmospheric releases of ¹²⁹I from nuclear plants become contaminated with ¹²⁹I; however, the relative importance of each land-surface ¹²⁹I-transfer pathway in the process of the contamination is not well understood. In this study, transfers of ¹²⁹I in an atmosphere-vegetation-soil system are modeled and incorporated into an existing land-surface model (SOLVEG-II). The model was also applied to the observed transfer of ¹²⁹I at a vegetated field impacted by atmospheric releases of ¹²⁹I (as gaseous I₂ and CH₃I) from the Rokkasho reprocessing plant, Japan, during 2007. Results from the model calculation and inter-comparison of the results with the measured environmental samples provide insights into the relative importance of each ¹²⁹I-transfer pathway in the processes of ¹²⁹I contamination of leaves and soil. The model calculation revealed that contamination of leaves of wild bamboo grasses was mostly caused by foliar adsorption of inorganic ¹²⁹I (81%) following wet deposition of ¹²⁹I. In contrast, accumulation of ¹²⁹I in the leaf due to foliar uptake of atmospheric ¹²⁹I₂ (2%) was lesser. Root uptake of soil ¹²⁹I was low, accounted for 17% of the ¹²⁹I of the leaf. The low root-uptake of ¹²⁹I in spite of the ¹²⁹I contained in the soil was ascribed to the fact that the most fraction (over 90%) of the soil ¹²⁹I existed in "soil-fixed" (not plant-available) form. Regarding the ¹²⁹I-transfer to the soil, wet deposition of ¹²⁹I was ten-fold more effective than dry deposition of atmospheric ¹²⁹I₂; however, the deposition of ¹²⁹I during the year represented only 2% of the model-assumed ¹²⁹I that pre-existed in the soil; indicating the importance of long-term accumulation of ¹²⁹I in terrestrial environments. The model calculation also revealed that root uptake of inorganic ¹²⁹I can be more influential than volatilization by methylation in exportation of ¹²⁹I from soil.

Stabilization/Solidification of Strontium Using Magnesium Silicate Hydrate Cement

Magnesium silicate hydrate (M–S–H) cement, formed by reacting MgO, SiO2, and H2O, was used to encapsulate strontium (Sr) radionuclide. Samples were prepared using light-burned magnesium oxide and silica fume, with sodium hexametaphosphate added to the mix water as a dispersant. The performance of the materials formed was evaluated by leach testing and the microstructure of the samples was also characterized. The stabilizing/solidifying effect on Sr radionuclide in the MgO–SiO2–H2O system with low alkalinity is demonstrated in the study. The leaching rate in a standard 42-day test was 2.53 × 10−4 cm/d, and the cumulative 42-day leaching fraction was 0.06 cm. This meets the relevant national standard performance for leaching requirements. Sr2+ was effectively incorporated into the M–S–H hydration products and new phase formation resulted in low Sr leaching being observed.

New predictions of 137Cs dynamics in forests after the Fukushima nuclear accident

Most of the area contaminated by the Fukushima Daiichi Nuclear Power Plant accident is covered by forest. In this paper, we updated model predictions of temporal changes in the ¹³⁷Cs dynamics using the latest observation data and newly provided maps of the predicted ¹³⁷Cs activity concentration for wood, which is the most commercially important part of the tree body. Overall, the previous prediction and latest observation data were in very good agreement. However, further validation revealed that the migration from the soil surface organic layer to the mineral soil was overestimated for evergreen needleleaf forests. The new prediction of the ¹³⁷Cs inventory showed that although the ¹³⁷Cs distribution within forests differed among forest types in the first 5 years, the difference diminished in the later phase. Besides, the prediction of the wood ¹³⁷Cs activity concentrations reproduced the different trends of the ¹³⁷Cs activity concentrations for cedar, oak, and pine trees. Our simulation suggests that the changes of the wood ¹³⁷Cs activity concentration over time will slow down after 5–10 years. Although the model uncertainty should be considered and monitoring and model updating must continue, the study provides helpful information on the ¹³⁷Cs dynamics within forest ecosystems and the changes in wood contamination.

Temporal dynamics of 137Cs distribution in soil and soil-to-crop transfer factor under different tillage systems after the Fukushima Daiichi Nuclear Power Plant accident in Japan

The accident at the Fukushima Daiichi Nuclear Power Plant (FDNPP) in Japan in 2011 released a large amount of radionuclides, primarily radiocesium-137 (¹³⁷Cs; half-life: 30 years), resulting in long-term contamination of soil and consequently crops. Tillage is a common agricultural management practice that alters the vertical distribution of nutrients in the soil. However, the effect of tillage on ¹³⁷Cs contamination in soil and crops over time remains unclear. In this study, we investigated the temporal changes in the vertical distribution of ¹³⁷Cs in the soil, concentration of ¹³⁷Cs in soybean and cover crops, and the transfer factor (TF) of ¹³⁷Cs from the soil to crops under three tillage systems (rotary cultivation [RC], moldboard plow [MP], and no tillage [NT]; main factors) using three cover crops (hairy vetch, winter rye, and fallow weeds; side factors). The amount of ¹³⁷Cs in the soil decreased exponentially with soil depth under the NT and RC treatments. By contrast, ¹³⁷Cs showed uniform distribution at each soil depth tested under the MP treatment since 2012. The exchangeable ¹³⁷Cs demonstrated a similar tendency as ¹³⁷Cs. The ¹³⁷Cs concentration in soybean (including grain and residue) and cover crops decreased exponentially with time. Consistently higher ¹³⁷Cs concentration was observed in soybean grains under the NT treatment, suggesting that tillage continuously reduced the concentration of ¹³⁷Cs in soybean over 7 years since the FDNPP accident. The TF of ¹³⁷Cs from soil to soybean and cover crops decreased continuously over time; however, ¹³⁷Cs concentration of soybean grain showed a positive linear correlation with its annual variation rate. Additionally, TF showed a positive logarithmic correlation with ¹³⁷Cs relaxation depth in the soil. These results enhance our understanding of the long-term behavior and radioecology of ¹³⁷Cs in agroecosystems in Japan since the radionuclide accident.

Six-year monitoring of the vertical distribution of radiocesium in three forest soils after the Fukushima Dai-ichi Nuclear Power Plant accident

After the Fukushima Dai-ichi Nuclear Power Plant accident on March 2011, several studies showed that the downward migration of ¹³⁷Cs from litter to mineral soil is more rapid in forests in Fukushima than in forests affected by the Chernobyl accident. Therefore, the downward migration within mineral soil layers is more important for predicting long-term dynamics of ¹³⁷Cs in forest ecosystems in Fukushima. In the present study, we monitored the detailed vertical distribution of ¹³⁷Cs in litter and soil layers for 6 y (2011-2017) following the previous study (2011-2012), and found that temporal changes in those distributions were different among mixed forest (MF), mature cedar (MC) and young cedar (YC) forests. The ¹³⁷Cs concentrations and inventories in the litter layer exponentially decreased with time for all sites, with more than 80-95% of the deposited ¹³⁷Cs on the forest floor distributed in mineral soil layers by 2017. The percentage of ¹³⁷Cs inventory in the litter layer to the total ¹³⁷Cs inventory in litter and mineral soil layers was well fitted by a single exponential equation with decreasing rate of 0.22-0.44 y^-1. The slower migration was observed in the YC site, probably because of higher initial interception of ¹³⁷Cs fallout by dense canopy. As the downward migration from litter to mineral soil progressed, the ¹³⁷Cs concentration in the first few cm of mineral soil surface gradually increased and became higher than the ¹³⁷Cs concentration in the litter within 2-3 y of the accident. The ¹³⁷Cs concentration in mineral soil layers exponentially decreased with depth throughout survey period, and an exponential equation fitted well. The relaxation depth of ¹³⁷Cs concentration in mineral soil layers estimated by the exponential equation were constantly increasing in the MC and YC sites with 0.08 cm y^-1. In contrast, there was no temporal increase in the relaxation depth in the MF site, indicating little migration to subsurface soil layer from not only litter layer but also surface soil layer. Further studies are necessary to identify the forests prone to the downward migration of ¹³⁷Cs and its factors regarding both forest and soil characteristics.

Characterizing vertical migration of 137Cs in organic layer and mineral soil in Japanese forests: Four-year observation and model analysis

Because of the Fukushima Dai-ichi Nuclear Power Plant accident, forest ecosystems in wide areas were contaminated with ¹³⁷Cs. It is important to characterize the behavior of ¹³⁷Cs after its deposition onto forest surface environments for evaluating and preventing long-term radiation risks. In the present study, ¹³⁷Cs vertical distributions in the soil profile were observed repeatedly at five forest sites with different vegetation types for 4.4 years after the accident in 2011, and ¹³⁷Cs migration in the organic layer and mineral soil was analyzed based on a comparison of models and observations. Cesium-137 migration from the organic layer to the underlying mineral soil was represented by a two-component exponential model. Cesium-137 migration from the organic layer was faster than that observed in European forests, suggesting that the mobility and bioavailability of ¹³⁷Cs could be suppressed rapidly in Japanese forests. At all sites, ¹³⁷Cs transfer in mineral soil could be reproduced by a simple diffusion equation model with continuous ¹³⁷Cs supply from the organic layer. The diffusion coefficients of ¹³⁷Cs in the mineral soil were estimated to be 0.042-0.55 cm² y^-1, which were roughly comparable with those of European forest soils affected by the Chernobyl Nuclear Power Plant accident. Model predictions using the determined model parameters indicated that 10 years after the accident, more than 70% of the deposited ¹³⁷Cs will migrate to the mineral soil but only less than 10% of the total ¹³⁷Cs inventory will penetrate deeper than 10 cm in the mineral soil across all sites. The results of the present study suggest that the ¹³⁷Cs deposited onto Japanese forest ecosystems will be retained in the surface layers of mineral soil for a long time.

Importance of root uptake of 14CO2 on 14C transfer to plants impacted by below-ground 14CH4 release

¹⁴C-labelled methane (¹⁴CH₄) released from deep underground radioactive waste disposal facilities can be a below-ground source of ¹⁴CO₂ owing to microbial oxidation of ¹⁴CH₄ to ¹⁴CO₂ in soils. Environmental ¹⁴C models assume that the transfer of ¹⁴CO₂ from soil to plant occurs via foliar uptake of ¹⁴CO₂. Nevertheless, the importance of ¹⁴CO₂ root uptake is not well understood. In the present study, below-ground transport and oxidation of ¹⁴CH₄ were modeled and incorporated into an existing land-surface ¹⁴CO₂ model (SOLVEG-II) to assess the relative importance of root uptake and foliar uptake on ¹⁴CO₂ transfer from soil to plants. Performance of the model in calculating the below-ground dynamics of ¹⁴CH₄ was validated by simulating a field experiment of ¹³CH₄ (as a substitute for ¹⁴CH₄) injection into subsoil in a wheat field in the UK. The proposed model simulation was then applied to ¹⁴C transfer in a hypothetical ecosystem impacted by continuous ¹⁴CH₄ input from the water table (bottom of 1-m thick soil), which simulated continuous release of ¹⁴CH₄ from a deep underground radioactive waste disposal facility. The contrast between the results obtained from the model calculation that assumed different distributions of roots (rooting depths of 11 cm, or 97 cm) and methane oxidation (characterized by e-folding depths of 5 cm, 20 cm, or 80 cm) in the soil provided insight into the relative importance of root uptake and foliar uptake pathways. In the shallowly rooted ecosystem with rooting depth of 11 cm, foliar uptake of ¹⁴CO₂ was significant, accounting for 80% of the ¹⁴C accumulation (as organic ¹⁴C) in the plant (leaf compartment). By contrast, in a deeply rooted ecosystem (rooting depth of 97 cm), where the root penetrated to depths close to the water-table, more than half (63%) the ¹⁴C accumulated in the plant was transferred via the root uptake pathway. We found that ¹⁴CO₂ root uptake (thus ¹⁴C accumulation in the plant) in this ecosystem depended on the distribution of methane oxidation in the soil; all ¹⁴C accumulated in the plant was transferred by the root uptake pathway when methane oxidation occurred at considerable depths (e-folding depths of 20 cm, or 80 cm) in the soil. The high level of ¹⁴CO₂ root uptake was ascribed to the oxidation of added ¹⁴CH₄ (i.e., production of ¹⁴CO₂) in the deep part of the soil and the subsequent high level of root uptake of the deep soil-water containing ¹⁴CO₂. These results indicate that ¹⁴CO₂ root uptake contributes significantly to ¹⁴CO₂ transfer to plants if ¹⁴CH₄ oxidation occurs at great depths and roots penetrate deeply into the soil. It is recommended that current environmental ¹⁴C models must be refined to consider the importance of the root uptake pathway to ensure that dose estimates of ¹⁴CH₄ release from deep underground waste disposal facilities are accurate.

Six-year monitoring of the vertical distribution of radiocesium in three forest soils after the Fukushima Dai-ichi Nuclear Power Plant accident

After the Fukushima Dai-ichi Nuclear Power Plant accident on March 2011, several studies showed that the downward migration of ¹³⁷Cs from litter to mineral soil is more rapid in forests in Fukushima than in forests affected by the Chernobyl accident. Therefore, the downward migration within mineral soil layers is more important for predicting long-term dynamics of ¹³⁷Cs in forest ecosystems in Fukushima. In the present study, we monitored the detailed vertical distribution of ¹³⁷Cs in litter and soil layers for 6 y (2011-2017) following the previous study (2011-2012), and found that temporal changes in those distributions were different among mixed forest (MF), mature cedar (MC) and young cedar (YC) forests. The ¹³⁷Cs concentrations and inventories in the litter layer exponentially decreased with time for all sites, with more than 80-95% of the deposited ¹³⁷Cs on the forest floor distributed in mineral soil layers by 2017. The percentage of ¹³⁷Cs inventory in the litter layer to the total ¹³⁷Cs inventory in litter and mineral soil layers was well fitted by a single exponential equation with decreasing rate of 0.22-0.44 y^-1. The slower migration was observed in the YC site, probably because of higher initial interception of ¹³⁷Cs fallout by dense canopy. As the downward migration from litter to mineral soil progressed, the ¹³⁷Cs concentration in the first few cm of mineral soil surface gradually increased and became higher than the ¹³⁷Cs concentration in the litter within 2-3 y of the accident. The ¹³⁷Cs concentration in mineral soil layers exponentially decreased with depth throughout survey period, and an exponential equation fitted well. The relaxation depth of ¹³⁷Cs concentration in mineral soil layers estimated by the exponential equation were constantly increasing in the MC and YC sites with 0.08 cm y^-1. In contrast, there was no temporal increase in the relaxation depth in the MF site, indicating little migration to subsurface soil layer from not only litter layer but also surface soil layer. Further studies are necessary to identify the forests prone to the downward migration of ¹³⁷Cs and its factors regarding both forest and soil characteristics.

Development and assessment of a simple ecological model (TRIPS) for forests contaminated by radiocesium fallout

The management of vast forested zones contaminated by radiocesium (rCs) following the Chernobyl and Fukushima fallout is of great social and economic concern in affected areas and requires appropriate dynamic models as predictive or questioning tools. Generally, the existing radio-ecological models need less fragmented data and more ecological realism in their quantitative description of the rCs cycling processes. The model TRIPS ("Transfer of Radionuclide In Perennial vegetation Systems") developed in this study privileged an integrated approach which makes the best use of mass balance studies and available explicit experimental data for Scots pine stands. A main challenge was the differentiation and calibration of foliar absorption as well as root uptake in order to well represent the rCs biocycling. The general dynamics of rCs partitioning was simulated with a relatively good precision against an independent series of observed values. In our scenario the rCs biological cycling enters a steady-state about 15 years after the atmospheric deposits. At that time, the simulations showed an equivalent contribution of foliage and root uptake to the tree contamination. But the root uptake seems not sufficient to compensate the activity decline in the tree. The initial foliar uptake and subsequent internal transfers were confirmed to have a great possible impact on the phasing of tree contamination. An extra finding concerns the roots system acting as a buffer in the early period. The TRIPS model is particularly useful in cases where site-specific integrated datasets are available, but it could also be used with adequate caution to generic sites. This development paves the way for simplification or integration of new modules, as well as for a larger number of other applications for the Chernobyl or Fukushima forests once the appropriate data become available. According to the sensitivity analysis that involves in particular reliable estimates of net foliar uptake as well as root uptake not disconnected from rCs exchange reactions in soil.

Calibration of forest 137Cs cycling model "FoRothCs" via approximate Bayesian computation based on 6-year observations from plantation forests in Fukushima

Predicting the environmental fate of ¹³⁷Cs in forest ecosystems along with the concentrations of ¹³⁷Cs in tree parts are important for the managements of radioactively contaminated forests. In this study, we calibrate the Forest RothC and Cs model (FoRothCs), a forest ecosystem ¹³⁷Cs dynamics model, using observational data obtained over six years from four forest sites with different levels of ¹³⁷Cs contamination from Fukushima Prefecture. To this end, we applied an approximate Bayesian computation (ABC) technique based on the observed ¹³⁷Cs concentrations (Bq kg^-1) of five compartments (leaf, branch, stem, litter, and soil) in a Japanese cedar plantation. The environmental decay (increment) constants of the five compartments were used as the summary statistics (i.e., the metric for model performance) to infer the five parameters related to ¹³⁷Cs transfer processes in FoRothCs. The ABC technique successfully reconciled the model outputs with the observed trends in ¹³⁷Cs concentrations at all sites during the study period. Furthermore, the estimated parameters are in agreement with the literature values (e.g., the root uptake rates of ¹³⁷Cs). Our study demonstrates that model calibration with ABC based on the trends in ¹³⁷Cs concentrations of multi compartments is useful for reducing the prediction uncertainty of ¹³⁷Cs dynamics in forest ecosystems.

Radiocesium distribution in aggregate-size fractions of cropland and forest soils affected by the Fukushima nuclear accident

The Fukushima Daiichi nuclear power plant accident caused serious radiocesium (¹³⁷Cs) contamination in soils in a range of terrestrial ecosystems. It is well documented that the interaction of ¹³⁷Cs with soil constituents, particularly clay minerals, in surface soil layers exerts strong control on the behavior of this radionuclide in the environment; however, there is little understanding of how soil aggregation-the binding of soil particles together into aggregates-can affect the mobility and bioavailability of ¹³⁷Cs in soils. To explore this, soil samples were collected at seven sites under different land-use conditions in Fukushima and were separated into four aggregate-size fractions: clay-sized (<2 μm); silt-sized (2-20 μm); sand-sized (20-212 μm); and macroaggregates (212-2000 μm). The fractions were then analyzed for ¹³⁷Cs content and extractability and mineral composition. In forest soils, aggregate formation was significant, and 69%-83% of ¹³⁷Cs was associated with macroaggregates and sand-sized aggregates. In contrast, there was less aggregation in agricultural field soils, and approximately 80% of ¹³⁷Cs was in the clay- and silt-sized fractions. Across all sites, the ¹³⁷Cs extractability was higher in the sand-sized aggregate fractions than in the clay-sized fractions. Mineralogical analysis showed that, in most soils, clay minerals (vermiculite and kaolinite) were present even in the larger-sized aggregate fractions. These results demonstrate that larger-sized aggregates are a significant reservoir of potentially mobile and bioavailable ¹³⁷Cs in organic-rich (forest and orchard) soils. Our study suggests that soil aggregation reduces the mobility of particle-associated ¹³⁷Cs through erosion and resuspension and also enhances the bioavailability of ¹³⁷Cs in soils.

The impact of radiocesium input forms on its extractability in Fukushima forest soils

The effects of ¹³⁷Cs deposit forms on its ageing in soil have not yet been reported. Soluble and Solid ¹³⁷Cs input forms were mixed with the mineral soils collected under Fukushima's coniferous and broadleaf forests, incubated under controlled laboratory, and examined the evolution of ¹³⁷Cs availability over time. Results show that the extracted ¹³⁷Cs fraction with water was less than 1% for the soluble input form and below detection limit for the solid input forms. Likewise, with an acetate reagent, the extracted ¹³⁷Cs fraction ranged from 46 to 56% for the soluble input and from 2 to 15% for the solid input, implying that the nature of the ¹³⁷Cs contamination strongly influences its extractability and mobility in soil. Although the degradation of organic materials was apparent, its impact on the ¹³⁷Cs extractability was found to be weak. Nevertheless, more Ac-available ¹³⁷Cs was obtained from broadleaf organic material mixes than the coniferous counterparts, suggesting that the lignified nature of latter tend to retain more ¹³⁷Cs. When extrapolated to a field context, more available ¹³⁷Cs fraction may be expected from wet-derived contaminated forest soils than contaminated via solid-derived inputs. Such information could be helpful for radioecological management schemes in contaminated forest environments.

Role of soil-to-leaf tritium transfer in controlling leaf tritium dynamics: Comparison of experimental garden and tritium-transfer model results

Environmental transfer models assume that organically-bound tritium (OBT) is formed directly from tissue-free water tritium (TFWT) in environmental compartments. Nevertheless, studies in the literature have shown that measured OBT/HTO ratios in environmental samples are variable and generally higher than expected. The importance of soil-to-leaf HTO transfer pathway in controlling the leaf tritium dynamics is not well understood. A model inter-comparison of two tritium transfer models (CTEM-CLASS-TT and SOLVEG-II) was carried out with measured environmental samples from an experimental garden plot set up next to a tritium-processing facility. The garden plot received one of three different irrigation treatments - no external irrigation, irrigation with low tritium water and irrigation with high tritium water. The contrast between the results obtained with the different irrigation treatments provided insights into the impact of soil-to-leaf HTO transfer on the leaf tritium dynamics. Concentrations of TFWT and OBT in the garden plots that were not irrigated or irrigated with low tritium water were variable, responding to the arrival of the HTO-plume from the tritium-processing facility. In contrast, for the plants irrigated with high tritium water, the TFWT concentration remained elevated during the entire experimental period due to a continuous source of high HTO in the soil. Calculated concentrations of OBT in the leaves showed an initial increase followed by quasi-equilibration with the TFWT concentration. In this quasi-equilibrium state, concentrations of OBT remained elevated and unchanged despite the arrivals of the plume. These results from the model inter-comparison demonstrate that soil-to-leaf HTO transfer significantly affects tritium dynamics in leaves and thereby OBT/HTO ratio in the leaf regardless of the atmospheric HTO concentration, only if there is elevated HTO concentrations in the soil. The results of this work indicate that assessment models should be refined to consider the importance of soil-to-leaf HTO transfer to ensure that dose estimates are accurate and conservative.

Approaches to modelling radioactive contaminations in forests - Overview and guidance

Modelling the radionuclide cycle in forests is important in case of contamination due to acute or chronic releases to the atmosphere and from underground waste repositories. This article describes the most important aspects to consider in forest model development. It intends to give an overview of the modelling approaches available and to provide guidance on how to address the quantification of radionuclide transport in forests. Furthermore, the most important gaps in modelling the radionuclide cycle in forests are discussed and suggestions are presented to address the variability of forest sites.

Forest type effects on the retention of radiocesium in organic layers of forest ecosystems affected by the Fukushima nuclear accident

The Fukushima Daiichi nuclear power plant disaster caused serious radiocesium (137Cs) contamination of forest ecosystems over a wide area. Forest-floor organic layers play a key role in controlling the overall bioavailability of 137Cs in forest ecosystems; however, there is still an insufficient understanding of how forest types influence the retention capability of 137Cs in organic layers in Japanese forest ecosystems. Here we conducted plot-scale investigations on the retention of 137Cs in organic layers at two contrasting forest sites in Fukushima. In a deciduous broad-leaved forest, approximately 80% of the deposited 137Cs migrated to mineral soil located below the organic layers within two years after the accident, with an ecological half-life of approximately one year. Conversely, in an evergreen coniferous forest, more than half of the deposited 137Cs remained in the organic layers, with an ecological half-life of 2.1 years. The observed retention behavior can be well explained by the tree phenology and accumulation of 137Cs associated with litter materials with different degrees of degradation in the organic layers. Spatial and temporal patterns of gamma-ray dose rates depended on the retention capability. Our results demonstrate that enhanced radiation risks last longer in evergreen coniferous forests than in deciduous broad-leaved forests.

Post-deposition early-phase migration and retention behavior of radiocesium in a litter-mineral soil system in a Japanese deciduous forest affected by the Fukushima nuclear accident

The fate of radiocesium (¹³⁷Cs) derived from the Fukushima nuclear accident and associated radiation risks are largely dependent on its migration and retention behavior in the litter-soil system of Japanese forest ecosystems. However, this behavior has not been well quantified. We established field lysimeters in a Japanese deciduous broad-leaved forest soon after the Fukushima nuclear accident to continuously monitor the downward transfer of ¹³⁷Cs at three depths: the litter-mineral soil boundary and depths of 5 cm and 10 cm in the mineral soil. Observations were conducted at two sites within the forest from May 2011 to May 2015. Results revealed similar temporal and depth-wise variations in ¹³⁷Cs downward fluxes for both sites. The ¹³⁷Cs downward fluxes generally decreased year by year at all depths, indicating that ¹³⁷Cs was rapidly leached from the forest-floor litter layer and was then immobilized in the upper (0-5 cm) mineral soil layer through its interaction with clay minerals. The ¹³⁷Cs fluxes also showed seasonal variation, which was in accordance with variations in the throughfall and soil temperature at the sites. There was no detectable ¹³⁷Cs flux at a depth of 10 cm in the mineral soil in the third and fourth years after the accident. The decreased inventory of mobile (or bioavailable) ¹³⁷Cs observed during early stages after deposition indicates that the litter-soil system in the Japanese deciduous forest provides only a temporary source for ¹³⁷Cs recycling in plants.