Difference of ecological half-life and transfer coefficient in aquatic invertebrates between high and low radiocesium contaminated streams

The Fukushima accident emitted radioactive substances into the environment, contaminating litter, algae, sand substrate, aquatic invertebrates, and fish in freshwater streams. Because these substances have substantial effects on stream ecology over many years, it is necessary to clarify the diffusion and decay mechanisms of radiocesium. The transfer coefficient differed among aquatic invertebrate groups, likely due to the differences in habitat. The ecological half-life of cesium was longer where the air dose rate was lower. The transfer coefficient was also higher in areas with lower air dose rate. The radiocesium concentration in algae was inversely related to stream current velocity in the radiocesium-contaminated area. However, this relationship was not observed in the lower air dose rate area: the radiocesium concentration in algae in the rapid-velocity areas tended to be higher than that in the slow-velocity areas. This reverse trend would lead to a longer period of freshwater contamination. The radiocesium concentration would continue to decrease in highly contaminated areas, but it would be difficult to reduce the radiocesium concentration in less-contaminated areas because different contamination mechanisms are at work. Controlling the water flow is key to regulating radiocesium concentration in freshwater ecosystems.


Scientific Reports
| (2020) 10:21819 | https://doi.org/10.1038/s41598-020-78844-8 www.nature.com/scientificreports/ large geographical areas 12,13 . However, there are few studies focusing on radioactive contamination in freshwater biota such as algae and aquatic invertebrates. Differences in the stream substrate structure play a large role in the degree of habitat contamination. Aquatic insects in stagnant pools generally exhibit higher radiocesium levels than those in flowing water, although the extent of cesium contamination in sand substrates do not significantly differ between pool and riffle sites 10 . Most of the radiocesium in ponds is in the sediment, and aquatic insects in the pond acquire radionuclides from the sediment through ingesting biotic components 11 . Aquatic insects inhabiting pools thus ingest high levels of radiocesium from the material settling at the bottom. Therefore, aquatic insects that can swim quickly and move between riffles and pools frequently, such as members of Ameletidae 12 , exhibit no differences in their radiocesium contents whether they are found in riffles or pools 10 . Considering the higher levels of contamination in aquatic insects in pools, the radiocesium concentration in algae attached to stones in the pool, which constitutes the diet of some aquatic insects, should therefore also be higher than that of algae in the riffles.
Chernobyl accident spread many radionuclides over a 2000 km from the pollution source, such as Finland, Sweden and Norway. Extensive cesium contamination in pools due to the Chernobyl fallout has also been detected in Norway 16 . Fish have been monitored for radioactivity in Øvre Heimdalsvatn, a Norwegian subalpine lake over the course of two decades 17,18 . Concentrations of cesium-137 in brown trout in the lake reached 8400 Bq/kg in 1987 and declined to 200-300 Bq/kg in 2008. Recently, the contamination level has remained low and is declining in an asymptotic manner 18 . However, the ecological consequences of this type of radioactive contamination in stream ecology are poorly understood.
Radiocesium contamination effects have differed depending on the habitat condition, and the habitat contamination continues for longer years. To precisely describe the mechanisms of diffusion and the export of cesium-134 and cesium-137 deposited in freshwater streams, the relationship between radiocesium concentration in sand, litter, algae and aquatic invertebrates and the days after diffusion as well as distance from FDNPP were considered in this study. The relationship between the concentration of cesium-134 and cesium-137 in periphytic algae and stream current velocity were also analysed. The effects of water flow on contamination levels in biota were then clarified, and the mechanism of decontamination particularly in areas with low levels of radiocesium contamination area was discussed.

