prilagođeno pretraživanje po punom tekstu

ŠUMARSKI LIST 7-8/2021 str. 21     <-- 21 -->        PDF

Distribution of 137Cs and 40K in the tissue of silver fir trees (Abies alba Mill.) from Lika (Croatia)
Distribucija 137Cs I 40K u tkivima stabala obične jele (Abies alba Mill.) iz Like (Hrvatska)
Marina Popijač
The research on activities of 137Cs and 40K, which was conducted on the silver fir (Abies alba Mill.) from Lika has included sampling of the trees in the field (rings of the bole from three different heights separated into bark, growth rings, roots, needles, shoots, and the soil surrounding the cut down trees), laboratory analysis of samples using the gamma-ray spectrometry and the statistical analysis of the collected data. The radial and vertical distribution of cesium (137Cs) in trees was investigated. 137Cs has contaminated forest ecosystems by remote atmospheric transport and radioactive precipitation as a result of nuclear test including the nuclear accident in Chernobyl. On a longer time scale, the variability of the 137Cs distribution determined in the organisms of the silver fir depended on the half-life, while the seasonal dynamics were influenced by the degree of physiological activity and the characteristics and functions of plant tissues. The highest activity of 137Cs was determined in the bark and the physiologically most active parts of the silver fir (shoots and needles). The highest activity concentration of the 137Cs in the growth rings was measured in the lowest parts of the silver fir trees. This research contributed to understanding the behavior of 137Cs, which entered the organisms of dominant tree species in the forest ecosystem, as well as its distribution in time and space.
Key words: forest ecosystem, radionuclides, distribution, bioindicators, silver fir
Forest ecosystems provide the basic requirements for maintenance of all life on Earth. Forests directly involved in the purification of groundwater and surface water, protect reservoirs of drinking water from pollution, play an important role in flood protection, protect the soil from erosion and sliding, and prevent the occurrence of avalanches and sliding snow weight (Prpić et al. 2005). The larger forest areas affecting the climate mitigation of major changes in temperature, increase humidity during summer, prevent changes in the microclimate of a particular area, protect against contaminated air, strong wind and noise, hold large amounts of dust, and provide shelter and food for numerous animal and plant species (EEA 2016). Anthropogenic pollution of forests, forest ecosystems and the overall environment is one of the largest problems of our time, which in many areas takes on catastrophic proportions. Many pollutions in forest ecosystems are monitored systematically, their appearance, performance, duration and disappearance are recorded. The radioactive cesium isotope l37Cs appeared in the environment in significant quantities after atmospheric nuclear testing during the 1950’s and 1960’s. (UNSCEAR 1993, FAS 2002), and a large amount of l37Cs

ŠUMARSKI LIST 7-8/2021 str. 22     <-- 22 -->        PDF

was discharged into the atmosphere during the nuclear accident at Chernobyl on April 26, 1986.
The deposition of radionuclides on the tree canopies is the most important contaminating form of 137Cs in forest ecosystems (Skoko et al. 2011). Given that the tree canopies are a major part of the biomass in the management of low aerodynamic resistance, they are effective “interceptor” of the radionuclide (Yamagata et al. 1969, Bunzl and Kracke 1988, Desmet and Myttenaere 1988, Sokolov et al. 1990). Generally, such “interception” in the trees (seen in the same unit volume of crown) is more effective for dry than wet deposition, especially for small particles and gases (Pröhl 2008). As the amount of 137Cs over time decreases exponentially, adoption root becomes the dominant source of intake of 137Cs in the trees long term (Ohashi et al. 2020). Accumulation of 137Cs in the trees is one of the most important problems in the wood use and forest pollution (Shinta et al. 2014). Wood as a raw material is permanently present in the environment of people, contributing to their radioactive exposure by increasing doses of radiation. However, compared to the annual equivalent dose per capita from various sources, the calculated annual equivalent doses do not pose a risk in beech, oak and silver fir wood panel samples tested from the territory of the Republic of Croatia (Hus et al. 2004).
The objective of the research was to determine the incidence and intensity of the dynamics of the 137Cs in the silver fir tissues (Abies alba Mill.), including their by-products entering the human and animal food chain, as well as in samples of soil directly under the selected trees. At study sites, on a longer time scale (2003 - 2017), measured variability of 137Cs in the tissue of silver fir trees depend on seasonal dynamics of 137Cs followed the level of physiological activity and decay of 137Cs during time. The highest activity of 137Cs was determined in the bark and the physiologically most active parts of the silver fir (shoots and needles) during the growing season. Activities of 40K were of approximately equal values throughout the observed period. This research will contribute to the understanding of the biogeochemical cycle of 137Cs within the forest ecosystem, and will also contribute to the monitoring of 137Cs wood raw material pollution on wood products. This research also tends to contribute to the study review discussion of the consequences of human activity has on the forest ecosystem.
Research area – Područje istraživanja
The research was conducted on the tissues of fir trees (Abies alba Mill.) in Lika (near Vrhovine), and the soils next to the trees. The chosen location is considered to be among the most contaminated during the Chernobyl disaster in Croatia (Barišić and Lulić 1990), and the sampling schedule is shown in Table 1.
Sample collection method – Metoda prikupljanja uzoraka
In the forest ecosystem at the site GJ Komarnica, subdivision 12a, studied the trees of silver fir (1, 2, 3) and recorded the coordinates of trees and soils samples. The measured trees were knocked down in the field by licensed workers of the company Hrvatske šume Ltd. after which the cut rings were used for further sectioning. Pits were dug next to fir trees in which root samples were collected for testing. The soil samples were collected along with the tree samples. The first sampling on a silver fir tree, 25 m high and 56 cm in (DBH) diameter at breast height, was performed during the vegetation dormancy on December 4, 2003. On a felled