Materials and methods
Study site. The study sites were located approximately 20-75 km from the Fukushima Daiichi Nuclear Power Plant in Fukushima Prefecture, Japan (Fig. 1). According to an aircraft radioactivity survey reported by the Ministry of Education, Culture, Sports, Sciences, and Technology of Japan 19 , the air dose rate in this region was 0.3-3.2 μSv/h, and the deposition of cesium-134 and cesium-137 ranged from less than 64,000 to 940,000 Bq/m 2 (Table 1) in June 2011. The study catchment area is mostly forested and dominated with deciduous trees. Other areas in the region are also forested as well, with Japanese cedar and cypress plantations used for timber production. A field survey was conducted at one headwater tributary (A) of the Nagase River and three headwater tributaries (B, C, and D) of the Kido River. The substrate of these sites was consisted with sand, cobble and rocks. Geological feature of the soil on all the sites was the same, biotite granite. Streams at sites B, C and D were covered with riparian forests and it was difficult for sunlight to penetrate directly. Stream width of site A was wider than sites B, C and D, so sunlight could penetrate through the forest cover and contact the stream surface only along the middle of the stream.   13 . The sand substrate was sampled in each riffle to a depth of 5-cm. When sand was not immediately visible in the stream substrate, stones were removed and the sand underneath the stones was sampled. Litter shed in the water was collected after gentle hand-rinsing. Leaf litter forms the base of stream food webs. Periphytic algae were collected by brushing the pebbles or rocks with a toothbrush. These algae are also primary producers at the base of stream food webs. Prior to brushing, we gently hand-rinsed the stone surface to remove other organic matter and aquatic invertebrates in the periphyton.
Aquatic invertebrates from thirteen groups (Perlidae Gen. spp., Nemouridae Gen. spp., Ephemera japonica, Ephemerellidae Gen. spp., Heptageniidae Gen. spp., Hydropsychidae Gen. spp., Stenopsychi spp., Rhyacophilidae Gen. spp., Epiophlebia superstes, Lanthus fujiacus, Tipulidae Gen. spp., and Corydalidae Gen. spp., Geothelphusa dehaani,) were qualitatively sampled from riffles at a depth of 10-15 cm at the four sites from July 2013 to April 2019. At each site, a D-frame net with a 1-mm mesh was placed downstream of the sampling area on the substrate in water. We then disturbed the substrate upstream of the net, allowing insects to drift into the D-frame net. The sampled aquatic invertebrates were identified to family level in the field and then frozen.
Three bricks (210 × 100 × 60 mm) were placed separately within the stream riffle at a depth of 10-20 cm on August 25, 2014 at each of the four sites. Then, periphytic algae growing on the bricks were collected by brushing the substrate with a toothbrush. Before brushing, we gently hand-rinsed the brick surface in running water to remove other organic matter from the periphytic algae. The sampling was carried out eight times: in October and December 2014; March, May, June, July and November 2015; and April 2016. Stream velocity of right side, upper reaches side and left side of each brick were measured and averaged. This averaged value was used as the stream velocity of each periphytic algae sample. Radiocesium analysis. Radiocesium was analysed according to the methods in previous studies 10,20 . Samples of sand substrate and litter were dried at 75 °C in an oven. Thereafter, samples of sand were placed in a sieve (mesh size 2 mm; Iida, Japan), and the sand that passed through the sieve was used, meaning that the sand substrate in this study included silt granules. Samples of algae were concentrated via evaporation and dried in an oven at 75 °C. Samples of aquatic invertebrates were also dried in an oven at 75 °C. All samples were homogenized and packed into 100-ml polystyrene containers (U-8). Gamma-ray spectrometric measurements were performed on each sample. The radioactive concentrations of cesium-134 (604 keV) and cesium-137 (662 keV) were measured using an HPGe coaxial detector system (GEM40P4-76, Seiko EG and G, Tokyo, Japan) at the Forestry and Forest Products Research Institute (FFPRI) with a time of 36,000 s or longer. Data with a standard error of < 10% were used for analysis. Gamma-ray peaks at 604.66 and 661.64 keV were used to identify cesium-134 and cesium-137, respectively. The measurement system was calibrated using a standard gamma-ray source (MX033U8PP; Japan Radioisotope Association, Tokyo, Japan), and a standard soil material (IAEA-444) was used to check the measurement accuracy (relative expanded uncertainty value: 4.6%).
The homogenized samples of aquatic invertebrates were further homogenized with Na 2 SO 4 for radiocesium analysis. The processed aquatic invertebrates were packed into 20-ml plastic tubes. Gamma-ray spectrometric measurements were performed on each tube. Total concentrations of cesium-134 and cesium-137 were measured using a gamma particle counter with a NaI crystal detector (2480WIZARD2, Perkin Elmer, Downers Grove, IL, USA) and a sampling time of 21,600 s at the FFPRI. The measurement system was calibrated with a standard gamma-ray source (E-265/107F16-3, Canberra Industries, Oak Ridge, TN, USA) and a standard soil material (IAEA-444) was used to check the measurement accuracy (relative expanded uncertainty value: 2.9%). A strong correlation between the data obtained with the HPGe coaxial detector system and the data obtained with the NaI crystal detector was observed (cesium-134: τ = 0.75, z = 5.61, n = 33, P < 0.001; cesium-137: τ = 0.76, z = 8.43, Table 1. Air dose rate and the deposition of Cs according to an aircraft radioactivity survey by MEXT (2011), averaged value of dose rate 1-m above the ground on the sampling date from 2013 to 2019 (n = 23) and five environmental factors on four sites on the sampling date from 2013 to 2019 (n = 23). ****P < 0.001. www.nature.com/scientificreports/ n = 63, P < 0.001; Kendall test). As the detection limit of the HPGe coaxial detector system was higher than that of the gamma particle counter with the NaI crystal detector, the data obtained with the gamma particle counter with the NaI crystal detector were used for the analysis of aquatic invertebrates because we could then use smaller numbers of aquatic invertebrates. Only samples with a dry mass > 50 mg were included in the analysis to avoid measurement errors due to the low quantity.
The relation between radiocesium concentration and the days after accident. The data measured for sand substrate, litter, algae and aquatic invertebrates were plotted on a logarithmic scale against the days after accident occurred, and correlation coefficients were calculated. The significance of the correlations was statistically assessed using the Pearson's product-moment correlation coefficients. Radiocesium concentrations that were below the detection level on a dry weight basis because of insufficient mass for radiocesium analysis were excluded from the analysis. Using the regression curves, the radiocesium decay was corrected to September 1, 2013 and the differences in the values among sites were statistically tested using the Friedman-test.
The ecological half-life and the transfer coefficient. Ecological half-life was estimated by fitting the measured data to an exponential curve with the equation: where the activity concentration in a sample group is plotted as a function of time. A is a constant, T ecol is the ecological half-life, and t is time after the accident. The ecological half-life of cesium-134 and cesium-137 was estimated through exponential curve fitting. A positive relationship between the radiocesium concentration and the days after the accident indicates the increase of radiocesium according to the time, thus the ecological half-life cannot be calculated under this scenario. It would be better to include the physical decay in calculations of the biological half-life using long-term data 21 . However the research of this study started 2 years after the accident and the research period was only 6 years. Moreover, we did not need to establish the biological half-life, but rather the ecological half-life. For this reason, in this study, we did not consider the physical decay of a radionuclide to calculate the ecological half-life, which enabled us to integrate the actual radionuclide transition in the field. Differences in ecological half-life among sites were statistically tested using the Friedman-test. The transfer coefficient (T) is an empirical coefficient that can be used in predictive models for natural ecosystems. It is defined by the following expression: The diet of some aquatic invertebrate species comprises litter rather than algae, whereas the diet of other species comprises both algae and litter. In this analysis, all the invertebrate groups were analysed under the assumption that their diets were dependent on algae or litter. Differences in the T value among sites were statistically tested using the Friedman-test.
The relation between radiocesium concentration and current velocity. The radiocesium concentrations based on dry weight were estimated for September 1, 2013. Radiocesium concentrations that were below the detection level on a dry weight basis, because of insufficient mass for radiocesium analysis, were excluded from statistical analysis. Nonparametric tests for statistical analysis were conducted. The Kendall rank correlation test was used to clarify the relationships between stream velocity and concentrations of radioactive cesium-134 and cesium-137. To examine the differences in the decline between sites A and D, the comparison of two regression slopes test was conducted.