ŠUMARSKI LIST 7-8/2021 str. 23     <-- 23 -->        PDF

tree silver fir rings were taken at a height of 0.1 m, 8 m and 16 m and then stratified on samples of biological material, that is, samples of bark and sapwood, and the growth rings with the time steps to the center of the trunk were set aside. On the sawn-off reel of a felled tree at a height of 0.1 m, rings were singled out in sections of two years for the period 1974 - 2003, in sections of four years for the period 1944 - 1973 and then one sample for the period from 1909 – 1943. The same collection dynamic is used for tree rings samples at altitudes of 8 and 16 meters, provided that at 4 meters the sample covered the period 1958 - 2003, and at a height of 16 meters the period 1978 - 2003. The main root and lateral roots (bark and wood) were sampled, and the group of samples were split according to diameter dimensions (<1 mm, 1-2 mm, 3-4 mm, 5-8 mm, 9-16 mm, 17 - 24 mm and 25 - 38 mm). Twigs (peak shoots) were sampled as individual annual samples for shoots that grew in the period 1994 - 2003 at heights of 8 m and 16 m, and at a height of 25 m in the period 1997 - 2003.
The second sampling on the silver fir tree was done during vegetation period on September 3, 2004. On the felled tree of 24.5 m height and 54 cm breast diameter, the rings were takes as in the first sampling, and then stratified on samples of biological material. On the reel of a felled tree, at a height of 0.1 m, rings were singled out in sections of five years each for the period 1945 - 2004, and in sections of ten years for the period 1903 - 1944. Samples were collected on the rings of a felled tree at heights of 8 and 16 meters at the same dynamics, with the period 1927 - 2004 being covered at eight meters, and the period 1960 - 2004 at a height of 16 meters. The main root and lateral roots (bark and wood), were sampled and separated in groups of samples as well as in the first sampling. Twigs (peak shoots) trees were sampled as individual annual samples for shoots that have grown in the period 1994 - 2004 at the height 8 m and 16 m, and a height of 25 m for the period 1996 – 2004. In the first and second sampling, along with peak shoots, samples of peak shoots needles were collected at the heights of 8 m, 16 m and 25 m.
Third sampling on the tree silver fir was conducted on March 27, 2017 for the purpose of comparing partial results after a long period since the last sampling (12 years). At Komarnica, the forest management unit, in subdivision 12a, felled the tree silver fir 26 m in height and 59 cm breast diameter, and the taken disc was then stratified on samples of biological material, that is, samples of the bark, bark and cambium, and annual rings with time intervals to the center of the trunk were singled out. The samples for this analysis were selected at 0.1 m height as the lowest level of the tree where the active life cycle takes place, and that is comparable to the first two samplings to determine the activity of 137Cs. In the period 1974 - 2016, the years were separated into sections of two years each, and for the period 1939 - 1973 in sections of five years, and one sample was separated for the period 1887 - 1938. Twigs (peak shoots) trees were sampled as individual annual samples for shoots that have grown at the height of 26 m in the period 2013 - 2016. Along with peak shoots, samples of peak shoots needles were collected. Composite soil samples were collected with felled trees in 2003, 2004 and 2017 for the analysis of 137Cs concentration in soils by gamma spectrometric method. For each analysis, approximately 1 kg of the average sample (taken from multiple points (5), homogenized, disrupted in a state) was collected.
Sample preparation and gamma spectrometric method – Priprema uzoraka i gamaspektrometrijska metoda
The activity of all samples was measured by gamma spectrometric method in the Laboratory for Radioecology, Division for Marine and Environmental Research, of the Ruđer Bošković Institute which is accredited for performed gamma spectrometric measurements since 2008, according to HRN EN ISO/IEC 17025. Before the measurements, all samples were homogenized, dried in an oven to constant mass, and then placed in the measuring geometry 125 cm3 volume, weighed and sealed. Canberra’s HPGe (High Purity Germanium) detector system was used for measurements, where each sample was measured for 80.000 seconds. The Genie 2000 program from the same manufacturer was used to analyze the recorded gamma spectra. The concentration of 137Cs activity in the samples was determined from the photo-peak at an energy of 661.6 keV, and the concentration of 40K activity from the photo-peak at an energy of 1460.7 keV. The reliability of the system was checked during intercomparison measurements regularly.
Since the analyzes showed that the composite uncertainty (U) is predominantly influenced by the uncertainty of the counting rate (Ubb), which also includes the uncertainty of the counting rate of the basic radiation, the uncertainty of determining the efficiency (Uef) and the uncertainty of determining the surface of photo-peak (Ufv - especially in the case of a small number of registered events) these parameters were taken for the calculation of the measurement uncertainty according to the relation:
                U2 = Ubb2 + Uef2 + Ufv2
Extended measurement uncertainties since 2008 were expressed with an overlap factor k = 2 (95% reliability of the results) as opposed to measurements performed until 2008, when measurement uncertainty calculations included only the uncertainty of determining the photo-peak area with an overlap factor k = 1 (68.27% reliability of the results). The measurement uncertainty was calculated as described above and shown in the results. At lower activities, the measurement uncertainty is high, but we consider the obtained results

ŠUMARSKI LIST 7-8/2021 str. 24     <-- 24 -->        PDF

to be statistically reliable. The reason for this is that the measurements are made in such a way to ensure the same measurement conditions each time (geometry of the sample, measurement time, the same detector and calibration).
Calculation elimination of radioactive decay impact – Računska eliminacija utjecaja radioaktivnog raspada
Given the large time gap between the first and second sampling and measurements (2003 and 2004) on the one hand, and the third sampling and measurements (2017) on the other hand, it could be assumed that the measured activity of 137Cs in the third measurement, due to the half-life of 137Cs, would be lower for about a quarter than in the samples from 2003 and 2004. The law of radioactive decay is well known and allows the calculation of radionuclide activity at any point in time. Therefore, it is possible to eliminate the influence by computation of radioactive decay for the measurement of 2017, in a way that all of the measured activity of 137Cs in the samples from the three measurements, recalculate for the date of July 1, 2003 allowing the comparison of the first and second measurements, and for July 1, 2016, when approximately one half-life of 137Cs has expired since the contamination (1986). All samples were computationally analyzed in the described manner, and the results for all samples were shown and interpreted in total (as an arranged pair of values for July 1, 2003 and July 1, 2016), ensuring that the interpretation of results of the difference between recalculated value of the third sampling and the first and second sampling can be  exclusively attributed to geochemical processes that occur in the environment over time (given that the impact of radioactive decay is eliminated by calculation).
Distribution of 137Cs and 40K in the tissues of silver fir – Distribucija 137Cs i 40K u tkivima obične jele
Distribution results of 137Cs in the tissues of silver fir (Abies alba Mill.) refer to the peak shoots and needles, the rings, the bark and cambium and root, and to reviewed comparative presentation of results in all silver fir sampled tissues.
Distribution of 137Cs in the peak shoots and needles of silver fir – Distribucija 137Cs u vršnim izbojcima i iglicama obične jele
Figure 1 shows the distribution of 137Cs activity in silver fir trees shoots of different heights and ages in all three sampling years, showing also the distribution according to the age of shoots (and not to the calendar year).
Comparison of age sequences of 137Cs activity in shoots of silver fir (Figure 1) for individual height between two years (2003 and 2004), using the non-parametric test for related samples (“Wilcoxon matched pairs test”), shows that the values of 2004 are statistically significantly (p = 0.05) higher at altitudes of 8 and 16 m than those in 2003, while at the 25 m difference between the two was not statistically significant. This could be interpreted as a result of an increase in the photosynthesis intensity, which is in the fir, as coniferous species, present throughout the year and of the trees metabolism level during the growing period, leading to a temporary accumulation of 137Cs in the shoots.
Figure 2 show the distribution of 137Cs in silver fir needles at different heights and different age in three sampling years.

ŠUMARSKI LIST 7-8/2021 str. 25     <-- 25 -->        PDF