Results
The dose rates in the air in June 2011 at sites C and D were higher than those at sites A and B. The deposition of radiocesium in June 2011 followed the same pattern. There were no differences in values for air or water temperature, pH, or DO among the four sites, whereas the EC and dose rate on the sampling dates differed significantly among sites ( Table 1).
Radiocesium concentrations in sand substrates at sites B, C and D significantly decreased with the years, but not at A. The radiocesium concentrations in litter significantly decreased with the years at the four sites. The radioccesium concentrations in algae decreased with the years and the relation was significant at site B. Radiocesium concentration in aquatic invertebrates at site D decreased with the years, and significant decreases were observed in Ephemera japonica, Lanthus fujiacus, and Tipulidae Gen. spp. for cesium-134 and Ephemera japonica, Rhyacophilidae Gen. spp., Lanthus fujiacus, and Tipulidae Gen. spp. for cesium-137. The radiocesium concentrations at site C also decreased with the years and the decreases were significant in Ephemera japonica, Heptageniidae Gen. spp., Stenopsyche spp., Epiophlebia superstes, and Lanthus fujiacus for cesium-134 and Nemouridae Gen. spp., Heptageniidae Gen. spp., Epiophlebia superstes, and Lanthus fujiacus for cesium-137, except for those in Rhyacophilidae Gen. spp. and Geothelphusa dehaani that increased. The radiocesium concentrations in aquatic invertebrates at site B decreased with the years and the decreases were significant in Nemouridae Gen. spp. and Ephemerellidae Gen. spp. for cesium-134 and Hydropsychidae Gen. spp. for cesium-137, except for those in  (Table 3). Because the relationship between the radiocesium concentration and the year was positive in some groups (Table 2), ecological half-life could not be calculated for those groups. The positive relationship between the radioesium concentration and the year indicates a longer ecological half-life, and positive relationship was often observed at site A. The ecological half-life was longer where the air dose rate was lower (cesium-134: χ 2 = 19.3, P < 0.001; cesium-137: χ 2 = 20.5, P < 0.001; Friedman-test).