Comparison of age sequences of 137Cs activity in shoots of silver fir with those in silver fir needles (Figure 2) for individual height between two years (2003 and 2004), also using non-parametric test for related samples (“Wilcoxon matched pairs test”), shows that the 2004 values are statistically significantly (p = 0.05) higher than in those in 2003, at all three heights. This result, very similar to that of the shoots, can be equally interpreted (as a result of increasing the rate of photosynthesis and metabolism levels during vegetation period, leading to a temporary accumulation of 137Cs in the needles). The regularity is shown even at a height of 25 m (not in the shoot case).
Comparison of age sequences of 137Cs activity in silver fir shoots with those in needles (Figures 2 and 3) was conducted separately for each height and year, also using the non-parametric test for related samples. The obtained results show that for the same year the values in the shoots are always statistically significantly (p = 0.05) higher than those in the needles of the same age at all three heights. This result can be associated with a smaller amount of flow in the peripheral tissues of the tree conduction system, which results in smaller amounts of 137Cs entering/appearing in these tissues. Figure 3 shows a transparent distribution of 137Cs through shoots and silver fir tree needles of different ages, separately for each height (8, 16 and 25 m). This figure shows what has already been pointed out: 1) higher values in 2004 during the vegetation period compared to those in the winter period of 2003 (at all heights, both in shoots and needles) and 2) generally higher values in shoots relative to those in needles (for both years, at all tree heights). Additionally, in this picture (looking only at the medians of age series related to the same tissue type in the same year connected by lines), it can qualitatively determine an overall increase in 137Cs activity in shoots and needles of the tree heights (regardless of age). Although due to the small number of sample point (three heights) it is not possible to test this increase for significantly positive correlation, it is appropriate to point out that such positive correlation is consistent with other results presented. That is, when viewed within the same type of tissue, the largest activity of 137Cs can be expected in areas of greatest physiological activities, which, if we talk about the shoots and needles, should be higher in the tissue sample from higher stem (having

ŠUMARSKI LIST 7-8/2021 str. 26     <-- 26 -->        PDF

less overshadowing impact of their own and other crowns, resulting in more radiated solar energy and higher intensity of photosynthesis).
Distribution of 137Cs in the annual rings of silver fir – Distribucija 137Cs u godovima obične jele
Figure 4 shows the distribution of 137Cs above the limit of detection in all sampled annual rings of silver fir trees at heights of 0.1, 8 and 16 m for samples from 2003. It is noted that in the annual rings of silver fir at all altitudes present much greater variability of 137Cs activity that is measurable way back in the past, reaching to the very center of the three rings (one on each level) from the fir tree felled for this survey in 2003, which was, as well as two trees felled in 2004 and 2017, more than hundred years old.
These long time series of radioactivity values can be divided into four intervals based on a graphic insight into the measured data (Figure 4):
1) In the period from the beginning of life trees from 1910 to 1943, constant value of 137Cs activity of 28.4 ± 0,9 Bq/kg (both in terms of the radius of the coil and in terms of the height of the tree) was measured.
2) The period 1944 to 1958, almost constant (in terms of the radius of the disc and the height of the tree) 137Cs activity was observed, varying from 22.7 ± 0.8 Bq/kg to 23.3 ± 0.9 Bq/kg,
3) For the period 1959 to approximately 1983, variation of 137Cs activity was also recorded in tree rings (without obvious regularity with regard to the year and the height of sampling), in the interval between 14.3 ± 0.7 Bq/kg, (16 m; 1982) to 36.3 ± 1.2 Bq/kg (8 m; 1974),
4) For the period from approximately 1984 to 2003, in which, a continuous increase in 137Cs activity was recorded in the youngest rings.
Based on the same results, it could even be assumed that the old xylem tissues that no longer participate in the flow from the root to the canopy represent a kind of accumulation reservoir of 137Cs, leaving the question of whether 137Cs in these (mostly dead) tissues is permanently excluded from migration within the plant and “captured” (but can disappear only by radioactive decay) or it is a temporary seasonal accumulation during the winter period (in which case, the constancy of the recorded values could explain the assumed constant capacity of dead wood). Additionally, given the existence of zone 3. in the dendrochronological sequence of fir, it can be assumed that the transition between the physiologically active and physiologically inactive part of the xylem in felled fir was very gradual (extending to more than thirty years of the dendrochronological sequence), and moreover, anatomically non-homogeneous (due to the variability of the recorded disc height sampling).
Figure 5 shows the distribution of 137Cs activity above the detection limit in all sampled annual rings of silver fir trees at heights of 0.1, 8 and 16 m for samples from 2004.
Graphic insight into the measured data shows the following:
1. Measurable 137Cs activity is present up to the oldest rings, as with the 2003 annual rings
2. The variability of 137Cs activity in the dendrochronological sequence from 2004 is lower than in those from 2003, which is mostly due to lower values in the youngest (and physiologically probably the most active) years (2004 compared to 2003), which could be attributed to the fact that 2004 samples were from the vegetation season (as a possible cause of increased migration of 137Cs from conductive tissues to the most physiologically active tissues in the canopy) and
3. Significantly weaker possibility of defining typical zones within the dendrochronological sequence in relation to the series from 2003, although the zone of increased 137Cs values in the youngest rings is also visible here, and,  on the other hand, the zone of almost constant 137Cs values in the oldest rings is well noticeable (on the basis of which it could be assumed that during the vegetation season there is even a migration of 137Cs from the inactive to the physiologically active xylem zone, and in the winter 137Cs is temporarily accumulated even in the physiologically inactive part of the xylem zone; compare above for dendrochronological sequence from 2003).
Divided time sequences can be observed at intervals for: a) the period from the beginning of tree life to 1958, in which a similar level (in terms of circle radius and tree height) of 137Cs activity was recorded, which varied from 15.2 ± 0.8 to 19.4 ± 0.9 Bq/kg at a height of 0.1 m, and from 20.0 ± 1.4 to 23.2 ± 1.5 Bq/kg at a height of 8 m, b) the period of 1959 to approximately 1982, in which a variation in 137Cs activity was recorded, both in terms of ring radius and in terms of tree height (without distinct regularity with respect to the  year formation), ranging between 10.4 ± 0.7 Bq/kg (0.1 m; 1980 - 1984 interval) to 25.1 ± 1.5 Bq/kg (16 m, the interval 1970 - 1974) and c) the period of about 1985 to 2004, in which substantially continuous increase of 137Cs is recorded in the most recent annual rings.
Figure 6 shows a comparison of 137Cs activity in the wood rings of silver fir trees (at a height of 0.1 m) between the average values of samples from 2003 - 2004, on one hand, and recalculated (on July 1, 2003) values from 2017 on the other.
Graphic insight into the measured data shows the following: a) lower values in recalculated samples of 2017 compared to the samples from the period 2003 - 2004, b) still measurable 137Cs activity along the dendrochronological sequence from 2017, and c) low variability in 137Cs activity

ŠUMARSKI LIST 7-8/2021 str. 27     <-- 27 -->        PDF

throughout the 2017 series (including the absence of a more evident increase in the youngest years).
Figure 7 shows the distribution of 137Cs in the annual rings of silver fir for different height trees and old growth rings in all three sampling years. These are the same values as in the charts of Figures 4 - 6, additionally displaying the distribution by age rings (and not according to the calendar year), and with a limited set of twenty early years within which the samples from 2003 and the 2004 recorded an increase in 137Cs activity with decreasing rings age.
Comparison of age series 137Cs activity in the annual rings of silver fir for different heights and year of sampling, was done using non-parametric test for related samples (“Wilcoxon matched pairs test”), Figure 8, and shows the following:
1) At the tree heights of 0.1 m to 8 m statistically significant (p = 0.05) increases were observed in 2003 (while at the height of 16 m this difference was not significant), which can be attributed to a vegetation dormancy in winter and

ŠUMARSKI LIST 7-8/2021 str. 28     <-- 28 -->        PDF

2) There was no statistically significant difference between series from the same sampling year, and from different heights, which can be interpreted as a homogeneous height distribution of 137Cs activity in the twenty youngest years of the tree rings.
Distribution of 137Cs in bark and cambium of silver fir – Distribucija 137Cs u kori i kambiju obične jele
In Figure 9, (comparing only the same tissue type in the same year) a general increase of 137Cs activity in the bark and cambium can be found correlating with the tree height (except at 0.1 m in 2003). In 2003, the highest concentration of 137Cs was measured in bark samples collected at a tree height of 16 meters, measuring 130.1 ± 1.5 Bq/kg, while the value of the measured sample at the same height in 2004 was 50.7 ± 1.2 Bq/kg. At the same height (16 m), the 137Cs activity in cambium was 124.1 ± 1.4 Bq/kg measured in 2003, and in 2004 it was 79.2 ± 1.6 Bq/kg. In the samples of the bark on tree height of 8 meters in 2003 measured 137Cs activity was 83.8 ± 1.4 Bq/kg, and in 2004 was 38.7 ± 1.0 Bq/kg. At the same height (8 m) in 2003 measured 137Cs activity in the cambium was 77.6 ± 1.4 Bq/kg, and in 2004 it was 65.2 ± 1.4 Bq/kg. At a height of 0.1 m, the measured 137Cs activity in the bark sample of 2003 was 70.6 ± 1.3 Bq/kg, and in 2004 sample it was 24.0 ± 0.3 Bq/kg. At the same height (0.1 m), the 137Cs activity in the cambium was 80.5 ± 1.5 Bq/kg measured in 2003, and in 2004 it was 27.0 ± 0.3 Bq/kg. The values of the measured 137Cs activity in 2017 and recalculated to day July 1, 2003 amounted to 13.8 ± 2.0 Bq/kg in bark and to 18.4 ± 2.67 Bq/kg, in cambium, recalculated to day July 1, 2016 amounted 10.2 ± 1.48 Bq/kg in bark and 13.6 ± 1.3 Bq/kg in cambium.
Distribution of 137Cs in the root of the silver fir – Distribucija 137Cs u korijenu obične jele
Distribution of 137Cs activity in the root thickness classes of silver fir of 2003 and 2004 samples is shown in Figure 10. The samples from the 2003 showed of 137Cs activity variation ranging from 62.1 ± 1.3 Bq/kg (root thickness 1-2 mm) to 98.4 ± 1.8 Bq/kg, while in samples from 2004, the range was from 31.4 ± 1.0 Bq/kg (root thickness 3-4 mm) to 71.2 ± 2.3 Bq/kg (tiniest roots). Meanwhile, the value nearest the maximum (recorded in the thinnest roots) for both years were recorded in the thickest roots (25-38 mm; 95.5 ± 2.7 Bq/kg in samples from 2003; 44.9 ± 1.5 Bq/kg in samples from 2004) which at the qualitative level indicate the absence of a link between root thickness and 137Cs activity. This conclusion is supported by the analysis of the association of 137Cs activity with the upper limit of the root thickness class, performed by the Spearman-rank correlation method, suitable for small samples, which did not result in a significant correlation (RSP = 0.1429 with p = 0.7599 for 2003). RSP = - 0.3214 with p = 0.4821 for 2004). Accordingly, it can be concluded that no significant variability in 137Cs activity has been observed in the fir root system, which could be due to the fact that these are plant tissues that overgrow the soil from which the tree draws 137Cs. By comparing the values of two sampling years separately for each diameter class roots (Figure 10), we can see that in 2003 (winter period) each diameter class recorded higher values than in 2004 (growing season). This result could be explained by the assumption that in winter, because of slower metabolism and less intensive flow through the vascular tissue of roots, roots temporarily accumulate 137Cs, whose

ŠUMARSKI LIST 7-8/2021 str. 29     <-- 29 -->        PDF

concentration during the vegetation period decreases due to migration of 137Cs in physiologically most active tissues (shoots and needles).
Comparison of 137Cs activities in different tissues of silver fir – Usporedba aktivnosti 137Cs u različitim tkivima obične jele
Figure 11 shows the distribution of 137Cs activity in rings, shoots, needles of silver fir tissues (the results from the youngest tissue samples generated in the sampling year for tissue types for which is possible to determine the age) for 2003 (winter sample), 2004 (vegetation season sample) and values from the respective samples collected in 2017 and recalculated (eliminated impact of radioactive decay) on July 1, 2003.In Figure 11, it can be primarily observed that in tissues (peak shoots, needles) from samples collected in 2004, values are significantly higher than in previously collected samples, which directly proves the presence of 137Cs redistribution in tree tissues over time (given that the absence of such redistribution over time should decrease the value due to radioactive decay). The existence of redistribution is confirmed by the results in the tissues in which the 137Cs activity was higher in 2003, because these differences (up to three times higher values for samples of bark and cambium near the soil, compared to 2004) were significantly higher than those that would result due only to radioactive decay during less than a year. Furthermore, it can be seen that in the winter period samples higher values of 137Cs were recorded in samples of dead bark, live bark with the cambium, and the young rings and roots, while in the growing season samples higher values were recorded in peak shoots and needles samples. It can be observed that for all tissue types the recalculated values from the samples collected in 2017 are

ŠUMARSKI LIST 7-8/2021 str. 30     <-- 30 -->        PDF