Discussion
The algae attached to brick surfaces in the low-velocity riffle had higher radiocesium concentrations than those in the high-velocity riffle at sites C and D. Therefore, in a contaminated area, more cesium-134 and cesium-137 would be taken up into the algal mat in a low-velocity area compared with in a high-velocity area. The radiocesium in stream water exists in both dissolved and suspended forms 22 Dissolved and suspended radiocesium can easily penetrate the algal mat, thus, the intake of radiocesium within the algal mat would occur regardless of the form of radiocesium. Consequently, the involvement of radiocesium in the algal mat would easily occur at low water velocities. Of course, cesium-134 and cesium-137 are also trapped in the litter mat and sand substrate. Considering this result, aquatic insects that consume the algae mat as part of their diet 15 and fish that consume the aquatic insects as part of their diet, which inhabit low-velocity areas, would also exhibit higher radiocesium concentrations via the food web.
Yoshimura and Akama 10 found no significant differences in the radiocesium levels in algae and aquatic insects between sites where the air dose rates differed (0.08-0.13 and 0.23-0.28 μSv/h) and also reported high variability in radiocesium levels in algae and aquatic insects (range: ca.300 Bq/kg). This may be due, in part, to the difference in the assemblages of the algae species attached to the stones at the sites. However, considering our results, the variability in the radiocesium concentration of algae might be due primarily to differences in the velocities of the streams where the algae were attached (ca: 500-5000 Bq/kg depend on the velocity at site D). The variability caused by the differences in stream velocity was sufficiently large to obscure the differences in contamination levels at sites with different deposition. Subsequently, the radiocesium concentration in aquatic insects that graze on algae would also show variability. Stream velocity plays a large role in habitat contamination, and might contribute to the variation in radiocesium concentration in the freshwater biota.
Particulate organic matter from the surrounding forest, which is contaminated with radiocesium, settles to the stream bottom when the flow velocity is low [23][24][25][26][27] and releases radiocesium at the stream bed when it decomposes [28][29][30] . The lack of flow in deep parts of a pool can also result in the accumulation of silt contaminated with radiocesium and free radiocesium in stagnant water. Most of the radionuclides in a pond were bound to the sediment, and aquatic insects in the pond acquire radionuclides from the sediment and ingested biotic components 14 . Consequently, fish living near the bottoms of lakes ingest material that is highly contaminated with radiocesium that had settled on the bottom. The higher concentrations of cesium-137 in brown trout (Salmo trutta) in lakes versus streams 20 suggest that the differences in contamination level are due to differences in the water flow rate. In freshwater charr (Salvelinus leucomaenis), there is substantial variation in the cesium-137 levels in equal-aged fish in a given stream 20 . Habitat selection by individual fish might cause variation in radiocesium concentrations in the diet, and levels of contamination in the diet can differ with stream velocity (habitat), even www.nature.com/scientificreports/ if the diet is the same. These factors would enhance the variability in radiocesium concentration in equal-aged fish inhabiting the same area. The transfer coefficient differed among aquatic invertebrate groups. This value was calculated on the assumption that all of the aquatic invertebrates ate algae or litter. However, the diet of all these aquatic invertebrates is not restricted only to algae or litter. Because the concentration in the litter was a little lower than that in the algae, the transfer coefficient based on litter would also be a little higher than that based on algae. Some groups such as Perlidae Gen. spp. are predators. However, the transfer coefficient for the predators did not largely differ from www.nature.com/scientificreports/ the other aquatic invertebrates, because the average radiocesium values in aquatic invertebrates were similar to those in algae and litter. Based on the differences in radiocesium values of the diet, the transfer coefficient calculated in this study also did not differ greatly from the actual value. The ecological half-life of cesium-137 ranges from 0 to more than 10 years depending on the sample type 31,32 . The value of ecological-half life in this study were largely within this range although the sample types differed. The transfer coefficient differed among aquatic invertebrate groups. Because the algae in this study were sampled from the stones in the riffle, the aquatic invertebrate groups with higher transfer coefficient would inhabit the place where radiocesium tends to be accumulated such as the low stream-velocity areas. Because the habitats of these groups with higher transfer coefficients were contaminated with radiocesium, they could not effectively exclude the radiocesium from the body. In fact, Ephemera japonica and Lanthus fujiacus inhabit low streamvelocity areas 15 . Estimated radiocesium concentrations also showed higher va in these groups. On the other hand, some other groups with lower transfer coefficients, such as Perlidae Gen. spp. and Corydalidae Gen. spp. inhabit the high stream-velocity area 15 where radiocesium does not tend to be accumulated.
The significantly higher cesium-134 and cesium-137 concentrations in algae in streams with low velocities was detected only at sites C and D, which were heavily contaminated areas. At sites A and B, which were less contaminated with radiocesium, no decline correlation between the radiocesium concentration and stream velocity was observed. Besides, at site A, the concentration tended to be higher at rapid stream velocities. Significant difference in the degree of slope coefficient between sites A and D would indicate that radiocesium dynamic are different between more and less contaminated areas.
In general, more sand and silt material is transported to high-velocity area than the low-velocity area of the stream because there is a higher volume of water in the former areas. However, in high-velocity area, most of the sand and silt passes through the stream without settling to the bottom of the stream. Because the air dose rate was higher at sites C and D, most of the silt grains might have been contaminated with radiocesium at sites C and D. Indeed, radiocesium amounts in the sand substrate (including silt) at sites C and D were higher than that at site A. Figure 3 shows the pattern diagram of sand settling to the stream bottom in highly radiocesium contaminated areas and less radiocesium contaminated areas. A rapid stream velocity would reduce the effect of contamination in algae (Fig. 3b) compared with the slow velocity (Fig. 3a). However, in the less contaminated areas, the phenomenon would result in a different outcome.
Generally, an increase in particle diameter leads to an increase in settling velocity in water 33,34 . An increase in weight also leads to an increase in settling velocity in water under the influence of gravity. Materials that are attached with substances such as radiocesium should be heavier and larger than unattached particles. Therefore, materials contaminated with radiocesium should sink to the bottom of the stream faster than uncontaminated materials 35,36 . In general, the percentage of silt and sand materials contaminated with radiocesium in the flow is not different between high-and low-velocity areas because they are connecting and form one stream. More silt materials are transported to the flow in the high-velocity areas than the low-velocity areas because of the higher volume of water in the high-velocity area. Thus, the number of silt grains with radiocesium that sink to the www.nature.com/scientificreports/ bottom would be greater in the high stream-velocity areas (Fig. 3c, d), because contaminated materials are heavier and larger than the uncontaminated materials. Consequently, the radiocesium concentrations in the algae mat in the high-velocity areas tended to be higher than those in the low-velocity areas in the less-contaminated area. Some of the aquatic invertebrate groups exhibited a positive relationship between radiocesium concentration and days, particularly at site A where the air dose rate was lower. As radiocesium concentration was low at site A, the variation in radiocesium concentration, which is dependent on environmental factors, was large enough to obscure any decrease in the concentration, and thus decreases were not detected for a few years. However, considering from the mechanism of radiocesium-contaminated silt grains accumulation at site A, a positive relationship between radiocesium concentration and the days after the accident would indicate that the ecological half-life should be far longer. Higher transfer coefficients at site A also indicate that radiocesium concentration in aquatic invertebrates were higher compared to those in algae. A continuous decrease in the radiocesium concentration would be obvious in more contaminated areas, but it is difficult to reduce the radiocesium concentration further when it has already reached a low level due to the alternation among contamination mechanisms.
Radiocesium is transported by the continuing input of radioactive substances from the catchment, and the radiocesium then would settle to the bottom of the stream. The radiocesium in stream water exists in dissolved and suspended forms 22 . The dissolved cesium-134 and cesium-137 in running water is largely not adsorbed by soil and is also readily taken up by microbes, algae and plants. This transport pathway of cesium-134 and cesium-137 eventually facilitates the uptake by aquatic invertebrates and freshwater fish at higher trophic levels in the food web. In this study, transport and transfer of Cs associated particles was discussed without considering the effects of irreversible and reversible sorption of Cs. These factors would also be essential for uptake in the biota. However, physiological experiments would be needed individually in order to fully understand these factors, and these will be the subject of future research. Difference in accumulation mechanisms between low and high velocity areas would lead to longer periods of contamination. Radiocesium is considered to be accumulated in some places such as pool or low stream-velocity area and it would also continuously affect the riffle stream substrate, especially in areas with lower concentration of radiocesiums. The transfer of radiocesium in aquatic ecosystems is also affected by other factors such as pH, temperature and dissolved organic carbon because the physiological function differs depend on the pH and water temperature [37][38][39] . Speciation of cesium and the concentration of other ions also affect the cesium transfer in the body [40][41][42] . Fishes are major functional components of forest stream ecosystems 43,44 . Radioactive contamination in fish should be avoided because they may be a source of food for human. A safety threshold of 100 Bq/kg of radiocesium was introduced in April 2012, but activity concentrations greater than this have been detected in fish hundreds of kilometres away from the FDNPP 20 . Clearly, it is critical to reduce radiocesium food contamination. The volume of water flow due to disturbances and stream velocity are key determinants of future contamination in freshwater ecosystems. Monitoring and continuous decontamination, even if the radiocesium concentration is low, are also necessary for the places where contamination levels have reached to a lower threshold.

Scientific Reports
| (2020) 10:21819 | https://doi.org/10.1038/s41598-020-78844-8 www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.