significantly less than the relative values from 2003 and 2004, but that these differences are not the same for all tissue types. These results can be explained by tree metabolism slowing down  and the intensity of photosynthesis during the vegetation dormancy, which results in poorer flow through the conducting tissues and temporary accumulation of 137Cs in the phloem, xylem and root, while on the other hand during the growing season, boosted metabolism and intensive circulation of substances in the relation root - xylem - needles - phloem – root,  the temporary accumulation of 137Cs is seen in the top shoots and needles.
Relationship of 40K and 137Cs activities in the tissues of silver fir – Odnos aktivnosti 40K i 137Cs u tkivima obične jele
Behavior and distribution of cesium in plants and organisms is similar to potassium because these elements are homologous (Shaw and Bell 1991, Robinson and Stone 1992). Research (Lovrenčić et al. 2008) confirms that 137Cs is involved in physiological processes in a similar way as potassium. For this purpose, a 40K distribution was measured simultaneously with that for 137Cs, in the same samples by gamma spectrometric measurement. The scatter plot in Figure 12 shows the ratio of the measured 40K and 137Cs activities in fir tissues (for twigs, needles, rings, bark and cambium). The obtained result is expressed collectively for all tissue types and shows the ratio of measured values for 40K and 137Cs, where each point on the graph represents an arranged pair of values for 40K and 137Cs measured in the same tissue type, sampled at the same time (the same year) and at the same height of an ordinary fir tree (when it was sampled from several heights for some type of tissue). We can observe a large scattering of data, with a linear correlation that is statistically significant (R = 0.4999; with p = 0.000), so this result is not inconsistent with the fact that 40K and 137Cs are homologous isotopes that the plant does not distinguish. The grouping of points from samples of different tissue types can also be observed on the same graph, so it was appropriate to perform the same analysis separately for individual tissue types.
Distribution of 40K and 137Cs in the soil – Distribucija 40K i 137Cs u tlu
While the focus of the research was the distribution of 137Cs activity in the plant tissue of silver fir, we collected composite samples of forest soil at the same site Vrhovine, on which we measured activity of 40K and 137Cs (Bq/kg) by gammaspectrometer in all three sampling years (Table 2). The level of 137Cs activity in the fir forest was much higher in the samples of soil collected in 2003 and 2004 in relation to the activity levels measured on 2017 samples (recalculated on day July 1, 2003 and July 1, 2016). The 40K activity level in 2017 was twice higher as in 2003 and 2004.
Activity of 137Cs in the tissues of silver fir trees – Aktivnost 137Cs u tkivima drveća obične jele
The results of this study showed that the distribution of 137Cs activity within each individual tissue was causally related to the degree of physiological activity in a particular tissue. In the radial distribution, this is reflected in an increased activity value with decreasing age (the maximum activity was observed in the youngest wood rings), while the vertical distribution is reflected in the increased activity correlated to the height of the tree. In the underground part (root samples) the similar dependence was measured in the vegetation period, while during the vegetation dormancy the level of 137Cs activity did not vary significantly with the reduction of the root dimensions. The results of this study show similar dynamics of 137Cs movement in fir tissues (Abies Alba Mill.) as in previous results by Popijač et al. (2004) fifteen years

ŠUMARSKI LIST 7-8/2021 str. 31     <-- 31 -->        PDF

after the Chernobyl accident on Medvednica (Sljeme) and Lovrenčić et al. (2008) in Gorski Kotar in which they monitored the dynamics of l37Cs movement on a monthly basis throughout the year and found higher 137Cs activities in the bark of fir and heartwood during the growing season, when the trees were physiologically more active. The previous results of the two studies suggest that silver fir is a good model species for 137Cs biomonitoring because they behave similarly in different habitats, although this claim would need to be extended to more localities.
Relationship between 40K and 137Cs activities in trees and soil – Odnos aktivnosti 40K i 137Cs u drveću i tlu
To understand the possible explanations of long-term presence of 137Cs in the forest ecosystem, this study, measuring concentration levels of 137Cs activity in the tissues of silver fir and soils near sampled trees, also covered the level of 40K activity, trying to explain the relationship between 40K and 137Cs activities in trees and soil, considering that the two isotopes are homologous. Research (Lovrenčić et al. 2008) also confirms that 137Cs is involved in physiological processes in a similar way as potassium, since these two elements are homologous. Their activity is inversely proportional to the age of the needles and twigs, that is, the highest is in the youngest tissue sections. The soil characteristics play a key role in the transfer of 137Cs, based on their texture, ability to exchange cations and organic matter content. Clay soils accumulate 137Cs and it was found that if the proportion of clay in the soil is higher, the 137Cs activity in the soil is also increased, because its sorption is conditioned by clay minerals (Kruyts and Delvaux 2002, Stauton et al. 2002). The content of K+ in the soil can cause the collapse of the extended interlayers (Rigol et al. 2002), because in this case the 137Cs is binding within the interlayer and is blocked, and unavailable for transfer processes. Potassium availability is strongly associated with sorption and desorption processes, as well as fixations that take place in the soil (Đurđević 2014). Previous research on the relationship of 40K and 137Cs in soils has been primarily related to measuring and determining the distribution and level of 137Cs activity in soils after the Chernobyl and Fukushima accidents (Smolders et al. 1997, Rochon et al. 1998, Gerzabek 1996, Zhu and Smolders 2000, Kruyts and Delvaux 2002; Kaunisto et al. 2002; Zibold et al. 2009). Their importance was highly related to dealing with the consequences that occurred after the entry of 137Cs in the environment and taking measures to secure food production and animal farming, but also protect forest ecosystems.
The study of distribution of 137Cs in the trees of silver fir (Abies alba Mill.) from Lika is based on three occasions within the interval 2003 - 2017 which included field ring samples of trees with three heights (separated on the bark and growth rings), roots, needles, peaks of shoots, and soils near felled trees, and after laboratory processing of collected samples in the gamma spectrometer, statistical analysis and interpretation of the data obtained, resulted in the following conclusions:
1. Following  the results of previous researches in a similar but wider area (which was limited to 137Cs activity in soil and honey, Barišić et al. 2018), this study confirmed long-term contamination of Lika forest ecosystems with radioactive precipitation that had entered the Vrhovine area by remote atmospheric transport after the nuclear tests and the Chernobyl nuclear power plant accident, which manifested itself in measurable 137Cs concentrations in plant tissues at least three decades after contamination.
2. Throughout the observed period (2003 - 2017), the variability of 137Cs distribution in the tissues of silver fir trees decreased continuously. The reason for this was on one hand radioactive decay (about 26% in all tissues during the observed period), while on the other hand it was due to the gradual elimination of 137Cs from tree tissue, mostly through dead bark, and needles, which caused a decrease in 137Cs activity during the observed period up to 96% in needles, peak shoots and fir rings (in relation to the recalculated values, after elimination of the influence of radioactive decay).
3. The highest 137Cs activity of silver fir tissues were measured in peak shoots and needles, and significantly lower activity in other tissues, regardless of the time of year. The conclusion is that the main elimination pathways of 137Cs are needles typically several years old and which the tree gradually sheds throughout the year.
4. In the soil habitat of the silver fir 137Cs migrated deeper, towards the root, and thus become biologically available.
5. The recorded concentration of 137Cs and 40K are statistically significant correlations for the total sample, and separately for the needles sample. These results (given that no statistically significant negative correlation was observed, while a significant positive correlation can still be interpreted by a random process) suggest that the tree does not distinguish between these two homologous elements.
6. The results of the research (the first of its kind in Croatia focused on the edificatory tree species tissue) contribute to the understanding of future of 137Cs that entered the tissue of edificatory tree species in the forest ecosystem, as well as its distribution in time and space, which (especially in comparison with rare similar research in the world) complements the emergence of the 137Cs biogeochemical cycle in the environment (especially in forest trees).
Many thanks Ph.D. Delko Barišić, Ph.D., Full Professor Oleg Antonić, Full Professor Nikola Kezić and Matija

ŠUMARSKI LIST 7-8/2021 str. 32     <-- 32 -->        PDF

Volner for all the advice and help in conducting this extensive research. I would like to thank my dear colleagues from Hrvatske šume Ltd., the Forest Administration Sisak, Petrinja forestry and the Forest Administration Gospić, Vrhovine forestry, who enabled me to collect samples for analysis. I would like to thank all the employees of the Laboratory for Radioecology at the “Ruđer Bošković” Institute for their support in the phase of conducting the laboratory analysis of the collected samples.
Barišić, D., S. Lulić, 1990: The contamination of ground surface layer in Republic Croatia as the consequence of Chernobyl accident. Proc. Int. Symp. Post-Chernobyl Environ. Radioac. Stud. in East European Countries, 20-25., Kazimierz, Poland.
Barišić, D., A. Vertačnik, S. Lulić: 1999: Caesium contamination and vertical distribution in undisturbed soils in Croatia. Journal of Environmental Radioactivity, 46: 361-374.
Barišić, D., N. Kezić, D. Bubalo, M. Dražić, L. Svečnjak, I. Seletković, Ž. Zgorelec, M. Popijač, D. Barišić, I. Tucaković, 2018: Mogućnost detekcije onečišćenja okoliša radiocezijem (137Cs) određivanjem aktivnosti 137Cs u uzorcima meda, Civitas Crisiensis: Radovi Zavoda za znanstvenostraživački i umjetnički rad Koprivničko-križevačke županije: 3, 71 – 86., Križevci
Bunzl, K., W. Kracke, 1988: Cumulative deposition of Cesium-137, Plutonium-238, Plutonium-239, Plutonium-240 and Americium-241 from global fallout in soils from forest, grassland and arable land in Bavaria, West Germany. Journal of Environmental Radioactivity, 8 (1): 1-14.
Canberra, 2000: Genie 2000 Basic Spectroscopy Software
Corine Land Cover 2000: izvor:
Desmet, GM., C. Myttenaere, 1988: Considerations on the role of natural ecosystems in the eventual contamination of man and his environment. Journal of Environmental Radioactivity, 6: 197-202.
Đurđević, B., 2014: Praktikum iz ishrane bilja, Sveučilište Josipa Jurja Strossmayera u Osijeku, Poljoprivredni fakultet, ISBN: 978-953-7871-30-7, Osijek
EEA (Europska agencija za okoliš): 2016
FAS (Federation of American Scientists): 2002
Gerzabek, M.H., 1996: Soil-to-plant transfer of Cs and Sr in Austria after the Chernobyl accident. Mitt. Österr. Bodenkundl. Ges., 53: 111-117., Vienna
Hus, M., K. Košutić, S. Lulić, 2004: 137Cs in wood, 11th IRPA (
IAEA (International Atomic Energy Agency), 2002: Modelling the migration and accumulation of radionuclides in forest ecosystems. Report of the Forest Working Group of the Biosphere Modelling and Assessment (BIOMASS) Theme 3. IAEA-BIOMASS-1, International Atomic Energy Agency, Vienna
IAEA (International Atomic Energy Agency), 2010: Handbook of Parameter Values for the Prediction of Radionuclide Transfer in Terrestrial and Freshwater Environments. Technical Report Series No. 472, Vienna
ICRP (The International Commission on Radiological Protection), 1989
Kaunisto, S., L. Aro, A. Rantavaara, 2002: Effect of fertilisation on the potassium and radiocaesium distribution in tree stands (Pinus sylvestris L.) and peat on a pine mire, Environmental Pollution, 117: 111-119.
Kruyts, N., B. Delvaux, 2002: Soil organic horizons as a major source for radiocesium biorecycling in forest ecosysstems. Journal of Environmental Radioactivity 58: 175-190.
Lovrenčić, I., M. Volner, D. Barišić, M. Popijač, N. Kezić, I. Seletković, S. Lulić, 2008: Distribution of Cs-137, K-40 and Be-7 in silver fir-tree (Abies alba L.) from Gorski kotar, Croatia. Journal of Radioanalytical and Nuclear Chemistry, 1: 71-79.
Ohashi, S., K. Kuroda, T. Fujiwara, T. Takano, 2020: Tracing radioactive cesium in stem wood of three Japanese conifer species 3 years after the Fukushima Dai-ichi Nuclear Power Plant accident, Journal of Wood Science, 66: online available,
Popijač, M., 2018: Distribucija 137Cs u pitomom kestenu (Castanea sativa Mill.) s Banovine i običnoj jeli (Abies alba Mill.) iz Like (Hrvatska), Doktorska disertacija, Sveučilište Josipa Jurja Strossmayera u Osijeku, Institut Ruđer Bošković, Zagreb
Popijač, M., I. Seletković, M. Volner, I. Lovrenčić, D. Barišić, N. Kezić, 2004: Dinamika kretanja 137Cs i 40K na stablima jele (Abies alba) na Sljemenu, Šumarski list. 5-6: 269-277.
Pröhl, G., 2008: Interception of dry and wet deposited radionuclides by vegetation. Journal of Environmental Radioactivity, online available,
Prpić, B., P. Jurjević, H. Jakovac, 2005: Procjena vrijednosti protuerozijske, hidrološke i vodozaštitne uloge šume, Šumarski list, 13:186-194.
Rigol, A., M. Vidal, G. Rauret, 2002: An overview of the effect of organic matter on soil – radiocaesium interaction: implications in root uptake, Journal of Environmental Radioactivity, 58: 191-216.
Robinson, W.L., E.L. Stone, 1992: The effect of potassium on the uptake of 137Cs in food crops grown on coral soils, Coconut at Bikini atoll. Hlth. Phys., 62: 496-511.
Rochon, P., D. Pare, C. Messie, 1998: Development of an improved model estimating the nutrient content of the bole for four boreal tree species, Canadian Journal of Forest Research, 28, 37–43.
Shaw, G., J.N.B. Bell, 1991: Competitive Effects of Potassium and Ammonium, on Caesium Uptake Kinetics in Wheat, Journal of Environmental Radioactivity, 13: 283-296.
Shinta, O., O. Naok, T. Atsushi, N. Wataru, T. Shigeyoshi, 2014: Radial and vertical distributions of radiocesium in tree stems of Pinus densiflora and Quercus serrata 1.5y after the Fukushima nuclear disaster, Journal of Environmental Radioactivity, 134: 54-60.
Skoko, B., Marović, G., Babić, D., I. Vicković, 2011. 137Cs na znanstveno-istraživačkom poligonu „Šumbar“. Zbornik radova VIII. simpozija Hrvatskog društva za zaštitu od zračenja s međunarodnim sudjelovanjem. Zagreb, 443-448.
Smolders, E., L. Sweeck, R. Merckx, A. Cremers, 1997: Cationic interactions inradiocaesium uptake from solution by spinach, Journal of Environmental Radioactivity, 34: 161–170.
Sokolov, V.E., I.N. Ryabov, I.A. Ryabtsev, F.A. Tikomirov, V.A. Shevechcnko, A.I. Taskaev, 1990: Ecological and genetic consequences of the Chernobyl atomic power plant accident, ICSU, SCOPERADPATH, Lancaster Radpath Meeting, 26-30 March 1990, University of Leicester, UK.

ŠUMARSKI LIST 7-8/2021 str. 33     <-- 33 -->        PDF

Stauton, S., C. Dumat, A. Zsolnay, 2002: Possible role of organic matter in radio-caesium adsorption in soils, Journal of Environmental Radioactivity, 58:163-173.
UNSCEAR, 1993: Sources and effects of ionizing radiation. United Nations Scientific Committee on the Effects of Atomic Radiation, Report to the General Assembly, with scientific annexes.
Zhu, Y.G., E. Smolders, 2000: Plant uptake of radiocaesium: a review of mechanisms, regulation and application, Journal of Experimental Botany, 51: 1635–1645.
Zibold, G., E. Klemt, 2005, Ecological half-lives of 137Cs and 90Sr in forest and freshwater ecosystems, Radioprotection, 40: 497-502.
Zibold, G., E. Klemt, I. Konopleva, A. Konoplev, 2009: Influence of fertilizing on the 137Cs soil–plant transfer in a spruce forest of Southern Germany, Journal of Environmental Radioactivity, 100: 489–496.
Yamagata, N., S. Matsuda, M. Chiba, 1969: Radioecology of 137Cs and 90Sr in a forest. J. Rad. Res., 10: 107-112.
Windows, Statistica 7.1.,
Istraživanje aktivnosti 137Cs i 40K provedeno na stablima obične jele (Abies alba Mll.) iz Like uključilo je terensko prikupljanje uzoraka (kolutove stabala sa tri visine i razdvajanih na koru i godove, korijen, iglice, vršne izbojke i tla uz oborena stabla), laboratorijsko mjerenje aktivnosti 137Cs i 40K u prikupljenim uzorcima gamaspektrometrijskom metodom i statističku obradu dobivenih podataka. Ispitivane su radijalne i vertikalne raspodjele cezija (137Cs) u stablima koji je kontaminirao šumske ekosustave putem daljinskog atmosferskog transporta i oborina kao posljedice nuklearnih pokusa, kao i havarije u Černobilu. Na dužoj vremenskoj skali utvrđena je varijabilnost distribucije 137Cs u tkivima stabala obične jele koja je u ispitivanom razdoblju ovisila i o vremenu poluraspada, dok je sezonska dinamika zavisila o razini fiziološke aktivnosti, te od karakteristika i funkcija biljnih tkiva. Najveća koncentracija aktivnosti 137Cs utvrđena je u kori i fiziološki najaktivnijim dijelovima stabala obične jele (vršnim izbojcima i iglicama). Najveća aktivnost 137Cs u godovima obične jele izmjerena je na najnižoj visini stabala. Ovo istraživanje doprinosi razumijevanju ponašanja 137Cs koji je ušao u organizam edifikatorske vrste drveća u šumskom ekosustavu, kao i njegove distribucije u vremenu i prostoru.
Ključne riječi: ekosustav, radionuklidi, distribucija, bioindikatori, obična jela